Skip to main content

Existence of the weak solutions for the compressible magnetohydrodynamic flows driven by stochastic forcing

Abstract

In this paper, we study the existence of weak solutions for the compressible magnetohydrodynamics flows driven by stochastic external forcing. Our method is based on solving the system for each fixed representative of the random variable and applying an abstract result on the measurability of multi-valued maps.

1 Introduction and main results

This paper considers the following compressible magnetohydrodynamic flows driven by a stochastic external force in the isentropic case [13]:

$$\begin{aligned}& d\rho + \mathbf{div}_{x}(\rho\mathbf{u})\, dt=0, \end{aligned}$$
(1.1)
$$\begin{aligned}& d(\rho\mathbf{u}) + \bigl(\mathbf{div}_{x}(\rho\mathbf{ u}\otimes\mathbf{u})+\nabla P(\rho)-(\nabla\times\mathbf{H})\times \mathbf{H} \\ & \quad {}-\nu\triangle\mathbf{u}-(\nu+\lambda)\nabla(\mathbf{div}_{x} \, \mathbf {u})\bigr)\, dt=\rho\, dW, \end{aligned}$$
(1.2)
$$\begin{aligned}& d\mathbf{H} - \bigl(\nabla\times(\mathbf{u}\times\mathbf{H})+ \nabla \times\bigl(\nu'\nabla\times\mathbf{H}\bigr) \bigr)\, dt=0, \end{aligned}$$
(1.3)

where \(\rho(t,x,\omega)\), \(\mathbf{u}(t,x,\omega)\), \(\mathbf {H}(t,x,\omega)\) are functions of the time \(t\in(0,T)\), the spatial position \(x\in\Omega\), and \(\omega\in\mathcal{O}=\{\mathcal{O},\mathcal {B},\mu\}\), \(\mathcal{O}\) is a topological probability case, with the family of Borel sets \(\mathcal{B}\), and a regular probability measure μ. The symbol denotes the Kronecker tensor product. \(\mathbf{H}(t,x,\omega)\) satisfies

$$ \mathbf{div}_{x}\, \mathbf{H}=0, $$
(1.4)

\(P(\rho)\) is the pressure, the viscosity coefficients of the flow satisfy \(2\nu+3\lambda>0\) and \(\nu>0\); \(\nu'>0\) is the magnetic diffusivity acting as a magnetic diffusion coefficient of the magnetic field, and all these kinetic coefficients and the magnetic diffusivity are independent of the magnitude and direction of the magnetic field, the perturbation W is a random variable represented for a.a. ω by a bounded function, sufficiently regular with respect to the spatial variable \(x\in\Omega\). Moreover, we supplement this with the no-slip boundary condition,

$$ \mathbf{u}|_{\partial\Omega}=0, \qquad \mathbf{H}|_{\partial\Omega}=0, $$

and the initial conditions

$$\begin{aligned}& \rho(0,\cdot)=\rho_{0},\qquad \rho\mathbf{u}(0,\cdot)=(\rho \mathbf{u})_{0},\qquad (\rho\mathbf{u})_{0}=0, \quad \mbox{if }\rho_{0}=0, \\& \mathbf{H}(0,\cdot)=\mathbf{H}_{0},\qquad \mathbf{div}_{x} \, \mathbf{H}_{0}=0. \end{aligned}$$

Although the electric field does not appear in the MHD equations (1.1)-(1.3), from (1.3), the electric field E is induced by moving conductive flow in the magnetic field, having the following relationship with the magnetic field H and the velocity u:

$$ \mathbf{E}=\nu'\nabla\times\mathbf{H}-\mathbf{u}\times\mathbf{H}. $$

The stochastic incompressible and compressible fluid systems have attracted great attention by considering that a system in reality is usually affected by external perturbations which in many cases are of great uncertainty or random influence. These random effects are introduced not only to compensate for the defects in some deterministic models but also to reveal the intrinsic phenomena. In this direction, the surveys corresponding to the stochastic case are in [46], and much work has been done on the so-called weak solutions in the probabilistic sense (or martingale solutions) for the stochastic incompressible Navier-Stokes equations; for example, see [79]. For the one dimensional compressible stochastic Navier-Stokes system, we refer to [10, 11] for details. Most recently, Feireisl et al. [12] gave results on the existence of weak solutions for the three dimensional stochastic compressible Navier-Stokes system. Their results included a number of stochastic external forces, e.g. Lévy noise.

The main difficulty of the study of MHD driven by stochastic external forcing is the presence of the magnetic field and its interaction with the hydrodynamic motion in the MHD flow of large oscillation. Inspired by the work of [12], this paper is devoted to the study of the existence of weak solutions for three dimensional stochastic compressible magnetohydrodynamic flows. One of the difficulties is to obtain a suitable energy estimate when there is a magnetic field and its interaction with the hydrodynamic motion is driven by the stochastic external forcing. To the best of our knowledge, there are only results on the existence of weak solutions for the stochastic incompressible magnetohydrodynamic flows, cf. [13, 14].

Inspired by the definition of the weak solutions to the compressible Naiver-Stokes equation in [12], we give the definition of the weak solutions to the problem (1.1)-(1.4).

The random variables \(\rho(t,x,\omega)\), \(\mathbf{u}(t,x,\omega)\), \(\mathbf{H}(t,x,\omega)\) is said to be a weak solution to the problem (1.1)-(1.4) of the compressible MHD equations with stochastic external forces if the following conditions hold:

• The density \(\rho\geq0\):

$$\begin{aligned}& \rho(\cdot,\omega)\in\mathbf{C}\bigl([0,T];\mathbf{L}^{1}(\Omega) \bigr), \\& \rho\in\mathbf{L}^{\infty}\bigl(0,T:\mathbf{L}^{q}\bigr) \quad \mbox{for certain }q>1 \mbox{ for a.a. }\omega\in\mathcal{O}, \end{aligned}$$

and the family of integral identities

$$\begin{aligned}& \int_{\Omega}\bigl(\rho(\tau,x,\omega) + b\bigl( \rho(\tau,x,\omega)\bigr)\bigr)\varphi (x)\,dx-\int_{\Omega}\bigl( \rho_{0}(x,\omega)+b\bigl(\rho_{0}(x,\omega)\bigr)\bigr) \varphi (x)\,dx \\& \quad = \int_{0}^{\tau}\int_{\Omega} \bigl(\rho(t,x,\omega)\mathbf{u}(t,x,\omega )\cdot\nabla_{x}\varphi(x) \\& \qquad {}+\bigl(b\bigl(\rho(t,x,\omega)\bigr)-b'\bigl(\rho(t,x, \omega)\bigr)\rho(t,x,\omega)\mathbf {div}_{x}\, \mathbf{u}(t,x,\omega) \varphi(x)\bigr)\bigr)\,dx\,dt \end{aligned}$$
(1.5)

for any test function \(\varphi\in\mathbf{C}^{1}(\bar{\Omega})\), any \(b\in\mathbf{C}^{\infty}_{c}[0,\infty)\), and a.a. \(\omega\in\mathcal{O}\).

• The momentum equation:

$$\begin{aligned}& \mathbf{u}(\cdot,\omega)\in\mathbf{L}^{2}\bigl(0,T; \mathbf{W}^{1,2}_{0}\bigl(\Omega ;\mathbf{R}^{3} \bigr)\bigr) \quad \mbox{for a.a. }\omega\in\mathcal{O}, \\& \rho\mathbf{u}\in\mathbf{L}^{\infty}\bigl(0,T:\mathbf{L}^{\frac {2q}{q+1}} \bigl(\Omega;\mathbf{R}^{3}\bigr)\bigr), \\& \rho(\mathbf{u}-W)\in \mathbf {C}_{\mathrm{weak}}\bigl([0,T];\mathbf{L}^{\frac{2q}{q+1}}\bigl(\Omega; \mathbf{R}^{3}\bigr)\bigr)\quad \mbox{for a.a. }\omega\in\Omega, \\& P(\rho)\in\mathbf{L}^{\infty}\bigl(0,T:\mathbf{L}^{1}(\Omega) \bigr)\cap\mathbf {L}^{\gamma}\bigl((0,T)\times\Omega\bigr)\quad \mbox{for some }\gamma>1\mbox{ for a.a. }\omega\in \mathcal{O}, \end{aligned}$$

and the following weak formulations hold:

$$\begin{aligned}& \int_{\Omega}\bigl(\rho(\tau,x,\omega) \bigl( \mathbf{u}(\tau,x,\omega)-W(\tau ,x,\omega)\bigr)\cdot\varphi(x)\bigr)\,dx-\int _{\Omega}(\rho\mathbf{u})_{0}(x,\omega )\cdot\varphi(x) \\& \quad = \int_{0}^{\tau}\int_{\Omega} \bigl(\rho(t,x,\omega)\mathbf{u}(t,x,\omega )\otimes\mathbf{u}(t,x,\omega): \nabla_{x}\varphi(x)+P\bigl(\rho(t,x,\omega )\bigr) \mathbf{div}_{x}\, \varphi(x)\bigr)\,dx\,dt \\& \qquad {} -\int_{0}^{\tau}\int_{\Omega} \biggl(\frac{1}{2}\bigl|\mathbf{H}(t,x,\omega)\bigr|^{2}\cdot \nabla_{x}\varphi(x)+\bigl(\mathbf{H}(t,x,\omega)\cdot\nabla\bigr) \mathbf {H}(t,x,\omega)\cdot\varphi(x)\biggr)\,dx\,dt \\& \qquad {} -\int_{0}^{\tau}\int_{\Omega} \bigl(\nu\nabla_{x}\mathbf{u}(t,x,\omega)\cdot \nabla_{x} \varphi(x)+(\nu+\lambda) \bigl(\mathbf{div}_{x}\, \mathbf{u}(t,x,\omega )\bigr)\cdot\nabla_{x}\varphi(x) \\& \qquad {} +\rho(t,x,\omega)\mathbf{u}(t,x,\omega)\cdot\nabla_{x} \bigl(W(t,x,\omega )\varphi(x)\bigr)\bigr)\,dx\,dt \end{aligned}$$
(1.6)

for all \(\varphi\in\mathbf{C}^{\infty}_{c}(\Omega;\mathbf{R}^{3})\), and for a.a. \(\omega\in\mathcal{O}\).

• The induction equation of electromagnetism:

$$ \mathbf{H}(\cdot,\omega)\in\mathbf{L}^{2}\bigl(0,T; \mathbf{W}^{1,2}_{0}\bigl(\Omega ;\mathbf{R}^{3} \bigr)\bigr) \quad \mbox{for a.a. }\omega\in\mathcal{O}, $$

and the following weak formulations hold:

$$\begin{aligned}& \int_{0}^{T}\int _{\Omega}\bigl(\bigl(d\mathbf{H}(\tau,x,\omega),\varphi(x)\bigr)+ \nu '\bigl(\nabla_{x}\times\mathbf{H}(\tau,x,\omega), \nabla_{x}\times\varphi (x)\bigr) \\& \quad {}-\bigl(\mathbf{u}(\tau,x,\omega)\times\mathbf{H}(\tau,x,\omega),\nabla _{x}\times\varphi(x)\bigr)\bigr)\,dx\,dt=0 \end{aligned}$$
(1.7)

for all \(\varphi\in\mathbf{C}^{\infty}_{c}(\Omega;\mathbf{R}^{3})\), and for a.a. \(\omega\in\mathcal{O}\).

• The energy inequality:

$$\begin{aligned}& \int_{\Omega}\biggl(\frac{1}{2}\rho(\tau,x, \omega)\bigl\vert \mathbf{u}(\tau,x,\omega )-W(\tau,x,\omega)\bigr\vert ^{2}+P\bigl(\rho(\tau,x,\omega)\bigr)+\frac{1}{2}\bigl\vert \mathbf{H}(\tau ,x,\omega)\bigr\vert ^{2}\biggr)\,dx \\& \qquad {} + \int_{\Omega}\bigl(\nu\bigl\vert \nabla_{x}\mathbf{u}(t,x,\omega)\bigr\vert ^{2}+(\nu+ \lambda )\bigl\vert \mathbf{div}_{x}\, \mathbf{u}(t,x,\omega)\bigr\vert ^{2}+\nu'\vert \nabla_{x}\times\mathbf {H}\vert ^{2}\bigr)\,dx \\& \quad \leq \int_{\Omega}\biggl(\frac{1}{2}\rho(s,x,\omega) \bigl\vert \mathbf{u}(s,x,\omega )-W(s,x,\omega)\bigr\vert ^{2}+P \bigl(\rho(s,x,\omega)\bigr)+\frac{1}{2}\bigl\vert \mathbf{H}(s,x,\omega ) \bigr\vert ^{2}\biggr)\,dx \\& \qquad {} -\int_{0}^{\tau}\int_{\Omega} \biggl(\rho\mathbf{u}\otimes\mathbf{u}:\nabla _{x}W+P(\rho) \mathbf{div}_{x}\, W-\frac{1}{2}\rho\mathbf{u}\cdot\nabla _{x}\vert W\vert ^{2}\biggr) (t,x,\omega)\,dx\,dt \\& \qquad {} -\int_{0}^{\tau}\int_{\Omega} \bigl(\nu\nabla_{x}\mathbf{u}\cdot\nabla_{x}W+(\nu +\lambda) \mathbf{div}_{x}\, \mathbf{u}\cdot\nabla_{x}W\bigr) (t,x, \omega )\,dx\,dt \\& \qquad {} +\int_{0}^{\tau}\int_{\Omega} \biggl(\frac{1}{2}\vert \mathbf{H}\vert ^{2}\cdot\nabla _{x}W+(\mathbf{H}\cdot\nabla)\mathbf{H}\cdot W\biggr) (t,x,\omega)\,dx \,dt \\& \qquad {} +\int_{0}^{\tau}\int_{\Omega}( \mathbf{u}\times\mathbf{H},\nabla_{x}\times W) (\tau,x,\omega)\,dx\,dt \end{aligned}$$
(1.8)

for a.a. \(\tau\geq s\geq0\) including \(s=0\), and a.a. \(\omega\in\mathcal {O}\), with

$$ \mathcal{P}(\rho)=\rho\int_{0}^{\rho} \frac{P(z)}{z^{2}}\, dz. $$
(1.9)

We remark that the family of integral identities (1.5) is taken from [12], which is a weak formulation of the renormalized equation of continuity introduced by [15]. In order to derive the family of integral identities (1.6)-(1.7), we use the following well-known identities:

$$\begin{aligned}& d(\rho\mathbf{u})-\rho \, dW=d\bigl[\rho(\mathbf{u}-W)\bigr]- \mathbf{div}_{x}(\rho \mathbf{u})W\, dt, \\& (\nabla\times\mathbf{H})\times\mathbf{H}=\frac{1}{2}\nabla\bigl(|\mathbf {H}|^{2}\bigr)-(\mathbf{H}\cdot\nabla)\mathbf{H}, \\& \nabla\times(\nabla\times\mathbf{H})=\nabla \mathbf{div}\, \mathbf{H}-\triangle \mathbf{H}, \\& \nabla\times(\mathbf{u}\times\mathbf{H})=(\mathbf{div}\, \mathbf {H}+\mathbf{H} \cdot\nabla)\mathbf{u}-(\mathbf{div}\, \mathbf{u}+\mathbf {u}\cdot\nabla) \mathbf{H}. \end{aligned}$$

Before stating our main result, we introduce the function space \(\mathcal{ID}\) for the data as

$$\begin{aligned} \mathcal{ID} =&\biggl\{ \bigl[\rho_{0},(\rho\mathbf{u})_{0}, \mathbf{H}_{0}\bigr]\Big|\rho_{0}\in \mathbf{L}^{q}( \Omega), (\rho\mathbf{u}_{0})\in\mathbf{L}^{\frac {2q}{q+1}}\bigl(\Omega; \mathbf{R}^{3}\bigr), \mathbf{H}_{0}\in\mathbf{L}^{2}( \Omega), \rho_{0}\geq0, \\ &\int_{\Omega}\rho_{0}\,dx=M>0, \int _{\Omega}\frac{1}{2}\frac{|(\rho\mathbf {u})_{0}|^{2}}{\rho_{0}}+P( \rho_{0})+\frac{1}{2}|\mathbf{H}_{0}|^{2} \,dx\leq E_{0}\biggr\} , \end{aligned}$$

which is a closed convex subset of the separable Banach space \(\mathbf {L}^{q}(\Omega)\times\mathbf{L}^{\frac{2q}{q+1}}(\Omega;\mathbf {R}^{3})\times\mathbf{W}_{0}^{1,2}(\Omega)\). In particular, it is a Suslin space.

Theorem 1.1

Let \(\Omega\in\mathbf{R}^{3}\) be a bounded Lipschitz domain. Assume that the pressure \(P\in\mathbf{C}[0,\infty)\cap\mathbf {C}^{2}[0,\infty)\) and satisfies

$$ P(0)=0, \qquad P'(0)>0,\qquad \lim _{\rho\rightarrow\infty}\frac{P'(\rho)}{\rho ^{p-1}}=P_{\infty}>0,\qquad p> \frac{3}{2}, $$
(1.10)

and \(\mathcal{W}\) is a complete separable metric space satisfying:

  1. (i)
    $$ \mathcal{W}\subset\mathbf{L}_{\mathrm{weak}}^{\infty}\bigl(0,T; \mathbf{W}_{0}^{1,\infty }\bigl(\Omega;\mathbf{R}^{3} \bigr)\bigr); $$
  2. (ii)

    if \(W_{n}\rightarrow W\) in \(\mathcal{W}\), then

    $$\begin{aligned}& \operatorname{ess}\sup_{t\in(0,T)}\bigl\Vert W_{n}(t, \cdot)\bigr\Vert _{\mathbf{W}_{0}^{1,\infty}(\Omega ;\mathbf{R}^{3})}\leq c, \\& W_{n}\rightarrow W\quad \textit{weakly } (*)\textit{ in } \mathbf{L}^{\infty}\bigl(0,T;\mathbf {W}^{1,\infty}(\Omega)\bigr). \end{aligned}$$

Let

$$ \omega\in\mathcal{O}\mapsto\bigl\{ \bigl[\rho_{0}(\cdot,\omega),(\rho \mathbf {u})_{0}(\cdot,\omega),\mathbf{H}_{0}(\cdot,\omega) \bigr],W\bigr\} \in\mathcal {ID}\times\mathcal{W} $$

be a random variable.

Then there exists a random variable \(\rho(\cdot,\omega)\), \(\mathbf {u}(\cdot,\omega)\), and \(\mathbf{H}(\cdot,\omega)\),

$$\begin{aligned}& \rho(\cdot,\omega)\in\mathbf{C}\bigl([0,T];\mathbf{L}^{1}(\Omega) \bigr), \\& \mathbf{u}(\cdot,\omega)\in\mathbf{L}^{2}\bigl(0,T; \mathbf{W}_{0}^{1,2}(\Omega )\bigr), \\& \mathbf{H}(\cdot,\omega)\in\mathbf{L}^{2}\bigl(0,T; \mathbf{W}_{0}^{1,2}(\Omega)\bigr), \end{aligned}$$

satisfying the compressible magnetohydrodynamic equations (1.5)-(1.8) for a.a. \(\omega\in\mathcal{O}\).

A simple example of \(P(\rho)\) satisfying condition (1.10) is \(P(\rho)=\frac{P_{\infty}}{p}\rho^{p}+a\rho\), which is not physical example. Here \(p>\frac{3}{2}\), \(P_{\infty}\), and a are positive constants. Comparing with the existence of weak solution to the MHD equations without any stochastic effect (see e.g. [16]), the main difference of the weak solution obtained in this paper is the use of a random variable \((\rho(\cdot,\omega),\mathbf{u}(\cdot,\omega), \mathbf{H}(\cdot,\omega ))\), for a.a. \(\omega\in\mathcal{O}\).

This paper is organized as follows. In Section 2, we first recall the standard definition of weak solutions to the MHD system and give a stability results for the low regularity of the driving force, then we exploit it and show the existence of weak solutions for system (1.5)-(1.8) with a fixed (irregular) force. Finally, in the last section, we give the proof of Theorem 1.1.

2 Existence of solutions for problems with irregular forces

In this section, we give the existence of solutions for problems with irregular forces. Firstly, we prove the weak sequential stability result, which is at the heart of the proof of the main result. The stability means that a sequence of weak solutions to (1.5)-(1.8), with pre-compact data \(\rho_{0}\), \((\rho\mathbf {u})_{0}\), \(\mathbf{H}_{0}\), and W, admits a subsequence that converges weakly to another solution of the same problem. This proof is based on the work of [12, 1618]. Before giving the stability result, we rewrite the weak formulation of the problem with dropping the parameter ω and fixing a function W.

• Equation of continuity:

$$\begin{aligned}& \int_{0}^{T}\int _{\Omega}\bigl(\rho + b(\rho)\bigr)\partial_{t}\varphi+ \bigl(\rho+b(\rho )\bigr)\mathbf{u}\cdot\nabla_{x}\varphi \,dx\,dt \\& \quad = \int_{0}^{\tau}\int_{\Omega} \bigl(\bigl(b'(\rho)\rho-b(\rho)\bigr)\mathbf {div}_{x}\, \mathbf{u}\bigr)\varphi \,dx\,dt-\int_{\Omega}\bigl( \rho_{0}+b(\rho_{0})\bigr)\varphi (0,\cdot)\,dx \end{aligned}$$
(2.1)

for any test function \(\varphi\in\mathbf{C}^{\infty}_{c}([0,T)\times\bar {\Omega})\), any \(b\in\mathbf{C}^{\infty}_{c}[0,\infty)\).

• The momentum equation:

$$\begin{aligned}& \int_{0}^{T}\int _{\Omega}\bigl(\rho(\mathbf{u}-W)\cdot\partial_{t} \varphi+\rho \mathbf{u}\otimes\mathbf{u}:\nabla_{x}\varphi+P(\rho) \mathbf {div}_{x}\, \varphi\bigr) \,dx\,dt \\& \quad = \int_{0}^{\tau}\int_{\Omega} \biggl(\frac{1}{2}|\mathbf{H}|^{2}\cdot\nabla _{x} \varphi+(\mathbf{H}\cdot\nabla)\mathbf{H}\cdot\varphi\biggr)\,dx\,dt \\& \qquad {}+\int _{0}^{\tau}\int_{\Omega}\bigl(\nu \nabla_{x}\mathbf{u}\cdot\nabla_{x}\varphi +(\nu+\lambda) ( \mathbf{div}_{x}\, \mathbf{u})\cdot\nabla_{x}\varphi \\& \qquad {} +\rho(t,x,\omega)\mathbf{u}\cdot\nabla_{x}(W\varphi)\bigr) \,dx\,dt-\int_{\Omega }(\rho\mathbf{u})_{0}\cdot \varphi(0,\cdot)\,dx \end{aligned}$$
(2.2)

for all \(\varphi\in\mathbf{C}^{\infty}_{c}([0,T)\times\Omega;\mathbf{R}^{3})\).

• The induction equation of electromagnetism:

$$\begin{aligned}& \int_{0}^{T}\int _{\Omega}\bigl(\mathbf{H}\cdot\partial_{t}\varphi + \nu(\nabla _{x}\times\mathbf{H},\nabla_{x}\times\varphi) \bigr)\,dx\,dt \\& \quad =\int_{0}^{T}\int_{\Omega}( \mathbf{u}\times\mathbf{H},\nabla_{x}\times \varphi)\,dx\,dt-\int _{\Omega}\mathbf{H}_{0}\cdot\varphi(0,\cdot)\,dx \end{aligned}$$
(2.3)

for all \(\varphi\in\mathbf{C}^{\infty}_{c}([0,T)\times\Omega;\mathbf{R}^{3})\).

• The energy inequality:

$$\begin{aligned}& -\int_{0}^{T}\int _{\Omega}\biggl(\frac{1}{2}\rho|\mathbf{u}-W|^{2} + P(\rho)+\frac {1}{2}|\mathbf{H}|^{2}\biggr)\,dx\,dt \\& \qquad {}+ \int_{0}^{T}\int_{\Omega} \bigl(\nu|\nabla_{x}\mathbf{u}|^{2}+(\nu+\lambda)|\mathbf {div}_{x}\, \mathbf{u}|^{2}+\nu|\nabla_{x}\times \mathbf{H}|^{2}\bigr)\,dx \\& \quad \leq \varphi(0)\int_{\Omega}\biggl(\frac{1}{2} \frac{|(\rho\mathbf {u})_{0}|^{2}}{\rho_{0}}+P(\rho_{0})+\frac{1}{2}|\mathbf{H}_{0}|^{2} \biggr)\,dx \\& \qquad {}-\int_{0}^{\tau}\int_{\Omega} \biggl(\rho\mathbf{u}\otimes\mathbf{u}:\nabla _{x}W+P(\rho) \mathbf{div}_{x}\, W-\frac{1}{2}\rho\mathbf{u}\cdot\nabla _{x}|W|^{2}\biggr)\,dx\,dt \\& \qquad {} -\int_{0}^{\tau}\int_{\Omega} \bigl(\nu\nabla_{x}\mathbf{u}\cdot\nabla_{x}W+(\nu +\lambda) \mathbf{div}_{x}\, \mathbf{u}\cdot\nabla_{x}W\bigr)\,dx\,dt \\& \qquad {} +\int_{0}^{\tau}\int_{\Omega} \biggl(\frac{1}{2}|\mathbf{H}|^{2}\cdot\nabla _{x}W+( \mathbf{H}\cdot\nabla)\mathbf{H}\cdot W\biggr)\,dx\,dt \\& \qquad {}+\int_{0}^{\tau}\int_{\Omega}( \mathbf{u}\times\mathbf{H},\nabla_{x}\times W)\,dx\,dt \end{aligned}$$
(2.4)

for any \(\varphi\in\mathbf{C}^{\infty}_{c}[0,T)\), \(\varphi\geq0\).

Proposition 2.1

Let \(\Omega\subset\mathbf{R}^{3}\) be a bounded Lipschitz domain. Assume that the pressure P satisfies (1.10). Let \(\{W_{n}\}_{n=1}^{\infty }\) be a sequence of functions which satisfies

$$ \operatorname{ess}\sup_{t\in(0,T)}\bigl\Vert W_{n}(t,\cdot)\bigr\Vert _{\mathbf{W}_{0}^{1,\infty}(\Omega ;\mathbf{R}^{3})}\leq C_{W},\quad W_{n}\rightarrow W\textit{ in }\mathbf{L}^{1}\bigl(0,T; \mathbf{W}_{0}^{1,1}\bigl(\Omega;\mathbf{R}^{3} \bigr)\bigr). $$
(2.5)

Let \(\{\rho_{n},\mathbf{u}_{n},\mathbf{H}_{n}\}\) be a sequence of weak solutions of the MHD system driven by \(W_{n}\) in \((0,T)\times\Omega\), with the initial data \(\rho_{0,n}\), \((\rho\mathbf{H})_{0,n}\), \(\mathbf {H}_{0,n}\) such that

$$\begin{aligned}& \rho_{0,n}\geq0,\qquad \int_{\Omega}\rho_{0,n} \,dx=M>0,\qquad \|\rho_{0,n}\| _{\mathbf{L}^{p}(\Omega)}\leq E_{0}, \qquad \int_{\Omega}\frac{|(\rho\mathbf {u})_{0,n}|^{2}}{\rho_{0,n}}\,dx\leq E_{0}, \\& \rho_{0,n}\rightarrow\rho_{0}\quad \textit{in } \mathbf{L}^{p}(\Omega),\qquad (\rho \mathbf{u})_{0,n} \rightarrow(\rho\mathbf{u})_{0}\quad \textit{weakly in }\mathbf {L}^{1}(\Omega), \\& \int_{\Omega}\frac{|(\rho\mathbf{u})_{0,n}|^{2}}{\rho _{0,n}}\,dx\rightarrow\int _{\Omega}\frac{|(\rho\mathbf {u})_{0}|^{2}}{\rho_{0}}\,dx. \end{aligned}$$

Then, at least for suitable subsequences,

$$\begin{aligned}& \rho_{n}\rightarrow\rho\quad \textit{in }\mathbf{C}_{\mathrm{weak}} \bigl(0,T;\mathbf {L}^{p}(\Omega)\bigr)\textit{ and in } \mathbf{L}^{1}\bigl((0,T)\times\Omega\bigr), \\& \mathbf{u}_{n}\rightarrow\mathbf{u}\quad \textit{weakly in }\mathbf {L}^{2}\bigl(0,T;\mathbf{W}_{0}^{1,2}\bigl(\Omega, \mathbf{R}^{3}\bigr)\bigr), \\& \mathbf{H}_{n}\rightarrow\mathbf{H}\quad \textit{weakly in }\mathbf {L}^{2}\bigl(0,T;\mathbf{W}_{0}^{1,2}\bigl(\Omega, \mathbf{R}^{3}\bigr)\bigr), \\& \rho_{n}(\mathbf{u}_{n}-W_{n})\rightarrow\rho( \mathbf{u}-W) \quad \textit{in }\mathbf {C}_{\mathrm{weak}}\bigl([0,T]; \mathbf{L}^{\frac{2p}{p+1}}\bigl(\Omega;\mathbf{R}^{3}\bigr)\bigr), \end{aligned}$$

where ρ, u, and H is a weak solution of the same problem driven by the force W, and having the initial data \(\rho _{0}\), \((\rho\mathbf{u})_{0}\), and \(\mathbf{H}_{0}\).

Proof

We first derive some uniform bounds estimates, which are exactly the same as the a priori bounded estimates available for the compressible MHD, cf. [16], and the stochastic compressible Navier-Stokes expression (see [12]). Applying the standard Gronwall inequality to the energy inequality (2.4), we get

$$\begin{aligned}& \operatorname{ess}\sup_{t\in(0,T)}\bigl\Vert \rho_{n}(t,\cdot)\bigr\Vert _{\mathbf{L}^{p}(\Omega)}\leq c(M,C_{w},E_{0},T), \end{aligned}$$
(2.6)
$$\begin{aligned}& \operatorname{ess}\sup_{t\in(0,T)}\bigl\Vert \sqrt{ \rho_{n}}\mathbf{u}_{n}(t,\cdot)\bigr\Vert _{\mathbf {L}^{2}(\Omega;\mathbf{R}^{3})} \leq c(M,C_{w},E_{0},T), \end{aligned}$$
(2.7)
$$\begin{aligned}& \operatorname{ess}\sup_{t\in(0,T)}\bigl\Vert \mathbf{H}_{n}(t,\cdot)\bigr\Vert _{\mathbf{L}^{2}(\Omega ;\mathbf{R}^{3})}\leq c(M,C_{w},E_{0},T), \end{aligned}$$
(2.8)

where M is the total mass, which is a constant. Moreover, one has

$$ \bigl\Vert \rho_{n}(t,\cdot)\bigr\Vert _{\mathbf{L}^{1}(\Omega)}=M \quad \forall t\in[0,T]. $$

Moreover, by assumption (2.5),

$$ \operatorname{ess}\sup_{t\in(0,T)}\int_{\Omega} \rho|W_{n}|^{2}\,dx\leq c(M,C_{w}). $$

We use Korn’s inequality,

$$\begin{aligned} \int_{0}^{T}\bigl\Vert \mathbf{u}_{n}(t, \cdot)\bigr\Vert ^{2}_{\mathbf{W}_{0}^{1,2}(\Omega;\mathbf {R}^{3})}\,dt\leq c(M,C_{w},E_{0},T), \\ \int_{0}^{T}\bigl\Vert \mathbf{H}_{n}(t, \cdot)\bigr\Vert ^{2}_{\mathbf{W}_{0}^{1,2}(\Omega;\mathbf {R}^{3})}\,dt\leq c(M,C_{w},E_{0},T). \end{aligned}$$

Consequently, there exists a suitable subsequence such that

$$\begin{aligned}& \rho_{n}\rightarrow\rho \quad \mbox{in } \mathbf{C}_{\mathrm{weak}}\bigl(0,T;\mathbf {L}^{p}(\Omega)\bigr), \end{aligned}$$
(2.9)
$$\begin{aligned}& \mathbf{u}_{n}\rightarrow\mathbf{u}\quad \mbox{weakly in }\mathbf {L}^{2}\bigl(0,T;\mathbf{W}_{0}^{1,2} \bigl(\Omega,\mathbf{R}^{3}\bigr)\bigr), \end{aligned}$$
(2.10)
$$\begin{aligned}& \mathbf{H}_{n}\rightarrow\mathbf{H}\quad \mbox{weakly in }\mathbf {L}^{2}\bigl(0,T;\mathbf{W}_{0}^{1,2} \bigl(\Omega,\mathbf{R}^{3}\bigr)\bigr), \end{aligned}$$
(2.11)
$$\begin{aligned}& \rho_{n}\mathbf{u}_{n}\rightarrow\rho \mathbf{u} \quad \mbox{weakly }(*)\mbox{ in } \mathbf{L}^{\infty}\bigl(0,T; \mathbf{L}^{\frac{2p}{p+1}}\bigl(\Omega,\mathbf{R}^{3}\bigr)\bigr). \end{aligned}$$
(2.12)

By the momentum equation (2.2), we derive

$$ \rho_{n}(\mathbf{u}_{n}-W_{n})\rightarrow\rho( \mathbf{u}-W)\quad \mbox{in }\mathbf {C}_{\mathrm{weak}}\bigl([0,T]; \mathbf{L}^{\frac{2p}{p+1}}\bigl(\Omega;\mathbf{R}^{3}\bigr)\bigr), $$

which combining with (2.5) and (2.12) gives

$$ \begin{aligned} &(\rho_{n}\mathbf{u}_{n}) (t,\cdot) \rightarrow(\rho\mathbf{u}) (t,\cdot)\quad \mbox{weakly in }\mathbf{L}^{\frac{2p}{p+1}} \bigl(\Omega;\mathbf{R}^{3}\bigr)\mbox{ for a.a. }t\in (0,T), \\ &(\rho_{n}\mathbf{u}_{n}) (t,\cdot) \rightarrow(\rho\mathbf{u}) (t,\cdot)\quad \mbox{in }\mathbf{L}^{q} \bigl(0,T;\mathbf{W}^{-1,2}\bigl(\Omega;\mathbf{R}^{3}\bigr) \bigr), \forall 1\leq q< \infty. \end{aligned} $$
(2.13)

By (2.6)-(2.7), (2.10), and (2.13), we have

$$ \rho_{n}\mathbf{u}_{n}\otimes \mathbf{u}_{n}\rightarrow\rho\mathbf {u}\otimes\mathbf{u}\quad \mbox{weakly in }\mathbf{L}^{2}\bigl(0,T;\mathbf{L}^{\frac {6p}{4p+3}}\bigl( \Omega;\mathbf{R}^{3\times3}\bigr)\bigr). $$
(2.14)

Using (2.10)-(2.11), we get

$$\begin{aligned}& (\nabla\mathbf{H}_{n}\times\mathbf{H}_{n}) \times\mathbf {H}_{n}\rightarrow(\nabla\mathbf{H}\times\mathbf{H}) \times\mathbf{H}\quad \mbox{in }\mathcal{D}\bigl((0,T)\times\Omega\bigr), \end{aligned}$$
(2.15)
$$\begin{aligned}& \nabla\times(\mathbf{u}_{n}\times\mathbf{H}_{n}) \rightarrow\nabla \times(\mathbf{u}\times\mathbf{H})\quad \mbox{in }\mathcal{D} \bigl((0,T)\times\Omega\bigr), \end{aligned}$$
(2.16)

where we use \(\mathcal{D}'\) to denote the sense of distributions.

Next, we derive the pressure estimate and the pointwise convergence of the densities. We denote the Bogovskii operator as \(\mathcal{B}\approx \mathbf{div}^{-1}\), and we choose the test function

$$ \varphi=\psi(t)\mathcal{B}\biggl[b(\rho_{n})-\frac{1}{|\Omega|}\int _{\Omega }b(\rho_{n})\,dx\biggr], \quad \psi\in \mathbf{C}_{c}^{\infty}(0,T), b(\rho)\approx\rho ^{\kappa}. $$

Following [12] and [16] step by step, we deduce that

$$ \int_{0}^{T}\int_{\Omega}P( \rho_{n})\rho_{n}^{\beta}\leq c(M,C_{w},E_{0},T) $$

for a certain positive constant β which depends on M, \(C_{w}\), and \(E_{0}\).

To derive the pointwise convergence of the densities

$$ \rho_{n}\rightarrow\rho\quad \mbox{a.a. in }(0,T)\times \Omega, $$
(2.17)

we take the following test function:

$$ \varphi=\psi\eta(x)\nabla_{x}\triangle^{-1}\bigl[\xi \overline{b(\rho)}\bigr] $$

in the limit equation

$$\begin{aligned}& \int_{0}^{T}\int _{\Omega}\bigl(\rho(\mathbf{u}-W)\cdot\partial_{t} \varphi+\rho \mathbf{u}\otimes\mathbf{u}:\nabla_{x}\varphi+\overline{P( \rho)}\mathbf {div}_{x}\, \varphi\bigr) \,dx\,dt \\& \quad = \int_{0}^{\tau}\int_{\Omega} \biggl(\frac{1}{2}|\mathbf{H}|^{2}\cdot\nabla _{x} \varphi+(\mathbf{H}\cdot\nabla)\mathbf{H}\cdot\varphi\biggr)\,dx\,dt \\& \qquad {}+\int _{0}^{\tau}\int_{\Omega}\bigl(\nu \nabla_{x}\mathbf{u}\cdot\nabla_{x}\varphi +(\nu+\lambda) ( \mathbf{div}_{x}\, \mathbf{u})\cdot\nabla_{x}\varphi \\& \qquad {} +\rho\mathbf{u}\cdot\nabla_{x}(W\varphi)\bigr)\,dx\,dt-\int _{\Omega}(\rho\mathbf {u})_{0}\cdot\varphi(0,\cdot) \,dx, \end{aligned}$$
(2.18)

where \(\overline{b(\rho)}\) or \(\overline{P(\rho)}\) denotes a weak limit of the compositions \(\overline{b(\rho_{n})}\) or \(\overline{P(\rho_{n})}\).

Using the same arguments as in [12, 19], we have

$$\begin{aligned}& \lim_{n\rightarrow\infty}\int_{0}^{T} \int_{\mathbf{R}^{3}}\psi\mathbf {u}_{n}\cdot\bigl( \nabla_{x}\triangle^{-1}\nabla_{x}\bigl[\xi b( \rho_{n})\bigr]\cdot\eta\rho _{n}\mathbf{u}_{n} - \xi b(\rho_{n})\nabla_{x}\triangle^{-1} \nabla_{x}[\eta\rho_{n}\mathbf {u}_{n}]\bigr)\,dx \,dt \\& \quad = \int_{0}^{T}\int_{\mathbf{R}^{3}} \psi\mathbf{u}\cdot\bigl(\nabla_{x}\triangle ^{-1} \nabla_{x}\bigl[\xi b(\rho)\bigr]\cdot\eta\rho\mathbf{u} -\xi \overline{b(\rho)}\nabla_{x}\triangle^{-1} \nabla_{x}[\eta\rho\mathbf{u}]\bigr)\,dx\,dt. \end{aligned}$$
(2.19)

By (2.5)-(2.9), we derive

$$ \rho_{n}W_{n}\rightarrow\rho W\quad \mbox{in } \mathbf{L}^{2}\bigl(0,T;\mathbf {W}^{-1,q'}\bigl(\Omega; \mathbf{R}^{3}\bigr)\bigr), $$

which implies that

$$\begin{aligned}& \lim_{n\rightarrow\infty}\int_{0}^{T} \int_{\mathbf{R}^{3}}\biggl[\psi\rho_{n}\mathbf {u}_{n} \cdot\nabla_{x}\bigl(\eta W_{n}\cdot\nabla_{x} \triangle^{-1}\bigl[\xi b(\rho_{n})\bigr] \bigr) \\& \qquad {}+ \rho_{n} W_{n}\cdot\partial_{t}\bigl(\psi\eta \nabla_{x}\triangle^{-1}\bigl[\xi b(\rho _{n}) \bigr]\bigr)\biggr]\,dx\,dt \\& \quad =\lim_{n\rightarrow\infty}\int_{0}^{T} \int_{\mathbf{R}^{3}}\biggl[\psi\rho\mathbf {u}\cdot\nabla_{x}(\eta W \cdot\nabla_{x}\triangle^{-1}\bigl[\xi\overline{b(\rho)} \bigr]\bigr) \\& \qquad {}+\rho W\cdot\partial_{t}\bigl(\psi\eta\nabla_{x} \triangle^{-1}\bigl[\xi\overline {b(\rho)}\bigr]\bigr)\biggr]\,dx\,dt. \end{aligned}$$
(2.20)

By (2.19)-(2.20), we obtain

$$\begin{aligned}& \lim_{n\rightarrow\infty}\int_{0}^{T}\int _{\mathbf{R}^{3}}\xi\eta\bigl(\xi P(\rho _{n})b( \rho_{n})-\bigl(\nu\nabla_{x}\mathbf{u}_{n}+(\nu+ \lambda) (\mathbf {div}_{x}\, \mathbf{u}_{n})\bigr)\cdot \nabla_{x}\triangle^{-1}\nabla_{x}\bigl[\xi b(\rho _{n})\bigr]\bigr)\,dx\,dt \\& \quad = \int_{0}^{T}\int_{\mathbf{R}^{3}} \xi\eta\bigl(\xi\overline{P(\rho)} \overline {b(\rho)}-\bigl(\nu \nabla_{x}\mathbf{u}+(\nu+\lambda) (\mathbf{div}_{x}\, \mathbf {u})\bigr)\cdot\nabla_{x}\triangle^{-1} \nabla_{x}\bigl[\xi\overline{b(\rho)}\bigr]\bigr)\,dx\,dt. \end{aligned}$$

In particular, we have the so-called weak continuity of the effective viscous flux

$$ \overline{P(\rho)b(\rho)}-\biggl(\frac{4}{3}\nu+\lambda\biggr) \overline{b(\rho)} (\mathbf{div}_{x}\, \mathbf{u}) = \overline{P(\rho)} \overline{b(\rho)}-\biggl(\frac{4}{3}\nu+\lambda\biggr) \overline {b(\rho)}\mathbf{div}_{x}\, \mathbf{u}, $$

which first was established in [20] and recovered in [12] for stochastic compressible Navier-Stokes system.

Now the rest of deriving (2.21) is by the same arguments as in [12, 19]. Moreover, for any \(1\leq q\leq p\), we can obtain

$$ \rho_{n}\rightarrow\rho \quad \mbox{in }\mathbf{C} \bigl([0,T];\mathbf{L}^{q}(\Omega)\bigr). $$
(2.21)

 □

Assume that W is regular and satisfies

$$ W\in\mathbf{C}^{1}\bigl([0,T];W^{1,\infty}\bigl(\Omega; \mathbf{R}^{3}\bigr)\bigr),\qquad W(0,\cdot)=0. $$

As shown in [12, 16], by a small modification of the proof for the compressible MHD system and using Proposition 2.1, we obtain the following result.

Proposition 2.2

Let \(\Omega\subset\mathbf{R}^{3}\) be a bounded Lipschitz domain. Assume that the pressure P satisfies (1.10),

$$ W\in\mathbf{L}^{\infty}_{\mathrm{weak}}\bigl(0,T;W^{1,\infty}_{0} \bigl(\Omega;\mathbf{R}^{3}\bigr)\bigr) $$

and

$$\begin{aligned}& \rho_{0}\in\mathbf{L}^{p}(\Omega),\qquad \rho_{0}\geq0, \qquad \int_{\Omega} \rho_{0}=M>0, \\& (\rho\mathbf{u})_{0}\in\mathbf{L}^{1}\bigl(\Omega, \mathbf{R}^{3}\bigr),\qquad \int_{\Omega } \frac{|\rho\mathbf{u}_{0}|^{2}}{\rho_{0}}\,dx< \infty, \qquad \mathbf{H}_{0}\in\mathbf {L}^{2}(\Omega). \end{aligned}$$

Then the compressible MHD system admits a weak solution in \((0,T)\times \Omega\) in the sense specified in (2.1)-(2.4).

Remark 2.1

We can use a sequence of approximation smooth function \(W_{n}\) for any regular W to obtain solutions for the problem driven by any regular W. The symbol \(\mathbf{L}^{\infty}_{\mathrm{weak}}\) denotes the space of weakly measurable functions.

3 Proof of Theorem 1.1

Assume that

$$ \omega\in\mathcal{O}\mapsto W(\cdot,\omega)\in\mathcal{W} $$

is a random variable, and the mapping

$$ \omega\in\mathcal{O}\mapsto\bigl[\rho_{0}(\cdot,\omega),(\rho\mathbf {u})_{0}(\cdot,\omega),\mathbf{H}_{0}(\cdot,\omega)\bigr] \in\mathcal{ID} $$

is a random variable on \(\mathcal{O}=\{\mathcal{O},\mathcal{B},\mu\}\), with a regular probability measure μ. Then for any \(\varepsilon>0\), there exists a compact set \(\mathcal {K}_{\varepsilon\subset\mathcal{O}}\) such that

$$ \mu(\mathcal{O}\backslash\mathcal{K}_{\varepsilon})< \varepsilon, $$

and the mapping

$$ \omega\mapsto\bigl[\rho_{0}(\cdot,\omega),(\rho\mathbf{u})_{0}( \cdot,\omega ),\mathbf{H}_{0}(\cdot,\omega)\bigr] $$

restricted to \(\mathcal{K}_{\varepsilon}\) is continuous.

Consider a multi-valued mapping

$$\begin{aligned} \mathcal{L}:\bigl\{ \bigl[\rho_{0},(\rho\mathbf{u})_{0}, \mathbf{H}_{0}\bigr],W\bigr\} \in&\mathcal {ID}\times\mathcal{W}\mapsto[ \rho,\mathbf{u},\mathbf{H}] \\ \in& \mathbf{C}\bigl([0,T];\mathbf{L}^{1}( \Omega)\bigr)\times\mathbf{L}^{2}\bigl(0,T;\mathbf {W}_{0}^{1,2} \bigl(\Omega,\mathbf{R}^{3}\bigr)\bigr). \end{aligned}$$

The existence of \(\mathcal{L}\) is guaranteed by Proposition 2.2. Proposition 2.1 implies that the value of \(\mathcal{L}\) are non-empty closed subsets of the Banach space \(\mathbf{C}([0,T];\mathbf{L}^{1}(\Omega ))\times\mathbf{L}^{2}(0,T;\mathbf{W}_{0}^{1,2}(\Omega; \mathbf{R}^{3}))\), and \(\mathcal{L}\) has a closed graph.

Finally, using the abstract result of Theorem 3.1 and Lemma 3.1 in [21], we conclude that there exists a weak solution of the MHD driven by the stochastic external force W in the sense of (1.5)-(1.8) such that the mapping

$$ \omega\mapsto\bigl[\rho(\cdot,\omega),\mathbf{u}(\cdot,\omega),\mathbf {H}( \cdot,\omega)\bigr]\in\mathbf{C}\bigl([0,T];\mathbf{L}^{1}(\Omega)\bigr) \times\mathbf {L}^{2}\bigl(0,T;\mathbf{W}_{0}^{1,2} \bigl(\Omega;\mathbf{R}^{3}\bigr)\bigr) $$

is a random variable, for \(\omega\in\mathcal{O}\). This completes the proof of Theorem 1.1.

References

  1. Kulikovskiy, AG, Lifshitz, EM: Magnetohydrodynamics. Addison-Wesley, Reading (1965)

    Google Scholar 

  2. Laudau, LD, Lifshitz, EM: Electrodynamics of Continuous Media, 2nd edn. Pergamon, New York (1984)

    Google Scholar 

  3. Polovin, RV, Demutskii, VP: Fundamentals of Magnetohydrodynamics. Consultants Bureau, New York (1990)

    Google Scholar 

  4. Bensoussan, A: Stochastic Navier-Stokes equations. Acta Appl. Math. 38, 268-304 (1995)

    Article  MathSciNet  Google Scholar 

  5. Flandoli, F: An Introduction to 3D Stochastic Fluid Dynamics. CIME Lecture Notes (2005)

    Google Scholar 

  6. Mattingly, J: On recent progress for the stochastic Navier-Stokes equations with non-Lipschitz random forcing. Rev. Math. Phys. 22, 669-697 (2010)

    Article  MathSciNet  Google Scholar 

  7. Flandoli, F, Gatarek, D: Martingale and stationary solutions for stochastic Navier-Stokes equations. Probab. Theory Relat. Fields 102, 367-391 (1995)

    Article  MATH  MathSciNet  Google Scholar 

  8. Flandoli, F, Romito, M: Markov selections for 3D stochastic Navier-Stokes equation. Probab. Theory Relat. Fields 140, 407-458 (2008)

    Article  MATH  MathSciNet  Google Scholar 

  9. Goldys, B, Röckner, M, Zhang, XC: Martingale solutions and Markov selections for stochastic partial differential equations. Stoch. Process. Appl. 119, 1725-1764 (2009)

    Article  MATH  Google Scholar 

  10. Tornatore, E, Fujita Yashima, H: One-dimensional equations of a barotropic viscous gas with a not very regular perturbation. Ann. Univ. Ferrara, Sez. VII 40, 137-168 (1994)

    MathSciNet  Google Scholar 

  11. Tornatore, E, Fujita Yashima, H: One-dimensional stochastic equations for a viscous barotropic gas. Ric. Mat. 46(2), 255-283 (1997)

    MATH  MathSciNet  Google Scholar 

  12. Feireisl, E, Maslowski, B, Novotny, A: Compressible fluid flows driven by stochastic forcing. J. Differ. Equ. 254, 1342-1358 (2013)

    Article  MATH  MathSciNet  Google Scholar 

  13. Sritharan, SS, Sundar, P: The stochastic magneto-hydrodynamic system. Infin. Dimens. Anal. Quantum Probab. Relat. Top. 2, 241-265 (1999)

    Article  MATH  MathSciNet  Google Scholar 

  14. Sunda, P: Stochastic magnetohydrodynamic system perturbed by general noise. Commun. Stoch. Anal. 4, 253-269 (2010)

    MathSciNet  Google Scholar 

  15. Diperna, R-J, Lions, P-L: Ordinary differential equations, transport theory and Sobolev spaces. Invent. Math. 98, 511-547 (1989)

    Article  MATH  MathSciNet  Google Scholar 

  16. Hu, X, Wang, D: Global existence and large time behavior of solutions to the three dimensional equations of compressible magnetohydrodynamic flows. Arch. Ration. Mech. Anal. 197, 203-238 (2010)

    Article  MATH  MathSciNet  Google Scholar 

  17. Coifman, R, Meyer, Y: On commutators of singular integrals and bilinear singular integrals. Trans. Am. Math. Soc. 212, 315-331 (1975)

    Article  MATH  MathSciNet  Google Scholar 

  18. Yan, WP: Existence of weak solutions to the three-dimensional density-dependent generalized incompressible magnetohydrodynamic flows. Discrete Contin. Dyn. Syst., Ser. A 35(3), 1359-1385 (2015)

    Article  MATH  Google Scholar 

  19. Feireisl, E, Novotny, A, Petzeltova, H: On the existence of globally defined weak solutions to the Navier-Stokes equations of compressible isentropic fluids. J. Math. Fluid Mech. 3, 358-392 (2001)

    Article  MATH  MathSciNet  Google Scholar 

  20. Lions, P-L: Mathematical Topics in Fluid Dynamics. Volume 2: Compressible Models. Oxford University Press, Oxford (1998)

    Google Scholar 

  21. Bensoussan, A, Temam, R: Équations stochastiques du type Navier-Stokes. J. Funct. Anal. 13, 195-222 (1973)

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors express their sincere thanks to the anonymous referees for the comments.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Junling Sun.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sun, J., Li, H. Existence of the weak solutions for the compressible magnetohydrodynamic flows driven by stochastic forcing. Bound Value Probl 2015, 74 (2015). https://doi.org/10.1186/s13661-015-0335-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0335-3

Keywords