JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 258528 10.1155/2013/258528 258528 Research Article Permeability Models for Magma Flow through the Earth's Mantle: A Lie Group Analysis Mindu N. Mason D. P. Mahomed Fazal M. Centre for Differential Equations Continuum Mechanics and Applications and School of Computational and Applied Mathematics University of the Witwatersrand Johannesburg Private Bag 3 Wits 2050 South Africa wits.ac.za 2013 31 3 2013 2013 29 12 2012 01 03 2013 2013 Copyright © 2013 N. Mindu and D. P. Mason. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The migration of melt through the mantle of the Earth is governed by a third-order nonlinear partial differential equation for the voidage or volume fraction of melt. The partial differential equation depends on the permeability of the medium which is assumed to be a function of the voidage. It is shown that the partial differential equation admits, as well as translations in time and space, other Lie point symmetries provided the permeability is either a power law or an exponential law of the voidage or is a constant. A rarefactive solitary wave solution of the partial differential equation is derived in the form of a quadrature for the exponential law for the permeability.

1. Introduction

The one-dimensional migration of melt upwards through the mantle of the Earth is governed by the third-order nonlinear partial differential equation (1)ϕt+z[K(ϕ)(1-2ϕtz)]=0, where ϕ(t,z) is the voidage or volume fraction of melt, t is time, z is the vertical spatial coordinate, and K is the permeability of the medium. The special case in which K(ϕ) is a power law, (2)K(ϕ)=ϕn, has been studied extensively, and solitary wave solutions have been derived . In this paper K(ϕ) will initially not be specified. For arbitrary forms of K(ϕ), (1) does not depend explicitly on t and z, and therefore it admits the Lie point symmetries (3)X1=t,X2=z. We will determine the forms of K(ϕ) for (1) to admit other Lie point symmetries besides the Lie point symmetries (3). This would be a significant property for (1) to posses because invariant solutions could then be constructed. One of the forms obtained for K(ϕ) is an exponential law relating the permeability to the voidage. We will derive a new rarefactive solitary wave solution of (1) with the exponential law for the permeability.

The variables ϕ, t, z, and K(ϕ) in (1) are dimensionless. The voidage ϕ(t,z) is scaled by the background voidage ϕ0. The background state is therefore defined by ϕ=1. The characteristic length in the z-direction, which is vertically upwards, is the compaction length δc defined by (4)δc=[K(ϕ0)(ξ+(4/3)η)μ]1/2, where μ is the coefficient of shear viscosity of the melt and ξ and η are the bulk and shear viscosity of the solid matrix, respectively. We will assume that ξ and η are constants as did Barcilon and Richter . Scott and Stevenson  assume that ξ and η are power laws of the voidage ϕ. The characteristic time is t0 defined by (5)t0=ϕ0gΔρ[μ(ξ+(4/3))ηK(ϕ0)]1/2, where g is the acceleration due to gravity and Δρ is the difference between the density of the solid matrix and the density of the melt. The permeability is scaled by K(ϕ0) and therefore (6)K(1)=1. When the voidage is zero, the permeability must also be zero and therefore (7)K(0)=0.

In the derivation of (1), it is assumed that the background voidage satisfies ϕ01. An outline of the derivation of (1) when K(ϕ) satisfies the power law (2) is given by Nakayama and Mason . The derivation is readily extended to the general case in which K=K(ϕ).

An outline of the paper is as follows. In Section 2 the Lie point symmetries of (1) are investigated and the forms of K(ϕ) are determined for (1) to admit, as well as the Lie point symmetries (3), other Lie point symmetries. In Section 3 a new rarefactive solitary wave solution of (1) is obtained when K(ϕ) depends on ϕ through an exponential law. Finally, the conclusions are summarized in Section 4.

2. Lie Point Symmetries

In this section K(ϕ) will not be specified initially. We will investigate the Lie point symmetries of (1) and the forms of K(ϕ) for these symmetries to exist.

Equation (1) is as follows: (8)ϕt+dKdϕϕz-dKdϕϕzϕtz-Kϕtzz=0, where a subscript denotes partial differentiation. We look for Lie point symmetries of the form (9)X=ξ1(t,z,ϕ)t+ξ2(t,z,ϕ)z+η(t,z,ϕ)ϕ. The coefficients ξ and η in the Lie point symmetry (9) should not be confused with the bulk and shear viscosities in the characteristic quantities (4) and (5). The invariance criterion is (10)X(ϕt+dKdϕϕz-dKdϕϕzϕtz-Kϕtzz)|(8)=0, where the third prolongation of X is of the form (11)X=ξ1t+ξ2z+ηϕ+ζ1ϕt+ζ2ϕz+ζ12ϕtz+ζ122ϕtzz+. The remaining terms in X are not required in (10) and (12)ζi=Di(η)-ϕsDi(ξs),i=1,2,ζij=Dj(ζi)-ϕisDj(ξs),i=1,2,ζijk=Dk(ζij)-ϕijsDk(ξs),i=1,2, with summation over the repeated index, s, from 1 to 2. The total derivatives D1 and D2 are defined by (13)D1=Dt=t+ϕtϕ+ϕttϕt+ϕztϕz+,D1=Dz=z+ϕzϕ+ϕtzϕt+ϕzzϕz+.

We replace ϕtzz in the determining equation (10) using the partial differential equation (8) and then separate the determining equation according to the following partial derivatives of ϕ: (14)ϕttz:ξ1z=0,ϕzϕttz:  ξ1ϕ=0,ϕzzz:ξ2t=0,ϕtϕzzz:ξ2ϕ=0,ϕzz:2ηtϕ=0,ϕtϕzz:2ηϕ2=0, It follows directly from (14) that (15)ξ1=f(t),ξ2=g(z),η=ϕA(z)+B(t,z), where f(t), g(z), A(z), and B(t,z) still have to be determined. Using (15), the invariance criterion separates further into the following system of equations: (16)ϕzϕtz:[ϕd2Kdϕ2-ϕK(dKdϕ)2+dKdϕ]A(z)+[d2Kdϕ2-1K(dKdϕ)2]B(t,z)=0,(17)ϕtz:[ϕdKdϕ+2K(ϕ)]dAdz+dKdϕBz-K(ϕ)d2gdz2=0,(18)ϕtϕz:dKdϕdAdz=0,(19)ϕz:[ϕd2Kdϕ2-ϕK(dKdϕ)2]A(z)+[d2Kdϕ2-1K(dKdϕ)2]B(t,z)-dKdϕ2Btz+dKdϕ(dfdt+dgdz)=0,(20)ϕt:K(ϕ)d2Adz2+ϕKdKdϕA(z)+1KdKdϕB(t,z)-2dgdz=0,(21)Remainder:ϕdKdϕdAdz-K(ϕ)3Btz2+dKdϕBz+Bt=0.

For arbitrary forms of K(ϕ), (16) to (21) are satisfied by (22)A(z)=0,B(t,z)=0,f(t)=c1,g(z)=c2, where c1 and c2 are constants and (1) admits the Lie point symmetries (3). We will seek possible forms of K(ϕ) for which (1) admits, as well as (3), other Lie point symmetries. From (18) we see that there are two general cases depending on whether the permeability, K, depends on ϕ or is constant.

2.1. Permeability Depends on Voidage

Consider first the case in which K(ϕ) depends on the voidage, ϕ; that is, the permeability is not constant, so that (23)dKdϕ0. Then, from (18), A(z)=A0, where A0 is a constant. By differentiating (20) with respect to t we find that B=B(z), and from (21) it follows that B(z)=B0, where B0 is a constant. Hence, from (15), (24)η(ϕ)=A0ϕ+B0. From (17), (25)ξ2=g(z)=c1z+c2, where c1 and c2 are constants. Also from (16) and (19), (26)dfdt=A0-c1, and therefore (27)ξ1=f(t)=(A0-c1)t+c3, where c3 is a constant. Equations (16) to (21) reduce to (28)(A0ϕ+B0)(d2Kdϕ2-1K(dKdϕ)2)+A0dKdϕ=0,(29)(A0ϕ+B0)dKdϕ-2c1K=0. It is readily verified that if K(ϕ) satisfies (29), then K(ϕ) satisfies (28) identically. Equation (28) therefore does not need to be considered further.

When (23) is satisfied Lie point symmetries of (1) exist provided that K(ϕ) satisfies (29) and are given by (24), (25), and (27). This case separates into two subcases depending on whether A00 or A0=0.

2.1.1. The Case <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M99"><mml:msub><mml:mrow><mml:mi>A</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo mathvariant="bold">≠</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula>

Consider first A00. The general solution of the ordinary differential equation (29) for K(ϕ) is (30)K(ϕ)=K0(A0ϕ+B0)2c1/A0, where K0 is a constant. But since K(0)=0, it follows that B0=0, and since K(1)=1 we obtain (31)K(ϕ)=ϕn,n=2c1A0. Since K is not a constant, n0. Equations (27), (25), and (24) become (32)ξ1(t)=c1n(2-n)t+c3,ξ2(z)=c1z+c2,η(ϕ)=2c1nϕ. The three Lie point symmetries of (1) with the power law (30) for K(ϕ) are presented in Table 1. The results agree with those derived by Maluleke and Mason [7, 9] for the generalized magma equation with m=0, where m is the exponent in the power law relating the bulk and shear viscosities of the solid matrix to the voidage.

The permeability K(ϕ) and the corresponding Lie point symmetries of (1).

Permeability Lie point symmetry
X 1 = ( 2 - n ) t t + n z z + 2 ϕ ϕ
K ( ϕ ) = ϕ n X 2 = z
X 3 = t

X 1 = - n t t + n z z + 2 ϕ
K ( ϕ ) = exp [ n ( ϕ - 1 ) ] X 2 = z
X 3 = t

X 1 = f ( t ) t
X 2 = ϕ ϕ
X 3 = z
K ( ϕ ) = 1 X 4 = cosh ( 2 z ) z + ϕ sinh ( 2 z ) ϕ
X 5 = sinh ( 2 z ) z + ϕ cosh ( 2 z ) ϕ
X 6 = B ( t , z ) ϕ , where B(t,z) satisfies
B t - 3 B t z 2 = 0

2.1.2. The Case <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M132"><mml:msub><mml:mrow><mml:mi>A</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo mathvariant="bold">=</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula>, <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M133"><mml:msub><mml:mrow><mml:mi>B</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo mathvariant="bold">≠</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula>

When A0=0 but B00, the general solution of (29) is (33)K(ϕ)=K0exp(nϕ),n=2c1B0, where K0 is a constant. Since K(ϕ) is not constant, n0, and because the permeability increases as the voidage increases, n>0. Also, K(1)=1 and therefore (34)K(ϕ)=exp(n(ϕ-1)),n>0. When ϕ=0, then (35)K(ϕ)=exp(-n)0. If n is large, K(0) is small. However, the exponential law (34) for K(ϕ) does not satisfy the condition K(0)=0. It is not a suitable model for physical phenomena with small values of the voidage. For instance, it would not be suitable to describe compressive solitary waves which contain ϕ=0 [8, 1012]. It is suitable for describing rarefactive solitary waves which satisfy ϕ1 and this will be considered in Section 3. Equations (27), (25), and (24) become (36)ξ1(t)=-c1t+c3,ξ2(z)=c1z+c2,η(ϕ)=2c1n. The three Lie point symmetries of (1) with the exponential law (34) for K(ϕ) are presented in Table 1.

2.2. Constant Permeability

Finally, consider constant permeability, K(ϕ)=K0. Since K(1)=1, it follows that K0=1. The model does not satisfy K(0)=0 and it cannot be used to describe physical phenomena in which the permeability depends on voidage.

When K(ϕ)=1, (16) to (21) are reduced to (37)2dAdz-d2gdz2=0,(38)d2Adz2-2dgdz=0,(39)Bt-3Btz2=0. Integrating (37) and (38) once with respect to z gives (40)2A(z)-dgdz=c1,(41)dAdz-2g(z)=c2, where c1 and c2 are constants. By eliminating A(z), we obtain (42)d2gdz2-4g=2c2 and therefore (43)g(z)=c3cosh(2z)+c4sinh(2z)-c22, where c3 and c4 are constants. From (40), (44)A(z)=c3sinh(2z)+c4cosh(2z)+c12. Thus, from (15) (45)ξ1(t)=f(t),ξ2(z)=c3cosh(2z)+c4sinh(2z)-c22,η(t,z,ϕ)=ϕ(c3sinh(2z)+c4cosh(2z)+c12)+B(t,z), where f(t) is arbitrary and B(t,z) satisfies (39) which is the partial differential equation (1) with K(ϕ)=1. The Lie point symmetries of (1) with constant permeability are presented in Table 1. The results agree with those derived by Maluleke and Mason [7, 9] for n=0 and m=0.

There are therefore three forms of K(ϕ) for which (1) has Lie point symmetries in addition to (3), namely, the power law (31), the exponential law (34), and constant permeability.The Lie point symmetries of (1) with the three forms of K(ϕ) are presented in Table 1. Equation (1) with the power law has been studied in detail. In Section 3 we will consider the exponential law and investigate rarefactive solitary wave solutions of (1) with K(ϕ) given by (34).

3. Rarefactive Solitary Wave

When the permeability satisfies the exponential law (34), the partial differential equation (1) becomes (46)ϕt+z[exp[n(ϕ-1)](1-2ϕtz)]=0. We now derive a rarefactive solitary wave solution, with ϕ1, of the partial differential equation (46).

The solution ϕ=Φ(t,z) is an invariant solution of (46) provided that (47)X(ϕ-Φ(t,z))|ϕ=Φ=0, where X is a linear combination of the Lie point symmetries of (46). The Lie point symmetries of (46) are given in Table 1. Consider the invariant solution generated by the linear combination (48)X=c1t+c2z, (47) becomes (49)c1Φt+c2Φz=0, where c1 and c2 are constants. The general solution of (49) is readily derived, and since ϕ=Φ(t,z) we obtain (50)ϕ=ψ(ζ),ζ=z-ct, where c=c2/c1. The group invariant solution (50) is a travelling wave solution and the constant c is the dimensionless speed of the wave.

Substitute (50) into (46). This gives the third-order ordinary differential equation (51)cdψdζ-ddζ[exp[n(ψ-1)](1+cd2ψdζ2)]=0, and integrating (51) once with respect to ξ we obtain (52)cd2ψdζ2=Aexp[n(ψ-1)]+cψexp[n(ψ-1)]-1, where A is a constant of integration. Since the right hand side of (52) depends only on ψ, we integrate (52) with respect to ψ. Now (53)d2ψdζ2=12ddψ((dψdζ)2) and (52) becomes (54)c2ddψ((dψdζ)2)=Aexp[n(ψ-1)]+cψexp[n(ψ-1)]-1, but (55)ψdψexp[n(ψ-1)]=-1nexp[n(ψ-1)],ψdψexp[n(ψ-1)]=-(1+nψ)n2exp[n(ψ-1)], and integrating (54) once with respect to ψ, we obtain (56)c2(dψdζ)2=-Anexp[n(ψ-1)]-c(1+nψ)n2exp[n(ψ-1)]-ψ+B, where B is a constant.

In order to obtain the three arbitrary constants, A, B, and c, we impose three boundary conditions suitable for a rarefactive solitary wave. The background state is ψ=1. Three boundary conditions for a rarefactive solitary wave are (57)ψ=1:dψdζ=0,d2ψdζ2=0,ψ=Ψ:dψdζ=0.

The amplitude of the solitary wave is Ψ-1 and ψ(ζ) has a local maximum when ψ=Ψ. Using (52) and (56), the three boundary conditions give (58)An+(1+n)cn2+1-B=0,A+c-1=0,Anexp[n(Ψ-1)]+c(1+nΨ)n2exp[n(Ψ-1)]+Ψ-B=0. Hence, (59)c=n[(n(Ψ-1)-1)exp[n(Ψ-1)]+1exp[n(Ψ-1)]-n(Ψ-1)-1] and, expressed in terms of c, (60)A=1-c,B=1n2(c+n+n2). Equation (56) becomes (61)(dψdζ)2=f(ψ), where (62)f(ψ)=2cn2[(n+c+nc(ψ-1))n+c-n2(ψ-1)-(n+c+nc(ψ-1))exp[-n(ψ-1)]ψn2]. Now from (53) and the boundary conditions (57), (63)f(1)=0,dfdψ(1)=0,f(Ψ)=0, and since the left hand side of (61) is nonnegative, for a solitary wave to exist it is necessary that f(ψ)>0 for 1<ψ<Ψ.

Before proceeding further with the solution, we first investigate the properties of f(ψ). In order to do that we define (64)F(ψ)=n[(n(ψ-1)-1)exp[n(ψ-1)]+1exp[n(ψ-1)]-n(ψ-1)-1]. From (59), it follows that (65)F(Ψ)=c. Equation (62) can be written in terms of F(ψ) as follows (66)f(ψ)=2n2[1-(1+n(ψ-1))exp[-n(ψ-1)]]×[F(Ψ)-F(ψ)F(Ψ)]. Now, (67)dFdψ=4n2exp[n(ψ-1)][exp[n(ψ-1)]-n(ψ-1)-1]2×[sinh2(n(ψ-1)2)-(n(ψ-1)2)2] and since sinh2x>x2 for x0, it follows that for n0 and ψ1, (68)dFdψ>0. Hence, F(ψ) is an increasing function of ψ. Also, (69)F(1)=n. Consider first n>0. Since F(ψ) is an increasing function of ψ, it follows that F(ψ)>n>0 and F(Ψ)-F(ψ)>0 for 1<ψ<Ψ. Also, it can be verified that (70)1-(1+n(ψ-1))exp[n(ψ-1)]>0, for ψ>1. Hence from (66), f(ψ)>0 for 1<ψ<Ψ and a rarefactive solitary wave solution exists. Consider next n<0. Then from (64), F(ψ)0 as ψ, and since F(1)=n<0 and F(ψ) is an increasing function of ψ, it follows that F(Ψ)<0. We still have F(Ψ)-F(ψ)>0 by (68) and (70) is still satisfied for n<0. Hence from (66), f(ψ)<0 for 1<ψ<Ψ and a rarefactive solitary wave solution does not exist.

The difference between F(ψ) for n>0 and n<0 is illustrated in Figure 1 for n=-0.1,0.5,1, and 2. For n<0 the permeability decreases as the voidage, ϕ, increases, which is generally not observed physically. When the permeability satisfies the power law (30), a rarefactive solitary wave solution exists for n>1 and does not exist for 0n1 .

The function f(ψ) defined by (62) plotted against ψ for n=-0.1,0.5,1, and 2 for Ψ=5.

Now from (61) and (62), the rarefactive solitary wave solution is given by the following: (71)ζ=±n(c2)2×ψΨ(×exp(-n(ψ-1))]1/2)-1(dψ)×(exp(-n(ψ-1))]1/2[(n+c+nc(ψ-1))n+c-n2(ψ-1)-(n+c+nc(ψ-1))×exp(-n(ψ-1))n2(n+c+nc(ψ-1))]1/2)-1), where the wave speed c is given by (59). In Figure 2, the solitary wave solution (71) is plotted against ζ for a range of values of n. In order to compare the solutions for different values of n, the length ζ is scaled in all cases with the same characteristic length, namely, δc calculated for n=1. We see that the width of the solitary wave increases as n increases. The increase in the permeability has the effect of spreading the solitary wave.

Rarefactive solitary wave for n=1,2, and 3 with Ψ=5. The length ζ is scaled with the characteristic length δc for n=1.

Consider now the dependence of the wave velocity, c, on the amplitude of the solitary wave, Ψ-1. Since F(Ψ)=c, it follows from (68) that (72)dcdΨ>0. The velocity increases as the amplitude of the solitary wave increases, and therefore larger amplitude waves travel faster. This property also holds for the solitary waves described by the power law (31) for n>1 . Equation (72) is a special case of the general result that F(ψ) is an increasing function of ψ for 1ψΨ which was central to the proof of the existence of solitary wave solutions for n>0. There is therefore a close connection between the property that larger amplitude solitary waves travel faster and the existence of solitary wave solutions. We also determine from (69) the limiting value (73)Ψ-1=0:c=n. This limiting value also holds for the dimensionless velocity of rarefactive solitary waves for the power law (31) with n>1.

In Figure 3, the dimensionless wave velocity, c, of the solitary wave, given by (59), is plotted against the amplitude Ψ-1 for n=1 and 2. In order to compare c for different values of n, the wave velocity is scaled for each value of n by the same characteristic velocity v0=δc/t0 calculated for n=1. We see that c increases steadily with Ψ-1 for both cases in agreement with (72). We also see that as n increases the velocity c increases. The increase in the permeability as n increases allows the melt to propagate at a greater speed through the solid matrix.

Velocity, c, of a rarefactive solitary wave for n=1 and 2 plotted against amplitude Ψ-1. Both velocities are scaled by the characteristic velocity, δc/t0, for n=1.

4. Conclusions

One of the functional forms which the permeability K(ϕ) must satisfy for the partial differential equation (1) to possess Lie point symmetries besides translation in time and space is the power law. This relation between the permeability and the voidage has been studied extensively to model the migration of melt through the mantle of the Earth. Mainly, travelling wave solutions have been considered, but more general group invariant solutions which include the third Lie point symmetry have been investigated by Harris and Clarkson .

Unlike the power law, the exponential law relating the permeability to the voidage is not valid at ϕ=0 because, when the voidage vanishes, the permeability should also vanish. It can be used when the migration of melt through the Earth’s mantle is modelled as a rarefactive solitary wave with ϕ1. For both the power law and exponential law, larger amplitude waves travel faster. There is a close connection between this property and the existence of solitary wave solutions.

The conservation laws for the partial differential equation (1) with the power law for K(ϕ) have been investigated by several authors [3, 6, 9]. The number of conservation laws is finite, except possibly for n=-1 which is not physical. This indicates that the solitary waves are not solitons. The conservation laws for the partial differential equation (46) with the exponential law for K(ϕ) still need to be investigated. If it is found that the number of conservation laws is finite it would indicate that the rarefactive solitary waves with the exponential law are also not solitons.

We have assumed that the bulk viscosity, ξ, and shear viscosity, η, of the solid matrix are constant. A large amount of research has been performed on models in which ξ+(4/3)η is related to the voidage by a power law with exponent m [1, 412]. The Lie group analysis for the permeability considered here could be extended to include the viscosity of the solid matrix.

Acknowledgments

N. Mindu thanks Shanduka and Tomorrow Trust for financial support. D. P. Mason thanks the National Research Foundation, Pretoria, South Africa, for financial support. The authors thank a referee for going through the paper thoroughly and suggesting improvements to the text.

Scott D. R. Stevenson D. J. Magma solitons Geophysical Research Letters 1984 4 1161 1164 Mckenzie D. The generation and compaction of partially molten rock Journal of Petrology 1984 25 3 713 765 2-s2.0-0021360880 10.1093/petrology/25.3.713 Barcilon V. Richter F. M. Nonlinear waves in compacting media Journal of Fluid Mechanics 1986 164 429 448 10.1017/S0022112086002628 MR844680 ZBL0587.76165 Takahashi D. Satsuma J. Explicit solutions of magma equation Journal of the Physical Society of Japan 1988 57 2 417 421 10.1143/JPSJ.57.417 MR936112 Nakayama M. Mason D. P. Rarefactive solitary waves in two-phase fluid flow of compacting media Wave Motion 1992 15 4 357 392 10.1016/0165-2125(92)90054-6 MR1166186 ZBL0761.76097 Harris S. E. Conservation laws for a nonlinear wave equation Nonlinearity 1996 9 1 187 208 10.1088/0951-7715/9/1/006 MR1374003 ZBL0892.35098 Maluleke G. H. Mason D. P. Optimal system and group invariant solutions for a nonlinear wave equation Communications in Nonlinear Science and Numerical Simulation 2004 9 1 93 104 10.1016/S1007-5704(03)00018-2 MR2036559 ZBL1036.35010 Harris S. E. Clarkson P. A. Painlevé analysis and similarity reductions for the magma equation Symmetry, Integrability and Geometry 2006 2 paper 068 17 10.3842/SIGMA.2006.068 MR2264884 ZBL1132.35328 Maluleke G. H. Mason D. P. Derivation of conservation laws for a nonlinear wave equation modelling melt migration using Lie point symmetry generators Communications in Nonlinear Science and Numerical Simulation 2007 12 4 423 433 10.1016/j.cnsns.2005.05.010 MR2283344 ZBL1110.35076 Nakayama M. Mason D. P. Compressive solitary waves in compacting media International Journal of Non-Linear Mechanics 1991 26 5 631 640 10.1016/0020-7462(91)90015-L MR1144416 ZBL0754.76085 Nakayama M. Mason D. P. On the existence of compressive solitary waves in compacting media Journal of Physics A 1994 27 13 4589 4599 MR1294961 10.1088/0305-4470/27/13/032 ZBL0841.76087 Nakayama M. Mason D. P. On the effect of background voidage on compressive solitary waves in compacting media Journal of Physics A 1995 28 24 7243 7261 MR1381417 10.1088/0305-4470/28/24/021 ZBL0877.35110