ترغب بنشر مسار تعليمي؟ اضغط هنا

Divergence-Free WENO Reconstruction-Based Finite Volume Scheme for Solving Ideal MHD Equations on Triangular Meshes

273   0   0.0 ( 0 )
 نشر من قبل Zhiliang Xu
 تاريخ النشر 2011
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

In this paper, we introduce a high-order accurate constrained transport type finite volume method to solve ideal magnetohydrodynamic equations on two-dimensional triangular meshes. A new divergence-free WENO-based reconstruction method is developed to maintain exactly divergence-free evolution of the numerical magnetic field. A new weighted flux interpolation approach is also developed to compute the z-component of the electric field at vertices of grid cells. We also present numerical examples to demonstrate the accuracy and robustness of the proposed scheme.



قيم البحث

اقرأ أيضاً

97 - Chunjae Park 2016
Recently, the $P_1$-nonconforming finite element space over square meshes has been proved stable to solve Stokes equations with the piecewise constant space for velocity and pressure, respectively. In this paper, we will introduce its locally diverge nce-free subspace to solve the elliptic problem for the velocity only decoupled from the Stokes equation. The concerning system of linear equations is much smaller compared to the Stokes equations. Furthermore, it is split into two smaller ones. After solving the velocity first, the pressure in the Stokes problem can be obtained by an explicit method very rapidly.
We present a new method for evolving the equations of magnetohydrodynamics (both Newtonian and relativistic) that is capable of maintaining a divergence-free magnetic field ($ abla cdot mathbf{B} = 0$) on adaptively refined, conformally moving meshes . The method relies on evolving the magnetic vector potential and then using it to reconstruct the magnetic fields. The advantage of this approach is that the vector potential is not subject to a constraint equation in the same way the magnetic field is, and so can be refined and moved in a straightforward way. We test this new method against a wide array of problems from simple Alfven waves on a uniform grid to general relativistic MHD simulations of black hole accretion on a nested, spherical-polar grid. We find that the code produces accurate results and in all cases maintains a divergence-free magnetic field to machine precision.
This paper presents stability and convergence analysis of a finite volume scheme (FVS) for solving aggregation, breakage and the combined processes by showing Lipschitz continuity of the numerical fluxes. It is shown that the FVS is second order conv ergent independently of the meshes for pure breakage problem while for pure aggregation and coupled equations, it shows second order convergent on uniform and non-uniform smooth meshes. Furthermore, it gives only first order convergence on non-uniform grids. The mathematical results of convergence analysis are also demonstrated numerically for several test problems.
A discontinuous Galerkin (DG) method suitable for large-scale astrophysical simulations on Cartesian meshes as well as arbitrary static and moving Voronoi meshes is presented. Most major astrophysical fluid dynamics codes use a finite volume (FV) app roach. We demonstrate that the DG technique offers distinct advantages over FV formulations on both static and moving meshes. The DG method is also easily generalized to higher than second-order accuracy without requiring the use of extended stencils to estimate derivatives (thereby making the scheme highly parallelizable). We implement the technique in the AREPO code for solving the fluid and the magnetohydrodynamic (MHD) equations. By examining various test problems, we show that our new formulation provides improved accuracy over FV approaches of the same order, and reduces post-shock oscillations and artificial diffusion of angular momentum. In addition, the DG method makes it possible to represent magnetic fields in a locally divergence-free way, improving the stability of MHD simulations and moderating global divergence errors, and is a viable alternative for solving the MHD equations on meshes where Constrained-Transport (CT) cannot be applied. We find that the DG procedure on a moving mesh is more sensitive to the choice of slope limiter than is its FV method counterpart. Therefore, future work to improve the performance of the DG scheme even further will likely involve the design of optimal slope limiters. As presently constructed, our technique offers the potential of improved accuracy in astrophysical simulations using the moving mesh AREPO code as well as those employing adaptive mesh refinement (AMR).
We present a practical approach for constructing meshes of general rough surfaces with given autocorrelation functions based on the unstructured meshes of nominally smooth surfaces. The approach builds on a well-known method to construct correlated r andom numbers from white noise using a decomposition of the autocorrelation matrix. We discuss important details arising in practical applications to the physicalmodeling of surface roughness and provide a software implementation to enable use of the approach with a broad range of numerical methods in various fields of science and engineering.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا