Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA mass conservative scheme for fluid-structure interaction problems by the staggered discontinuous Galerkin method
135
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this paper, we develop a new mass conservative numerical scheme for the simulations of a class of fluid-structure interaction problems. We will use the immersed boundary method to model the fluid-structure interaction, while the fluid flow is governed by the incompressible Navier-Stokes equations. The immersed boundary method is proven to be a successful scheme to model fluid-structure interactions. To ensure mass conservation, we will use the staggered discontinuous Galerkin method to discretize the incompressible Navier-Stokes equations. The staggered discontinuous Galerkin method is able to preserve the skew-symmetry of the convection term. In addition, by using a local postprocessing technique, the weakly divergence free velocity can be used to compute a new postprocessed velocity, which is exactly divergence free and has a superconvergence property. This strongly divergence free velocity field is the key to the mass conservation. Furthermore, energy stability is improved by the skew-symmetric discretization of the convection term. We will present several numerical results to show the performance of the method.
rate research
Read More
A mass-conservative Lagrange--Galerkin scheme of second order in time for convection-diffusion problems is presented, and convergence with optimal error estimates is proved in the framework of $L^2$-theory. The introduced scheme maintains the advantages of the Lagrange--Galerkin method, i.e., CFL-free robustness for convection-dominated problems and a symmetric and positive coefficient matrix resulting from the discretization. In addition, the scheme conserves the mass on the discrete level. Unconditional stability and error estimates of second order in time are proved by employing two new key lemmas on the truncation error of the material derivative in conservative form and on a discrete Gronwall inequality for multistep methods. The mass-conservation property is achieved by the Jacobian multiplication technique introduced by Rui and Tabata in 2010, and the accuracy of second order in time is obtained based on the idea of the multistep Galerkin method along characteristics originally introduced by Ewing and Russel in 1981. For the first time step, the mass-conservative scheme of first order in time by Rui and Tabata in 2010 is employed, which is efficient and does not cause any loss of convergence order in the $ell^infty(L^2)$- and $ell^2(H^1_0)$-norms. For the time increment $Delta t$, the mesh size $h$ and a conforming finite element space of polynomial degree $k$, the convergence order is of $O(Delta t^2 + h^k)$ in the $ell^infty(L^2)cap ell^2(H^1_0)$-norm and of $O(Delta t^2 + h^{k+1})$ in the $ell^infty(L^2)$-norm if the duality argument can be employed. Error estimates of $O(Delta t^{3/2}+h^k)$ in discre
In this paper, we construct an efficient numerical scheme for full-potential electronic structure calculations of periodic systems. In this scheme, the computational domain is decomposed into a set of atomic spheres and an interstitial region, and different basis functions are used in different regions: radial basis functions times spherical harmonics in the atomic spheres and plane waves in the interstitial region. These parts are then patched together by discontinuous Galerkin (DG) method. Our scheme has the same philosophy as the widely used (L)APW methods in materials science, but possesses systematically spectral convergence rate. We provide a rigorous a priori error analysis of the DG approximations for the linear eigenvalue problems, and present some numerical simulations in electronic structure calculations.
We consider a class of time dependent second order partial differential equations governed by a decaying entropy. The solution usually corresponds to a density distribution, hence positivity (non-negativity) is expected. This class of problems covers important cases such as Fokker-Planck type equations and aggregation models, which have been studied intensively in the past decades. In this paper, we design a high order discontinuous Galerkin method for such problems. If the interaction potential is not involved, or the interaction is defined by a smooth kernel, our semi-discrete scheme admits an entropy inequality on the discrete level. Furthermore, by applying the positivity-preserving limiter, our fully discretized scheme produces non-negative solutions for all cases under a time step constraint. Our method also applies to two dimensional problems on Cartesian meshes. Numerical examples are given to confirm the high order accuracy for smooth test cases and to demonstrate the effectiveness for preserving long time asymptotics.
In this article, using the weighted discrete least-squares, we propose a patch reconstruction finite element space with only one degree of freedom per element. As the approximation space, it is applied to the discontinuous Galerkin methods with the upwind scheme for the steady-state convection-diffusion-reaction problems over polytopic meshes. The optimal error estimates are provided in both diffusion-dominated and convection-dominated regimes. Furthermore, several numerical experiments are presented to verify the theoretical error estimates, and to well approximate boundary layers and/or internal layers.
We propose a new discontinuous Galerkin method based on the least-squares patch reconstruction for the biharmonic problem. We prove the optimal error estimate of the proposed method. The two-dimensional and three-dimensional numerical examples are presented to confirm the accuracy and efficiency of the method with several boundary conditions and several types of polygon meshes and polyhedral meshes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
