Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHigh order numerical simulations of the Richtmyer Meshkov instability in a relativistic fluid
156
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study the Richtmyer--Meshkov (RM) instability of a relativistic perfect fluid by means of high order numerical simulations with adaptive mesh refinement (AMR). The numerical scheme adopts a finite volume Weighted Essentially Non-Oscillatory (WENO) reconstruction to increase accuracy in space, a local space-time discontinuous Galerkin predictor method to obtain high order of accuracy in time and a high order one-step time update scheme together with a cell-by-cell space-time AMR strategy with time-accurate local time stepping. In this way, third order accurate (both in space and in time) numerical simulations of the RM instability are performed, spanning a wide parameter space. We present results both for the case in which a light fluid penetrates into a higher density one (Atwood number $A>0$), and for the case in which a heavy fluid penetrates into a lower density one (Atwood number $A<0$). We find that, for large Lorentz factors gamma_s of the incident shock wave, the relativistic RM instability is substantially weakened and ultimately suppressed. More specifically, the growth rate of the RM instability in the linear phase has a local maximum which occurs at a critical value of gamma_s ~ [1.2,2]. Moreover, we have also revealed a genuine relativistic effect, absent in Newtonian hydrodynamics, which arises in three dimensional configurations with a non-zero velocity component tangent to the incident shock front. In this case, the RM instability is strongly affected, typically resulting in less efficient mixing of the fluid.
rate research
Read More
We study the growth of Richtmyer-Meshkov mixing layers from an initial surface with spatially localized perturbations. We use two symmetric forms of the initial patch, which allow simulation data to be averaged to generate a two-dimensional statistical representation of the three dimensional turbulent flow. We find that as the mixing layer grows, the turbulent structures tend to form into discrete packets separated from the surface, with material entrainment into them dominated by a laminar entrainment flow inward from the surrounding regions where the surface was originally smooth. The entrainment appears to be controlled by the propagation of vortex pairs which appear at the boundary of the region of initial perturbations. This suggests that the growth of RM mixing from isolated features, as may be found in manufactured Inertial Confinement Fusion capsules, has a rather different mechanism than the growth of an RM mixing layer when the perturbations are uniform. This may be a challenge for some existing engineering models.
The two-fluid (ions and electrons) plasma Richtmyer-Meshkov instability of a cylindrical light/heavy density interface is numerically investigated without an initial magnetic field. Varying the Debye length scale, we examine the effects of the coupling between the electron and ion fluids. When the coupling becomes strong, the electrons are restricted to co-move with the ions and the resulting evolution is similar to the hydrodynamic neutral fluid case. The charge separation that occurs between the electrons and ions results in self-generated electromagnetic fields. We show that the Biermann battery effect dominates the generation of magnetic field when the coupling between the electrons and ions is weak. In addition to the Rayleigh-Tayler stabilization effect during flow deceleration, the interfaces are accelerated by the induced spatio-temporally varying Lorentz force. As a consequence, the perturbations develop into the Rayleigh-Taylor instability, leading to an enhancement of the perturbation amplitude compared with the hydrodynamic case.
We provide numerical simulations of an incompressible pressure-thickening and shear-thinning lubricant flowing in a plane slider bearing. We study the influence of several parameters, namely the ratio of the characteristic lengths $varepsilon>0$ (with $varepsilonsearrow0$ representing the Reynolds lubrication approximation); the coefficient of the exponential pressure--viscosity relation $alpha^*geq0$; the parameter $G^*geq0$ related to the Carreau--Yasuda shear-thinning model and the modified Reynolds number $mathrm{Re}_varepsilongeq0$. The finite element approximations to the steady isothermal flows are computed without resorting to the lubrication approximation. We obtain the numerical solutions as long as the variation of the viscous stress $mathbf{S}=2eta(p,mathrm{tr}mathbf{D}^2)mathbf{D}$ with the pressure is limited, say $|partialmathbf{S}/partial p|leq1$. We show conclusively that the existing practice of avoiding the numerical difficulties by cutting the viscosity off for large pressures leads to results that depend sorely on the artificial cut-off parameter. We observe that the piezoviscous rheology generates pressure differences across the fluid film.
We provide a road towards obtaining gravitational waveforms from inspiraling material binaries with an accuracy viable for third-generation gravitational wave detectors, without necessarily advancing computational hardware or massively-parallel software infrastructure. We demonstrate a proof-of-principle 1+1-dimensional numerical implementation that exhibits up to 7th-order convergence for highly dynamic barotropic stars in curved spacetime, and numerical errors up to 6 orders of magnitude smaller than a standard method. Aside from high-order interpolation errors (Runges phenomenon), there are no obvious fundamental obstacles to obtaining convergence of even higher order. The implementation uses a novel surface-tracking method, where the surface is evolved and high-order accurate boundary conditions are imposed there. Computational memory does not need to be allocated to fluid variables in the vacuum region of spacetime. We anticipate the application of this new method to full $3! +! 1$-dimensional simulations of the inspiral phase of compact binary systems with at least one material body. The additional challenge of a deformable surface must be addressed in multiple spatial dimensions, but it is also an opportunity to input more precise surface tension physics.
Computational fluid dynamics is a direct modeling of physical laws in a discretized space. The basic physical laws include the mass, momentum and energy conservations, physically consistent transport process, and similar domain of dependence and influence between the physical reality and the numerical representation. Therefore, a physically soundable numerical scheme must be a compact one which involves the closest neighboring cells within the domain of dependence for the solution update under a CFL number $(sim 1 )$. In the construction of explicit high-order compact scheme, subcell flow distributions or the equivalent degree of freedoms beyond the cell averaged flow variables must be evolved and updated, such as the gradients of the flow variables inside each control volume. The direct modeling of flow evolution under generalized initial condition will be developed in this paper. The direct modeling will provide the updates of flow variables differently on both sides of a cell interface and limit high-order time derivatives of the flux function nonlinearly in case of discontinuity in time, such as a shock wave moving across a cell interface within a time step. The direct modeling unifies the nonlinear limiters in both space for the data reconstruction and time for the time-dependent flux transport. Under the direct modeling framework, as an example, the high-order compact gas-kinetic scheme (GKS) will be constructed. The scheme shows significant improvement in terms of robustness, accuracy, and efficiency in comparison with the previous high-order compact GKS.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
