Supersymmetric Version of the Euler System and Its Invariant Solutions

Abstract

In this paper, we formulate a supersymmetric extension of the Euler system of equations. We compute a superalgebra of Lie symmetries of the supersymmetric system. Next, we classify the one-dimensional subalgebras of this superalgebra into 49 equivalence conjugation classes. For some of the subalgebras, the invariants have a non-standard structure. For nine selected subalgebras, we use the symmetry reduction method to find invariants, orbits and reduced systems. Through the solutions of these reduced systems, we obtain solutions of the supersymmetric Euler system. The obtained solutions include bumps, kinks, multiple wave solutions and solutions expressed in terms of arbitrary functions.

1. Introduction

Supersymmetric generalizations of fluid dynamical equations have recently generated much interest among mathematicians and physicists. Such extensions have been constructed recently for certain hydrodynamic-type models (see, e.g., [1,2,3,4]). A study of polytropic supersymmetric models has also been performed by Das and Popowicz [5] in the specific case when $p=k{\rho }^k$. Euler-type systems have also been analyzed from the Grassmann variable point of view by Fatyga, Kostelecky and Truax [6]. Certain solutions were found in terms of infinite series. In addition, supersymmetric versions of the Chaplygin gas in $(1+1)$ and $(2+1)$ dimensions, derived from the action for a superstring and a Nambu-Goto supermembrane, respectively, were proposed by Jackiw et al. (see [7] and references therein). It was suggested that non-Abelian fluid mechanics may describe a quark-gluon plasma [7]. In this context, we propose to investigate the Euler equations for a one-dimensional compressible fluid flow.

The purpose of this paper is to formulate a supersymmetric version of the Euler system of equations describing the nonstationary one-dimensional ideal (non-viscous) compressible fluid flow and construct their invariant solutions. To our knowledge, a supersymmetric extension of the Euler system has not been formulated before. The symmetry group transformations of the classical version of the Euler system of equations:

$$\begin{array}{cc}& \rho u_t+\rho u{u}_x+p_x=0\hfill \\ & {\rho }_t+u{\rho }_x+\rho u_x=0\hfill \\ & p_t+u{p}_x+kp{u}_x=0\hfill \end{array}$$

(1)

has been investigated by Ovsiannikov [8], where u is the flow velocity, p and ρ are the pressure and density of the fluid flow, respectively, and k is the polytropic exponent. Under the above assumptions, the examined system of Equations (1) forms a quasilinear hyperbolic system. It was proven [8,9] that the basic system (1) is invariant under the Galilean similitude group, sim$\left(2\right)$ (when $k\ne 3$), and under the Galilean-projective group, simPr$\left(2\right)$ (i.e., the Galilean-similitude group extended by the projective transformation when $k=3$). The Lie algebra of infinitesimal symmetries of fluid dynamics Equations (1) is spanned by:

$$\begin{array}{cc}& P_0={\partial }_t,P_1={\partial }_x,B=t{\partial }_x+{\partial }_u,D=t{\partial }_t+x{\partial }_x,\hfill \\ & F=-t{\partial }_t+u{\partial }_u-2\rho {\partial }_{\rho },G=\rho {\partial }_{\rho }+p{\partial }_p,\hfill \\ & C=t^2{\partial }_t+tx{\partial }_x+(x-tu){\partial }_u-t\rho {\partial }_{\rho }-3tp{\partial }_p,whenk=3\hfill \end{array}$$

(2)

The infinitesimal operators, $P_0$ and $P_1$, are translation operators in the directions of t and x, respectively. The operator, B, corresponds to the Galilean transformation, and operators D, F and G are dilation operators. Operator C corresponds to the projective conformal transformation, when the specific value of the polytropic exponent is $k=3$. The Lie algebra, $\mathfrak{L}$, generated by (2) can be decomposed in the form:

$$\mathfrak{L}=\{C,P_0,F,D,G\}.15ex+\supset \{P_1,B\}$$

(3)

In the work of one of the authors of this paper [10], all subalgebras of the Lie algebra, $\mathfrak{L}$, were classified, and the symmetry reduction method (SRM) was used to generate invariant solutions of the Euler system (1), including, among other solutions, travelling waves, bumps, kinks, double-periodic solutions, etc. Many of these solutions had previously been found by other methods, such as the method of characteristics (see, for example, [8,9,11,12] and references therein).

The SRM, on the other hand, certainly has a broader range of application, since it leads to different sets of solutions obtained from the Lie algebra of symmetries (2) and its subalgebras classified into different conjugacy classes. This approach is of great utility in the investigation of the classical Euler Equations (1). The results obtained for the classical system (1) were so promising that it seemed worthwhile to try to apply the method to the case of a supersymmetric extension of a hydrodynamic-type system and construct invariant solutions through the use of a superspace and superfield formalism.

This paper is organized as follows. In Section 2, we construct a supersymmetric extension of the Euler-Equations (1) through a superspace and superfield formalism involving two independent fermionic variables. In Section 3, we discuss in detail a number of symmetries of the proposed supersymmetric system. For this purpose, we use a generalization of the prolongation formalism of vector fields extended in such a way as to encompass both even and odd Grassmann variables. In this case, we adapt the symmetry criterion in order to determine a superalgebra associated with the supersymmetric system. Next, we classify the one-dimensional subalgebras of the superalgebra into conjugation classes under the action of the associated supergroup. The subalgebras are classified in such a way that each representative subalgebra corresponding to a different conjugation class generates a different type of solution [13,14,15,16]. This analysis allows us to perform the symmetry reductions systematically. Section 4 contains the reduced systems corresponding to those subalgebras with its standard invariant structures. A new method is proposed for solving the reduced systems of the equations. Several new invariant solutions of the supersymmetric Euler Equation are constructed, some of which contain the freedom of arbitrary functions of one or two arguments involving expressions in terms of bosonic or fermionic variables. In the case when the arbitrary functions are expressed only in the independent fermionic variables, we are able to obtain the solution based on a truncated Taylor expansion. We give some physical interpretation of the obtained results in the context of non-linear field models. Finally, Section 5 contains some final remarks and future perspectives.

2. Supersymmetric Extension

In order to supersymmetrize the Euler Equation (1), we extend the space of independent variables, $\{x,t\}$, to the superspace, $\{x,t,{\theta }_1,{\theta }_2\}$, where ${\theta }_1$ and ${\theta }_2$ are independent Grassmannian-odd variables. The variables, x and t, represent the bosonic (even Grassmannian) coordinates on two-dimensional Euclidean space. The quantities, ${\theta }_1$ and ${\theta }_2$, are anticommuting fermionic (odd Grassmann) variables. We replace the real bosonic-valued fields of velocity, $u(x,t)$, density, $\rho (x,t)$, and pressure, $p(x,t)$, by the Grassmann-even superfields defined by:

$$\begin{array}{cc}& \mathcal{U}(x,t,{\theta }_1,{\theta }_2)=F(x,t)+{\theta }_1{\varphi }_1(x,t)+{\theta }_2{\varphi }_2(x,t)+{\theta }_1{\theta }_{2}u(x,t)\hfill \\ & \mathcal{R}(x,t,{\theta }_1,{\theta }_2)=G(x,t)+{\theta }_1{\psi }_1(x,t)+{\theta }_2{\psi }_2(x,t)+{\theta }_1{\theta }_2\rho (x,t)\hfill \\ & \mathcal{P}(x,t,{\theta }_1,{\theta }_2)=H(x,t)+{\theta }_1{\omega }_1(x,t)+{\theta }_2{\omega }_2(x,t)+{\theta }_1{\theta }_{2}p(x,t)\hfill \end{array}$$

(4)

respectively. Here, the quantities, ${\varphi }_i$, ${\psi }_i$ and ${\omega }_i$ $i=1,2$, are fermionic fields, while F, G and H are bosonic fields. The supersymmetric extension of the Euler Equations (1) is constructed in such a way that it remains invariant under the supersymmetry transformations:

$$x\to x-\underline{{\eta }_1}{\theta }_1,{\theta }_1\to {\theta }_1+\underline{{\eta }_1}andt\to t-\underline{{\eta }_2}{\theta }_2,{\theta }_2\to {\theta }_2+\underline{{\eta }_2}$$

(5)

where $\underline{{\eta }_1}$ and $\underline{{\eta }_2}$ are arbitrary constant odd parameters. In what follows, we use the convention that underlined constants represent Grassmann-odd parameters. These transformations are generated by the infinitesimal supersymmetry operators:

$$Q_1=-{\theta }_1{\partial }_x+{\partial }_{{\theta }_1}and{Q}_2=-{\theta }_2{\partial }_t+{\partial }_{{\theta }_2}$$

(6)

satisfying the anticommutation relations:

$$\{Q_i,Q_i\}=-2{\partial }_{x^i},i=1,2$$

(7)

where $x^1$ and $x^2$ stand for x and t, respectively. In order to make the extended model invariant under the supersymmetry generators, $Q_1$ and $Q_2$, we introduce the covariant derivative operators:

$$D_1={\theta }_1{\partial }_x+{\partial }_{{\theta }_1}and{D}_2={\theta }_2{\partial }_t+{\partial }_{{\theta }_2}$$

(8)

which possess the property that each of the derivatives, $D_i$, $i=1,2$, anticommutes with both supersymmetry operators, $Q_i$, $i=1,2$:

$$[D_i,D_i]=2{\partial }_{x^i},\{D_1,D_2\}=\{D_i,Q_j\}=0fori,j=1,2$$

(9)

Thus, if we write our supersymmetry equations in terms of the superfields, $\mathcal{U}$, $\mathcal{R}$ and $\mathcal{P}$, and their covariant derivatives of various orders, then they will be invariant under the transformations generated by $Q_1$ and $Q_2$. The most general supersymmetric extension of the Euler Equation (1) invariant under the transformations generated by $Q_1$ and $Q_2$ is constructed by considering linear combinations of the products of the various covariant derivatives of the superfields, $\mathcal{U}$, $\mathcal{R}$ and $\mathcal{P}$. These multiply together to produce the given terms as coefficients, whose components reproduce each term of the classical Equation (1). The result of this analysis provides the following form of the supersymmetric Euler Equations:

$$\begin{array}{cc}& -a_1\mathcal{R}\left(D_1{D}_2{}^3\mathcal{U}\right)+a_2\left(D_1\mathcal{R}\right)\left(D_2{}^3\mathcal{U}\right)-a_3\left(D_2\mathcal{R}\right)\left(D_1{D}_2{}^2\mathcal{U}\right)+\hfill \\ & (a_1+a_2+a_3-1)\left(D_1{D}_2\mathcal{R}\right)\left(D_2{}^2\mathcal{U}\right)+a_4\mathcal{R}\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^3{D}_2\mathcal{U}\right)\hfill \\ & -a_5\left(D_1\mathcal{R}\right)\left(D_2\mathcal{U}\right)\left(D_1{}^3{D}_2\mathcal{U}\right)-a_6\left(D_1\mathcal{R}\right)\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^2{D}_2\mathcal{U}\right)\hfill \\ & +a_7\left(D_2\mathcal{R}\right)\left(D_1\mathcal{U}\right)\left(D_1{}^3{D}_2\mathcal{U}\right)+a_8\left(D_2\mathcal{R}\right)\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^3\mathcal{U}\right)+a_9\left(D_1{D}_2\mathcal{R}\right)\mathcal{U}\left(D_1{}^3{D}_2\mathcal{U}\right)\hfill \\ & -a_{10}\left(D_1{D}_2\mathcal{R}\right)\left(D_1\mathcal{U}\right)\left(D_1{}^2{D}_2\mathcal{U}\right)+a_{11}\left(D_1{D}_2\mathcal{R}\right)\left(D_2\mathcal{U}\right)\left(D_1{}^3\mathcal{U}\right)\hfill \\ & +\left(1-\sum _{i=4}^{11}{a}_i\right)\left(D_1{D}_2\mathcal{R}\right)\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^2\mathcal{U}\right)+D_1{}^2\mathcal{P}=0\hfill \\ & \\ & D_2{}^2\mathcal{R}-a_{12}\mathcal{U}\left(D_1{}^3{D}_2\mathcal{R}\right)+a_{13}\left(D_1\mathcal{U}\right)\left(D_1{}^2{D}_2\mathcal{R}\right)-a_{14}\left(D_2\mathcal{U}\right)\left(D_1{}^3\mathcal{R}\right)\hfill \\ & +(a_{12}+a_{13}+a_{14}-1)\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^2\mathcal{R}\right)-a_{15}\mathcal{R}\left(D_1{}^3{D}_2\mathcal{U}\right)+a_{16}\left(D_1\mathcal{R}\right)\left(D_1{}^2{D}_2\mathcal{U}\right)\hfill \\ & -a_{17}\left(D_2\mathcal{R}\right)\left(D_1{}^3\mathcal{U}\right)+(a_{15}+a_{16}+a_{17}-1)\left(D_1{D}_2\mathcal{R}\right)\left(D_1{}^2\mathcal{U}\right)=0\hfill \\ & \\ \mathrm{and}\hfill \\ & \\ & D_2{}^2\mathcal{P}-a_{18}\mathcal{U}\left(D_1{}^3{D}_2\mathcal{P}\right)+a_{19}\left(D_1\mathcal{U}\right)\left(D_1{}^2{D}_2\mathcal{P}\right)-a_{20}\left(D_2\mathcal{U}\right)\left(D_1{}^3\mathcal{P}\right)\hfill \\ & +(a_{18}+a_{19}+a_{20}-1)\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^2\mathcal{P}\right)-a_{21}k\mathcal{P}\left(D_1{}^3{D}_2\mathcal{U}\right)+a_{22}k\left(D_1\mathcal{P}\right)\left(D_1{}^2{D}_2\mathcal{U}\right)\hfill \\ & -a_{23}k\left(D_2\mathcal{P}\right)\left(D_1{}^3\mathcal{U}\right)+(a_{21}+a_{22}+a_{23}-1)k\left(D_1{D}_2\mathcal{P}\right)\left(D_1{}^2\mathcal{U}\right)=0\hfill \end{array}$$

(10)

where $a_1,a_2,\dots ,a_{23}$ are twenty-three arbitrary bosonic parameters. In this paper, we focus our analysis on the simplest case, where all of the parameters, $a_i$, $i=1,2,\dots ,23$, vanish:

$$\begin{array}{cc}\hfill \left(D_1{D}_2\mathcal{R}\right)\left(D_2{}^2\mathcal{U}\right)-\left(D_1{D}_2\mathcal{R}\right)\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^2\mathcal{U}\right)-D_1{}^2\mathcal{P}& =0\hfill \\ \hfill D_2{}^2\mathcal{R}-\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^2\mathcal{R}\right)-\left(D_1{D}_2\mathcal{R}\right)\left(D_1{}^2\mathcal{U}\right)& =0\hfill \\ \hfill D_2{}^2\mathcal{P}-\left(D_1{D}_2\mathcal{U}\right)\left(D_1{}^2\mathcal{P}\right)-k\left(D_1{D}_2\mathcal{P}\right)\left(D_1{}^2\mathcal{U}\right)& =0\hfill \end{array}$$

(11)

In what follows, we will refer to system (11) as the supersymmetric extension of the Euler Equations, since it reduces to the classical Euler Equation (1) when ${\theta }_1$ and ${\theta }_2$ tend to zero.

In this paper, we use the convention that partial derivatives involving odd variables satisfy the Leibniz rule:

$${\partial }_{{\theta }_i}\left(fg\right)=\left({\partial }_{{\theta }_i}f\right)g+{(-1)}^{\mathrm{deg}\left(f\right)}f\left({\partial }_{{\theta }^i}g\right)$$

(12)

where:

$$\mathrm{deg}\left(f\right)=\left\{\begin{array}{cc}0\hfill & \mathrm{if}f \mathrm{is} \mathrm{even},\hfill \\ 1\hfill & \mathrm{if}ff \mathrm{is} \mathrm{odd},\hfill \end{array}\right.$$

(13)

and the notation:

$$f_{{\theta }_1{\theta }_2}=\frac{\partial }{\partial {\theta }_2}\left(\frac{\partial f}{\partial {\theta }_1}\right)$$

(14)

The partial derivatives with respect to odd coordinates satisfy ${\partial }_{{\theta }_i}{\theta }_j={\delta }_j^i$. The operators, ${\partial }_{{\theta }_i}$, $Q_1$, $Q_2$, $D_1$ and $D_2$, change the parity of the function acted on in the sense that it converts a bosonic function to a fermionic function and vice versa. For example, ${\partial }_{{\theta }_1}\mathcal{U}$ is an odd superfield, while ${\partial }_{{\theta }_1{\theta }_2}\mathcal{U}$ is an even superfield, and so on. For further details, see the book by Cornwell [17] and the reference by DeWitt [18]. The chain rule for a Grassmann-valued composite function, $f\left(g\right(x\left)\right)$, is [18,19]

$$\frac{\partial f}{\partial x}=\frac{\partial g}{\partial x}·\frac{\partial f}{\partial g}$$

(15)

The interchangeability of mixed derivatives (with proper respect to the ordering of odd variables) is assumed throughout.

The even supernumbers, variables, fields, etc., are assumed to be elements of the even part, ${\Lambda }_{even}$, of the underlying abstract real Grassmann ring, Λ. The odd supernumbers, variables, fields, etc., lie in its odd part, ${\Lambda }_{odd}$. We shall assume throughout the paper that the functions, $u(x,t)$, $\rho (x,t)$ and $p(x,t)$, have values in the invertible subset of ${\Lambda }_{even}$ plus $\left\{0\right\}$, i.e., nonvanishing nilpotent values are ruled out. This technical assumption allows us to perform the necessary simplifications in our calculations without splitting off of singular subcases.

Re-writing Equation (11) in terms of the independent variables, x, t, ${\theta }_1$ and ${\theta }_2$, we obtain the form:

$$\begin{array}{cc}& {\theta }_1{\theta }_2{\mathcal{R}}_{xt}{\mathcal{U}}_t+{\theta }_1{\mathcal{R}}_{x{\theta }_2}{\mathcal{U}}_t-{\theta }_2{\mathcal{R}}_{t{\theta }_1}{\mathcal{U}}_t-{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_t+{\theta }_1{\theta }_2{\mathcal{R}}_{xt}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x-{\theta }_1{\theta }_2{\mathcal{R}}_{x{\theta }_2}{\mathcal{U}}_{t{\theta }_1}{\mathcal{U}}_x\hfill \\ & +{\theta }_1{\mathcal{R}}_{x{\theta }_2}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x+{\theta }_1{\theta }_2{\mathcal{R}}_{t{\theta }_1}{\mathcal{U}}_{x{\theta }_2}{\mathcal{U}}_x-{\theta }_2{\mathcal{R}}_{t{\theta }_1}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x+{\theta }_1{\theta }_2{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{xt}{\mathcal{U}}_x\hfill \\ & +{\theta }_1{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{x{\theta }_2}{\mathcal{U}}_x-{\theta }_2{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{t{\theta }_1}{\mathcal{U}}_x-{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x-{\mathcal{P}}_x=0\hfill \end{array}$$

(16)

$$\begin{array}{cc}& {\mathcal{R}}_t-{\theta }_1{\theta }_2{\mathcal{U}}_{xt}{\mathcal{R}}_x-{\theta }_1{\mathcal{U}}_{x{\theta }_2}{\mathcal{R}}_x+{\theta }_2{\mathcal{U}}_{t{\theta }_1}{\mathcal{R}}_x+{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{R}}_x-{\theta }_1{\theta }_2{\mathcal{R}}_{xt}{\mathcal{U}}_x-{\theta }_1{\mathcal{R}}_{x{\theta }_2}{\mathcal{U}}_x\hfill \\ & +{\theta }_2{\mathcal{R}}_{t{\theta }_1}{\mathcal{U}}_x+{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x=0\hfill \end{array}$$

(17)

$$\begin{array}{cc}& {\mathcal{P}}_t-{\theta }_1{\theta }_2{\mathcal{U}}_{xt}{\mathcal{P}}_x-{\theta }_1{\mathcal{U}}_{x{\theta }_2}{\mathcal{P}}_x+{\theta }_2{\mathcal{U}}_{t{\theta }_1}{\mathcal{P}}_x+{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{P}}_x-k{\theta }_1{\theta }_2{\mathcal{P}}_{xt}{\mathcal{U}}_x-k{\theta }_1{\mathcal{P}}_{x{\theta }_2}{\mathcal{U}}_x\hfill \\ & +k{\theta }_2{\mathcal{P}}_{t{\theta }_1}{\mathcal{U}}_x+k{\mathcal{P}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x=0\hfill \end{array}$$

(18)

The system (16)–(18) is a second-order nonlinear system of partial differential Equations, whose coefficients depend on the fermionic variables.

3. Symmetries of the Supersymmetric Euler Equations

A symmetry supergroup, G, of a supersymmetric system is a (local) supergroup of transformations acting on a Cartesian product of supermanifolds, $X\times U$, where X is the space of independent variables, $(x,t,{\theta }_1,{\theta }_2)$, and U is the space of dependent superfields, $(\mathcal{U},\mathcal{R},\mathcal{P})$. The action of the group, G, on the function, $\mathcal{U}(x,t,{\theta }_1,{\theta }_2)$, $\mathcal{R}(x,t,{\theta }_1,{\theta }_2)$ and $\mathcal{P}(x,t,{\theta }_1,{\theta }_2)$, maps solutions of (11) to solutions of (11). If we assume that G is a Lie supergroup, as described in [20,21], one can associate with it its Lie superalgebra, $\mathfrak{G}$, whose elements are infinitesimal symmetries of (11).

The supersymmetric Euler Equation (11) are invariant under the Lie superalgebra, $\mathfrak{S}$, generated by the following seven infinitesimal vector fields:

$$\begin{array}{cc}& L_1=2x{\partial }_x+{\theta }_1{\partial }_{{\theta }_1}+3\mathcal{U}{\partial }_{\mathcal{U}}+4\mathcal{P}{\partial }_{\mathcal{P}},L_2=2t{\partial }_t+{\theta }_2{\partial }_{{\theta }_2}-\mathcal{U}{\partial }_{\mathcal{U}}-4\mathcal{P}{\partial }_{\mathcal{P}},\hfill \\ & L_3=\mathcal{R}{\partial }_{\mathcal{R}}+\mathcal{P}{\partial }_{\mathcal{P}},P_1={\partial }_x,P_2={\partial }_t,Q_1=-{\theta }_1{\partial }_x+{\partial }_{{\theta }_1}\hfill \\ & Q_2={\theta }_2{\partial }_t+{\partial }_{{\theta }_2}\hfill \end{array}$$

(19)

The generators, $P_1$ and $P_2$, represent translations in space and time, respectively, while $L_1$, $L_2$ and $L_3$ generate dilations in both even and odd variables. In addition, we recover the supersymmetry transformations, $Q_1$ and $Q_2$, which we identified previously in (6). In order to determine this superalgebra of infinitesimal symmetries, we have made use of the theory described in the book by Olver [22]. The commutation (anticommutation in the case of two odd operators) relations of the superalgebra, $\mathfrak{S}$, of the supersymmetric Euler Equations are given in Table 1.

Table 1. Supercommutation table for the Lie superalgebra, $\mathfrak{S}$, spanned by the vector fields (19).

	${\mathbf{L}}_1$	${\mathbf{P}}_1$	${\mathbf{Q}}_1$	${\mathbf{L}}_2$	${\mathbf{P}}_2$	${\mathbf{Q}}_2$	${\mathbf{L}}_3$
${\mathbf{L}}_{\mathbf{1}}$	0	$-2P_1$	$-Q_1$	0	0	0	0
${\mathbf{P}}_{\mathbf{1}}$	$2P_1$	0	0	0	0	0	0
${\mathbf{Q}}_{\mathbf{1}}$	$Q_1$	0	$-2P_1$	0	0	0	0
${\mathbf{L}}_{\mathbf{2}}$	0	0	0	0	$-2P_2$	$-Q_2$	0
${\mathbf{P}}_{\mathbf{2}}$	0	0	0	$2P_2$	0	0	0
${\mathbf{Q}}_{\mathbf{2}}$	0	0	0	$Q_2$	0	$-2P_2$	0
${\mathbf{L}}_{\mathbf{3}}$	0	0	0	0	0	0	0

Table 1. Supercommutation table for the Lie superalgebra, $\mathfrak{S}$, spanned by the vector fields (19).

	${\mathbf{L}}_1$	${\mathbf{P}}_1$	${\mathbf{Q}}_1$	${\mathbf{L}}_2$	${\mathbf{P}}_2$	${\mathbf{Q}}_2$	${\mathbf{L}}_3$
${\mathbf{L}}_{\mathbf{1}}$	0	$-2P_1$	$-Q_1$	0	0	0	0
${\mathbf{P}}_{\mathbf{1}}$	$2P_1$	0	0	0	0	0	0
${\mathbf{Q}}_{\mathbf{1}}$	$Q_1$	0	$-2P_1$	0	0	0	0
${\mathbf{L}}_{\mathbf{2}}$	0	0	0	0	$-2P_2$	$-Q_2$	0
${\mathbf{P}}_{\mathbf{2}}$	0	0	0	$2P_2$	0	0	0
${\mathbf{Q}}_{\mathbf{2}}$	0	0	0	$Q_2$	0	$-2P_2$	0
${\mathbf{L}}_{\mathbf{3}}$	0	0	0	0	0	0	0

It should be noted that the generators of the superalgebra, $\mathfrak{S}$, listed in (19), contain three dilations, $L_1$, $L_2$ and $L_3$. This is analogous to the Lie algebra, $\mathfrak{L}$, of the classical version, whose generators (in Equation (2)) also contain three dilations, D, F and G. The decomposition of the superalgebra, $\mathfrak{S}$, of the supersymmetric system is very different from the decomposition of the Lie algebra of the classical system, $\mathfrak{L}$.

The Lie superalgebra, $\mathfrak{S}$, can be decomposed into the following combination of direct and semi-direct sums:

$$\mathfrak{S}=\left(\left\{L_1\right\}.15ex+\supset \{P_1,Q_1\}\right)\oplus \left(\left\{L_2\right\}.15ex+\supset \{P_2,Q_2\}\right)\oplus \left\{L_3\right\}$$

(20)

The one-dimensional subalgebras of $\mathfrak{S}$ can be classified into conjugacy classes under the action of the group generated by $\mathfrak{S}$. We obtain the following list of subalgebras:

$$\begin{array}{cc}& {\mathcal{L}}_1=\left\{P_2\right\},{\mathcal{L}}_2=\left\{\underline{\nu }{Q}_2\right\},{\mathcal{L}}_3=\{P_2+\underline{\nu }{Q}_2\},{\mathcal{L}}_4=\left\{L_2\right\},\hfill \\ & {\mathcal{L}}_5=\left\{P_1\right\},{\mathcal{L}}_6=\left\{\underline{\mu }{Q}_1\right\},{\mathcal{L}}_7=\{P_1+\underline{\mu }{Q}_1\},{\mathcal{L}}_8=\left\{L_1\right\},\hfill \\ & {\mathcal{L}}_9=\{P_1+\epsilon P_2\},{\mathcal{L}}_{10}=\{P_1+\underline{\nu }{Q}_2\},{\mathcal{L}}_{11}=\{P_1+\epsilon P_2+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{12}=\{L_2+\epsilon P_1\},{\mathcal{L}}_{13}=\{P_2+\underline{\mu }{Q}_1\},{\mathcal{L}}_{14}=\{\underline{\mu }{Q}_1+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{15}=\{P_2+\underline{\mu }{Q}_1+\underline{\nu }{Q}_2\},{\mathcal{L}}_{16}=\{L_2+\underline{\mu }{Q}_1\},{\mathcal{L}}_{17}=\{P_1+\epsilon P_2+\underline{\mu }{Q}_1\},\hfill \\ & {\mathcal{L}}_{18}=\{P_1+\underline{\mu }{Q}_1+\underline{\nu }{Q}_2\},{\mathcal{L}}_{19}=\{P_1+\epsilon P_2+\underline{\mu }{Q}_1+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{20}=\{L_2+\epsilon P_1+\underline{\mu }{Q}_1\}{\mathcal{L}}_{21}=\{L_1+\epsilon P_2\},{\mathcal{L}}_{22}=\{L_1+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{23}=\{L_1+\epsilon P_2+\underline{\nu }{Q}_2\},{\mathcal{L}}_{24}=\{L_1+k{L}_2\},{\mathcal{L}}_{25}=\left\{L_3\right\},\hfill \\ & {\mathcal{L}}_{26}=\{L_3+\epsilon P_2\},{\mathcal{L}}_{27}=\{L_3+\underline{\nu }{Q}_2\},{\mathcal{L}}_{28}=\{L_3+\epsilon P_2+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{29}=\{L_3+k{L}_2\}{\mathcal{L}}_{30}=\{L_3+\epsilon P_1\},{\mathcal{L}}_{31}=\{L_3+\underline{\mu }{Q}_1\},\hfill \\ & {\mathcal{L}}_{32}=\{L_3+\epsilon P_1+\underline{\mu }{Q}_1\},{\mathcal{L}}_{33}=\{L_3+k{L}_1\}{\mathcal{L}}_{34}=\{L_3+{\epsilon }_1{P}_1+{\epsilon }_2{P}_2\},\hfill \\ & {\mathcal{L}}_{35}=\{L_3+\epsilon P_1+\underline{\nu }{Q}_2\},{\mathcal{L}}_{36}=\{L_3+{\epsilon }_1{P}_1+{\epsilon }_2{P}_2+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{37}=\{L_3+k{L}_2+\epsilon P_1\},{\mathcal{L}}_{38}=\{L_3+\epsilon P_2+\underline{\mu }{Q}_1\},\hfill \\ & {\mathcal{L}}_{39}=\{L_3+\underline{\mu }{Q}_1+\underline{\nu }{Q}_2\},{\mathcal{L}}_{40}=\{L_3+\epsilon P_2+\underline{\mu }{Q}_1+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{41}=\{L_3+k{L}_2+\underline{\mu }{Q}_1\},{\mathcal{L}}_{42}=\{L_3+{\epsilon }_1{P}_1+{\epsilon }_2{P}_2+\underline{\mu }{Q}_1\},\hfill \\ & {\mathcal{L}}_{43}=\{L_3+\epsilon P_1+\underline{\mu }{Q}_1+\underline{\nu }{Q}_2\},{\mathcal{L}}_{44}=\{L_3+{\epsilon }_1{P}_1+{\epsilon }_2{P}_2+\underline{\mu }{Q}_1+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{45}=\{L_3+k{L}_2+\epsilon P_1+\underline{\mu }{Q}_1\}{\mathcal{L}}_{46}=\{L_3+k{L}_1+\epsilon P_2\},\hfill \\ & {\mathcal{L}}_{47}=\{L_3+k{L}_1+\underline{\nu }{Q}_2\},{\mathcal{L}}_{48}=\{L_3+k{L}_1+\epsilon P_2+\underline{\nu }{Q}_2\},\hfill \\ & {\mathcal{L}}_{49}=\{L_3+k{L}_1+\ell L_2\}\hfill \end{array}$$

(21)

where ε, ${\epsilon }_1$ and ${\epsilon }_2$ are either one or $-1$, the parameters, k and ℓ, are non-zero, and $\underline{\mu }$ and $\underline{\nu }$ are fermionic constants. These representative subalgebras allow us to determine invariant solutions of the supersymmetric Euler Equations (11) using the symmetry reduction method. It should be noted that some of the subalgebras of the Lie symmetry superalgebra have non-standard invariant structures in the sense that they do not lead to symmetry reductions in the usual sense. For instance, the subalgebra, ${\mathcal{L}}_2=\left\{\underline{\nu }{Q}_2\right\}$, has the nonstandard invariant, $\underline{\nu }f(x,t,{\theta }_1,{\theta }_2,\mathcal{U},\mathcal{R},\mathcal{P})$, where f is an arbitrary function of its arguments. Such non-standard invariants were found by the authors for several other supersymmetric hydrodynamic-type systems, including the supersymmetric sinh-Gordon Equation [23], the supersymmetric Klein-Gordon polynomial Equation [23] and supersymmetric polytropic gas dynamics [4].

4. Invariant Solutions of the Supersymmetric Euler Equations

We now proceed to use the symmetry reduction method to determine the invariants and orbits of nine different subalgebras. Consequently, we obtain the reduced system specific to the given subalgebra.

(1) For subalgebra ${\mathcal{L}}_1=\left\{P_2\right\}$, the invariants are

$$\sigma =x,{\theta }_1,{\theta }_2,\mathcal{U},\mathcal{R},\mathcal{P}.$$

Therefore, the superfields are of the form:

$$\mathcal{U}(x,{\theta }_1,{\theta }_2),\mathcal{R}(x,{\theta }_1,{\theta }_2),\mathcal{P}(x,{\theta }_1,{\theta }_2),$$

and the reduced system is:

$$\begin{array}{cc}& {\theta }_1{\mathcal{R}}_{x{\theta }_2}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x+{\theta }_1{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{x{\theta }_2}{\mathcal{U}}_x-{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x-{\mathcal{P}}_x=0,\hfill \\ & \\ & -{\theta }_1{\mathcal{U}}_{x{\theta }_2}{\mathcal{R}}_x+{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{R}}_x-{\theta }_1{\mathcal{R}}_{x{\theta }_2}{\mathcal{U}}_x+{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x=0\hfill \\ & \\ & -{\theta }_1{\mathcal{U}}_{x{\theta }_2}{\mathcal{P}}_x+{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{P}}_x-k{\theta }_1{\mathcal{P}}_{x{\theta }_2}{\mathcal{U}}_x+k{\mathcal{P}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_x=0\hfill \end{array}$$

(22)

One solution of this system is:

$$\begin{array}{cc}& \mathcal{U}=k_1-\underline{k_3}\underline{\mu }x+{\theta }_1\underline{k_3}+{\theta }_2\underline{k_2}\hfill \\ & \mathcal{R}=\gamma \left(x\right)+{\theta }_1{\psi }_1\left(x\right)+{\theta }_2\underline{k_4}\hfill \\ & \mathcal{P}=k_5+{\theta }_1\underline{k_6}+{\theta }_2\underline{k_7}\hfill \end{array}$$

(23)

where $k_1$ and $k_5$ are arbitrary bosonic constants, $\underline{k_2}$, $\underline{k_3}$, $\underline{k_4}$, $\underline{k_6}$ and $\underline{k_7}$ are arbitrary fermionic constants, γ is an arbitrary bosonic function of x and ${\psi }_1$ is an arbitrary fermionic function of x. This solution is stationary. The choice of arbitrary functions in the density superfield, $\mathcal{R}$, does not influence the velocity and pressure superfields ($\mathcal{U}$ and $\mathcal{P}$, respectively), which are linear in x, ${\theta }_1$ and ${\theta }_2$.

(2) For subalgebra, ${\mathcal{L}}_3=\{P_2+\underline{\nu }{Q}_2\}$, the invariants are:

$$\sigma =x,{\theta }_1,{\tau }_2=\theta -\underline{\nu }t,\mathcal{U},\mathcal{R},\mathcal{P},$$

Therefore, the superfields are of the form:

$$\mathcal{U}(x,{\theta }_1,{\tau }_2),\mathcal{R}(x,{\theta }_1,{\tau }_2),\mathcal{P}(x,{\theta }_1,{\tau }_2),$$

and the reduced system is:

$$\begin{array}{cc}& {\theta }_1\underline{\nu }{\mathcal{R}}_{x{\tau }_2}{\mathcal{U}}_{{\tau }_2}+\underline{\nu }{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\tau }_2}-{\theta }_1{\tau }_2\underline{\nu }{\mathcal{R}}_{x{\tau }_2}{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x+{\theta }_1{\tau }_2{\mathcal{R}}_{x{\tau }_2}\underline{\nu }{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x+{\theta }_1{\mathcal{R}}_{x{\tau }_2}{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x\hfill \\ & -{\theta }_1{\tau }_2\underline{\nu }{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{x{\tau }_2}{\mathcal{U}}_x+{\tau }_2\underline{\nu }{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x-{\theta }_1{\tau }_2\underline{\nu }{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{x{\tau }_2}{\mathcal{U}}_x+{\theta }_1{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{x{\tau }_2}{\mathcal{U}}_x\hfill \\ & +{\tau }_2\underline{\nu }{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x-{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x-{\mathcal{P}}_x=0\hfill \\ & \\ & -\underline{\nu }{\mathcal{R}}_{{\tau }_2}+{\theta }_1{\tau }_2\underline{\nu }{\mathcal{U}}_{x{\tau }_2}{\mathcal{R}}_x-{\theta }_1{\mathcal{U}}_{x{\tau }_2}{\mathcal{R}}_x-{\tau }_2\underline{\nu }{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{R}}_x+{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{R}}_x+{\theta }_1{\tau }_2\underline{\nu }{\mathcal{R}}_{x{\tau }_2}{\mathcal{U}}_x\hfill \\ & -{\theta }_1{\mathcal{R}}_{x{\tau }_2}{\mathcal{U}}_x-{\tau }_2\underline{\nu }{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x+{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x=0\hfill \\ & \\ & -\underline{\nu }{\mathcal{P}}_{{\tau }_2}+{\theta }_1{\tau }_2\underline{\nu }{\mathcal{U}}_{x{\tau }_2}{\mathcal{P}}_x-{\theta }_1{\mathcal{U}}_{x{\tau }_2}{\mathcal{P}}_x-{\tau }_2\underline{\nu }{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{P}}_x+{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{P}}_x+k{\theta }_1{\tau }_2\underline{\nu }{\mathcal{P}}_{x{\tau }_2}{\mathcal{U}}_x\hfill \\ & -k{\theta }_1{\mathcal{P}}_{x{\tau }_2}{\mathcal{U}}_x-k{\tau }_2\underline{\nu }{\mathcal{P}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x+k{\mathcal{P}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_x=0\hfill \end{array}$$

(24)

We have two solutions. The first solution is:

$$\begin{array}{cc}& \mathcal{U}=\underline{\nu }h\left(x\right)+{\theta }_1{\varphi }_1\left(x\right)+{\theta }_2\underline{\nu }\ell \left(x\right)+{\theta }_1{\theta }_2\underline{\nu }f\left(x\right)\hfill \\ & \mathcal{R}=\underline{\nu }j\left(x\right)+{\theta }_2\underline{\nu }g\left(x\right)\hfill \\ & \mathcal{P}=C_1+{\theta }_2\underline{\nu }{C}_2\hfill \end{array}$$

(25)

where $C_1$ and $C_2$ are arbitrary bosonic constants, ℓ, and g are arbitrary bosonic functions of x and f and h, j and ${\varphi }_1$ are arbitrary fermionic functions of x. This solution represents a stationary flow. The choice of arbitrary functions in $\mathcal{U}$ and $\mathcal{R}$ does not influence the pressure superfield, $\mathcal{P}$, which is linear in ${\theta }_2$. The arbitrary functions can be chosen completely arbitrarily. In particular, we can choose bumps, kinks and multiple wave solutions. The postulated form of the arbitrary functions will change as the solution evolves in time. For example, as an illustration, if we choose the arbitrary function:

$$g\left(x\right)=\mathrm{cn}\left[1+cosh{\left(arctan\left(ax\right)\right)}^{-1},\kappa \right],a\in \mathbb{R}$$

(26)

and the modulus, κ, of the elliptic function is such that $0\le {\kappa }^2\le 1$, then the solution has one purely real and one purely imaginary period. For a real argument of the elliptic function, we have the relation:

$$-1\le \mathrm{cn}\left[1+cosh{\left(arctan\left(ax\right)\right)}^{-1},\kappa \right]\le 1$$

(27)

and so, the solution represents a cnoidal wave bump. Asymptotically, this solution approaches a constant. The second solution is:

$$\begin{array}{cc}& \mathcal{U}=\underline{h}(\underline{\mu }x+C_3{\theta }_1+\underline{\mu }\underline{\nu }{\theta }_2)\hfill \\ & \mathcal{R}=\gamma \left(x\right)+{\theta }_1{\psi }_1\left(x\right)+{\theta }_2(\underline{\nu }x+\underline{\nu })+{\theta }_1{\theta }_2\left(\underline{h}\underline{\mu }\right)-{\theta }_1\underline{\nu }t\left(\underline{h}\underline{\mu }\right)\hfill \\ & \mathcal{P}=C_1+{\theta }_1\underline{C_4}+{\theta }_2\underline{C_2}-\underline{\nu }t\underline{C_2}+{\theta }_1{\theta }_2\underline{\nu }\underline{C_5}\hfill \end{array}$$

(28)

where $C_1$ and $C_3$ are arbitrary bosonic constants, $\underline{h}$, $\underline{C_2}$, $\underline{C_4}$ and $\underline{C_5}$ are arbitrary fermionic constants, γ is an arbitrary bosonic function of x and ${\psi }_1$ is an arbitrary fermionic function of x. This solution is non-stationary, and the choice of $\gamma \left(x\right)$ and ${\psi }_1\left(x\right)$ in the density superfield, $\mathcal{R}$, does not influence $\mathcal{U}$ and $\mathcal{P}$. The time dependence in the pressure superfield, $\mathcal{P}$, does not affect $\mathcal{U}$ and $\mathcal{R}$.

(3) For subalgebra, ${\mathcal{L}}_4=\left\{L_2\right\}$, the invariants are:

$$\sigma =x,{\theta }_1,{\tau }_2=t^{-1/2}{\theta }_2,\mathcal{A}=t^{1/2}\mathcal{U},\mathcal{R},\mathcal{C}=t^2\mathcal{P}.$$

Therefore, we have the superfields:

$$\mathcal{U}=t^{-1/2}\mathcal{A}(x,{\theta }_1,{\tau }_2),\mathcal{R}(x,{\theta }_1,{\tau }_2),\mathcal{P}=t^{-2}\mathcal{C}(x,{\theta }_1,{\tau }_2).$$

The corresponding system is given by:

$$\begin{array}{cc}& -\frac{1}{2}{\theta }_1{\mathcal{R}}_{x{\tau }_2}\mathcal{A}+\frac{1}{2}{\theta }_1{\tau }_2{\mathcal{R}}_{x{\tau }_2}{\mathcal{A}}_{{\tau }_2}+\frac{1}{2}{\mathcal{R}}_{{\theta }_1{\tau }_2}\mathcal{A}+\frac{1}{2}{\tau }_2{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{A}}_{{\tau }_2}+\frac{1}{2}{\theta }_1{\tau }_2{\mathcal{R}}_{x{\tau }_2}{\mathcal{A}}_{{\theta }_1}{\mathcal{A}}_x\hfill \\ & +{\theta }_1{\mathcal{R}}_{x{\tau }_2}{\mathcal{A}}_{{\theta }_1{\tau }_2}{\mathcal{A}}_x-\frac{1}{2}{\theta }_1{\tau }_2{\mathcal{R}}_{{\theta }_1{\tau }_2}{\left({\mathcal{A}}_x\right)}^2+\frac{1}{2}{\tau }_2{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{A}}_{{\theta }_1}{\mathcal{A}}_x+{\theta }_1{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{A}}_{x{\tau }_2}{\mathcal{A}}_x\hfill \\ & -{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{A}}_{{\theta }_1{\tau }_2}{\mathcal{A}}_x-{\mathcal{C}}_x=0\hfill \end{array}$$

(29)

$$\begin{array}{cc}& \frac{1}{2}{\theta }_1{\tau }_2{\mathcal{A}}_x{\mathcal{R}}_x-\frac{1}{2}{\tau }_2{\mathcal{A}}_{{\theta }_1}{\mathcal{R}}_x-{\theta }_1{\mathcal{A}}_{x{\tau }_2}{\mathcal{R}}_x+{\mathcal{A}}_{{\theta }_1{\tau }_2}{\mathcal{R}}_x-{\theta }_1{\mathcal{R}}_{x{\tau }_2}{\mathcal{A}}_x+{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{A}}_x-\frac{1}{2}{\tau }_2{\mathcal{R}}_{{\tau }_2}=0\hfill \end{array}$$

(30)

$$\begin{array}{cc}& -2\mathcal{C}-\frac{1}{2}{\tau }_2{\mathcal{C}}_{{\tau }_2}+\frac{1}{2}{\theta }_1{\tau }_2{\mathcal{A}}_x{\mathcal{C}}_x-\frac{1}{2}{\tau }_2{\mathcal{A}}_{{\theta }_1}{\mathcal{C}}_x-{\theta }_1{\mathcal{A}}_{x{\tau }_2}{\mathcal{C}}_x+{\mathcal{A}}_{{\theta }_1{\tau }_2}{\mathcal{C}}_x+2k{\theta }_1{\tau }_2{\mathcal{C}}_x{\mathcal{A}}_x\hfill \\ & -2k{\tau }_2{\mathcal{C}}_{{\theta }_1}{\mathcal{A}}_x-k{\theta }_1{\mathcal{C}}_{x{\tau }_2}{\mathcal{A}}_x+k{\mathcal{C}}_{{\theta }_1{\tau }_2}{\mathcal{A}}_x=0\hfill \end{array}$$

(31)

We obtain the solution:

$$\begin{array}{cc}& \mathcal{U}=t^{-1/2}{\theta }_1\underline{\mu }f\left(x\right)\hfill \\ & \mathcal{R}={\theta }_1\underline{\mu }g\left(x\right)\hfill \\ & \mathcal{P}=0\hfill \end{array}$$

(32)

where f and g are arbitrary bosonic functions of x. This solution is non-stationary. If the functions, $f\left(x\right)$ and $g\left(x\right)$, are chosen to be bounded or periodic, then they remain bounded, and the velocity superfield, $\mathcal{U}$, is dampened with time, t.

(4) For subalgebra, ${\mathcal{L}}_7=\{P_1+\underline{\mu }{Q}_1\}$, the invariants are:

$$\sigma =t,{\tau }_1={\theta }_1-\underline{\mu }x,{\theta }_2,\mathcal{U},\mathcal{R},\mathcal{P}.$$

Therefore, we have the superfields:

$$\mathcal{U}=\mathcal{U}(t,{\tau }_1,{\theta }_2),\mathcal{R}=\mathcal{R}(t,{\tau }_1,{\theta }_2),\mathcal{P}=\mathcal{P}(t,{\tau }_1,{\theta }_2).$$

The corresponding system is given by:

$$\begin{array}{cc}& -\underline{\mu }{\tau }_1{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\mathcal{U}}_t-\underline{\mu }{\tau }_1{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_t-{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\mathcal{U}}_t-{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_t+\underline{\mu }{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{{\tau }_1}\hfill \\ & +\underline{\mu }{\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{t{\tau }_1}{\mathcal{U}}_{{\tau }_1}+\underline{\mu }{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{{\tau }_1}+\underline{\mu }{\mathcal{P}}_{{\tau }_1}=0\hfill \\ & \\ & {\mathcal{R}}_t-\underline{\mu }{\theta }_2{\mathcal{U}}_{t{\tau }_1}{\mathcal{R}}_{{\tau }_1}-\underline{\mu }{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{R}}_{{\tau }_1}-\underline{\mu }{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\mathcal{U}}_{{\tau }_1}-\underline{\mu }{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{{\tau }_1}=0\hfill \\ & \\ & {\mathcal{P}}_t-\underline{\mu }{\theta }_2{\mathcal{U}}_{t{\tau }_1}{\mathcal{P}}_{{\tau }_1}-\underline{\mu }{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{P}}_{{\tau }_1}-k\underline{\mu }{\theta }_2{\mathcal{P}}_{t{\tau }_1}{\mathcal{U}}_{{\tau }_1}-k\underline{\mu }{\mathcal{P}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{{\tau }_1}=0\hfill \end{array}$$

(33)

We have two solutions. The first solution is:

$$\begin{array}{cc}& \mathcal{U}=\alpha \left(t\right)+({\theta }_1-\underline{\mu }x){\varphi }_1\left(t\right)+{\theta }_2{\varphi }_2\left(t\right)+({\theta }_1-\underline{\mu }x){\theta }_2\beta \left(t\right)\hfill \\ & \mathcal{R}=C_4+{\theta }_1\underline{\mu }\ell \hfill \\ & \mathcal{P}=C_3+{\theta }_1\underline{\mu }{C}_2\hfill \end{array}$$

(34)

where $C_2$, $C_3$, $C_4$ and ℓ are arbitrary bosonic constants, α and β are arbitrary bosonic functions of t and ${\varphi }_1$ and ${\varphi }_2$ are arbitrary fermionic functions of t. This solution is non-stationary, and the velocity, $\mathcal{U}$, does not influence the other superfields, $\mathcal{R}$ and $\mathcal{P}$.

It should be noted that both supersymmetries, $Q_1$ and $Q_2$, are broken by the solution (34). Indeed, the supersymmetric transformation generated by $Q_1$ is: $x\to x-\underline{{\eta }_1}{\theta }_1,{\theta }_1\to {\theta }_1+\underline{{\eta }_1}$, so that the superfield, $\mathcal{U}(x,t,{\theta }_1,{\theta }_2)$, is transformed to $\mathcal{U}(x+\underline{{\eta }_1}{\theta }_1,t,{\theta }_1-\underline{{\eta }_1},{\theta }_2)$. Therefore, if the superfield, $\mathcal{U}$, is to be preserved by the supersymmetry, $Q_1$, we must have:

$$(-\underline{{\eta }_1}-\underline{\mu }\underline{{\eta }_1}{\theta }_1){\varphi }_1(t)=0,and(-\underline{{\eta }_1}-\underline{\mu }\underline{{\eta }_1}{\theta }_1){\theta }_2\beta \left(t\right)=0$$

Similarly, if the superfields, $\mathcal{R}$ and $\mathcal{P}$, are to be preserved, we must have $\underline{{\eta }_1}\underline{\mu }\ell =0$ and $\underline{{\eta }_1}\underline{\mu }{C}_2=0,$ respectively. Therefore, the solution (34) is only preserved if $\underline{{\eta }_1}=0$, so the supersymmetry, $Q_1$, is broken by the solution (34). In a similar way, it can be shown that the supersymmetry, $Q_2$, is also broken by the solution (34).

The second solution of system (33) is:

$$\begin{array}{cc}& \mathcal{U}=(\underline{\mu }t\underline{C_1}+C_5+\underline{C_7}\underline{\mu }t)+({\theta }_1-\underline{\mu }x)(-\underline{\mu }t{C}_2-\underline{C_6}+C_8\underline{\mu }t)+{\theta }_2(-\underline{\mu }t{C}_3-\underline{C_7})\hfill \\ & +({\theta }_1-\underline{\mu }x){\theta }_2(\underline{\mu }t\underline{C_4}+C_8)\hfill \\ & \mathcal{R}=(\underline{\mu }{C}_8\underline{C_9}t-\underline{\mu }\underline{C_6}{C}_{10}t+C_{13})+({\theta }_1-\underline{\mu }x)\underline{C_9}+{\theta }_2(-2\underline{\mu }{C}_8{C}_{10}t+\underline{C_{14}})+C_{10}\hfill \\ & \mathcal{P}=(\underline{\mu }{C}_8\underline{C_{11}}t-k\underline{\mu }\underline{C_6}{C}_{12}t+C_{15})+({\theta }_1-\underline{\mu }x)\underline{C_{11}}+{\theta }_2(-(k+1)\underline{\mu }{C}_8{C}_{12}t+\underline{C_{16}})\hfill \\ & +C_{12}\hfill \end{array}$$

(35)

where $C_2$, $C_3$, $C_5$, $C_8$, $C_{10}$, $C_{12}$, $C_{13}$ and $C_{15}$ are arbitrary bosonic constants and $\underline{C_1}$, $\underline{C_4}$, $\underline{C_6}$, $\underline{C_7}$, $\underline{C_9}$, $\underline{C_{11}}$, $\underline{C_{14}}$ and $\underline{C_{16}}$ are arbitrary fermionic constants. The solution is non-stationary. All of the superfields are linear in x and t. In this case, they influence each other, since the superfields have parameters in common.

(5) For subalgebra, ${\mathcal{L}}_8=\left\{L_1\right\}$, the invariants are:

$$\sigma =t,{\tau }_1=x^{-1/2}{\theta }_1,{\theta }_2,\mathcal{A}=x^{-3/2}\mathcal{U},\mathcal{R},\mathcal{C}=x^{-2}\mathcal{P}.$$

Therefore, we have the superfields:

$$\mathcal{U}=x^{3/2}\mathcal{A}(t,{\tau }_1,{\theta }_2),\mathcal{R}(t,{\tau }_1,{\theta }_2),\mathcal{P}=x^2\mathcal{C}(t,{\tau }_1,{\theta }_2).$$

The corresponding system is given by:

$$\begin{array}{cc}& -{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\mathcal{A}}_t-{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{A}}_t+\frac{9}{4}{\tau }_1{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\mathcal{A}}_{{\theta }_2}\mathcal{A}-\frac{3}{2}{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\mathcal{A}}_{{\tau }_1{\theta }_2}\mathcal{A}+\frac{1}{2}{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\mathcal{A}}_{{\tau }_1{\theta }_2}{\tau }_1{\mathcal{A}}_{{\tau }_1}\hfill \\ & +\frac{9}{4}{\tau }_1{\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{A}}_t\mathcal{A}+\frac{9}{4}{\tau }_1{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{A}}_{{\theta }_2}\mathcal{A}-\frac{3}{2}{\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{A}}_{t{\tau }_1}\mathcal{A}+\frac{1}{2}{\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{A}}_{t{\tau }_1}{\tau }_1{\mathcal{A}}_{{\tau }_1}\hfill \\ & -\frac{3}{2}{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{A}}_{{\tau }_1{\theta }_2}\mathcal{A}+\frac{1}{2}{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{A}}_{{\tau }_1{\theta }_2}{\tau }_1{\mathcal{A}}_{{\tau }_1}-2\mathcal{C}+\frac{1}{2}{\tau }_1{\mathcal{C}}_{{\tau }_1}=0\hfill \\ & \end{array}$$

(36)

$$\begin{array}{cc}& {\theta }_2{\mathcal{A}}_{t{\tau }_1}{\tau }_1{\mathcal{R}}_{{\tau }_1}+{\mathcal{A}}_{{\tau }_1{\theta }_2}{\tau }_1{\mathcal{R}}_{{\tau }_1}-3{\theta }_2{\mathcal{R}}_{t{\tau }_1}\mathcal{A}+{\theta }_2{\mathcal{R}}_{t{\tau }_1}{\tau }_1{\mathcal{A}}_{{\tau }_1}-3{\mathcal{R}}_{{\tau }_1{\theta }_2}\mathcal{A}+{\mathcal{R}}_{{\tau }_1{\theta }_2}{\tau }_1{\mathcal{A}}_{{\tau }_1}\hfill \\ & -2{\mathcal{R}}_t=0\hfill \\ & \end{array}$$

(37)

$$\begin{array}{cc}& {\mathcal{C}}_t-3{\tau }_1{\theta }_2{\mathcal{A}}_t\mathcal{C}-3{\tau }_1{\mathcal{A}}_{{\theta }_2}\mathcal{C}+2{\theta }_2{\mathcal{A}}_{t{\tau }_1}\mathcal{C}-\frac{1}{2}{\theta }_2{\mathcal{A}}_{t{\tau }_1}{\tau }_1{\mathcal{C}}_{{\tau }_1}+2{\mathcal{A}}_{{\tau }_1{\theta }_2}\mathcal{C}-\frac{1}{2}{\mathcal{A}}_{{\tau }_1{\theta }_2}{\tau }_1{\mathcal{C}}_{{\tau }_1}\hfill \\ & -3k{\tau }_1{\theta }_2{\mathcal{C}}_t\mathcal{A}-3k{\tau }_1{\mathcal{C}}_{{\theta }_2}\mathcal{A}+\frac{3}{2}k{\mathcal{C}}_{t{\tau }_1}\mathcal{A}-\frac{1}{2}k{\mathcal{C}}_{t{\tau }_1}{\tau }_1{\mathcal{A}}_{{\tau }_1}+\frac{3}{2}k{\mathcal{C}}_{{\tau }_1{\theta }_2}\mathcal{A}\hfill \\ & -\frac{1}{2}k{\mathcal{C}}_{{\tau }_1{\theta }_2}{\tau }_1{\mathcal{A}}_{{\tau }_1}=0\hfill \end{array}$$

(38)

We obtain the solution:

$$\begin{array}{cc}& \mathcal{U}=x^{3/2}\left(\underline{\mu }f\left(t\right)+x^{-1/2}{\theta }_1\underline{\nu }q\left(t\right)+{\theta }_2\underline{\mu }h\left(t\right)+x^{-1/2}{\theta }_1{\theta }_2\underline{\mu }g\left(t\right)\right),\hfill \\ & \mathcal{R}=C_1+{\theta }_2\underline{C_2},\hfill \\ & \mathcal{P}=x^2{\theta }_2\underline{\mu }\hfill \end{array}$$

(39)

where $C_1$ and $C_2$ are arbitrary bosonic constants, h and q are arbitrary bosonic functions of t and f and g are arbitrary fermionic functions of t. This solution is non-stationary. The velocity and pressure superfields are subject to infinite growth in the variable, x. The arbitrary functions, $f\left(t\right)$, $q\left(t\right)$, $h\left(t\right)$ and $g\left(t\right)$, do not influence the superfields, $\mathcal{R}$ and $\mathcal{P}$.

(6) For subalgebra, ${\mathcal{L}}_9=\{P_1+\epsilon P_2\}$, the invariants are:

$$\sigma =x-\epsilon t,{\theta }_1,{\theta }_2,\mathcal{U},\mathcal{R},\mathcal{P}.$$

Obtaining the superfields:

$$\mathcal{U}(\sigma ,{\theta }_1,{\theta }_2),\mathcal{R}(\sigma ,{\theta }_1,{\theta }_2),\mathcal{P}(\sigma ,{\theta }_1,{\theta }_2),$$

we are led to the reduced system:

$$\begin{array}{cc}& {\theta }_1{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }-\epsilon {\theta }_1{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-{\theta }_2{\mathcal{R}}_{\sigma {\theta }_1}{\mathcal{U}}_{\sigma }+\epsilon {\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon {\theta }_1{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & +\epsilon {\theta }_1{\theta }_2{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma {\theta }_1}{\mathcal{U}}_{\sigma }+{\theta }_1{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon {\theta }_1{\theta }_2{\mathcal{R}}_{\sigma {\theta }_1}{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+\epsilon {\theta }_2{\mathcal{R}}_{\sigma {\theta }_1}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & -\epsilon {\theta }_1{\theta }_2{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }+{\theta }_1{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+\epsilon {\theta }_2{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma {\theta }_1}{\mathcal{U}}_{\sigma }-{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma }-{\mathcal{P}}_{\sigma }=0\hfill \end{array}$$

(40)

$$\begin{array}{cc}& -\epsilon {\mathcal{R}}_{\sigma }+\epsilon {\theta }_1{\theta }_2{\mathcal{U}}_{\sigma \sigma }{\mathcal{R}}_{\sigma }-{\theta }_1{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{R}}_{\sigma }-\epsilon {\theta }_2{\mathcal{U}}_{\sigma {\theta }_1}{\mathcal{R}}_{\sigma }+{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{R}}_{\sigma }+\epsilon {\theta }_1{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }\hfill \\ & -{\theta }_1{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon {\theta }_2{\mathcal{R}}_{\sigma {\theta }_1}{\mathcal{U}}_{\sigma }+{\mathcal{R}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma }=0\hfill \end{array}$$

(41)

$$\begin{array}{cc}& -\epsilon {\mathcal{P}}_{\sigma }+\epsilon {\theta }_1{\theta }_2{\mathcal{U}}_{\sigma \sigma }{\mathcal{P}}_{\sigma }-{\theta }_1{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{P}}_{\sigma }-\epsilon {\theta }_2{\mathcal{U}}_{\sigma {\theta }_1}{\mathcal{P}}_{\sigma }+{\mathcal{U}}_{{\theta }_1{\theta }_2}{\mathcal{P}}_{\sigma }+\epsilon k{\theta }_1{\theta }_2{\mathcal{P}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }\hfill \\ & -k{\theta }_1{\mathcal{P}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon k{\theta }_2{\mathcal{P}}_{\sigma {\theta }_1}{\mathcal{U}}_{\sigma }+k{\mathcal{P}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_{\sigma }=0\hfill \end{array}$$

(42)

We have two solutions. The first solution is:

$$\begin{array}{cc}& \mathcal{U}=C_0+{\theta }_1{C}_1\underline{\mu }+{\theta }_2{C}_1\underline{\nu }+{\theta }_1{\theta }_2{C}_2\underline{\mu }\underline{\nu }\hfill \\ & \mathcal{R}=C_6+{\theta }_1\underline{C_7}+{\theta }_2\underline{C_8}+{\theta }_1{\theta }_2{C}_9\hfill \\ & \mathcal{P}=C_{10}+{\theta }_1\underline{C_3}+{\theta }_2\underline{C_4}+{\theta }_1{\theta }_2{C}_5\hfill \end{array}$$

(43)

where $C_0$, $C_1$, $C_2$, $C_5$, $C_6$, $C_9$ and $C_{10}$ are arbitrary bosonic constants and $\underline{C_3}$, $\underline{C_4}$, $\underline{C_7}$ and $\underline{C_8}$ are arbitrary fermionic constants. This solution does not depend on the symmetry variable, σ. However, the superfields are expressible in terms of the Grassmann variables, ${\theta }_1$ and ${\theta }_2$. The second solution is:

$$\begin{array}{cc}& \mathcal{U}=\underline{\mu }\underline{\nu }\ell (x-\epsilon t)+{\theta }_1\underline{\mu }\underline{\nu }q(x-\epsilon t)+{\theta }_2\underline{\mu }\underline{\nu }q(x-\epsilon t)+{\theta }_1{\theta }_2\underline{\mu }j(x-\epsilon t)\hfill \\ & \mathcal{R}=C_3+{\theta }_1\underline{C_4}+{\theta }_2\underline{C_4}\hfill \\ & \mathcal{P}=C_0+{\theta }_1{C}_1\underline{\mu }+{\theta }_2{C}_1\underline{\nu }+{\theta }_1{\theta }_2{C}_2\underline{\mu }\underline{\nu }\hfill \end{array}$$

(44)

where $C_0$, $C_1$, $C_2$ and $C_3$ are arbitrary bosonic constants, $\underline{C_4}$ is an arbitrary fermionic constant, ℓ is an arbitrary bosonic function of its argument and q and j are arbitrary fermionic functions of their arguments. This solution represents a travelling wave. The choice of the functions, ℓ, q and j, in $\mathcal{U}$ does not affect the superfields, $\mathcal{R}$ and $\mathcal{P}$, which depend only on ${\theta }_1$ and ${\theta }_2$.

(7) For subalgebra, ${\mathcal{L}}_{10}=\{P_1+\underline{\nu }{Q}_2\}$, the invariants are:

$$\sigma =t+\underline{\nu }{\theta }_{2}x,{\theta }_1,{\tau }_2={\theta }_2-\underline{\nu }x,\mathcal{U},\mathcal{R},\mathcal{P},$$

and the superfields are:

$$\mathcal{U}(\sigma ,{\theta }_1,{\tau }_2),\mathcal{R}(\sigma ,{\theta }_1,{\tau }_2),\mathcal{P}(\sigma ,{\theta }_1,{\tau }_2),$$

which allows us to reduce to the system:

$$\begin{array}{cc}& -\nu {\theta }_1{\theta }_2{\mathcal{R}}_{\sigma {\tau }_2}{\mathcal{U}}_{\sigma }+\nu {\theta }_1{\mathcal{R}}_{\sigma }{\mathcal{U}}_{\sigma }-\nu {\theta }_1{\theta }_2{\mathcal{R}}_{\sigma {\tau }_2}{\mathcal{U}}_{\sigma }-{\theta }_2{\mathcal{R}}_{\sigma {\theta }_1}{\mathcal{U}}_{\sigma }+\nu x{\mathcal{R}}_{\sigma {\theta }_1}{\mathcal{U}}_{\sigma }-{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{\sigma }\hfill \\ & +\nu {\theta }_2{\mathcal{R}}_{\sigma {\theta }_1}{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\tau }_2}+\nu {\theta }_2{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{\sigma {\theta }_1}{\mathcal{U}}_{{\tau }_2}-\nu {\theta }_2{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{\sigma }+\nu {\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\tau }_2}-\nu {\theta }_2{\mathcal{P}}_{\sigma }\hfill \\ & +\nu {\mathcal{P}}_{{\tau }_2}=0\hfill \end{array}$$

(45)

$$\begin{array}{cc}& {\mathcal{R}}_{\sigma }-\nu {\theta }_2{\mathcal{U}}_{\sigma {\theta }_1}{\mathcal{R}}_{{\tau }_2}+\nu {\theta }_2{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{R}}_{\sigma }-\nu {\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{R}}_{{\tau }_2}-\nu {\theta }_2{\mathcal{R}}_{\sigma {\theta }_1}{\mathcal{U}}_{{\tau }_2}+\nu {\theta }_2{\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{\sigma }\hfill \\ & -\nu {\mathcal{R}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\tau }_2}=0\hfill \end{array}$$

(46)

$$\begin{array}{cc}& {\mathcal{P}}_{\sigma }-\nu {\theta }_2{\mathcal{U}}_{\sigma {\theta }_1}{\mathcal{P}}_{{\tau }_2}+\nu {\theta }_2{\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{P}}_{\sigma }-\nu {\mathcal{U}}_{{\theta }_1{\tau }_2}{\mathcal{P}}_{{\tau }_2}-k\nu {\theta }_2{\mathcal{P}}_{\sigma {\theta }_1}{\mathcal{U}}_{{\tau }_2}+k\nu {\theta }_2{\mathcal{P}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{\sigma }\hfill \\ & -k\nu {\mathcal{P}}_{{\theta }_1{\tau }_2}{\mathcal{U}}_{{\tau }_2}=0\hfill \end{array}$$

(47)

We obtain the solution:

$$\begin{array}{cc}& \mathcal{U}=\alpha (t+\underline{\nu }{\theta }_{2}x)+{\theta }_1{\varphi }_1(t+\underline{\nu }{\theta }_{2}x)+({\theta }_2-\underline{\nu }x){\varphi }_2(t+\underline{\nu }{\theta }_{2}x)+{\theta }_1({\theta }_2-\underline{\nu }x)\beta (t+\underline{\nu }{\theta }_{2}x)\hfill \\ & \mathcal{R}=\underline{C_1}\underline{\nu }+{\theta }_2{C}_2\underline{\nu }\hfill \\ & \mathcal{P}=k_1\underline{C_1}\underline{\nu }+{\theta }_2{k}_1{C}_2\underline{\nu }\hfill \end{array}$$

(48)

where $k_1$ and $C_2$ are arbitrary bosonic constants, $\underline{C_1}$ is an arbitrary fermionic constant, α and β are arbitrary bosonic function of their arguments and ${\varphi }_1$ and ${\varphi }_2$ are arbitrary fermionic functions of their arguments. This solution is a travelling wave involving the Grassmann variable, ${\theta }_2$. The superfields, $\mathcal{R}$ and $\mathcal{P}$, depend linearly on ${\theta }_2$.

(8) For subalgebra, ${\mathcal{L}}_{17}=\{P_1+\epsilon P_2\underline{\mu }{Q}_1\}$, the invariants:

$$\sigma =x-\epsilon t+\epsilon \underline{\mu }t{\theta }_1,{\tau }_1={\theta }_1-\epsilon \underline{\mu }t,{\theta }_2,\mathcal{U},\mathcal{R},\mathcal{P},$$

give us superfields of the form:

$$\mathcal{U}(\sigma ,{\tau }_1,{\theta }_2),\mathcal{R}(\sigma ,{\tau }_1,{\theta }_2),\mathcal{P}(\sigma ,{\tau }_1,{\theta }_2),$$

which lead us to the following reduced system:

$$\begin{array}{cc}& {\tau }_1{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }+{\tau }_1{\theta }_2\underline{\mu }{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{{\tau }_1}+{\tau }_1{\theta }_2\underline{\mu }{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }+\epsilon \underline{\mu }t{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }-\epsilon {\tau }_1{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & -\epsilon {\tau }_1{\mathcal{R}}_{\sigma {\theta }_2}\underline{\mu }{\mathcal{U}}_{{\tau }_1}-\underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+\epsilon {\theta }_2\underline{\mu }t{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }-{\theta }_2\underline{\mu }{\mathcal{R}}_{\sigma }{\mathcal{U}}_{\sigma }+{\theta }_2\underline{\mu }{\tau }_1{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }\hfill \\ & +{\theta }_2{\mathcal{R}}_{\sigma {\tau }_1}\underline{\mu }{\tau }_1{\mathcal{U}}_{\sigma }-{\theta }_2{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }-{\theta }_2{\mathcal{R}}_{\sigma {\tau }_1}\underline{\mu }{\mathcal{U}}_{{\tau }_1}+\underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon \underline{\mu }{\tau }_1{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }+\epsilon {\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & +\epsilon \underline{\mu }{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{{\tau }_1}-{\tau }_1{\theta }_2{\mathcal{R}}_{\sigma \sigma }\underline{\mu }t{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon {\tau }_1{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon {\tau }_1{\theta }_2\underline{\mu }{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & -\underline{\mu }t{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }-{\tau }_1{\theta }_2{\mathcal{R}}_{\sigma {\theta }_2}\underline{\mu }t{\mathcal{U}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }+\epsilon {\tau }_1{\theta }_2{\mathcal{R}}_{\sigma {\theta }_2}\underline{\mu }{\left({\mathcal{U}}_{\sigma }\right)}^2+\epsilon {\tau }_1{\theta }_2{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }\hfill \\ & +\underline{\mu }t{\theta }_2{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }+\epsilon {\tau }_1{\mathcal{R}}_{\sigma {\theta }_2}\underline{\mu }t{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+{\tau }_1{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }+\epsilon \underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & +{\tau }_1{\theta }_2\underline{\mu }t{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon {\tau }_1{\theta }_2\underline{\mu }{\mathcal{R}}_{\sigma }{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon {\tau }_1{\theta }_2{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-\underline{\mu }t{\theta }_2{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & -{\theta }_2\underline{\mu }t{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }+\epsilon {\theta }_2\underline{\mu }{\mathcal{R}}_{\sigma }{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon {\theta }_2\underline{\mu }{\tau }_1{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }+{\theta }_2{\mathcal{R}}_{\sigma {\tau }_1}\underline{\mu }t{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & +\epsilon {\theta }_2{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }-{\tau }_1{\theta }_2\underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }-\epsilon {\tau }_1{\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }-\epsilon {\tau }_1{\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}\underline{\mu }{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }\hfill \\ & -\underline{\mu }t{\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }+\epsilon {\tau }_1\underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+{\tau }_1{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+\epsilon \underline{\mu }t{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & +{\theta }_2\underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }-{\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}\underline{\mu }t{\mathcal{U}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }+\epsilon {\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}\underline{\mu }{\left({\mathcal{U}}_{\sigma }\right)}^2-\epsilon {\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}\underline{\mu }{\tau }_1{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }\hfill \\ & +\epsilon {\theta }_2{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }-\epsilon \underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon \underline{\mu }t{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }-{\mathcal{P}}_{\sigma }=0\hfill \end{array}$$

(49)

$$\begin{array}{cc}& \epsilon \underline{\mu }{\tau }_1{\mathcal{R}}_{\sigma }-\epsilon {\mathcal{R}}_{\sigma }-\epsilon \underline{\mu }{\mathcal{R}}_{{\tau }_1}+\epsilon {\tau }_1{\theta }_2{\mathcal{U}}_{\sigma \sigma }{\mathcal{R}}_{\sigma }+\epsilon {\tau }_1{\theta }_2\underline{\mu }{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{R}}_{\sigma }+\underline{\mu }t{\theta }_2{\mathcal{U}}_{\sigma \sigma }{\mathcal{R}}_{\sigma }-{\tau }_1{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{R}}_{\sigma }\hfill \\ & -\epsilon \underline{\mu }t{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{R}}_{\sigma }+{\theta }_2\underline{\mu }t{\mathcal{U}}_{\sigma \sigma }{\mathcal{R}}_{\sigma }-\epsilon {\theta }_2\underline{\mu }{\mathcal{U}}_{\sigma }{\mathcal{R}}_{\sigma }+\epsilon {\theta }_2\underline{\mu }{\tau }_1{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{R}}_{\sigma }-\epsilon {\theta }_2{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{R}}_{\sigma }\hfill \\ & +\epsilon \underline{\mu }t{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{R}}_{\sigma }+{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{R}}_{\sigma }+\epsilon {\tau }_1{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }+\epsilon {\tau }_1{\theta }_2\underline{\mu }{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }+\underline{\mu }t{\theta }_2{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }-{\tau }_1{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }\hfill \\ & -\epsilon \underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+{\theta }_2\underline{\mu }t{\mathcal{R}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }-\epsilon {\theta }_2\underline{\mu }{\mathcal{R}}_{\sigma }{\mathcal{U}}_{\sigma }+\epsilon {\theta }_2\underline{\mu }{\tau }_1{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }-\epsilon {\theta }_2{\mathcal{R}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }\hfill \\ & +\epsilon \underline{\mu }t{\mathcal{R}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+{\mathcal{R}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }=0\hfill \end{array}$$

(50)

$$\begin{array}{cc}& \epsilon \underline{\mu }{\tau }_1{\mathcal{P}}_{\sigma }-\epsilon {\mathcal{P}}_{\sigma }-\epsilon \underline{\mu }{\mathcal{P}}_{{\tau }_1}+\epsilon {\tau }_1{\theta }_2{\mathcal{U}}_{\sigma \sigma }{\mathcal{P}}_{\sigma }+\epsilon {\tau }_1{\theta }_2\underline{\mu }{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{P}}_{\sigma }+\underline{\mu }t{\theta }_2{\mathcal{U}}_{\sigma \sigma }{\mathcal{P}}_{\sigma }-{\tau }_1{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{P}}_{\sigma }\hfill \\ & -\epsilon \underline{\mu }t{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{P}}_{\sigma }+{\theta }_2\underline{\mu }t{\mathcal{U}}_{\sigma \sigma }{\mathcal{P}}_{\sigma }-\epsilon {\theta }_2\underline{\mu }{\mathcal{U}}_{\sigma }{\mathcal{P}}_{\sigma }+\epsilon {\theta }_2\underline{\mu }{\tau }_1{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{P}}_{\sigma }-\epsilon {\theta }_2{\mathcal{U}}_{\sigma {\tau }_1}{\mathcal{P}}_{\sigma }\hfill \\ & +\epsilon \underline{\mu }t{\mathcal{U}}_{\sigma {\theta }_2}{\mathcal{P}}_{\sigma }+{\mathcal{U}}_{{\tau }_1{\theta }_2}{\mathcal{P}}_{\sigma }+\epsilon k{\tau }_1{\theta }_2{\mathcal{P}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }+\epsilon k{\tau }_1{\theta }_2\underline{\mu }{\mathcal{P}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }+\underline{\mu }kt{\theta }_2{\mathcal{P}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }\hfill \\ & -k{\tau }_1{\mathcal{P}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }-\epsilon \underline{\mu }kt{\mathcal{P}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+k{\theta }_2\underline{\mu }t{\mathcal{P}}_{\sigma \sigma }{\mathcal{U}}_{\sigma }-\epsilon k{\theta }_2\underline{\mu }{\mathcal{P}}_{\sigma }{\mathcal{U}}_{\sigma }+\epsilon k{\theta }_2\underline{\mu }{\tau }_1{\mathcal{P}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }\hfill \\ & -\epsilon k{\theta }_2{\mathcal{P}}_{\sigma {\tau }_1}{\mathcal{U}}_{\sigma }+\epsilon k\underline{\mu }t{\mathcal{P}}_{\sigma {\theta }_2}{\mathcal{U}}_{\sigma }+k{\mathcal{P}}_{{\tau }_1{\theta }_2}{\mathcal{U}}_{\sigma }=0\hfill \end{array}$$

(51)

We obtain the solution:

$$\begin{array}{cc}& \mathcal{U}=C_0+{\theta }_1\underline{\mu }+{\theta }_2\underline{C_1},\hfill \\ & \mathcal{R}=\underline{\mu }\underline{C_8}-\underline{\mu }\underline{C_7}(x-\epsilon t)+({\theta }_1-\epsilon \underline{\mu }t)\underline{C_7}+{\theta }_2\underline{\mu }{C}_{10}+{\theta }_2\underline{\mu }{C}_9(x-\epsilon t)+({\theta }_1-\epsilon \underline{\mu }t){\theta }_2{C}_9,\hfill \\ & \mathcal{P}=C_3+{\theta }_1\underline{\mu }\ell +{\theta }_2\underline{C_5}+{\theta }_1{\theta }_2\underline{C_2}\underline{\mu }\hfill \end{array}$$

(52)

where $C_0$, $C_3$, $C_9$, $C_{10}$ and ℓ are arbitrary bosonic constants and $\underline{C_1}$, $\underline{C_2}$, $\underline{C_5}$, $\underline{C_7}$ and $\underline{C_8}$ are arbitrary fermionic constants. This solution is non-stationary. The density, $\mathcal{R}$, is a travelling wave. The superfields, $\mathcal{U}$ and $\mathcal{P}$, depend only on ${\theta }_1$ and ${\theta }_2$.

(9) For subalgebra, ${\mathcal{L}}_{30}=\{L_3+\epsilon P_1\}$, the invariants are:

$$t,{\theta }_1,{\theta }_2,\mathcal{U},\mathcal{B}=e^{-\epsilon x}\mathcal{R},\mathcal{C}{e}^{-\epsilon x}\mathcal{P}.$$

Therefore, we have the superfields:

$$\mathcal{U}(t,{\theta }_1,{\theta }_2),\mathcal{R}=e^{\epsilon x}\mathcal{B}(t,{\theta }_1,{\theta }_2),\mathcal{P}=e^{\epsilon x}\mathcal{C}(t,{\theta }_1,{\theta }_2),$$

and the reduced system is:

$$\begin{array}{cc}& {\theta }_1{\theta }_2\epsilon {\mathcal{B}}_t{\mathcal{U}}_t-{\theta }_2{\mathcal{B}}_{t{\theta }_1}{\mathcal{U}}_t+{\theta }_1\epsilon {\mathcal{B}}_{{\theta }_2}{\mathcal{U}}_t-{\mathcal{B}}_{{\theta }_1{\theta }_2}{\mathcal{U}}_t-\epsilon \mathcal{C}=0,\hfill \\ & \\ & {\mathcal{B}}_t+{\theta }_2{\mathcal{U}}_{t{\theta }_1}\epsilon \mathcal{B}+{\mathcal{U}}_{{\theta }_1{\theta }_2}\epsilon \mathcal{B}=0,\hfill \\ & \\ & {\mathcal{C}}_t+{\theta }_2{\mathcal{U}}_{t{\theta }_1}\epsilon \mathcal{C}+{\mathcal{U}}_{{\theta }_1{\theta }_2}\epsilon \mathcal{C}=0\hfill \end{array}$$

(53)

We have two solutions. The first solution is:

$$\begin{array}{cc}& \mathcal{U}=C_1+{\theta }_1\underline{C_2}+{\theta }_2\underline{C_3}+{\theta }_1{\theta }_2{C}_4,\hfill \\ & \mathcal{R}=e^{\epsilon (x-C_{4}t)}\left[k_1+{\theta }_1\underline{k_3}+{\theta }_2\underline{k_4}+{\theta }_1{\theta }_2{k}_5\right],\hfill \\ & \mathcal{P}=0\hfill \end{array}$$

(54)

where $C_1$, $C_4$, $k_1$ and $k_5$ are arbitrary bosonic constants and $\underline{C_2}$, $\underline{C_3}$, $\underline{k_3}$ and $\underline{k_4}$ are arbitrary fermionic constants. This is a non-stationary solution. The density superfield, $\mathcal{R}$, contains a factor of $e^{\epsilon (x-C_{4}t)}$. Depending on whether ε is $-1$ or one, at a fixed position, x, we have either damping or growth. The velocity is not influenced, however. The second solution is:

$$\begin{array}{cc}& \mathcal{U}=At+B+{\theta }_1\underline{C_6}+{\theta }_1{\theta }_2\underline{\mu }\underline{C_5},\hfill \\ & \mathcal{R}=e^{\epsilon x}\left[{\theta }_1\underline{\mu }+{\theta }_2\underline{\mu }+{\theta }_1{\theta }_2\underline{\mu }\underline{C_5}\right],\hfill \\ & \mathcal{P}=A{e}^{\epsilon x}\left[-\epsilon \underline{\mu }\underline{C_5}+{\theta }_1\underline{\mu }\right]\hfill \end{array}$$

(55)

where A and B are arbitrary bosonic constants and $\underline{C_5}$ and $\underline{C_6}$ are arbitrary fermionic constants. This is a non-stationary solution. The density and pressure contain a damping (or growth) factor. The velocity superfield is linear in time t. The velocity and pressure are linked through the constant, A. Depending on the sign of A, the velocity, $\mathcal{U}$, is either increasing or decreasing.

5. Final Remarks and Open Problems

In this paper, we discuss first the construction of a supersymmetric extension of the Euler system (1) through a superspace and superfield formalism. This analysis includes a supersymmetric extension of a one-dimensional ideal compressible non-viscous fluid flow. This allows us to determine a Lie superalgebra of infinitesimal symmetries, which generate the Lie point symmetries of the Euler system of Equations. We observe that, in analogy with the classical Euler system, the Lie symmetry superalgebra of the supersymmetric Euler system contains three independent dilations. Next, we discuss the classification of its subalgebras. Through the use of a generalized version of the symmetry reduction method, we have demonstrated how to find exact invariant solutions of the supersymmetric model. A systematic use of the structure of the invariance supergroup of the Euler Equations allows us to generate all symmetry variables. For certain subalgebras, the invariants have a nonstandard structure and, therefore, do not lead to invariant solutions. This phenomenon of nonstandard invariants has also been observed, among other places, in the analysis of the symmetries of the supersymmetric sine-Gordon Equation [24] and of supersymmetric polytropic gas dynamics [4]. The SRM enables us to reduce, after some transformations, the basic system of PDEs to many possible ODEs. In each case, the three superfields can be decomposed in terms of its independent fermionic variables, ${\theta }_1$ and ${\theta }_2$. This allows us to determine the invariant solutions of the supersymmetric Euler system component-wise. For nine specific subalgebras involving different types of generators, a number of invariant solutions were found. These included solutions involving several arbitrary functions of one variable. In particular, this includes bumps, kinks and multiple wave solutions. If these arbitrary functions depend on x, the postulated form will change as the solution evolves in time. Such a large number of arbitrary function degrees of freedom were not found in the previous analyses by the authors of supersymmetric hydrodynamic systems (i.e., polytropic gas [4], integrable models (sine-Gordon, sinh-Gordon, polynomial Klein-Gordon [23,24]). Some of the obtained solutions involved damping and growth. However, no blow-up phenomenon (gradient catastrophe) was observed for any of our solutions.

In the future, it would be interesting to expand our analysis in several directions. One such possibility would be to apply the above supersymmetric extension methods to the Euler system in higher dimensions. Due to the complexity of the calculations involved, this would require the development of a computer algebra Lie symmetry package capable of handling odd and even Grassmann variables. To the best of our knowledge, such a package does not presently exist. Conservation laws are well-established for the classical Euler Equations, but it has been observed that, for many hydrodynamic-type systems, such conservation laws are broken in their corresponding supersymmetric extensions (see, e.g., [4] and references therein). The question of which quantities are conserved by the supersymmetric model still remains an open question for hydrodynamic-type Equations. These topics will be investigated in our future work.