Abstract and Applied Analysis

Volume 2014, Article ID 585167, 13 pages

http://dx.doi.org/10.1155/2014/585167

## Derivation of Conservation Laws for the Magma Equation Using the Multiplier Method: Power Law and Exponential Law for Permeability and Viscosity

Programme for Differential Equations, Continuum Mechanics and Applications, School of Computational and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050, South Africa

Received 24 February 2014; Accepted 1 April 2014; Published 4 May 2014

Academic Editor: Chaudry M. Khalique

Copyright © 2014 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.

#### Abstract

The derivation of conservation laws for the magma equation using the multiplier method for both the power law and exponential law relating the permeability and matrix viscosity to the voidage is considered. It is found that all known conserved vectors for the magma equation and the new conserved vectors for the exponential laws can be derived using multipliers which depend on the voidage and spatial derivatives of the voidage. It is also found that the conserved vectors are associated with the Lie point symmetry of the magma equation which generates travelling wave solutions which may explain by the double reduction theorem for associated Lie point symmetries why many of the known analytical solutions are travelling waves.

#### 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 where is the voidage or volume fraction of melt, is time, is the vertical spatial coordinate, is the permeability of the medium, and is the viscosity of the matrix phase. The variables , , and and the physical quantities and in (1) are dimensionless. The voidage is scaled by the background voidage . The background state is therefore defined by . The characteristic length in the -direction, which is vertically upwards, is the compaction length defined by where is the coefficient of shear viscosity of the melt. The characteristic time is defined by where 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 and therefore When the voidage is zero the permeability must also be zero and therefore The viscosity is scaled by so that and will be infinite because the matrix viscosity is infinite when the voidage vanishes. In the derivation of (1) it is assumed that the background voidage satisfies .

The partially molten medium consists of a solid matrix and a fluid melt which are modelled as two immiscible fully connected fluids of constant but different densities. The density of the melt is less than the density of the solid matrix and the melt migrates through the compacting medium by the buoyancy force due to the difference in density between the melt and the solid matrix. Changes of phase are not included in the model. It is assumed that the melting has occurred and only migration of the melt under gravity is described by (1) [1].

In the model proposed by Scott and Stevenson [2], consider where and . Harris and Clarkson [3] have investigated this model using Painleve analysis. Mindu and Mason [4] showed that the magma equation also admits Lie point symmetries other than translations in time and space if the permeability is in the form of an exponential law: Conservation laws for (1) when the permeability and matrix viscosity satisfy the power laws (7) have been obtained using the direct method by Barcilon and Richter [5] and Harris [6] and using Lie point symmetry generators by Maluleke and Mason [7].

In this paper we will derive the conservation laws for the partial differential equation (1) using the multiplier method. We will consider power laws given by (7) and also the exponential laws where and , relating the permeability and matrix viscosity to the voidage. The permeability increases as the voidage increases while the viscosity of the matrix decreases as the voidage increases. The exponential laws are not suitable models when the voidage is small because and . They are suitable for describing rarefaction for which .

An outline of the paper is as follows. In Section 2 we present the formulae and theory that we will use in the paper. In Section 3 conservation laws for the magma equation, with power laws relating the permeability and viscosity to the voidage, are derived using the multiplier method. Further in Section 4 conservation laws for the magma equation, with exponential laws relating the permeability and viscosity to the voidage, are derived using the multiplier method. Finally the conclusions are summarized in Section 5.

#### 2. Formulae and Theory

Consider an th order partial differential equation in the variables , where denotes the collection of th-order partial derivatives of . The equation evaluated on the surface given by (10), where runs from 1 to and is the total derivative defined by is called a conservation law for the differential equation (10). The vector is a conserved vector for the partial differential equation and are its components. Thus, a conserved vector gives rise to a conservation law. A Lie point symmetry generator where runs from 1 to , is said to be associated with the conserved vector for the partial differential equation (10) if [8, 9] The association of a Lie point symmetry with a conserved vector can be used to integrate the partial differential equation twice by the double reduction theorem of Sjöberg [10].

Conserved vectors for a partial differential equation can be generated from known conserved vectors and Lie point symmetries of the partial differential equation. For where runs from 1 to , is a conserved vector for the partial differential equation although it may be a linear combination of known conserved vectors [8, 9].

We now present the multiplier method for the derivation of conservation laws for partial differential equations. We will outline its application to the partial differential equation (1) in two independent variables.

(1) Multiply the partial differential equation (1) by the multiplier, , to obtain the conservation law where is the partial differential equation (1), and , and The multiplier depends on , , , and the partial derivatives of . The more derivatives included in the multiplier the wider the range of conserved vectors that can be derived.

(2) The determining equation for the multiplier is obtained by operating on (16) by the Euler operator defined by [11] Since the Euler operator annihilates divergence expressions this gives [11]

(3) The determining equation (19) is separated by equating the coefficients of like powers and products of the derivatives of because is an arbitrary function.

(4) When is a solution of the partial differential equation, , (16) becomes a conservation law. The condition is imposed on (16). The product of the multiplier and the partial differential equation is then written in conserved form by elementary manipulations. This yields the conserved vectors by setting all the constants equal to zero except one in turn.

#### 3. Conservation Laws for the Magma Equation with Power Law Permeability and Viscosity by the Multiplier Method

When the permeability and viscosity are related to the voidage by the power laws (7) the magma equation becomes

##### 3.1. Lower Order Conservation Laws

In order to derive conservation laws for (20) consider first a multiplier of the form A multiplier for the partial differential equation has the property where The determining equation for the multiplier is where is defined by (18). Separating (24) with respect to products and powers of the partial derivatives of we obtain the following system of equations: Equation (26) is the same as (25). It is readily verified that every solution of (25) is a solution of (27). We therefore need to consider only (25). The general solution of (25) is

There are several cases to consider depending on the values of and . The special cases are illustrated as lines and points in the plane in Figure 1.

*(i)**, **, **.* From (22) and (28),
Equation (30) is satisfied for arbitrary functions . When is a solution of the partial differential equation (20), then
Hence, any conserved vector of the partial differential equation (20) with and satisfying the conditions of this case and with multiplier of the form is a linear combination of the two conserved vectors
The conserved vector (32) is the elementary conserved vector.

*(ii)**, **.* Proceeding as before we obtain
The conserved vector (34) is the elementary conserved vector with . The multiplier for (35) is, from (29),

*(iii)**, **.* We find that
The conserved vector (37) is the elementary conserved vector with . The multiplier for (38) is, from (28),

*(iv)**, **.* We obtain
The conserved vector (40) is the elementary conserved vector with , . The multiplier for (41) is given by (29).

*(v)**.* We obtain
The conserved vector (42) is the elementary conserved vector with . The multiplier for (43) is (39).

*(vi)**, **, **.* We obtain
The conserved vector (44) is the elementary conserved vector with . The multiplier for (45) is

##### 3.2. The Search for Higher Order Conservation Laws

We now consider a multiplier of the form As before the determining equation for the multiplier is where is given by (23). By equating the coefficient of the highest order derivative term, , to zero in (48) we have and therefore Hence, (47) does not give a new multiplier or a new conserved vector.

Consider next the multiplier The determining equation for the multiplier is where is given by (23). By Equating the coefficients of , , and to zero in (52), we obtain the following system of equations: From (55) it follows that Substituting (56) into (53) we find that and therefore Thus, (51) becomes Now substitute (59) into (54) which gives The solution to (60) is where is a constant. Equation (59) becomes

Now substitute (62) into (48) and then equate the coefficient of in (48) to zero. This gives and integrating (63) we have Thus, (62) becomes Lastly, substitute (65) into (48) and equate the coefficient of to zero, which gives It follows from (66) that there are two cases.

*Case 1 (). *If , then from (66) we have . Therefore,
Thus, does not give a new multiplier and therefore new conservation laws will not be derived.

*Case 2 (). *If , then from (66), we have . Thus, (65) becomes
Now substitute (68) into (48) and equate the coefficient of in (48) to zero, which gives
Solving (69) we have
provided that and . Since , these two special cases correspond to the points and on the plane in Figure 1.

Consider first the general case excluding the points and . Equation (68) becomes Substituting (71) into (48) we find that . Hence, Since the multiplier (72) contains two constants, and , it leads to two conserved vectors. The conserved vector corresponding to , is the elementary conserved vector (32). The constants , lead to the conserved vector

The case with and has already been considered. The multiplier is given by (29) and the conserved vectors by (40) and (41).

Consider with and . The differential equation (69) becomes The general solution of (75) is When the multiplier (68) becomes On substituting (77) into the determining equation (48) we find that and the multiplier reduces to The multiplier (78) again contains two arbitrary constants, and . Setting the constants , gives the elementary conserved vector (32). Setting , leads to the conserved vector

Consider next multipliers of the form The determining equation for the multiplier is where is given by (23). Equating the coefficient of in (82) to zero, we find that and therefore which has already been considered.

Harris [6] proved that, except possibly for the two special cases with and with , there are no more independent conserved vectors. She proved this result using the direct method for conservation laws.

The multipliers and the corresponding conserved vectors for the partial differential equation (20) are listed in Table 1. This table was presented by Maluleke and Mason [7] without the multipliers. These conserved vectors agree with the results obtained by Barcilon and Richter [5] and Harris [6].

##### 3.3. Association of Lie Point Symmetries with Conserved Vectors

The Lie point symmetries for the partial differential equation (20) are listed in Table 2. These Lie point symmetries were derived by Maluleke and Mason [7, 12]. Using (14) we will investigate which of the Lie point symmetries are associated with the conserved vectors for the Magma equation (20).

*(i)**, **.* Consider first the Lie symmetry generator
and the elementary conserved vector (32). Applying (14) we find that (85) is associated with the conserved vector (32) provided that , that is, provided that

*(ii)**, **, **.* Consider next the Lie point symmetry generator (85), with the conserved vector (33). Applying (14) we find that (85) is associated with (33) provided that , that is, provided is given by (86).

*(iii)**, **, **.* Now consider
and the conserved vector (45). Applying (14) we find that (87) is associated with (45) provided that , that is, provided that is given by (86).

*(iv)**.* Consider next
and the conserved vector (42). Applying (14) we find that (88) is associated with (42) provided that , that is, provided that is (86).

*(v)**, **.* Consider
and the conserved vector (41). Applying (14) we find that (89) is associated with (41) provided that , that is, provided that

*(vi)**, **.* Consider next
and the conserved vector (35). Applying (14) we find that (91) is associated with (35) provided that , that is, provided that is given by (86).

*(vii)**, **.* Consider
with the conserved vector (38). Applying (14) we find that (92) is associated with (38) provided that , that is, provided that is given by (86).

*(viii)**, **, **.* Consider
with the conserved vector given by (73) and (74). Applying (14) we find that (93) is associated with the conserved vector with components (73) and (74) provided that , that is, provided that is given by (86).

*(ix)**, **, **.* Finally consider the Lie point symmetry (93) and the conserved vector with components (79) and (80). Applying (14) we find that (93) is associated with this conserved vector provided that , that is, provided that is given by (86).

Except for the conserved vector (41) the conserved vectors are all associated with the Lie point symmetry which generates travelling wave solutions. The Lie point symmetry with which the conserved vector (41) is associated contains (86) as a special case. In all cases new conserved vectors are not generated by (15).

Next we derive the conservation laws for the magma equation with an exponential law for the permeability and viscosity using the multiplier method.

#### 4. Conservation Laws for the Magma Equation with an Exponential Law for the Permeability and Viscosity by the Multiplier Method

When the permeability and viscosity are related to the voidage by exponential laws the magma equation becomes

##### 4.1. Lower Order Conservation Laws

In order to derive conservation laws for (94) consider a multiplier of the form (21). A multiplier for the partial differential equation has the property (22), where now The determining equation for the multiplier is given by (24), where is given by (18). Separating (24) with respect to products and powers of the partial derivatives of we obtain the following system of equations: Equation (96) is the same as (97). It is readily verified that every solution of (96) is a solution of (98). We therefore need to consider only (96). The general solution of (96) is

We are considering and and therefore only if and . Proceeding as before we have for various combinations of and different conserved vectors.

*(i)**, **.* This gives the conserved vectors
The conserved vector (101) is the elementary conserved vector. The multiplier for (102) is, from (99),

*(ii)**.* We obtain two conserved vectors
The conserved vector (104) is the elementary conserved vector with multiplier and . The multiplier for (105) is, from (100),

*(iii)** and **.* We again obtain two conserved vectors
The conserved vector (107) is the elementary conserved vector with multiplier and . The multiplier of the conserved vector (108) is by (99)

*(iv)**, **.* Finally we obtain the conserved vectors
The conserved vector (110) is the elementary conserved vector with multiplier and , while the multiplier for (111) is by (99):

##### 4.2. The Search for Higher Order Conservation Laws

We now consider a multiplier of the form The determining equation for the multiplier is (48), where is given by (95). By equating the coefficient of the highest order derivative term, , to zero in (48) we obtain again (49), so that . The multiplier therefore reduces to that of the previous case and new conserved vectors are not derived.

Consider next the multiplier As before, the determining equation for the multiplier is (48), where is given by (95). By equating the coefficients of , , and to zero in (48), the following system of equations is obtained: By using (115) and (117), it is readily shown that (56) again holds. Substituting (56) into (116) gives and therefore Equation (56) now becomes Substituting (120) into the determining equation (48) and then equating the coefficients of in (48) to zero gives and hence The multiplier becomes Finally we substitute (123) back into (48) and equate the coefficient of to zero. This yields There are three cases to consider, , , and . The conserved vectors for , are given by (107) and (108). We now consider the two remaining cases.

*Case 1 (, ). *Then, and (123) reduces to
The multiplier is therefore a function of only which does not yield new conserved vectors.

*Case 2 (, ). *Then, and the multiplier (123) becomes
Equation (126) is substituted back into the determining equation (48) and by equating the coefficient of to zero we obtain
The general solution to (127) is
where since and . Thus, (123) becomes
Finally substituting (129) into the determining equation (48) gives and therefore
Two conserved vectors are obtained since the multiplier (130) contains two arbitrary constants. The constant gives the elementary conserved vector (100) while the constant gives the new conserved vector
which exists if and , .

Consider next multipliers of the form (81). The determining equation for the multiplier is (82), where is given by (95). By equating to zero the coefficient of in (82) we again derive (83) and the multiplier therefore reduces to the form (84) which has already been considered.

The multipliers and the corresponding conserved vectors for the partial differential equation (94) are listed in Table 3. The plane is illustrated in Figure 2.

##### 4.3. Association of Lie Point Symmetries with Conserved Vectors

The Lie point symmetries of the partial differential equation (94) are given in Table 4. We use (14) to investigate which Lie point symmetries of (94) are associated with the conserved vectors for (94).

*(i)**, **, **.* Consider first the Lie point symmetry generator
and the elementary conserved vector (101). We find that (133) is associated with the elementary conserved vector provided that , that is, provided that is given by (86).

*(ii)**, **, *. Consider next the Lie point symmetry generator (133) and the conserved vector (102). It can be verified that (102) is associated with (133) provided that , that is, provided that is given by (86).

*(iii)**, *. Now consider the Lie point symmetry
and the conserved vector (108). Using again (14) we find that (134) is associated with (108) provided that , that is, provided that is given by (86).

*(iv)**, *. Consider the Lie point symmetry
and the conserved vector (111). We find that (135) is associated with (111) provided that , that is, provided that

We see that, except for (135) , the conserved vector is associated with the Lie point symmetry (86) which generates a travelling wave solution. The conserved vector (111) is associated with (136) which includes (86) as a special case. In all cases, (15) does not yield a new conserved vector.

#### 5. Conclusion

In this paper the multiplier method was used to derive the conservation laws for the magma equation for the case in which the permeability and viscosity satisfy a power law. The results agree with those of Harris [6], who derived the conserved vectors using the direct method. Unlike the direct method the functional form of the conserved vector does not need to be assumed with the multiplier method. Instead the variables on which the multiplier depend have to be chosen but this can be done by starting with a simple form and including higher order partial derivatives later to derive higher order conservation laws. The determining equation for the multiplier is readily obtained with the aid of the Euler operator.

Conserved vectors for the magma equation when the permeability and matrix viscosity depend on the voidage by exponential laws were derived using the multiplier method. Their properties are similar to the properties of the conserved vectors for the power law relations.

We investigated the association of Lie point symmetries of the magma equation with the conserved vectors. For all conserved vectors considered except two the associated Lie point symmetry was the Lie point symmetry which generates travelling wave solutions [4, 5].

We were not able to derive new conservation laws for the partial differential equation (20) or determine if the number of conservation laws for (20) is finite or infinite. Harris [6] has proved that except possibly for the two special cases, with and with , there are no more independent conserved vectors. Our results derived using multipliers are consistent with the results of Harris. All known conserved vectors of (20) and also the new conserved vectors for (94) can be derived from multipliers which depend only on and the partial derivatives of with respect to . We find that the multipliers and whose variables ended in odd order partial derivatives of with respect to did not generate new conserved vectors but instead reduced to the multipliers and , respectively. This also applies to the multipliers for the conserved vectors for (94).

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

N. Mindu thanks JCSE Coachlab and the National Research Foundation, Pretoria, South Africa, for financial support. D. P. Mason thanks the National Research Foundation, Pretoria, South Africa, for financial support.

#### References

- D. Mckenzie, “The generation and compaction of partially molten rock,”
*Journal of Petrology*, vol. 25, no. 3, pp. 713–765, 1984. View at Publisher · View at Google Scholar · View at Scopus - D. R. Scott and D. J. Stevenson, “Magma solitons,”
*Geophysical Research Letters*, vol. 4, pp. 1161–1164, 1984. View at Google Scholar - S. E. Harris and P. A. Clarkson, “Painlevé analysis and similarity reductions for the magma equation,”
*SIGMA. Symmetry, Integrability and Geometry. Methods and Applications*, vol. 2, paper 068, 17 pages, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. Mindu and D. P. Mason, “Permeability models for magma flow through the Earth's mantle: a Lie group analysis,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 258528, 8 pages, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. Barcilon and F. M. Richter, “Nonlinear waves in compacting media,”
*Journal of Fluid Mechanics*, vol. 164, pp. 429–448, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. E. Harris, “Conservation laws for a nonlinear wave equation,”
*Nonlinearity*, vol. 9, no. 1, pp. 187–208, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. H. Maluleke and D. P. Mason, “Derivation of conservation laws for a nonlinear wave equation modelling melt migration using Lie point symmetry generators,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 12, no. 4, pp. 423–433, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. H. Kara and F. M. Mahomed, “Relationship between symmetries and conservation laws,”
*International Journal of Theoretical Physics*, vol. 39, no. 1, pp. 23–40, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. H. Kara and F. M. Mahomed, “A basis of conservation laws for partial differential equations,”
*Journal of Nonlinear Mathematical Physics*, vol. 9, supplement 2, pp. 60–72, 2002. View at Publisher · View at Google Scholar · View at MathSciNet - A. Sjöberg, “Double reduction of PDEs from the association of symmetries with conservation laws with applications,”
*Applied Mathematics and Computation*, vol. 184, no. 2, pp. 608–616, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Naz, D. P. Mason, and F. M. Mahomed, “Conservation laws and conserved quantities for laminar two-dimensional and radial jets,”
*Nonlinear Analysis: Real World Applications*, vol. 10, no. 5, pp. 2641–2651, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. H. Maluleke and D. P. Mason, “Optimal system and group invariant solutions for a nonlinear wave equation,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 9, no. 1, pp. 93–104, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet