In this paper,we consider a model of compressible isentropic two-fluid magneto-hydrodynamics without resistivity in a strip domain in three dimensional space.By exploiting the two-tier energy method developed in[Anal ...In this paper,we consider a model of compressible isentropic two-fluid magneto-hydrodynamics without resistivity in a strip domain in three dimensional space.By exploiting the two-tier energy method developed in[Anal PDE,2013,6:1429–1533],we prove the global well-posedness of the governing model around a uniform magnetic field which is non-parallel to the horizontal boundary.Moreover,we show that the solution converges to the steady state at an almost exponential rate as time goes to infinity.Compared to the work of Tan and Wang[SIAM J Math Anal,2018,50:1432–1470],we need to overcome the difficulties caused by particles.展开更多
We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscal...We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.展开更多
We report on the magnetohydrodynamic impact on the axisymmetric flow of Al_(2)O_(3)/Cu nanoparticles suspended in H_(2)O past a stretched/shrinked sheet.With the use of partial differential equations and the correspon...We report on the magnetohydrodynamic impact on the axisymmetric flow of Al_(2)O_(3)/Cu nanoparticles suspended in H_(2)O past a stretched/shrinked sheet.With the use of partial differential equations and the corresponding thermophysical characteristics of nanoparticles,the physical flow process is illustrated.The resultant nonlinear system of partial differential equations is converted into a system of ordinary differential equations using the suitable similarity transformations.The transformed differential equations are solved analytically.Impacts of the magnetic parameter,solid volume fraction and stretching/shrinking parameter on momentum and temperature distribution have been analyzed and interpreted graphically.The skin friction and Nusselt number were also evaluated.In addition,existence of dual solution was deduced for the shrinking sheet and unique solution for the stretching one.Further,Al_(2)O_(3)/H_(2)O nanofluid flow has better thermal conductivity on comparing with Cu/H_(2)O nanofluid.Furthermore,it was found that the first solutions of the stream are stable and physically realizable,whereas those of the second ones are unstable.展开更多
Observations of transmission spectra reveal that hot Jupiters and Neptunes are likely to possess escaping atmospheres driven by stellar radiation.Numerous models predict that magnetic fields may exert significant infl...Observations of transmission spectra reveal that hot Jupiters and Neptunes are likely to possess escaping atmospheres driven by stellar radiation.Numerous models predict that magnetic fields may exert significant influences on the atmospheres of hot planets.Generally,the escaping atmospheres are not entirely ionized,and magnetic fields only directly affect the escape of ionized components within them.Considering the chemical reactions between ionized components and neutral atoms,as well as collision processes,magnetic fields indirectly impact the escape of neutral atoms,thereby influencing the detection signals of planetary atmospheres in transmission spectra.In order to simulate this process,we developed a magnetohydrodynamic multi-fluid model based on MHD code PLUTO.As an initial exploration,we investigated the impact of magnetic fields on the decoupling of H^(+)and H in the escaping atmosphere of the hot Neptune GJ436b.Due to the strong resonant interactions between H and H^(+),the coupling between them is tight even if the magnetic field is strong.Of course,alternatively,our work also suggests that merging H and H^(+)into a single flow can be a reasonable assumption in MHD simulations of escaping atmospheres.However,our simulation results indicate that under the influence of magnetic fields,there are noticeable regional differences in the decoupling of H^(+)and H.With the increase of magnetic field strength,the degree of decoupling also increases.For heavier particles such as O,the decoupling between O and H^(+)is more pronounced.Our findings provide important insights for future studies on the decoupling processes of heavy atoms in the escaping atmospheres of hot Jupiters and hot Neptunes under the influence of magnetic fields.展开更多
The paper aims to estimate the thickness of the boundary layer for the planar MHD system with vanishing shear viscosity μ. Under some conditions on the initial and boundary data, we show that the thickness is of the ...The paper aims to estimate the thickness of the boundary layer for the planar MHD system with vanishing shear viscosity μ. Under some conditions on the initial and boundary data, we show that the thickness is of the order √μ|lnμ|. Note that this estimate holds also for the Navier-Stokes system so that it extends the previous works even without the magnetic effect.展开更多
We consider a nonlinear hyperbolic system of two conservation laws which arises in ideal magnetohydrodynamics and includes second-order terms accounting for magnetic resistivity and Hall effect. We show that the initi...We consider a nonlinear hyperbolic system of two conservation laws which arises in ideal magnetohydrodynamics and includes second-order terms accounting for magnetic resistivity and Hall effect. We show that the initial value problem for this model may lead to solutions exhibiting complex wave structures, including undercompressive nonclassical shock waves. We investigate numerically the subtle competition that takes place between the hyperbolic, diffusive, and dispersive parts of the system. Following Abeyratne, Knowles, LeFloch, and Truskinovsky, who studied similar questions arising in fluid and solid flows, we determine the associated kinetic function which characterizes the dynamics of undereompressive shocks driven by resistivity and Hall effect. To this end, we design a new class of "schemes with eontroled dissipation", following recent work by LeFloch and Mohammadian. It is now recognized that the equivalent equation associated with a scheme provides a guideline to design schemes that capture physically relevant, nonclassical shocks. We propose a new class of schemes based on high-order entropy conservative, finite differences for the hyperbolic flux, and high-order central differences for the resistivity and Hall terms. These schemes are tested for several regimes of (co-planar or not) initial data and parameter values, and allow us to analyze the properties of nonclassical shocks and establish the existence of monotone kinetic functions in magnetohydrodynamics.展开更多
We conduct an electron magnetohydrodynamics magnetic reconnection experiment with guide-field in our Keda linear magnetized plasma device, in which two pulsed currents with the same direction are conducted in parallel...We conduct an electron magnetohydrodynamics magnetic reconnection experiment with guide-field in our Keda linear magnetized plasma device, in which two pulsed currents with the same direction are conducted in parallel with the axial direction of the main chamber of the device using two long aluminum sticks. After approximately 5μs, an X-type magnetic field line topology is formed at the center of the chamber. With the formation of the X-type topology of magnetic field lines, we can also find the rapid increase of the current and ratio of the common flux to the private flux in this area. Additionally, a reduction in the plasma density and the plasma density concentration along one pair of separatrices can also be found.展开更多
The interaction between a converging cylindrical shock and double density interfaces in the presence of a saddle magnetic field is numerically investigated within the framework of ideal magnetohydrodynamics.Three flui...The interaction between a converging cylindrical shock and double density interfaces in the presence of a saddle magnetic field is numerically investigated within the framework of ideal magnetohydrodynamics.Three fluids of differing densities are initially separated by the two perturbed cylindrical interfaces.The initial incident converging shock is generated from a Riemann problem upstream of the first interface.The effect of the magnetic field on the instabilities is studied through varying the field strength.It shows that the Richtmyer-Meshkov and Rayleigh-Taylor instabilities are mitigated by the field,however,the extent of the suppression varies on the interface which leads to non-axisymmetric growth of the perturbations.The degree of asymmetry of the interfacial growth rate is increased when the seed field strength is increased.展开更多
Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phas...Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phases:find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution,and then locally solve linearized residual subproblems by one of three iterations(Stokes-type,Newton,and Oseen-type)on subdomains with fine grid in parallel to approximate the high frequency component.Optimal error estimates with regard to two mesh sizes and iterative steps of the proposed algorithms are given.Some numerical examples are implemented to verify the algorithm.展开更多
We investigate the local existence of smooth solutions of a 3D ideal magneto-hydrodynamics (MHD) equations in a bounded domain and give a blow-up criteria to thisequations with respect to vorticists.
A nonlinear Galerkin mixed element (NGME) method for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of the NGME solution are derived.
The Galerkin-Petrov least squares method is combined with the mixed finite element method to deal with the stationary, incompressible magnetohydrodynamics system of equations with viscosity. A Galerkin-Petrov least sq...The Galerkin-Petrov least squares method is combined with the mixed finite element method to deal with the stationary, incompressible magnetohydrodynamics system of equations with viscosity. A Galerkin-Petrov least squares mixed finite element format for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of its solution are derived. Through this method, the combination among the mixed finite element spaces does not demand the discrete Babuska-Brezzi stability conditions so that the mixed finite element spaces could be chosen arbitrartily and the error estimates with optimal order could be obtained.展开更多
Lattice Boltzmann Method is recently developed within numerical schemes for simulating a variety of physical systems. In this paper a new lattice Bhatnagar-Gross-Krook (LBGK) model for two-dimensional incompressible m...Lattice Boltzmann Method is recently developed within numerical schemes for simulating a variety of physical systems. In this paper a new lattice Bhatnagar-Gross-Krook (LBGK) model for two-dimensional incompressible magnetohydrodynamics (IMHD) is presented. The model is an extension of a hydrodynamics lattice BGK model with 9 velocities on a square lattice, resulting in a model with 17 velocities. Most of the existing LBGK models for MHD can be viewed as compressible schemes to simulate incompressible flows. The compressible effect might lead to some undesirable errors in numerical simulations. In our model the compressible effect has been overcome successfully. The model is then applied to the Hartmann flow, giving reasonable results.展开更多
The gas-kinetic theory based flux splitting method has been successfully proposed for solving one-and two-dimensional ideal magnetohydrodynamics by Xu et al. [J.Comput.Phys.,1999;2000],respectively.This paper extends ...The gas-kinetic theory based flux splitting method has been successfully proposed for solving one-and two-dimensional ideal magnetohydrodynamics by Xu et al. [J.Comput.Phys.,1999;2000],respectively.This paper extends the kinetic method to solve three-dimensional ideal magnetohydrodynamics equations,where an adaptive parameter 17 is used to control the numerical dissipation in the flux splitting method. Several numerical examples are given to demonstrate that the proposed method can achieve high numerical accuracy and resolve strong discontinuous waves in three dimensional ideal MHD problems.展开更多
The solutions of incompressible ideal Hall magnetohydrodynamics are obtained by using the traveling wave method. It is shown that the velocity and magnetic field parallel to the wave vector can be arbitrary constants....The solutions of incompressible ideal Hall magnetohydrodynamics are obtained by using the traveling wave method. It is shown that the velocity and magnetic field parallel to the wave vector can be arbitrary constants. The velocity and magnetic field perpendicular to the wave vector are both helical waves. Moreover, the amplitude of the velocity perpendicular to the wave vector is related to the wave number and the circular frequency. In addition, further studies indicate that, no matter whether the uniform ambient magnetic field exists or not, the forms of the travelling wave solutions do not change.展开更多
Using the Fourier helical decomposition, we obtain the absolute statistical equilibrium spectra of left- and right-handed helical modes in the incompressible ideal Hall magnetohydrodynamics (MHD). It is shown that the...Using the Fourier helical decomposition, we obtain the absolute statistical equilibrium spectra of left- and right-handed helical modes in the incompressible ideal Hall magnetohydrodynamics (MHD). It is shown that the left-handed helical modes play a major role on the spectral transfer properties of turbulence when the generalized helicity and magnetic helicity are both positive. In contrast, the right-handed helical modes will play a major role when both are negative. Furthermore, we also find that if the generalized helicity and magnetic helicity have opposite signs, the tendency of equilibrium spectra to condense at the large or small wave numbers will be presented in different helical sectors. This indicates that the generalized helicity dominates the forward cascade and the magnetic helicity dominates the inverse cascade properties of the Hall MHD turbulence.展开更多
An analysis is presented with magnetohydrodynamics natural convective flow of a viscous Newtonian fluid saturated porous medium in a vertical slot. The flow in the porous media has been modeled using the Brinkman mode...An analysis is presented with magnetohydrodynamics natural convective flow of a viscous Newtonian fluid saturated porous medium in a vertical slot. The flow in the porous media has been modeled using the Brinkman model. The fully-developed two-dimensional flow from capped to open ends is considered for which a continuum of solutions is obtained. The influence of pertinent parame-ters on the flow is delineated and appropriate conclusions are drawn. The asymptotic behaviour and the volume flux are analyzed and incorporated graphically for the three-parameter family of solution.展开更多
A numerical study has been carried out to investigate the temperature distribution and the natural convection heat transfer in axisymmetric two-dimensional vertical saturated porous cylinder with steady state laminar ...A numerical study has been carried out to investigate the temperature distribution and the natural convection heat transfer in axisymmetric two-dimensional vertical saturated porous cylinder with steady state laminar flow. A comparison between two situations is done under the effect of MHD (magnetohydrodynamics) and radiation. In the two situations, the vertical walls of the cylinder are cooled with constant wall temperature and a constant heat generation subjected along the centerline of the cylinder. The first case for cylinder with insulated upper surface and cooled bottom surface while the second case for cylinder with cooled upper surface and insulated bottom surface. The governing equations used are continuity, momentum and energy equations which are transformed to dimensionless equations. The finite difference approach is used to obtain all the computational results using the MATLAB-7 programming. The parameters affected the system are Rayleigh number ranging within (102≤ Ra ≤104), radiation parameter (0≤ Rd ≤ 2) and magnetohydrodynamics MHD (Mn) (0 ≤ Mn≤ 2).The results show that the temperature of Case 1 is more than that in Case 2 at constant Ra, Mn and Rd while the value of the stream in Case 2 is greater than that in Case 1. Nu increase with the increase of Rd and increasing Mn caused the temperature to increase and the streamline dropped while Nu decreased. A correlation has been set up to give the average Nusselt number variation with Ra, Rd and Mn for which the results are found to be in good agreement with previously published researches.展开更多
基金supported by the National Natural Science Foundation of China(12101095)the Natural Science Foundation of Chongqing(CSTB2022NSCQ-MSX0949,2022NSCQ-MSX2878,CSTC2021jcyj-msxmX0224)+2 种基金the Science and Technology Research Program of Chongqing Municipal Education Commission(KJQN202100517,KJQN202300542,KJQN202100511)the Research Project of Chongqing Education Commission(CXQT21014)the grant of Chongqing Young Experts’Workshop.
文摘In this paper,we consider a model of compressible isentropic two-fluid magneto-hydrodynamics without resistivity in a strip domain in three dimensional space.By exploiting the two-tier energy method developed in[Anal PDE,2013,6:1429–1533],we prove the global well-posedness of the governing model around a uniform magnetic field which is non-parallel to the horizontal boundary.Moreover,we show that the solution converges to the steady state at an almost exponential rate as time goes to infinity.Compared to the work of Tan and Wang[SIAM J Math Anal,2018,50:1432–1470],we need to overcome the difficulties caused by particles.
基金support from the European Space Agency(ESA)PRODEX(PROgramme de Développement d’Expériences scientifiques)Project mission(No.PEA4000134960)Partial funding was provided by the Romanian Ministry of Research,Innovation and Digitalization under Romanian National Core Program LAPLAS VII(Contract No.30N/2023)+2 种基金the Belgian Solar-Terrestrial Centre of Excellencesupported by the project Belgian Research Action through Interdisciplinary Networks(BRAIN-BE)2.0(Grant No.B2/223/P1/PLATINUM)funded by the Belgian Office for Research(BELSPO)partially supported by a grant from the Romanian Ministry of Education and Research(CNCS-UEFISCDI,Project No.PN-III-P1-1.1TE-2021-0102)。
文摘We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.
基金LMP acknowledges financial support from ANID through Convocatoria Nacional Subvención a Instalación en la Academia Convocatoria Año 2021,Grant SA77210040。
文摘We report on the magnetohydrodynamic impact on the axisymmetric flow of Al_(2)O_(3)/Cu nanoparticles suspended in H_(2)O past a stretched/shrinked sheet.With the use of partial differential equations and the corresponding thermophysical characteristics of nanoparticles,the physical flow process is illustrated.The resultant nonlinear system of partial differential equations is converted into a system of ordinary differential equations using the suitable similarity transformations.The transformed differential equations are solved analytically.Impacts of the magnetic parameter,solid volume fraction and stretching/shrinking parameter on momentum and temperature distribution have been analyzed and interpreted graphically.The skin friction and Nusselt number were also evaluated.In addition,existence of dual solution was deduced for the shrinking sheet and unique solution for the stretching one.Further,Al_(2)O_(3)/H_(2)O nanofluid flow has better thermal conductivity on comparing with Cu/H_(2)O nanofluid.Furthermore,it was found that the first solutions of the stream are stable and physically realizable,whereas those of the second ones are unstable.
基金supported by the Strategic Priority Research Program of Chinese Academy of Sciences,grant No.XDB 41000000National Natural Science Foundation of China(NSFC,Grant No.12288102)+4 种基金support of the National Natural Science Foundation of China(NSFC,Grant No.11973082)support of the National Natural Science Foundation of China(NSFC,Grant No.42305136)supported by the National Key R&D Program of China(Grant No.2021YFA1600400/2021YFA1600402)Natural Science Foundation of Yunnan Province(No.202201AT070158)the International Centre of Supernovae,Yunnan Key Laboratory(No.202302AN360001)。
文摘Observations of transmission spectra reveal that hot Jupiters and Neptunes are likely to possess escaping atmospheres driven by stellar radiation.Numerous models predict that magnetic fields may exert significant influences on the atmospheres of hot planets.Generally,the escaping atmospheres are not entirely ionized,and magnetic fields only directly affect the escape of ionized components within them.Considering the chemical reactions between ionized components and neutral atoms,as well as collision processes,magnetic fields indirectly impact the escape of neutral atoms,thereby influencing the detection signals of planetary atmospheres in transmission spectra.In order to simulate this process,we developed a magnetohydrodynamic multi-fluid model based on MHD code PLUTO.As an initial exploration,we investigated the impact of magnetic fields on the decoupling of H^(+)and H in the escaping atmosphere of the hot Neptune GJ436b.Due to the strong resonant interactions between H and H^(+),the coupling between them is tight even if the magnetic field is strong.Of course,alternatively,our work also suggests that merging H and H^(+)into a single flow can be a reasonable assumption in MHD simulations of escaping atmospheres.However,our simulation results indicate that under the influence of magnetic fields,there are noticeable regional differences in the decoupling of H^(+)and H.With the increase of magnetic field strength,the degree of decoupling also increases.For heavier particles such as O,the decoupling between O and H^(+)is more pronounced.Our findings provide important insights for future studies on the decoupling processes of heavy atoms in the escaping atmospheres of hot Jupiters and hot Neptunes under the influence of magnetic fields.
基金supported in part by NNSFC(11271381 and 11431015)supported in part by the Joint NSFC-RGC Research Fund,N-CityU 102/12+1 种基金supported in part by the Program for Liaoning Excellent Talents in University(LJQ2013124)the Fundamental Research Fund for the Central Universities
文摘The paper aims to estimate the thickness of the boundary layer for the planar MHD system with vanishing shear viscosity μ. Under some conditions on the initial and boundary data, we show that the thickness is of the order √μ|lnμ|. Note that this estimate holds also for the Navier-Stokes system so that it extends the previous works even without the magnetic effect.
基金The first author (PLF) was partially supported by the Centre National de la Recherche Scientifique (CNRS) the Agence Nationale de la Recherche (ANR)
文摘We consider a nonlinear hyperbolic system of two conservation laws which arises in ideal magnetohydrodynamics and includes second-order terms accounting for magnetic resistivity and Hall effect. We show that the initial value problem for this model may lead to solutions exhibiting complex wave structures, including undercompressive nonclassical shock waves. We investigate numerically the subtle competition that takes place between the hyperbolic, diffusive, and dispersive parts of the system. Following Abeyratne, Knowles, LeFloch, and Truskinovsky, who studied similar questions arising in fluid and solid flows, we determine the associated kinetic function which characterizes the dynamics of undereompressive shocks driven by resistivity and Hall effect. To this end, we design a new class of "schemes with eontroled dissipation", following recent work by LeFloch and Mohammadian. It is now recognized that the equivalent equation associated with a scheme provides a guideline to design schemes that capture physically relevant, nonclassical shocks. We propose a new class of schemes based on high-order entropy conservative, finite differences for the hyperbolic flux, and high-order central differences for the resistivity and Hall terms. These schemes are tested for several regimes of (co-planar or not) initial data and parameter values, and allow us to analyze the properties of nonclassical shocks and establish the existence of monotone kinetic functions in magnetohydrodynamics.
基金Supported by the National Natural Science Foundation of China under Grant Nos 41331067 and 41527804the Key Research Program of Frontier Sciences of Chinese Academy of Sciences under Grant No QYZDJ-SSW-DQC010the Fundamental Research Funds for the Central Universities
文摘We conduct an electron magnetohydrodynamics magnetic reconnection experiment with guide-field in our Keda linear magnetized plasma device, in which two pulsed currents with the same direction are conducted in parallel with the axial direction of the main chamber of the device using two long aluminum sticks. After approximately 5μs, an X-type magnetic field line topology is formed at the center of the chamber. With the formation of the X-type topology of magnetic field lines, we can also find the rapid increase of the current and ratio of the common flux to the private flux in this area. Additionally, a reduction in the plasma density and the plasma density concentration along one pair of separatrices can also be found.
基金This work was supported by the KAUST Office of Spon-sored Research under Award No.URF/1/2162-01.
文摘The interaction between a converging cylindrical shock and double density interfaces in the presence of a saddle magnetic field is numerically investigated within the framework of ideal magnetohydrodynamics.Three fluids of differing densities are initially separated by the two perturbed cylindrical interfaces.The initial incident converging shock is generated from a Riemann problem upstream of the first interface.The effect of the magnetic field on the instabilities is studied through varying the field strength.It shows that the Richtmyer-Meshkov and Rayleigh-Taylor instabilities are mitigated by the field,however,the extent of the suppression varies on the interface which leads to non-axisymmetric growth of the perturbations.The degree of asymmetry of the interfacial growth rate is increased when the seed field strength is increased.
基金Project supported by the National Natural Science Foundation of China(Nos.11971410 and12071404)the Natural Science Foundation of Hunan Province of China(No.2019JJ40279)+2 种基金the Excellent Youth Program of Scientific Research Project of Hunan Provincial Department of Education(Nos.18B064 and 20B564)the China Postdoctoral Science Foundation(Nos.2018T110073 and 2018M631402)the International Scientific and Technological Innovation Cooperation Base of Hunan Province for Computational Science(No.2018WK4006)。
文摘Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phases:find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution,and then locally solve linearized residual subproblems by one of three iterations(Stokes-type,Newton,and Oseen-type)on subdomains with fine grid in parallel to approximate the high frequency component.Optimal error estimates with regard to two mesh sizes and iterative steps of the proposed algorithms are given.Some numerical examples are implemented to verify the algorithm.
基金supported by NRF-2015R1A5A1009350the National Research Foundation of Korea Grant funded by the Korean Government(NRF-2016R1D1A1B03930422)
文摘We investigate the local existence of smooth solutions of a 3D ideal magneto-hydrodynamics (MHD) equations in a bounded domain and give a blow-up criteria to thisequations with respect to vorticists.
基金Project supported by the National Natural Science Foundation of China (Nos.10471100 and 40437017)
文摘A nonlinear Galerkin mixed element (NGME) method for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of the NGME solution are derived.
基金Project supported by the National Natural Science Foundation of China (Nos.10471100 and 40437017)the Science and Technology Foundation of Beijing Jiaotong University
文摘The Galerkin-Petrov least squares method is combined with the mixed finite element method to deal with the stationary, incompressible magnetohydrodynamics system of equations with viscosity. A Galerkin-Petrov least squares mixed finite element format for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of its solution are derived. Through this method, the combination among the mixed finite element spaces does not demand the discrete Babuska-Brezzi stability conditions so that the mixed finite element spaces could be chosen arbitrartily and the error estimates with optimal order could be obtained.
文摘Lattice Boltzmann Method is recently developed within numerical schemes for simulating a variety of physical systems. In this paper a new lattice Bhatnagar-Gross-Krook (LBGK) model for two-dimensional incompressible magnetohydrodynamics (IMHD) is presented. The model is an extension of a hydrodynamics lattice BGK model with 9 velocities on a square lattice, resulting in a model with 17 velocities. Most of the existing LBGK models for MHD can be viewed as compressible schemes to simulate incompressible flows. The compressible effect might lead to some undesirable errors in numerical simulations. In our model the compressible effect has been overcome successfully. The model is then applied to the Hartmann flow, giving reasonable results.
基金supported by the National Basic Research Program under the Grant 2005CB321703the National Natural Science Foundation of China(No.10925101 and No.10828101)+4 种基金the Doctoral Program of the Education Ministry of China(No.20070001036)the Program for New Century Excellent Talents in University(No.NCET-07-0022)supported by Hong Kong Research Grant Council 621709National Natural Science Foundation of China(Project No.10928205)National Key Basic Research Program(2009CB724101)
文摘The gas-kinetic theory based flux splitting method has been successfully proposed for solving one-and two-dimensional ideal magnetohydrodynamics by Xu et al. [J.Comput.Phys.,1999;2000],respectively.This paper extends the kinetic method to solve three-dimensional ideal magnetohydrodynamics equations,where an adaptive parameter 17 is used to control the numerical dissipation in the flux splitting method. Several numerical examples are given to demonstrate that the proposed method can achieve high numerical accuracy and resolve strong discontinuous waves in three dimensional ideal MHD problems.
基金Supported by the National Natural Science Foundation of China under Grant No 11375190
文摘The solutions of incompressible ideal Hall magnetohydrodynamics are obtained by using the traveling wave method. It is shown that the velocity and magnetic field parallel to the wave vector can be arbitrary constants. The velocity and magnetic field perpendicular to the wave vector are both helical waves. Moreover, the amplitude of the velocity perpendicular to the wave vector is related to the wave number and the circular frequency. In addition, further studies indicate that, no matter whether the uniform ambient magnetic field exists or not, the forms of the travelling wave solutions do not change.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11375190 and 11547137
文摘Using the Fourier helical decomposition, we obtain the absolute statistical equilibrium spectra of left- and right-handed helical modes in the incompressible ideal Hall magnetohydrodynamics (MHD). It is shown that the left-handed helical modes play a major role on the spectral transfer properties of turbulence when the generalized helicity and magnetic helicity are both positive. In contrast, the right-handed helical modes will play a major role when both are negative. Furthermore, we also find that if the generalized helicity and magnetic helicity have opposite signs, the tendency of equilibrium spectra to condense at the large or small wave numbers will be presented in different helical sectors. This indicates that the generalized helicity dominates the forward cascade and the magnetic helicity dominates the inverse cascade properties of the Hall MHD turbulence.
文摘An analysis is presented with magnetohydrodynamics natural convective flow of a viscous Newtonian fluid saturated porous medium in a vertical slot. The flow in the porous media has been modeled using the Brinkman model. The fully-developed two-dimensional flow from capped to open ends is considered for which a continuum of solutions is obtained. The influence of pertinent parame-ters on the flow is delineated and appropriate conclusions are drawn. The asymptotic behaviour and the volume flux are analyzed and incorporated graphically for the three-parameter family of solution.
文摘A numerical study has been carried out to investigate the temperature distribution and the natural convection heat transfer in axisymmetric two-dimensional vertical saturated porous cylinder with steady state laminar flow. A comparison between two situations is done under the effect of MHD (magnetohydrodynamics) and radiation. In the two situations, the vertical walls of the cylinder are cooled with constant wall temperature and a constant heat generation subjected along the centerline of the cylinder. The first case for cylinder with insulated upper surface and cooled bottom surface while the second case for cylinder with cooled upper surface and insulated bottom surface. The governing equations used are continuity, momentum and energy equations which are transformed to dimensionless equations. The finite difference approach is used to obtain all the computational results using the MATLAB-7 programming. The parameters affected the system are Rayleigh number ranging within (102≤ Ra ≤104), radiation parameter (0≤ Rd ≤ 2) and magnetohydrodynamics MHD (Mn) (0 ≤ Mn≤ 2).The results show that the temperature of Case 1 is more than that in Case 2 at constant Ra, Mn and Rd while the value of the stream in Case 2 is greater than that in Case 1. Nu increase with the increase of Rd and increasing Mn caused the temperature to increase and the streamline dropped while Nu decreased. A correlation has been set up to give the average Nusselt number variation with Ra, Rd and Mn for which the results are found to be in good agreement with previously published researches.