Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn extended finite element method with smooth nodal stress
108
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The enrichment formulation of double-interpolation finite element method (DFEM) is developed in this paper. DFEM is first proposed by Zheng emph{et al} (2011) and it requires two stages of interpolation to construct the trial function. The first stage of interpolation is the same as the standard finite element interpolation. Then the interpolation is reproduced by an additional procedure using the nodal values and nodal gradients which are derived from the first stage as interpolants. The re-constructed trial functions are now able to produce continuous nodal gradients, smooth nodal stress without post-processing and higher order basis without increasing the total degrees of freedom. Several benchmark numerical examples are performed to investigate accuracy and efficiency of DFEM and enriched DFEM. When compared with standard FEM, super-convergence rate and better accuracy are obtained by DFEM. For the numerical simulation of crack propagation, better accuracy is obtained in the evaluation of displacement norm, energy norm and the stress intensity factor.
rate research
Read More
For the Hodge--Laplace equation in finite element exterior calculus, we introduce several families of discontinuous Galerkin methods in the extended Galerkin framework. For contractible domains, this framework utilizes seven fields and provides a unifying inf-sup analysis with respect to all discretization and penalty parameters. It is shown that the proposed methods can be hybridized as a reduced two-field formulation.
In this paper, a stabilized extended finite element method is proposed for Stokes interface problems on unfitted triangulation elements which do not require the interface align with the triangulation. The velocity solution and pressure solution on each side of the interface are separately expanded in the standard nonconforming piecewise linear polynomials and the piecewise constant polynomials, respectively. Harmonic weighted fluxes and arithmetic fluxes are used across the interface and cut edges (segment of the edges cut by the interface), respectively. Extra stabilization terms involving velocity and pressure are added to ensure the stable inf-sup condition. We show a priori error estimates under additional regularity hypothesis. Moreover, the errors {in energy and $L^2$ norms for velocity and the error in $L^2$ norm for pressure} are robust with respect to the viscosity {and independent of the location of the interface}. Results of numerical experiments are presented to {support} the theoretical analysis.
We introduce a new method for the numerical approximation of time-harmonic acoustic scattering problems stemming from material inhomogeneities. The method works for any frequency $omega$, but is especially efficient for high-frequency problems. It is based on a time-domain approach and consists of three steps: emph{i)} computation of a suitable incoming plane wavelet with compact support in the propagation direction; emph{ii)} solving a scattering problem in the time domain for the incoming plane wavelet; emph{iii)} reconstruction of the time-harmonic solution from the time-domain solution via a Fourier transform in time. An essential ingredient of the new method is a front-tracking mesh adaptation algorithm for solving the problem in emph{ii)}. By exploiting the limited support of the wave front, this allows us to make the number of the required degrees of freedom to reach a given accuracy significantly less dependent on the frequency $omega$, as shown in the numerical experiments. We also present a new algorithm for computing the Fourier transform in emph{iii)} that exploits the reduced number of degrees of freedom corresponding to the adapted meshes.
We develop a general framework for construction and analysis of discrete extension operators with application to unfitted finite element approximation of partial differential equations. In unfitted methods so called cut elements intersected by the boundary occur and these elements must in general by stabilized in some way. Discrete extension operators provides such a stabilization by modification of the finite element space close to the boundary. More precisely, the finite element space is extended from the stable interior elements over the boundary in a stable way which also guarantees optimal approximation properties. Our framework is applicable to all standard nodal based finite elements of various order and regularity. We develop an abstract theory for elliptic problems and associated parabolic time dependent partial differential equations and derive a priori error estimates. We finally apply this to some examples of partial differential equations of different order including the interface problems, the biharmonic operator and the sixth order triharmonic operator.
We present a stabilized finite element method for the numerical solution of cavitation in lubrication, modeled as an inequality-constrained Reynolds equation. The cavitation model is written as a variable coefficient saddle-point problem and approximated by a residual-based stabilized method. Based on our recent results on the classical obstacle problem, we present optimal a priori estimates and derive novel a posteriori error estimators. The method is implemented as a Nitsche-type finite element technique and shown in numerical computations to be superior to the usually applied penalty methods.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
