# Wave Equations in Bianchi Space-Times

Authors: S. Jamal; A. H. Kara; R Narain
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

We investigate the wave equation in Bianchi type III space-time. We construct a Lagrangian of the model, calculate and classify the Noether symmetry generators, and construct corresponding conserved forms. A reduction of the underlying equations is performed to obtain invariant solutions.


---

## Body

## 1. Introduction

The study of partial differential equations (PDEs) in terms of Lie point symmetries is well known and well established [1–5], where these symmetries can be used to obtain, inter alia, exact analytic solutions of the PDEs. In addition, Noether symmetries are also widely investigated and are associated with PDEs that possess a Lagrangian. Noether [6] discovered the interesting   link between symmetries and conservation laws showing that for every infinitesimal transformation admitted by the action integral of a system there exists a conservation law. Investigations have been devoted to understand Noether symmetries of Lagrangians that arise from certain pseudo-Riemannian metrics of interest [7, 8]. Recently, a study was aimed at understanding the effect of gravity on the solutions of the wave equation by solving the wave equation in various space-time geometries [9].In [10], the Bianchi universes were investigated using Noether symmetries. The authors of [11] studied the Noether symmetries of Bianchi type I and III space-times in scalar coupled theories. Therein, they obtained the exact solutions for potential functions, scalar field, and the scale factors, see also [12].We pursue an investigation of the symmetries of the wave equation in Bianchi III space-time. We construct solutions of these equations and find conservation laws associated with Noether symmetries. The plan of the paper is as follows.In Section2, we discuss the procedure to obtain an expression representing Noether symmetries and conservation laws. In Section 3, we derive and classify strict Noether symmetries of the Bianchi III space-time. Also in Section 3, we briefly describe the relation of Noether symmetries to conservation laws. We then illustrate the reduction of the wave equation and obtain invariant solutions.


## 2. Definitions and Notation

We briefly outline the notation and pertinent results used in this work. In this regard, the reader is referred to [13].The convention that repeated indices imply summation is used. Letx=(x1,x2,…,xn)∈Rn be independent variable with xi, and let u=(u1,u2,…,um)∈Rm be the dependent variable with coordinates uα. Furthermore, let π:Rn+m→Rn be the projection map π(x,u)=x. Also, suppose that s:χ⊂Rn→𝒰⊂Rn+m is a smooth map such that π∘s=1χ, where 1χ is the identity map on χ. The r-jet bundle 𝒥r(𝒰) is given by the equivalence classes of sections of 𝒰. The coordinates on 𝒥r(𝒰) are denoted by (xi,uα,…,ui1⋯irα), where 1≤i1≤⋯≤ir≤n and ui1⋯irα corresponds to the partial derivatives of uα with respect to xi1,…,xir. The partial derivatives of u with respect to x are connected by the operator of total differentiation
(2.1)Di=∂∂xi+uiα∂∂uα+uijα∂∂ujα+⋯,i=1,…,n,
as
(2.2)uiα=Di(uα),uijα=DjDi(uα),…
The collection of all first-order derivatives uia will be denoted by u(1). Similarly, the collections of all higher order derivatives will be denoted by u(2),u(3),….Ther-jet bundle on U will be written as 𝒥r(𝒰)={(x,u,u(1),…,u(r))/(x,u)∈𝒰}. We now review the space of differential forms on 𝒥r(𝒰). To this end, let Ωkr(𝒰) be the vector space of differential k-forms on 𝒥r(𝒰) with differential d. A smooth differential k-form on 𝒥r(𝒰) is given by
(2.3)ω=fi1,i2,…,ikdxi1∧dxi2∧⋯∧dxik,
where each component fi1,i2,…,ik∈Ω0r(𝒰). Note that for differential functions f∈Ω0r(𝒰),
(2.4)Df=Djfdxj,
where D is the total differential or the total exterior derivative. Moreover, the total exterior derivative of ω is
(2.5)Dω=Dfi1,i2,…,ikdxi1∧dxi2∧⋯∧dxik,
and by invoking (2.4) one has
(2.6)Dω=Dfi1,i2,…,ikdxj∧dxi1∧dxi2∧⋯∧dxik.
The total differential D has properties analogous to the algebraic properties of the usual exterior derivative d:
(2.7)D(ω∧υ)=Dω∧υ+(-1)kω∧Dυ,
for ω a k-form and υ an l-form and D(Dω)=0. Also, it is known that if D(Dω)=0, then ω is a locally exact k-form, that is, ω=Dυ for some (k-1)-form υ, [14].


### 2.1. Action of Symmetries

Consider anrth-order system of partial differential equations of n independent variables and m dependent variables:
(2.8)Eβ(x,u,u(1),…,u(r))=0,β=1,…,m~.Definition 2.1.
A conserved form of (2.8) is a differential (n-1)-form,
(2.9)ω=Ti(x,u,u(1),…,u(r-1))(∂∂xi⌟(dx1∧⋯∧dxn)),
defined on 𝒥r-1(𝒰) if
(2.10)Dω=0
is satisfied on the surface given by (2.8).Remark 2.2.
When Definition2.1 is satisfied, (2.10) is called a conservation law for (2.8).
It is clear that (2.10) evaluated on the surface (2.8) implies that
(2.11)DiTi=0
on the surface given by (2.8), which is also referred to as a conservation law of (2.8). The tuple T=(T1,…,Tn), Tj∈Ω0r-1(𝒰), j=1,…,n, is called a conserved vector of (2.8).We now review some definitions and results relating to Euler-Lagrange, Lie, and Noether operators ([15, 16] and references therein).Let𝒜=⋃r=0p‍Ω0r(𝒰) for some p<∞. Then 𝒜 is the universal space of differential functions of finite orders.Consider a symmetry operator given by the infinite formal sum:(2.12)X=ξi∂∂xi+ηα∂∂uα+∑s≥1‍ζi1⋯isα∂∂ui1⋯isα,
where ξi, ηα∈𝒜, and the additional coefficients are determined uniquely by the prolongation formulae
(2.13)ζi1⋯isα=Di1⋯Dis(Wα)+ξjuji1…isα,s>1.
In (2.13), Wα is the Lie characteristic function given by
(2.14)Wα=ηα-ξjujα.
In particular, a symmetry operator of the form X~=Qα∂/∂uα+⋯, where Qα∈𝒜, is called a canonical or evolutionary representation of X, and Qα is called its characteristic.An operatorX is said to be a Noether symmetry corresponding to a Lagrangian L∈𝒜, if there exists a vector Bi=(B1,…,Bn),Bi∈𝒜, such that
(2.15)X(L)+LDi(ξi)=Di(Bi).
If Bi=0,(i=1,…,n), then X is referred to as a strict Noether symmetry corresponding to a Lagrangian L∈𝒜. This case is also obtained by setting the Lie derivative on the n-form Ldx1∧⋯∧dxn in the direction of X to zero, that is,
(2.16)ℒXLdx1∧⋯∧dxn=X(Ldx1∧⋯∧dxn)=0,
where ℒ is the Lie derivative operator.In view of the above discussions and definitions, the Noether theorem [6] is formulated as follows.Noether’s Theorem
For any Noether symmetryX corresponding to a given Lagrangian L∈𝒜, there corresponds a vector Ti=(T1,…,Tn),Ti∈𝒜, defined by
(2.17)Ti=Bi-Ni(L),i=1,…,n,
which is a conserved current of the Euler-Lagrange equations δL/δuα=0, α=1,…,m, where δ/δuα is the Euler-Lagrange operator given by
(2.18)δδuα=∂∂uα+∑s≥1‍(-1)sDi1⋯Dis∂∂ui1⋯isα,α=1,…,m,
and the Noether operator associated with the operator X is given by
(2.19)Ni=ξi+Wαδδuiα+∑s≥1‍Di1⋯Dis(Wα)δδuii1⋯isα,i=1,…,n,
in which the Euler-Lagrange operators with respect to derivatives of uα are obtained from (2.18) by replacing uα by the corresponding derivatives, for example,
(2.20)δδuiα=∂∂uiα+∑s≥1‍(-1)sDj1⋯Djs∂∂uij1⋯jsα,i=1,…,n,α=1,…,m.


## 2.1. Action of Symmetries

Consider anrth-order system of partial differential equations of n independent variables and m dependent variables:
(2.8)Eβ(x,u,u(1),…,u(r))=0,β=1,…,m~.Definition 2.1.
A conserved form of (2.8) is a differential (n-1)-form,
(2.9)ω=Ti(x,u,u(1),…,u(r-1))(∂∂xi⌟(dx1∧⋯∧dxn)),
defined on 𝒥r-1(𝒰) if
(2.10)Dω=0
is satisfied on the surface given by (2.8).Remark 2.2.
When Definition2.1 is satisfied, (2.10) is called a conservation law for (2.8).
It is clear that (2.10) evaluated on the surface (2.8) implies that
(2.11)DiTi=0
on the surface given by (2.8), which is also referred to as a conservation law of (2.8). The tuple T=(T1,…,Tn), Tj∈Ω0r-1(𝒰), j=1,…,n, is called a conserved vector of (2.8).We now review some definitions and results relating to Euler-Lagrange, Lie, and Noether operators ([15, 16] and references therein).Let𝒜=⋃r=0p‍Ω0r(𝒰) for some p<∞. Then 𝒜 is the universal space of differential functions of finite orders.Consider a symmetry operator given by the infinite formal sum:(2.12)X=ξi∂∂xi+ηα∂∂uα+∑s≥1‍ζi1⋯isα∂∂ui1⋯isα,
where ξi, ηα∈𝒜, and the additional coefficients are determined uniquely by the prolongation formulae
(2.13)ζi1⋯isα=Di1⋯Dis(Wα)+ξjuji1…isα,s>1.
In (2.13), Wα is the Lie characteristic function given by
(2.14)Wα=ηα-ξjujα.
In particular, a symmetry operator of the form X~=Qα∂/∂uα+⋯, where Qα∈𝒜, is called a canonical or evolutionary representation of X, and Qα is called its characteristic.An operatorX is said to be a Noether symmetry corresponding to a Lagrangian L∈𝒜, if there exists a vector Bi=(B1,…,Bn),Bi∈𝒜, such that
(2.15)X(L)+LDi(ξi)=Di(Bi).
If Bi=0,(i=1,…,n), then X is referred to as a strict Noether symmetry corresponding to a Lagrangian L∈𝒜. This case is also obtained by setting the Lie derivative on the n-form Ldx1∧⋯∧dxn in the direction of X to zero, that is,
(2.16)ℒXLdx1∧⋯∧dxn=X(Ldx1∧⋯∧dxn)=0,
where ℒ is the Lie derivative operator.In view of the above discussions and definitions, the Noether theorem [6] is formulated as follows.Noether’s Theorem
For any Noether symmetryX corresponding to a given Lagrangian L∈𝒜, there corresponds a vector Ti=(T1,…,Tn),Ti∈𝒜, defined by
(2.17)Ti=Bi-Ni(L),i=1,…,n,
which is a conserved current of the Euler-Lagrange equations δL/δuα=0, α=1,…,m, where δ/δuα is the Euler-Lagrange operator given by
(2.18)δδuα=∂∂uα+∑s≥1‍(-1)sDi1⋯Dis∂∂ui1⋯isα,α=1,…,m,
and the Noether operator associated with the operator X is given by
(2.19)Ni=ξi+Wαδδuiα+∑s≥1‍Di1⋯Dis(Wα)δδuii1⋯isα,i=1,…,n,
in which the Euler-Lagrange operators with respect to derivatives of uα are obtained from (2.18) by replacing uα by the corresponding derivatives, for example,
(2.20)δδuiα=∂∂uiα+∑s≥1‍(-1)sDj1⋯Djs∂∂uij1⋯jsα,i=1,…,n,α=1,…,m.


## 3. Bianchi III Space-Time

Consider the Bianchi III metric:(3.1)ds2=-β2dt2+t2L(dx2+e-2ax/Ndy2)+t2L/mdz2.
The wave equation in (3.1) takes the form [17]
(3.2)-1β(2L+1)t2Le-(a/N)xut-1βt2L+1e-(a/N)xutt-aβNte-(a/N)xux+βte-(a/N)xuxx+βte(a/N)xuyy+βt2L+1-(2L/m)e-(a/N)xuzz=0.In [17], some aspects of the wave equation on the Bianchi metric were studied. The multiplier method [1] was adopted to determine some of the conserved densities. This lengthy procedure ultimately leads to the construction of only three symmetries and its associated conserved vectors.In this paper, we investigate the wave equation on the Bianchi III metric using Noether’s theorem and the method of differential forms. We obtain a wide range of results and also perform symmetry reductions of the wave equation for some cases to obtain invariant solutions. For the purposes of Sections3.1 and 3.2, we denote the Lagrangian by L.


### 3.1. The Strict Noether Symmetries of (3.2)

We classify the cases that yield strict Noether symmetries (gauge is zero) of (3.2), via the Lagrangian
(3.3)L=-β2te-(a/N)xux2-β2te(a/N)xuy2-β2t2L+1-(2L/m)e-(a/N)xuz2+12βt2L+1e-(a/N)xut2.
Many of the calculations have been left out as they are tedious—the details are available to the reader in a number of texts that have been cited here.The principle Noether algebra is(3.4)X1=∂y,X2=∂z,X3=∂u,X4=Na∂x+y∂y,X5=-2ayN∂x+1N2(-a2y2+N2e2ax/N)∂y.
Furthermore, specific cases of L and m give rise to the symmetries X1 to X5 from above, and some additional symmetries.Case 1 (L=1, m=1/3).
The additional symmetries are,(3.5)X6=-2z∂z+t∂t,X7=-(4z2t4+β24t4)∂z+zt∂t.Case 2 (L=1, m=-1).
The additional symmetries are,(3.6)X8=-z∂z-t2∂t+u∂u.Case 3 (L=1, m=1).
The additional symmetries are,(3.7)X9=-t∂t+u∂u.Table1 contains the conserved forms corresponding to each Noether symmetry Xi(i=1,…,9), that is, it lists the three form ω. The four form is Dω, and Dω vanishes on the solutions of (3.2). Thus
(3.8)ω=-Ttdx∧dy∧dz+Tzdx∧dy∧dt-Tydx∧dz∧dt+Txdy∧dz∧dt,
so that
(3.9)DtTt+DxTx+DyTy+DzTz=0
in (3.2).Table 1
Conserved Formω.
T
t
=
e
-
a
x
/
N
t
1
+
2
L
(
-
u
y
u
t
+
u
u
t
y
)
2
β
,
T
x
=
1
2
e
-
a
x
/
N
t
β
(
u
y
u
x
-
u
u
x
y
),
X
1
T
y
=
e
-
a
x
/
N
t
-
2
L
/
m
(
e
2
a
x
/
N
N
t
1
+
(
2
L
/
m
)
β
2
u
y
2
+
u
(
N
t
1
+
2
L
β
2
u
z
z
-
t
2
L
/
m
(
a
t
β
2
u
x
+
N
(
-
t
β
2
u
x
x
+
t
2
L
(
(
1
+
2
L
)
u
t
+
t
u
t
t
)
)
)
)
)
2
N
β,
T
z
=
1
2
e
-
a
x
/
N
t
1
+
(
2
L
(
-
1
+
m
)
/
m
)
β
(
u
z
u
y
-
u
u
y
z
)
T
t
=
e
-
a
x
/
N
t
1
+
2
L
(
-
u
z
u
t
+
u
u
t
z
)
2
β
,
T
x
=
1
2
e
-
a
x
/
N
t
β
(
u
z
u
x
-
u
u
x
z
),
X
2
T
y
=
1
2
e
a
x
/
N
t
β
(
u
z
u
y
-
u
u
y
z
),
T
z
=
e
-
a
x
/
N
t
-
2
L
/
m
(
N
t
1
+
2
L
β
2
u
z
2
+
t
2
L
/
m
u
(
e
2
a
x
/
N
N
t
β
2
u
y
y
-
a
t
β
2
u
x
+
N
(
t
β
2
u
x
x
-
t
2
L
(
(
1
+
2
L
)
u
t
+
t
u
t
t
)
)
)
)
2
N
β
X
3
T
t
=
-
e
-
a
x
/
N
t
1
+
2
L
u
t
β
,
T
x
=
e
-
a
x
/
N
t
β
u
x,
T
y
=
e
a
x
/
N
t
β
u
y
,
T
z
=
e
-
a
x
/
N
t
1
+
2
L
-
(
2
L
/
m
)
β
u
z
T
t
=
-
e
-
a
x
/
N
t
1
+
2
L
(
a
y
u
y
u
t
+
N
u
x
u
t
-
u
(
a
y
u
t
y
+
N
u
t
x
)
)
2
a
β,
T
x
=
1
2
a
β
(
e
-
a
x
/
N
t
-
2
L
/
m
(
t
1
+
(
2
L
/
m
)
β
2
u
x
(
a
y
u
y
+
N
u
x
)
X
4
+u(Nt1+2Lβ2uzz+t2L/m(e2ax/NNtβ2uyy-atβ2ux-atyβ2uxy-Nt2Lut-2LNt2Lut-Nt1+2Lutt)))),
T
y
=
1
2
a
N
β
e
-
a
x
/
N
t
-
2
L
/
m
(
e
2
a
x
/
N
N
t
1
+
(
2
L
/
m
)
β
2
u
y
(
a
y
u
y
+
N
u
x
)
T
y
=
1
2
a
N
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
a
N
t
1
+
2
L
y
β
2
u
z
z
-
t
2
L
/
m
(
a
e
2
a
x
/
N
N
t
β
2
u
y
+
a
2
t
y
β
2
u
x
+
N
(
e
2
a
x
/
N
N
t
β
2
u
x
y
+
a
y
(
-
t
β
2
u
x
x
+
t
2
L
(
(
1
+
2
L
)
u
t
+
t
u
t
t
)
)
)
)
)
),
T
z
=
e
-
a
x
/
N
t
1
+
(
2
L
(
-
1
+
m
)
/
m
)
β
(
u
z
(
a
y
u
y
+
N
u
x
)
-
u
(
a
y
u
y
z
+
N
u
x
z
)
)
2
a
T
t
=
1
2
N
2
β
(
e
-
a
x
/
N
t
1
+
2
L
(
-
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
u
t
+
2
a
N
y
u
x
u
t
+
u
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
t
y
-
2
a
N
y
u
t
x
)
)
),
T
x
=
1
2
N
2
β
e
-
a
x
/
N
t
-
2
L
/
m
(
t
1
+
(
2
L
/
m
)
β
2
u
x
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
-
2
a
N
y
u
x
)
T
x
=
1
2
N
2
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
-
2
a
N
t
1
+
2
L
y
β
2
u
z
z
+
t
2
L
/
m
(
-
2
a
e
2
a
x
/
N
N
t
β
2
u
y
-
2
a
e
2
a
x
/
N
N
t
y
β
2
u
y
y
+
2
a
2
t
y
β
2
u
x
-
e
2
a
x
/
N
N
2
t
β
2
u
x
y
+
a
2
t
y
2
β
2
u
x
y
X
5
T
x
=
1
2
N
2
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
-
2
a
N
t
1
+
2
L
y
β
2
u
z
z
+
t
2
L
/
m
+
2
a
N
t
2
L
y
u
t
+
4
a
L
N
t
2
L
y
u
t
+
2
a
N
t
1
+
2
L
y
u
t
t
)
)
),
T
y
=
1
2
N
3
β
e
-
a
x
/
N
t
-
2
L
/
m
(
e
2
a
x
/
N
N
t
1
+
(
2
L
/
m
)
β
2
u
y
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
-
2
a
N
y
u
x
)
T
y
=
1
2
N
3
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
N
t
1
+
2
L
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
β
2
u
z
z
+
t
2
L
/
m
(
2
a
2
e
2
a
x
/
N
N
t
y
β
2
u
y
+
a
t
(
e
2
a
x
/
N
N
2
+
a
2
y
2
)
β
2
u
x
T
y
=
1
2
N
3
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
N
t
1
+
2
L
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
β
2
u
z
z
+
t
2
L
/
m
+
N
(
2
a
e
2
a
x
/
N
N
t
y
β
2
u
x
y
+
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
(
t
β
2
u
x
x
-
t
2
L
(
(
1
+
2
L
)
u
t
+
t
u
t
t
)
)
)
)
)
),
T
z
=
e
-
a
x
/
N
t
1
+
(
2
L
(
-
1
+
m
)
/
m
)
β
(
u
z
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
-
2
a
N
y
u
x
)
-
u
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
z
-
2
a
N
y
u
x
z
)
)
2
N
2
T
t
=
1
2
N
t
2
β
(
e
-
a
x
/
N
(
N
t
5
u
t
(
2
z
u
z
-
t
u
t
)
+
u
(
N
β
2
u
z
z
+
t
4
(
e
2
a
x
/
N
N
β
2
u
y
y
-
a
β
2
u
x
+
N
β
2
u
x
x
-
2
N
t
u
t
-
2
N
t
z
u
t
z
)
)
)
),
X
6
T
x
=
1
2
e
-
a
x
/
N
t
β
(
-
2
z
u
z
u
x
+
t
u
x
u
t
+
u
(
2
z
u
x
z
-
t
u
t
x
)
),
T
y
=
1
2
e
a
x
/
N
t
β
(
-
2
z
u
z
u
y
+
t
u
y
u
t
+
u
(
2
z
u
y
z
-
t
u
t
y
)
),
T
z
=
1
2
N
t
3
β
(
e
-
a
x
/
N
(
N
β
2
u
z
(
-
2
z
u
z
+
t
u
t
)
+
u
(
2
N
β
2
u
z
-
2
e
2
a
x
/
N
N
t
4
z
β
2
u
y
y
+
2
a
t
4
z
β
2
u
x
-
2
N
t
4
z
β
2
u
x
x
+
6
N
t
5
z
u
t
-
N
t
β
2
u
t
z
+
2
N
t
6
z
u
t
t
)
)
)
T
t
=
1
8
N
t
2
β
(
e
-
a
x
/
N
(
N
t
u
t
(
(
4
t
4
z
2
+
β
2
)
u
z
-
4
t
5
z
u
t
)
+
u
(
4
N
β
2
u
z
+
4
N
z
β
2
u
z
z
+
4
e
2
a
x
/
N
N
t
4
z
β
2
u
y
y
-
4
a
t
4
z
β
2
u
x
T
t
=
1
8
N
t
2
β
(
e
-
a
x
/
N
(
N
t
u
t
(
(
4
t
4
z
2
+
β
2
)
u
z
-
4
t
5
z
u
t
)
+
u
+
4
N
t
4
z
β
2
u
x
x
-
8
N
t
5
z
u
t
-
4
N
t
5
z
2
u
t
z
-
N
t
β
2
u
t
z
)
)
),
X
7
T
x
=
e
-
a
x
/
N
β
(
-
(
4
t
4
z
2
+
β
2
)
u
z
u
x
+
4
t
5
z
u
x
u
t
+
(
(
4
t
4
z
2
+
β
2
)
u
x
z
-
4
t
5
z
u
t
x
)
)
8
t
3,
T
y
=
e
a
x
/
N
β
(
-
(
4
t
4
z
2
+
β
2
)
u
z
u
y
+
4
t
5
z
u
y
u
t
+
u
(
(
4
t
4
z
2
+
β
2
)
u
y
z
-
4
t
5
z
u
t
y
)
)
8
t
3,
T
z
=
1
8
N
t
7
β
e
-
a
x
/
N
(
-
N
β
2
u
z
(
(
4
t
4
z
2
+
β
2
)
u
z
-
4
t
5
z
u
t
)
+
t
4
u
(
8
N
z
β
2
u
z
-
e
2
a
x
/
N
N
β
2
(
4
t
4
z
2
+
β
2
)
u
y
y
+
4
a
t
4
z
2
β
2
u
x
+
a
β
4
u
x
-
4
N
t
4
z
2
β
2
u
x
x
-
N
β
4
u
x
x
T
z
=
1
8
N
t
7
β
e
-
a
x
N
(
-
N
β
2
u
z
(
(
4
t
4
z
2
+
β
2
)
u
z
-
4
t
5
z
u
t
)
+
t
4
u
+
12
N
t
5
z
2
u
t
-
N
t
β
2
u
t
-
4
N
t
z
β
2
u
t
z
+
4
N
t
6
z
2
u
t
t
+
N
t
2
β
2
u
t
t))
T
t
=
1
4
N
β
(
e
-
a
x
/
N
t
2
(
-
N
t
u
t
(
2
z
u
z
+
t
u
t
)
+
u
(
N
t
4
β
2
u
z
z
+
e
2
a
x
/
N
N
β
2
u
y
y
-
a
β
2
u
x
+
N
β
2
u
x
x
-
2
N
t
u
t
+
2
N
t
z
u
t
z
)
)
),
X
8
T
x
=
1
4
e
-
a
x
/
N
t
β
(
2
z
u
z
u
x
+
t
u
x
u
t
-
u
(
2
z
u
x
z
+
t
u
t
x
)
),
T
y
=
1
4
e
a
x
/
N
t
β
(
2
z
u
z
u
y
+
t
u
y
u
t
-
u
(
2
z
u
y
z
+
t
u
t
y
)
),
T
z
=
1
4
N
β
(
e
-
a
x
/
N
t
(
N
t
4
β
2
u
z
(
2
z
u
z
+
t
u
t
)
-
u
(
2
N
t
4
β
2
u
z
-
2
e
2
a
x
/
N
n
z
β
2
u
y
y
+
2
a
z
β
2
u
x
-
2
N
z
β
2
u
x
x
+
6
N
t
z
u
t
+
N
t
5
β
2
u
t
z
+
2
N
t
2
z
u
t
t
)
)
)
T
t
=
1
2
N
β
(
e
-
a
x
/
N
t
2
(
-
N
t
2
u
t
2
+
u
(
N
β
2
u
z
z
+
e
2
a
x
/
N
N
β
2
u
y
y
-
a
β
2
u
x
+
N
β
2
u
x
x
-
2
N
t
u
t
)
)
),
X
9
T
x
=
1
2
e
-
a
x
/
N
t
2
β
(
u
x
u
t
-
u
u
t
x
)
,
T
y
=
1
2
e
a
x
/
N
t
2
β
(
u
y
u
t
-
u
u
t
y
),
T
z
=
1
2
e
-
a
x
/
N
t
2
β
(
u
z
u
t
-
u
u
t
z
)


### 3.2. Symmetry Reduction and Invariant Solutions

We briefly show how the order of the (1+3) wave equation (3.2) can be reduced. Ultimately the equation with four independent variables is reduced to an ordinary differential equation.


#### 3.2.1. Reduction—Using the Principle Noether Algebra

We begin reducing (3.2) using X2 followed by X4. The characteristic equations are
(3.10)adxN=dt0=dyy=du0.
Integrating yields s=ye-(a/N)x and (3.2) is reduced to
(3.11)-1βt2Lutt-1β(2L+1)t2Lut+βt(1+a2N2s2)uss+2sa2N2βtus=0,
with u=u(s,t).A Lagrangian of (3.11) is
(3.12)L=-1βt2L+1ut22+βt(1+a2N2s2)us22.
It turns out that we if we let L=1, N=1 in (3.11), we can obtain its Noether symmetries, namely,
(3.13)X1*=t∂t-u∂u,X2*=∂u.
We reduce (3.11) with X1*, and the characteristic equations are
(3.14)dtt=ds0=du-u.
Integrating yields Y=tu and (3.11) is reduced to the ordinary differential equation:
(3.15)1βY+β(1+s2a2)Y′′+2a2βsY′=0,
with Y=Y(s), and which has a solution in terms of special functions, that is,
(3.16)Y(s)=C1LegendreP[-aβ+-4+a2β22aβ,ias]+C2LegendreQ[-aβ+-4+a2β22aβ,ias],
where C1, C2 are arbitrary constants, LegendreP[n,x] refers to the Legendre polynomial Pn(x), and Q[n,z] refers to the Legendre function of the second kind Qn(z).


#### 3.2.2. Reduction—Case1, L=1, m=1/3

We reduce (3.2) using X1 followed by X6 from Case 1. The characteristic equations are
(3.17)dz-2z=dtt=dx0=du0.
Integrating yields r=t2z and (3.2) is reduced to
(3.18)-8βre-axNur-aβNe-axNux+βe-axNuxx+e-axN(β2-4r2β)urr=0,
with u=u(r,x).A Lagrangian of (3.18) is
(3.19)L=βe-ax/Nux22+e-ax/N(β2-4r2β)ur22.
Hence, we obtain the Noether symmetries of (3.18), namely,
(3.20)X3*=∂x+a2Nu∂u,X4*=∂u.
We reduce (3.18) with X3*, and the characteristic equations are
(3.21)dx1=dr0=2Nduau.
Integrating yields Z=e-ax/2Nu and (3.18) is reduced to the ordinary differential equation:
(3.22)-8βrZ′+(β2-4r2β)Z′′-a2β4N2Z=0,
with Z=Z(r), and which has a solution in terms of special functions, that is,
(3.23)Z(r)=C1LegendreP[-2N+4N2-a2β24N,2rβ]+C2LegendreQ[-2N+4N2-a2β24N,2rβ],
where, as before, C1, C2 are arbitrary constants, LegendreP[n,x] refers to the Legendre polynomial Pn(x), and Q[n,z] refers to the Legendre function of the second kind Qn(z).


## 3.1. The Strict Noether Symmetries of (3.2)

We classify the cases that yield strict Noether symmetries (gauge is zero) of (3.2), via the Lagrangian
(3.3)L=-β2te-(a/N)xux2-β2te(a/N)xuy2-β2t2L+1-(2L/m)e-(a/N)xuz2+12βt2L+1e-(a/N)xut2.
Many of the calculations have been left out as they are tedious—the details are available to the reader in a number of texts that have been cited here.The principle Noether algebra is(3.4)X1=∂y,X2=∂z,X3=∂u,X4=Na∂x+y∂y,X5=-2ayN∂x+1N2(-a2y2+N2e2ax/N)∂y.
Furthermore, specific cases of L and m give rise to the symmetries X1 to X5 from above, and some additional symmetries.Case 1 (L=1, m=1/3).
The additional symmetries are,(3.5)X6=-2z∂z+t∂t,X7=-(4z2t4+β24t4)∂z+zt∂t.Case 2 (L=1, m=-1).
The additional symmetries are,(3.6)X8=-z∂z-t2∂t+u∂u.Case 3 (L=1, m=1).
The additional symmetries are,(3.7)X9=-t∂t+u∂u.Table1 contains the conserved forms corresponding to each Noether symmetry Xi(i=1,…,9), that is, it lists the three form ω. The four form is Dω, and Dω vanishes on the solutions of (3.2). Thus
(3.8)ω=-Ttdx∧dy∧dz+Tzdx∧dy∧dt-Tydx∧dz∧dt+Txdy∧dz∧dt,
so that
(3.9)DtTt+DxTx+DyTy+DzTz=0
in (3.2).Table 1
Conserved Formω.
T
t
=
e
-
a
x
/
N
t
1
+
2
L
(
-
u
y
u
t
+
u
u
t
y
)
2
β
,
T
x
=
1
2
e
-
a
x
/
N
t
β
(
u
y
u
x
-
u
u
x
y
),
X
1
T
y
=
e
-
a
x
/
N
t
-
2
L
/
m
(
e
2
a
x
/
N
N
t
1
+
(
2
L
/
m
)
β
2
u
y
2
+
u
(
N
t
1
+
2
L
β
2
u
z
z
-
t
2
L
/
m
(
a
t
β
2
u
x
+
N
(
-
t
β
2
u
x
x
+
t
2
L
(
(
1
+
2
L
)
u
t
+
t
u
t
t
)
)
)
)
)
2
N
β,
T
z
=
1
2
e
-
a
x
/
N
t
1
+
(
2
L
(
-
1
+
m
)
/
m
)
β
(
u
z
u
y
-
u
u
y
z
)
T
t
=
e
-
a
x
/
N
t
1
+
2
L
(
-
u
z
u
t
+
u
u
t
z
)
2
β
,
T
x
=
1
2
e
-
a
x
/
N
t
β
(
u
z
u
x
-
u
u
x
z
),
X
2
T
y
=
1
2
e
a
x
/
N
t
β
(
u
z
u
y
-
u
u
y
z
),
T
z
=
e
-
a
x
/
N
t
-
2
L
/
m
(
N
t
1
+
2
L
β
2
u
z
2
+
t
2
L
/
m
u
(
e
2
a
x
/
N
N
t
β
2
u
y
y
-
a
t
β
2
u
x
+
N
(
t
β
2
u
x
x
-
t
2
L
(
(
1
+
2
L
)
u
t
+
t
u
t
t
)
)
)
)
2
N
β
X
3
T
t
=
-
e
-
a
x
/
N
t
1
+
2
L
u
t
β
,
T
x
=
e
-
a
x
/
N
t
β
u
x,
T
y
=
e
a
x
/
N
t
β
u
y
,
T
z
=
e
-
a
x
/
N
t
1
+
2
L
-
(
2
L
/
m
)
β
u
z
T
t
=
-
e
-
a
x
/
N
t
1
+
2
L
(
a
y
u
y
u
t
+
N
u
x
u
t
-
u
(
a
y
u
t
y
+
N
u
t
x
)
)
2
a
β,
T
x
=
1
2
a
β
(
e
-
a
x
/
N
t
-
2
L
/
m
(
t
1
+
(
2
L
/
m
)
β
2
u
x
(
a
y
u
y
+
N
u
x
)
X
4
+u(Nt1+2Lβ2uzz+t2L/m(e2ax/NNtβ2uyy-atβ2ux-atyβ2uxy-Nt2Lut-2LNt2Lut-Nt1+2Lutt)))),
T
y
=
1
2
a
N
β
e
-
a
x
/
N
t
-
2
L
/
m
(
e
2
a
x
/
N
N
t
1
+
(
2
L
/
m
)
β
2
u
y
(
a
y
u
y
+
N
u
x
)
T
y
=
1
2
a
N
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
a
N
t
1
+
2
L
y
β
2
u
z
z
-
t
2
L
/
m
(
a
e
2
a
x
/
N
N
t
β
2
u
y
+
a
2
t
y
β
2
u
x
+
N
(
e
2
a
x
/
N
N
t
β
2
u
x
y
+
a
y
(
-
t
β
2
u
x
x
+
t
2
L
(
(
1
+
2
L
)
u
t
+
t
u
t
t
)
)
)
)
)
),
T
z
=
e
-
a
x
/
N
t
1
+
(
2
L
(
-
1
+
m
)
/
m
)
β
(
u
z
(
a
y
u
y
+
N
u
x
)
-
u
(
a
y
u
y
z
+
N
u
x
z
)
)
2
a
T
t
=
1
2
N
2
β
(
e
-
a
x
/
N
t
1
+
2
L
(
-
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
u
t
+
2
a
N
y
u
x
u
t
+
u
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
t
y
-
2
a
N
y
u
t
x
)
)
),
T
x
=
1
2
N
2
β
e
-
a
x
/
N
t
-
2
L
/
m
(
t
1
+
(
2
L
/
m
)
β
2
u
x
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
-
2
a
N
y
u
x
)
T
x
=
1
2
N
2
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
-
2
a
N
t
1
+
2
L
y
β
2
u
z
z
+
t
2
L
/
m
(
-
2
a
e
2
a
x
/
N
N
t
β
2
u
y
-
2
a
e
2
a
x
/
N
N
t
y
β
2
u
y
y
+
2
a
2
t
y
β
2
u
x
-
e
2
a
x
/
N
N
2
t
β
2
u
x
y
+
a
2
t
y
2
β
2
u
x
y
X
5
T
x
=
1
2
N
2
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
-
2
a
N
t
1
+
2
L
y
β
2
u
z
z
+
t
2
L
/
m
+
2
a
N
t
2
L
y
u
t
+
4
a
L
N
t
2
L
y
u
t
+
2
a
N
t
1
+
2
L
y
u
t
t
)
)
),
T
y
=
1
2
N
3
β
e
-
a
x
/
N
t
-
2
L
/
m
(
e
2
a
x
/
N
N
t
1
+
(
2
L
/
m
)
β
2
u
y
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
-
2
a
N
y
u
x
)
T
y
=
1
2
N
3
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
N
t
1
+
2
L
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
β
2
u
z
z
+
t
2
L
/
m
(
2
a
2
e
2
a
x
/
N
N
t
y
β
2
u
y
+
a
t
(
e
2
a
x
/
N
N
2
+
a
2
y
2
)
β
2
u
x
T
y
=
1
2
N
3
β
e
-
a
x
/
N
t
-
2
L
/
m
+
u
(
N
t
1
+
2
L
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
β
2
u
z
z
+
t
2
L
/
m
+
N
(
2
a
e
2
a
x
/
N
N
t
y
β
2
u
x
y
+
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
(
t
β
2
u
x
x
-
t
2
L
(
(
1
+
2
L
)
u
t
+
t
u
t
t
)
)
)
)
)
),
T
z
=
e
-
a
x
/
N
t
1
+
(
2
L
(
-
1
+
m
)
/
m
)
β
(
u
z
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
-
2
a
N
y
u
x
)
-
u
(
(
e
2
a
x
/
N
N
2
-
a
2
y
2
)
u
y
z
-
2
a
N
y
u
x
z
)
)
2
N
2
T
t
=
1
2
N
t
2
β
(
e
-
a
x
/
N
(
N
t
5
u
t
(
2
z
u
z
-
t
u
t
)
+
u
(
N
β
2
u
z
z
+
t
4
(
e
2
a
x
/
N
N
β
2
u
y
y
-
a
β
2
u
x
+
N
β
2
u
x
x
-
2
N
t
u
t
-
2
N
t
z
u
t
z
)
)
)
),
X
6
T
x
=
1
2
e
-
a
x
/
N
t
β
(
-
2
z
u
z
u
x
+
t
u
x
u
t
+
u
(
2
z
u
x
z
-
t
u
t
x
)
),
T
y
=
1
2
e
a
x
/
N
t
β
(
-
2
z
u
z
u
y
+
t
u
y
u
t
+
u
(
2
z
u
y
z
-
t
u
t
y
)
),
T
z
=
1
2
N
t
3
β
(
e
-
a
x
/
N
(
N
β
2
u
z
(
-
2
z
u
z
+
t
u
t
)
+
u
(
2
N
β
2
u
z
-
2
e
2
a
x
/
N
N
t
4
z
β
2
u
y
y
+
2
a
t
4
z
β
2
u
x
-
2
N
t
4
z
β
2
u
x
x
+
6
N
t
5
z
u
t
-
N
t
β
2
u
t
z
+
2
N
t
6
z
u
t
t
)
)
)
T
t
=
1
8
N
t
2
β
(
e
-
a
x
/
N
(
N
t
u
t
(
(
4
t
4
z
2
+
β
2
)
u
z
-
4
t
5
z
u
t
)
+
u
(
4
N
β
2
u
z
+
4
N
z
β
2
u
z
z
+
4
e
2
a
x
/
N
N
t
4
z
β
2
u
y
y
-
4
a
t
4
z
β
2
u
x
T
t
=
1
8
N
t
2
β
(
e
-
a
x
/
N
(
N
t
u
t
(
(
4
t
4
z
2
+
β
2
)
u
z
-
4
t
5
z
u
t
)
+
u
+
4
N
t
4
z
β
2
u
x
x
-
8
N
t
5
z
u
t
-
4
N
t
5
z
2
u
t
z
-
N
t
β
2
u
t
z
)
)
),
X
7
T
x
=
e
-
a
x
/
N
β
(
-
(
4
t
4
z
2
+
β
2
)
u
z
u
x
+
4
t
5
z
u
x
u
t
+
(
(
4
t
4
z
2
+
β
2
)
u
x
z
-
4
t
5
z
u
t
x
)
)
8
t
3,
T
y
=
e
a
x
/
N
β
(
-
(
4
t
4
z
2
+
β
2
)
u
z
u
y
+
4
t
5
z
u
y
u
t
+
u
(
(
4
t
4
z
2
+
β
2
)
u
y
z
-
4
t
5
z
u
t
y
)
)
8
t
3,
T
z
=
1
8
N
t
7
β
e
-
a
x
/
N
(
-
N
β
2
u
z
(
(
4
t
4
z
2
+
β
2
)
u
z
-
4
t
5
z
u
t
)
+
t
4
u
(
8
N
z
β
2
u
z
-
e
2
a
x
/
N
N
β
2
(
4
t
4
z
2
+
β
2
)
u
y
y
+
4
a
t
4
z
2
β
2
u
x
+
a
β
4
u
x
-
4
N
t
4
z
2
β
2
u
x
x
-
N
β
4
u
x
x
T
z
=
1
8
N
t
7
β
e
-
a
x
N
(
-
N
β
2
u
z
(
(
4
t
4
z
2
+
β
2
)
u
z
-
4
t
5
z
u
t
)
+
t
4
u
+
12
N
t
5
z
2
u
t
-
N
t
β
2
u
t
-
4
N
t
z
β
2
u
t
z
+
4
N
t
6
z
2
u
t
t
+
N
t
2
β
2
u
t
t))
T
t
=
1
4
N
β
(
e
-
a
x
/
N
t
2
(
-
N
t
u
t
(
2
z
u
z
+
t
u
t
)
+
u
(
N
t
4
β
2
u
z
z
+
e
2
a
x
/
N
N
β
2
u
y
y
-
a
β
2
u
x
+
N
β
2
u
x
x
-
2
N
t
u
t
+
2
N
t
z
u
t
z
)
)
),
X
8
T
x
=
1
4
e
-
a
x
/
N
t
β
(
2
z
u
z
u
x
+
t
u
x
u
t
-
u
(
2
z
u
x
z
+
t
u
t
x
)
),
T
y
=
1
4
e
a
x
/
N
t
β
(
2
z
u
z
u
y
+
t
u
y
u
t
-
u
(
2
z
u
y
z
+
t
u
t
y
)
),
T
z
=
1
4
N
β
(
e
-
a
x
/
N
t
(
N
t
4
β
2
u
z
(
2
z
u
z
+
t
u
t
)
-
u
(
2
N
t
4
β
2
u
z
-
2
e
2
a
x
/
N
n
z
β
2
u
y
y
+
2
a
z
β
2
u
x
-
2
N
z
β
2
u
x
x
+
6
N
t
z
u
t
+
N
t
5
β
2
u
t
z
+
2
N
t
2
z
u
t
t
)
)
)
T
t
=
1
2
N
β
(
e
-
a
x
/
N
t
2
(
-
N
t
2
u
t
2
+
u
(
N
β
2
u
z
z
+
e
2
a
x
/
N
N
β
2
u
y
y
-
a
β
2
u
x
+
N
β
2
u
x
x
-
2
N
t
u
t
)
)
),
X
9
T
x
=
1
2
e
-
a
x
/
N
t
2
β
(
u
x
u
t
-
u
u
t
x
)
,
T
y
=
1
2
e
a
x
/
N
t
2
β
(
u
y
u
t
-
u
u
t
y
),
T
z
=
1
2
e
-
a
x
/
N
t
2
β
(
u
z
u
t
-
u
u
t
z
)


## 3.2. Symmetry Reduction and Invariant Solutions

We briefly show how the order of the (1+3) wave equation (3.2) can be reduced. Ultimately the equation with four independent variables is reduced to an ordinary differential equation.


### 3.2.1. Reduction—Using the Principle Noether Algebra

We begin reducing (3.2) using X2 followed by X4. The characteristic equations are
(3.10)adxN=dt0=dyy=du0.
Integrating yields s=ye-(a/N)x and (3.2) is reduced to
(3.11)-1βt2Lutt-1β(2L+1)t2Lut+βt(1+a2N2s2)uss+2sa2N2βtus=0,
with u=u(s,t).A Lagrangian of (3.11) is
(3.12)L=-1βt2L+1ut22+βt(1+a2N2s2)us22.
It turns out that we if we let L=1, N=1 in (3.11), we can obtain its Noether symmetries, namely,
(3.13)X1*=t∂t-u∂u,X2*=∂u.
We reduce (3.11) with X1*, and the characteristic equations are
(3.14)dtt=ds0=du-u.
Integrating yields Y=tu and (3.11) is reduced to the ordinary differential equation:
(3.15)1βY+β(1+s2a2)Y′′+2a2βsY′=0,
with Y=Y(s), and which has a solution in terms of special functions, that is,
(3.16)Y(s)=C1LegendreP[-aβ+-4+a2β22aβ,ias]+C2LegendreQ[-aβ+-4+a2β22aβ,ias],
where C1, C2 are arbitrary constants, LegendreP[n,x] refers to the Legendre polynomial Pn(x), and Q[n,z] refers to the Legendre function of the second kind Qn(z).


### 3.2.2. Reduction—Case1, L=1, m=1/3

We reduce (3.2) using X1 followed by X6 from Case 1. The characteristic equations are
(3.17)dz-2z=dtt=dx0=du0.
Integrating yields r=t2z and (3.2) is reduced to
(3.18)-8βre-axNur-aβNe-axNux+βe-axNuxx+e-axN(β2-4r2β)urr=0,
with u=u(r,x).A Lagrangian of (3.18) is
(3.19)L=βe-ax/Nux22+e-ax/N(β2-4r2β)ur22.
Hence, we obtain the Noether symmetries of (3.18), namely,
(3.20)X3*=∂x+a2Nu∂u,X4*=∂u.
We reduce (3.18) with X3*, and the characteristic equations are
(3.21)dx1=dr0=2Nduau.
Integrating yields Z=e-ax/2Nu and (3.18) is reduced to the ordinary differential equation:
(3.22)-8βrZ′+(β2-4r2β)Z′′-a2β4N2Z=0,
with Z=Z(r), and which has a solution in terms of special functions, that is,
(3.23)Z(r)=C1LegendreP[-2N+4N2-a2β24N,2rβ]+C2LegendreQ[-2N+4N2-a2β24N,2rβ],
where, as before, C1, C2 are arbitrary constants, LegendreP[n,x] refers to the Legendre polynomial Pn(x), and Q[n,z] refers to the Legendre function of the second kind Qn(z).


## 3.2.1. Reduction—Using the Principle Noether Algebra

We begin reducing (3.2) using X2 followed by X4. The characteristic equations are
(3.10)adxN=dt0=dyy=du0.
Integrating yields s=ye-(a/N)x and (3.2) is reduced to
(3.11)-1βt2Lutt-1β(2L+1)t2Lut+βt(1+a2N2s2)uss+2sa2N2βtus=0,
with u=u(s,t).A Lagrangian of (3.11) is
(3.12)L=-1βt2L+1ut22+βt(1+a2N2s2)us22.
It turns out that we if we let L=1, N=1 in (3.11), we can obtain its Noether symmetries, namely,
(3.13)X1*=t∂t-u∂u,X2*=∂u.
We reduce (3.11) with X1*, and the characteristic equations are
(3.14)dtt=ds0=du-u.
Integrating yields Y=tu and (3.11) is reduced to the ordinary differential equation:
(3.15)1βY+β(1+s2a2)Y′′+2a2βsY′=0,
with Y=Y(s), and which has a solution in terms of special functions, that is,
(3.16)Y(s)=C1LegendreP[-aβ+-4+a2β22aβ,ias]+C2LegendreQ[-aβ+-4+a2β22aβ,ias],
where C1, C2 are arbitrary constants, LegendreP[n,x] refers to the Legendre polynomial Pn(x), and Q[n,z] refers to the Legendre function of the second kind Qn(z).


## 3.2.2. Reduction—Case1, L=1, m=1/3

We reduce (3.2) using X1 followed by X6 from Case 1. The characteristic equations are
(3.17)dz-2z=dtt=dx0=du0.
Integrating yields r=t2z and (3.2) is reduced to
(3.18)-8βre-axNur-aβNe-axNux+βe-axNuxx+e-axN(β2-4r2β)urr=0,
with u=u(r,x).A Lagrangian of (3.18) is
(3.19)L=βe-ax/Nux22+e-ax/N(β2-4r2β)ur22.
Hence, we obtain the Noether symmetries of (3.18), namely,
(3.20)X3*=∂x+a2Nu∂u,X4*=∂u.
We reduce (3.18) with X3*, and the characteristic equations are
(3.21)dx1=dr0=2Nduau.
Integrating yields Z=e-ax/2Nu and (3.18) is reduced to the ordinary differential equation:
(3.22)-8βrZ′+(β2-4r2β)Z′′-a2β4N2Z=0,
with Z=Z(r), and which has a solution in terms of special functions, that is,
(3.23)Z(r)=C1LegendreP[-2N+4N2-a2β24N,2rβ]+C2LegendreQ[-2N+4N2-a2β24N,2rβ],
where, as before, C1, C2 are arbitrary constants, LegendreP[n,x] refers to the Legendre polynomial Pn(x), and Q[n,z] refers to the Legendre function of the second kind Qn(z).


## 4. Conclusion

We classified the Noether symmetry generators, determined some conserved forms, and reduced some cases of the underlying equations associated with the wave equation on the Bianchi III manifold. The first reduction done above involved the principle Noether algebra, whilst the second dealt with a particular case. To obtain other reductions, one needs to conclude a three-dimensional subalgebra of symmetries to reduce to an ordinary differential equation whose solution would be an invariant solution invariant under the subalgebra. Alternatively, a lower dimensional subalgebra can be used to reduce to a partial differential equation which may be tackled using other methods. The final solution in this case will be invariant only under the lower dimensional algebra. In general, the procedure performed above is the most convenient.



---

*Source: 765361-2012-12-04.xml*