Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDivergence-conforming methods for transient doubly-diffusive flows: A priori and a posteriori error analysis
124
0
0.0
(
0
)
Added by
Ricardo Ruiz Baier
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The analysis of the double-diffusion model and $mathbf{H}(mathrm{div})$-conforming method introduced in [Burger, Mendez, Ruiz-Baier, SINUM (2019), 57:1318--1343] is extended to the time-dependent case. In addition, the efficiency and reliability analysis of residual-based {it a posteriori} error estimators for the steady, semi-discrete, and fully discrete problems is established. The resulting methods are applied to simulate the sedimentation of small particles in salinity-driven flows. The method consists of Brezzi-Douglas-Marini approximations for velocity and compatible piecewise discontinuous pressures, whereas Lagrangian elements are used for concentration and salinity distribution. Numerical tests confirm the properties of the proposed family of schemes and of the adaptive strategy guided by the {it a posteriori} error indicators.
rate research
Read More
Under some regularity assumptions, we report an a priori error analysis of a dG scheme for the Poisson and Stokes flow problem in their dual mixed formulation. Both formulations satisfy a Babuv{s}ka-Brezzi type condition within the space H(div) x L2. It is well known that the lowest order Crouzeix-Raviart element paired with piecewise constants satisfies such a condition on (broken) H1 x L2 spaces. In the present article, we use this pair. The continuity of the normal component is weakly imposed by penalizing jumps of the broken H(div) component. For the resulting methods, we prove well-posedness and convergence with constants independent of data and mesh size. We report error estimates in the methods natural norms and optimal local error estimates for the divergence error. In fact, our finite element solution shares for each triangle one DOF with the CR interpolant and the divergence is locally the best-approximation for any regularity. Numerical experiments support the findings and suggest that the other errors converge optimally even for the lowest regularity solutions and a crack-problem, as long as the crack is resolved by the mesh.
A number of non-standard finite element methods have been proposed in recent years, each of which derives from a specific class of PDE-constrained norm minimization problems. The most notable examples are $mathcal{L}mathcal{L}^*$ methods. In this work, we argue that all high-order methods in this class should be expected to deliver substandard uniform h-refinement convergence rates. In fact, one may not even see rates proportional to the polynomial order $p > 1$ when the exact solution is a constant function. We show that the convergence rate is limited by the regularity of an extraneous Lagrange multiplier variable which naturally appears via a saddle-point analysis. In turn, limited convergence rates appear because the regularity of this Lagrange multiplier is determined, in part, by the geometry of the domain. Numerical experiments support our conclusions.
For the Stokes equation over 2D and 3D domains, explicit a posteriori and a priori error estimation are novelly developed for the finite element solution. The difficulty in handling the divergence-free condition of the Stokes equation is solved by utilizing the extended hypercircle method along with the Scott-Vogelius finite element scheme. Since all terms in the error estimation have explicit values, by further applying the interval arithmetic and verified computing algorithms, the computed results provide rigorous estimation for the approximation error. As an application of the proposed error estimation, the eigenvalue problem of the Stokes operator is considered and rigorous bounds for the eigenvalues are obtained. The efficiency of proposed error estimation is demonstrated by solving the Stokes equation on both convex and non-convex 3D domains.
We present and analyze an a posteriori error estimator based on mesh refinement for the solution of the hypersingular boundary integral equation governing the Laplacian in three dimensions. The discretization under consideration is a non-conforming domain decomposition method based on the Nitsche technique. Assuming a saturation property, we establish quasi-reliability and efficiency of the error estimator in comparison with the error in a natural (non-conforming) norm. Numerical experiments with uniform and adaptively refined meshes confirm our theoretical results.
Hybrid quantum/molecular mechanics models (QM/MM methods) are widely used in material and molecular simulations when MM models do not provide sufficient accuracy but pure QM models are computationally prohibitive. Adaptive QM/MM coupling methods feature on-the-fly classification of atoms during the simulation, allowing the QM and MM subsystems to be updated as needed. In this work, we propose such an adaptive QM/MM method for material defect simulations based on a new residual based it a posteriori error estimator, which provides both lower and upper bounds for the true error. We validate the analysis and illustrate the effectiveness of the new scheme on numerical simulations for material defects.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
