Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPoint Forces in Elasticity Equation and Their Alternatives in Multi Dimensions
71
0
0.0
(
0
)
Added by
Qiyao Peng
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We consider several mathematical issues regarding models that simulate forces exerted by cells. Since the size of cells is much smaller than the size of the domain of computation, one often considers point forces, modelled by Dirac Delta distributions on boundary segments of cells. In the current paper, we treat forces that are directed normal to the cell boundary and that are directed toward the cell centre. Since it can be shown that there exists no smooth solution, at least not in $H^1$ for solutions to the governing momentum balance equation, we analyse the convergence and quality of the approximation. Furthermore, the expected finite element problems that we get necessitate scrutinizing alternative model formulations, such as the use of smoothed Dirac distributions, or the so-called smoothed particle approach as well as the so-called hole approach where cellular forces are modelled through the use of (natural) boundary conditions. In this paper, we investigate and attempt to quantify the conditions for consistency between the various approaches. This has resulted in error analyses in the $H^1$-norm of the numerical solution based on Galerkin principles that entail Lagrangian basis functions. The paper also addresses well-posedness in terms of existence and uniqueness. The current analysis has been performed for the linear steady-state (hence neglecting inertia and damping) momentum equations under the assumption of Hookes law.
rate research
Read More
A finite element elasticity complex on tetrahedral meshes is devised. The $H^1$ conforming finite element is the smooth finite element developed by Neilan for the velocity field in a discrete Stokes complex. The symmetric div-conforming finite element is the Hu-Zhang element for stress tensors. The construction of an $H(textrm{inc})$-conforming finite element for symmetric tensors is the main focus of this paper. The key tools of the construction are the decomposition of polynomial tensor spaces and the characterization of the trace of the $textrm{inc}$ operator. The polynomial elasticity complex and Koszul elasticity complex are created to derive the decomposition of polynomial tensor spaces. The trace of the $textrm{inc}$ operator is induced from a Greens identity. Trace complexes and bubble complexes are also derived to facilitate the construction. Our construction appears to be the first $H(textrm{inc})$-conforming finite elements on tetrahedral meshes without further splits.
We propose a numerical method to approximate the scattering amplitudes for the elasticity system with a non-constant matrix potential in dimensions $d=2$ and $3$. This requires to approximate first the scattering field, for some incident waves, which can be written as the solution of a suitable Lippmann-Schwinger equation. In this work we adapt the method introduced by G. Vainikko in cite{V} to solve such equations when considering the Lame operator. Convergence is proved for sufficiently smooth potentials. Implementation details and numerical examples are also given.
In this paper, we design a novel class of arbitrarily high-order structure-preserving numerical schemes for the time-dependent Gross-Pitaevskii equation with angular momentum rotation in three dimensions. Based on the idea of the scalar auxiliary variable approach which is proposed in the recent papers [J. Comput. Phys., 416 (2018) 353-407 and SIAM Rev., 61(2019) 474-506] for developing energy stable schemes for gradient flow systems, we firstly reformulate the Gross-Pitaevskii equation into an equivalent system with a modified energy conservation law. The reformulated system is then discretized by the Gauss collocation method in time and the standard Fourier pseudo-spectral method in space, respectively. We show that the proposed schemes can preserve the discrete mass and modified energy exactly. Numerical results are addressed to verify the efficiency and high-order accuracy of the proposed schemes.
The derivation of nonlocal strong forms for many physical problems remains cumbersome in traditional methods. In this paper, we apply the variational principle/weighted residual method based on nonlocal operator method for the derivation of nonlocal forms for elasticity, thin plate, gradient elasticity, electro-magneto-elasticity and phase field fracture method. The nonlocal governing equations are expressed as integral form on support and dual-support. The first example shows that the nonlocal elasticity has the same form as dual-horizon non-ordinary state-based peridynamics. The derivation is simple and general and it can convert efficiently many local physical models into their corresponding nonlocal forms. In addition, a criterion based on the instability of the nonlocal gradient is proposed for the fracture modelling in linear elasticity. Several numerical examples are presented to validate nonlocal elasticity and the nonlocal thin plate .
We develop an algorithm for the concurrent (on-the-fly) estimation of parameters for a system of evolutionary dissipative partial differential equations in which the state is partially observed. The intuitive nature of the algorithm makes its extension to several different systems immediate, and it allows for recovery of multiple parameters simultaneously. We test this algorithm on the Kuramoto-Sivashinsky equation in one dimension and demonstrate its efficacy in this context.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
