Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSmoothed Particle Magnetohydrodynamics Simulations of Protostellar Jets and Turbulent Dynamos
138
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We presents results from Smoothed Particle Magnetohydrodynamics simulations of collapsing molecular cloud cores, and dynamo amplification of the magnetic field in the presence of Mach 10 magnetised turbulence. Our star formation simulations have produced, for the first time ever, highly collimated magnetised protostellar jets from the first hydrostatic core phase. Up to 40% of the initial core mass may be ejected through this outflow. The primary difficulty in performing these simulations is maintaining the divergence free constraint of the magnetic field, and to address this issue, we have developed a new divergence cleaning method which has allowed us to stably follow the evolution of these protostellar jets for long periods. The simulations performed of supersonic MHD turbulence are able to exponentially amplify magnetic energy by up to 10 orders of magnitude via turbulent dynamo. To reduce numerical dissipation, a new shock detection algorithm is utilised which is able to track magnetic shocks throughout a large range of magnetic field strengths.
rate research
Read More
Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small-scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of star-forming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning method may be used to achieve $ abla cdot {bf B}=0$ to machine precision, albeit at significant computational expense.
Small-scale dynamo action is often held responsible for the generation of quiet-Sun magnetic fields. We aim to determine the excitation conditions and saturation level of small-scale dynamos in non-rotating turbulent convection at low magnetic Prandtl numbers. We use high resolution direct numerical simulations of weakly stratified turbulent convection. We find that the critical magnetic Reynolds number for dynamo excitation increases as the magnetic Prandtl number is decreased, which might suggest that small-scale dynamo action is not automatically evident in bodies with small magnetic Prandtl numbers as the Sun. As a function of the magnetic Reynolds number (${rm Rm}$), the growth rate of the dynamo is consistent with an ${rm Rm}^{1/2}$ scaling. No evidence for a logarithmic increase of the growth rate with ${rm Rm}$ is found.
The magnetohydrodynamic (MHD) description of plasmas with relativistic particles necessarily includes an additional new field, the chiral chemical potential associated with the axial charge (i.e., the number difference between right- and left-handed relativistic fermions). This chiral chemical potential gives rise to a contribution to the electric current density of the plasma (emph{chiral magnetic effect}). We present a self-consistent treatment of the emph{chiral MHD equations}, which include the back-reaction of the magnetic field on a chiral chemical potential and its interaction with the plasma velocity field. A number of novel phenomena are exhibited. First, we show that the chiral magnetic effect decreases the frequency of the Alfv{e}n wave for incompressible flows, increases the frequencies of the Alfv{e}n wave and of the fast magnetosonic wave for compressible flows, and decreases the frequency of the slow magnetosonic wave. Second, we show that, in addition to the well-known laminar chiral dynamo effect, which is not related to fluid motions, there is a dynamo caused by the joint action of velocity shear and chiral magnetic effect. In the presence of turbulence with vanishing mean kinetic helicity, the derived mean-field chiral MHD equations describe turbulent large-scale dynamos caused by the chiral alpha effect, which is dominant for large fluid and magnetic Reynolds numbers. The chiral alpha effect is due to an interaction of the chiral magnetic effect and fluctuations of the small-scale current produced by tangling magnetic fluctuations (which are generated by tangling of the large-scale magnetic field by sheared velocity fluctuations). These dynamo effects may have interesting consequences in the dynamics of the early universe, neutron stars, and the quark--gluon plasma.
There have been some issues in the past in attempts to simulate magnetic fields using the Smoothed Particle Hydrodynamics (SPH) method. SPH is well suited to star formation problems because of its Lagrangian nature. We present new, stable and conservative methods for magnetohydrodynamics (MHD) in SPH and present numerical tests on both waves and shocks in one dimension to show that it gives robust and accurate results.
Numerical methods to improve the treatment of magnetic fields in smoothed field magnetohydrodynamics (SPMHD) are developed and tested. Chapter 2 is a review of SPMHD. In Chapter 3, a mixed hyperbolic/parabolic scheme is developed which cleans divergence error from the magnetic field. Average divergence error is an order of magnitude lower for all test cases considered, and allows for the stable simulation of the gravitational collapse of magnetised molecular cloud cores. The effectiveness of the cleaning may be improved by explicitly increasing the hyperbolic wave speed or by cycling the cleaning equations between timesteps. In the latter, it is possible to achieve DivB=0. Chapter 4 develops a switch to reduce dissipation of the magnetic field from artificial resistivity. Compared to the existing switch in the literature, this leads to sharper shock profiles in shocktube tests, lower overall dissipation of magnetic energy, and importantly, is able to capture magnetic shocks in the highly super-Alfvenic regime. Chapter 5 compares these numerical methods against grid-based MHD methods (using the Flash code) in simulations of the small-scale dynamo amplification of a magnetic field in driven, isothermal, supersonic turbulence. Both codes exponentially amplify the magnetic energy at a constant rate, though SPMHD shows a resolution dependence that arises from the scaling of the numerical dissipation terms. The time-averaged saturated magnetic spectra have similar shape, and both codes have PDFs of magnetic field strength that are log-normal, which become lopsided as the magnetic field saturates. We conclude that SPMHD is able to reliably simulate the small-scale dynamo amplification of magnetic fields. Chapter 6 concludes the thesis and presents some preliminary work demonstrating that SPMHD can activate the magneto-rotational instability in 2D shearing box tests.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
