Do you want to publish a course? Click here

A family of energy stable, skew-symmetric finite difference schemes on collocated grids

122   0   0.0 ( 0 )
 Added by Julius Reiss
 Publication date 2014
and research's language is English
 Authors Julius Reiss




Ask ChatGPT about the research

A simple scheme for incompressible, constant density flows is presented, which avoids odd-even decoupling for the Laplacian on a collocated grids. Energy stability is implied by maintaining strict energy conservation. Momentum is conserved. Arbitrary order in space and time can easily be obtained. The conservation properties hold on transformed grids.



rate research

Read More

We present a fully conservative, skew-symmetric finite difference scheme on transformed grids. The skew-symmetry preserves the kinetic energy by first principles, simultaneously avoiding a central instability mechanism and numerical damping. In contrast to other skew-symmetric schemes no special averaging procedures are needed. Instead, the scheme builds purely on point-wise operations and derivatives. Any explicit and central derivative can be used, permitting high order and great freedom to optimize the scheme otherwise. This also allows the simple adaption of existing finite difference schemes to improve their stability and damping properties.
In this paper, we develop a new free-stream preserving (FP) method for high-order upwind conservative finite-difference (FD) schemes on the curvilinear grids. This FP method is constrcuted by subtracting a reference cell-face flow state from each cell-center value in the local stencil of the original upwind conservative FD schemes, which effectively leads to a reformulated dissipation. It is convenient to implement this method, as it does not require to modify the original forms of the upwind schemes. In addition, the proposed method removes the constraint in the traditional FP conservative FD schemes that require a consistent discretization of the mesh metrics and the fluxes. With this, the proposed method is more flexible in simulating the engineering problems which usually require a low-order scheme for their low-quality mesh, while the high-order schemes can be applied to approximate the flow states to improve the resolution. After demonstrating the strict FP property and the order of accuracy by two simple test cases, we consider various validation cases, including the supersonic flow around the cylinder, the subsonic flow past the three-element airfoil, and the transonic flow around the ONERA M6 wing, etc., to show that the method is suitable for a wide range of fluid dynamic problems containing complex geometries. Moreover, these test cases also indicate that the discretization order of the metrics have no significant influences on the numerical results if the mesh resolution is not sufficiently large.
154 - R. Eymard 2020
In this paper, we present a class of finite volume schemes for incompressible flow problems. The unknowns are collocated at the center of the control volumes, and the stability of the schemes is obtained by adding to the mass balance stabilization terms involving the pressure jumps across the edges of the mesh.
This paper develops the high-order accurate entropy stable (ES) finite difference schemes for the shallow water magnetohydrodynamic (SWMHD) equations.They are built on the numerical approximation of the modified SWMHD equations with the Janhunen source term. First, the second-order accurate well-balanced semi-discrete entropy conservative (EC) schemes are constructed, satisfying the entropy identity for the given convex entropy function and preserving the steady states of the lake at rest (with zero magnetic field). The key is to match both discretizations for the fluxes and the non-flat river bed bottom and Janhunen source terms, and to find the affordable EC fluxes of the second-order EC schemes. Next, by using the second-order EC schemes as building block, high-order accurate well-balanced semi-discrete EC schemes are proposed. Then, the high-order accurate well-balanced semi-discrete ES schemes %satisfying the entropy inequality are derived by adding a suitable dissipation term to the EC scheme with the WENO reconstruction of the scaled entropy variables in order to suppress the numerical oscillations of the EC schemes. After that, the semi-discrete schemes are integrated in time by using the high-order strong stability preserving explicit Runge-Kutta schemes to obtain the fully-discrete high-order well-balanced schemes. The ES property of the Lax-Friedrichs flux is also proved and then the positivity-preserving ES schemes are studied by using the positivity-preserving flux limiter. Finally, extensive numerical tests are conducted to validate the accuracy, the well-balanced, ES and positivity-preserving properties, and the ability to capture discontinuities of our schemes.
We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are present. The proposed coupling technique requires minimal changes in the existing schemes while maintaining strict stability, accuracy, and energy conservation. Results are demonstrated on linear and nonlinear scalar conservation laws in two spatial dimensions.
comments
Fetching comments Fetching comments
mircosoft-partner

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