Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA Priori Analysis of a Higher Order Nonlinear Elasticity Model for an Atomistic Chain with Periodic Boundary Condition
56
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Nonlinear elastic models are widely used to describe the elastic response of crystalline solids, for example, the well-known Cauchy-Born model. While the Cauchy-Born model only depends on the strain, effects of higher order strain gradients are significant and higher order continuum models are preferred, in various applications such as defect dynamics and modeling of carbon nanotubes. In this paper, we rigorously derive a higher order nonlinear elasticity model for crystals from its atomistic description in one dimension. We show that, compared to the second order accuracy of the Cauchy-Born model, the higher order continuum model in this paper is of fourth oder accuracy with respect to the interatomic spacing in the thermal dynamic limit. In addition, we discuss the key issues for the derivation of higher order continuum models in more general cases. The theoretical convergence results are demonstrated by numerical experiments.
rate research
Read More
In this work we propose and analyze a novel Hybrid High-Order discretization of a class of (linear and) nonlinear elasticity models in the small deformation regime which are of common use in solid mechanics. The proposed method is valid in two and three space dimensions, it supports general meshes including polyhedral elements and nonmatching interfaces, enables arbitrary approximation order, and the resolution cost can be reduced by statically condensing a large subset of the unknowns for lineariz
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.
A hybrid surface integral equation partial differential equation (SIE-PDE) formulation without the boundary condition requirement is proposed to solve the electromagnetic problems. In the proposed formulation, the computational domain is decomposed into two emph{overlapping} domains: the SIE and PDE domains. In the SIE domain, complex structures with piecewise homogeneous media, e.g., highly conductive media, are included. An equivalent model for those structures is constructed through replacing them by the background medium and introducing a surface equivalent electric current density on an enclosed boundary to represent their electromagnetic effects. The remaining computational domain and homogeneous background medium replaced domain consist of the PDE domain, in which inhomogeneous or non-isotropic media are included. Through combining the surface equivalent electric current density and the inhomogeneous Helmholtz equation, a hybrid SIE-PDE formulation is derived. Unlike other hybrid formulations, where the transmission condition is usually used, no boundary conditions are required in the proposed SIE-PDE formulation, and it is mathematically equivalent to the original physical model. Through careful construction of basis functions to expand electric fields and the equivalent current density, the discretized formulation is compatible on the interface of the SIE and PDE domain. Finally, its accuracy and efficiency are validated through two numerical examples. Results show that the proposed SIE-PDE formulation can obtain accurate results including both near and far fields, and significant performance improvements in terms of CPU time and memory consumption compared with the FEM are achieved.
This paper is devoted to studying a type of contact problems modeled by hemivariational inequalities with small periodic coefficients appearing in PDEs, and the PDEs we considered are linear, second order and uniformly elliptic. Under the assumptions, it is proved that the original problem can be homogenized, and the solution weakly converges. We derive an $O(epsilon^{1/2})$ estimation which is pivotal in building the computational framework. We also show that Robin problems--- a special case of contact problems, it leads to an $O(epsilon)$ estimation in $L^2$ norm. Our computational framework is based on finite element methods, and the numerical analysis is given, together with experiments to convince the estimation.
This paper provides an a~priori error analysis of a localized orthogonal decomposition method (LOD) for the numerical stochastic homogenization of a model random diffusion problem. If the uniformly elliptic and bounded random coefficient field of the model problem is stationary and satisfies a quantitative decorrelation assumption in form of the spectral gap inequality, then the expected $L^2$ error of the method can be estimated, up to logarithmic factors, by $H+(varepsilon/H)^{d/2}$; $varepsilon$ being the small correlation length of the random coefficient and $H$ the width of the coarse finite element mesh that determines the spatial resolution. The proof bridges recent results of numerical homogenization and quantitative stochastic homogenization.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
