Elastic Equation Research Articles

Strain accumulation and relaxation on crack formation in epitaxial AlN film on Si (111) substrate

The formation of cracks is often observed in the epitaxial growth of ultrawide-bandgap aluminum nitride (AlN) semiconductor films on economical and versatile silicon (Si) substrates due to the significant differences in in-plane lattice parameters and thermal expansion coefficients between the film and the substrate, which hampers the development of template, buffer layer, and device structure with a relatively thick AlN layer for devices. The present study aims to elucidate the conditions of crack formation through a simple but comprehensive estimation of strain energy accumulation and relaxation by lattice strain, misfit dislocation density, and crack formation. Strain energy in the epitaxial film from lattice and thermal mismatches is evaluated by an elastic strain equation tailored to the epitaxy of the hexagonal crystal structure. The effects of temperature, thickness, and dislocation density on the lattice and dislocation strain energies of the film are also considered. Finally, the comparison in the changes in the total strain energy and cleavage energy with decreasing temperature shows that cleavage energy is higher than strain energy if the film is thinner than 400 nm but becomes lower than the strain energy if the film is thicker than 400 nm during cooldown, suggesting the crack formation, which matches well with experimental observations.

Asymptotic long-wave model for a high-contrast two-layered elastic plate

The paper is concerned with the derivation of asymptotically consistent equations governing the long-wave flexural response of a two-layered rectangular plate with high-contrast elastic properties. In the general case, the plate is under dynamic and variable surface, volume, and edge forces. Performing the asymptotic integration of the three-dimensional (3D) elasticity equations in the transverse direction and satisfying boundary conditions on the faces and interface, we derived the sequence of two-dimensional (2D) differential equations with respect to required functions in the first two approximations. The eight independent restraints for the generalized displacements and stress resultants are considered to formulate the 16 independent variants of boundary conditions. One of the main results of the paper is the Timoshenko–Reissner type equation capturing the effect of the softer layer and taking into account the in-plane deformation induced by the edge forces. Comparative calculations of natural frequencies were carried out based on our and alternative models.

Discovering interpretable physical models using symbolic regression and discrete exterior calculus

Computational modeling is a key resource to gather insight into physical systems in modern scientific research and engineering. While access to large amount of data has fueled the use of machine learning to recover physical models from experiments and increase the accuracy of physical simulations, purely data-driven models have limited generalization and interpretability. To overcome these limitations, we propose a framework that combines symbolic regression (SR) and discrete exterior calculus (DEC) for the automated discovery of physical models starting from experimental data. Since these models consist of mathematical expressions, they are interpretable and amenable to analysis, and the use of a natural, general-purpose discrete mathematical language for physics favors generalization with limited input data. Importantly, DEC provides building blocks for the discrete analog of field theories, which are beyond the state-of-the-art applications of SR to physical problems. Further, we show that DEC allows to implement a strongly-typed SR procedure that guarantees the mathematical consistency of the recovered models and reduces the search space of symbolic expressions. Finally, we prove the effectiveness of our methodology by re-discovering three models of continuum physics from synthetic experimental data: Poisson equation, the Euler’s elastica and the equations of linear elasticity. Thanks to their general-purpose nature, the methods developed in this paper may be applied to diverse contexts of physical modeling.

Efficient and reliable divergence-conforming methods for an elasticity-poroelasticity interface problem

We present a finite element discretization to model the interaction between a poroelastic structure and an elastic medium. The consolidation problem considers fully coupled deformations across an interface, ensuring continuity of displacement and total traction, as well as no-flux for the fluid phase. Our formulation of the poroelasticity equations incorporates displacement, fluid pressure, and total pressure, while the elasticity equations adopt a displacement-pressure formulation. Notably, the transmission conditions at the interface are enforced without the need for Lagrange multipliers. We demonstrate the stability and convergence of the divergence-conforming finite element method across various polynomial degrees. The a priori error bounds remain robust, even when considering large variations in intricate model parameters such as Lamé constants, permeability, and storativity coefficient. To enhance computational efficiency and reliability, we develop residual-based a posteriori error estimators that are independent of the aforementioned coefficients. Additionally, we devise parameter-robust and optimal block diagonal preconditioners. Through numerical examples, including adaptive scenarios, we illustrate the scheme's properties such as convergence and parameter robustness.

On the enriched mixed Galerkin method with θ scheme for the elastic wave equation

In this paper, we propose the application of an enriched Galerkin (EG) method to solve the elastic wave equation in the “displacement pseudo-pressure” mixed form, which is known for being locking-free. Unlike the continuous Galerkin (CG) method and the discontinuous Galerkin (DG) method, the EG method is based on the weak form of the DG method, and the enriched space combines a continuous piecewise linear vector-valued function space with some discontinuous piecewise linear functions. Specifically, we employ the enriched space to discretize the displacement, while the pseudo-pressure is discretized using the P0 element space. As a result, compared to the DG method, the EG method significantly reduces the number of degrees of freedom. Furthermore, we adopt a more general θ-scheme for time discretization, which behaves as an implicit scheme for 14≤θ≤1 and as an explicit scheme for 0≤θ<14. We establish error estimates for both the semi-discrete and fully-discrete schemes. The accuracy of the proposed method is validated through numerical examples, which not only confirm the theoretical analysis but also provide simulations of elastic wave propagation.

Time domain modeling of elastic waves using a stress-based unsplit-field perfectly matched layer with enhanced numerical stability

A stress-based unsplit-field perfectly matched layer (PML) method with enhanced numerical stability is proposed to simulate elastic waves in heterogeneous unbounded domains. The proposed PML method attenuates the waves by forcing the satisfaction of mixed-field equations constructed using a complex coordinate stretching in the frequency domain. Unlike other unsplit-field PMLs, which introduce coordinate stretching before combining the constitutive and compatibility equations, the proposed method first combines the two equations and then stretches the resulting equation. The stretched equation is converted back to the time domain using the inverse Fourier transform to yield transient elastic wave equations in the PML-truncated domain. Semi-discrete equations of motion for displacement and stresses were obtained using standard Lagrange-family finite elements. In a series of numerical examples involving semi-infinite elastic media, the accuracy of the wave solutions calculated using the proposed PML method was evaluated using a displacement L2-norm, a normalized L2-error, and an energy norm. The examples further validated the wave absorption performance of the PML and evaluated the long-term stability of solutions utilizing the Ricker pulse excitation with various central frequencies. The numerical experiments indicate that the solutions obtained using the proposed PML are significantly more stable than those obtained using the existing strain-based unsplit-field PMLs. The proposed method is particularly useful for applications requiring explicit time integration of semi-discrete equations of motion because the mass matrix can be diagonalized. The PML method and its associated solution strategies can be applied to various engineering problems, such as soil-structure interaction analysis, monitoring of structures, and geophysical probing.

Linearized waveform inversion for vertical transversely isotropic elastic media: Methodology and multi-parameter crosstalk analysis

Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs. Seismic migration aims at obtaining structural images of subsurface geologic discontinuities. More specifically, seismic migration estimates the reflectivity function (stacked average reflectivity or pre-stack angle-dependent reflectivity) from seismic reflection data. On the other hand, seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations. Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters. Although there exist many differences, pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually, seismic migration and inversion techniques assume an acoustic or isotropic elastic medium. Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property. We present a linearized waveform inversion (LWI) scheme for weakly anisotropic elastic media with vertical transversely isotropic (VTI) symmetry. It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations (ΔVp0/Vp0,ΔVs0/Vs0,Δϵ,Δδ) from the long-wavelength reference model. Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration. We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method. The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images. We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters. Especially, the Thomsen parameter δ is the most difficult to resolve. We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.

An analytic model for the flow induced in syringomyelia cavities.

A simple two-dimensional fluid-structure-interaction problem, involving viscous oscillatory flow in a channel separated by an elastic membrane from a fluid-filled slender cavity, is analyzed to shed light on the flow dynamics pertaining to syringomyelia, a neurological disorder characterized by the appearance of a large tubular cavity (syrinx) within the spinal cord. The focus is on configurations in which the velocity induced in the cavity, representing the syrinx, is comparable to that found in the channel, representing the subarachnoid space surrounding the spinal cord, both flows being coupled through a linear elastic equation describing the membrane deformation. An asymptotic analysis for small stroke lengths leads to closed-form expressions for the leading-order oscillatory flow, and also for the stationary flow associated with the first-order corrections, the latter involving a steady distribution of transmembrane pressure. The magnitude of the induced flow is found to depend strongly on the frequency, with the result that for channel flow rates of non-sinusoidal waveform, as those found in the spinal canal, higher harmonics can dominate the sloshing motion in the cavity, in agreement with previous in vivo observations. Under some conditions, the cycle-averaged transmembrane pressure, also showing a marked dependence on the frequency, changes sign on increasing the cavity transverse dimension (i.e. orthogonal to the cord axis), underscoring the importance of cavity size in connection with the underlying hydrodynamics. The analytic results presented here can be instrumental in guiding future numerical investigations, needed to clarify the pathogenesis of syringomyelia cavities.

Full-wave Image Reconstruction in Transcranial Photoacoustic Computed Tomography using a Finite Element Method.

Transcranial photoacoustic computed tomography presents challenges in human brain imaging due to skull-induced acoustic aberration. Existing full-wave image reconstruction methods rely on a unified elastic wave equation for skull shear and longitudinal wave propagation, therefore demanding substantial computational resources. We propose an efficient discrete imaging model based on finite element discretization. The elastic wave equation for solids is solely applied to the hard-tissue skull region, while the soft-tissue or coupling-medium region that dominates the simulation domain is modeled with the simpler acoustic wave equation for liquids. The solid-liquid interfaces are explicitly modeled with elastic-acoustic coupling. Furthermore, finite element discretization allows coarser, irregular meshes to conform to object geometry. These factors significantly reduce the linear system size by 20 times to facilitate accurate whole-brain simulations with improved speed. We derive a matched forward-adjoint operator pair based on the model to enable integration with various optimization algorithms. We validate the reconstruction framework through numerical simulations and phantom experiments.

РАСЧЁТ НА ПРОЧНОСТЬ И ПРОВЕРКА ЖЁСТКОСТИ СТАТИЧЕСКИ ОПРЕДЕЛИМОЙ БАЛКИ ДВУТАВРОВОГО СЕЧЕНИЯ

This paper presents a strength calculation for a statically determinate I-beam. Diagrams of internal force factors (transverse forces and bending moments) are constructed and the number of the I-beam is selected according to the given external load and the values of the material strength characteristics. The given material characteristics are permissible normal stress, elastic modulus and permissible shear stress. The strength of a beam with a selected cross-section for shear stresses was tested. To construct the I-beam elastic line, the deflections of cross-sections with a unit interval along the length are calculated. The beam deflections were calculated using the universal elastic line equation. The obtained value of the maximum beam deflection made it possible to draw a conclusion about the beam sufficient hardness and operability under given external loads.

Solving overhead transmission line engineering problems with elastic catenary equations

Solving overhead transmission line engineering problems with elastic catenary equations

Numerical solution of unsteady elastic equations with C-Bézier basis functions

&lt;abstract&gt;&lt;p&gt;In this paper, the finite element method is applied to solve the unsteady elastic equations, C-Bézier basis functions are used to construct the shape function spaces, the semi-discrete scheme of the unsteady elastic equations is obtained by Galerkin finite element method and then the fully discretized Galerkin method is obtained by further discretizing the time variable with $ \theta $-scheme finite difference. Furthermore, for several numerical examples, the accuracy of approximate solutions are improved by 1–3 order-of magnitudes compared with the Lagrange basis function in $ L^\infty $ norm, $ L^2 $ norm and $ H^1 $ semi-norm, and the numerical examples show that the method proposed possesses a faster convergence rate. It is fully demonstrated that the C-Bézier basis functions have a better approximation effect in simulating unsteady elastic equations.&lt;/p&gt;&lt;/abstract&gt;

Blow up and lifespan of solutions for elastic membrane equation with delay

The primary objective of this research is to examine a nonlinear elastic membrane equation incorporating delay and source terms within a bounded domain. We obtain sufficient conditions on the initial data and the involved functionals for which the energy of solutions with non positive initial energy as well as positive initial energy blow up in a finite-time. In addition, this research work provides estimates for the lifespan of these solutions.

Optimal surrogate boundary selection and scalability studies for the shifted boundary method on octree meshes

The accurate and efficient simulation of Partial Differential Equations (PDEs) in and around arbitrarily defined geometries is critical for many application domains. Immersed boundary methods (IBMs) alleviate the usually laborious and time-consuming process of creating body-fitted meshes around complex geometry models (described by CAD or other representations, e.g., STL, point clouds), especially when high levels of mesh adaptivity are required. In this work, we advance the field of IBM in the context of the recently developed Shifted Boundary Method (SBM). In the SBM, the location where boundary conditions are enforced is shifted from the actual boundary of the immersed object to a nearby surrogate boundary, and boundary conditions are corrected utilizing Taylor expansions. This approach allows choosing surrogate boundaries that conform to a Cartesian mesh without losing accuracy or stability. Our contributions in this work are as follows: (a) we show that the SBM numerical error can be greatly reduced by an optimal choice of the surrogate boundary, (b) we mathematically prove the optimal convergence of the SBM for this optimal choice of the surrogate boundary, (c) we deploy the SBM on massively parallel octree meshes, including algorithmic advances to handle incomplete octrees, and (d) we showcase the applicability of these approaches with a wide variety of simulations involving complex shapes, sharp corners, and different topologies. Specific emphasis is given to Poisson’s equation and the linear elasticity equations.

Modelling and Investigation of the Dynamic Behavior of a Penny-Shaped Interface Crack in Piezoelectric Bimaterials

In this section, the dynamic propagation behavior of a penny-shaped interface crack in piezoelectric bimaterials is analyzed. The objective of this paper is to use the boundary conditions of the penny-shaped interface crack to study the dynamic propagation of the crack under the action of load, so as to provide some valuable implications for the fracture mechanics of the piezoelectric bimaterials and simulate the interface crack between piezoelectric bimaterials, it is necessary to establish a suitable model and give appropriate boundary conditions according to the actual situation. The elastic displacement and potential equations are constructed according to the structural characteristics of the circular crack. In the case of a given displacement or stress, the Laplace transform and Hankel transform are used to simplify the problem into an integral equation with unknown functions. According to the boundary conditions, the corresponding unknowns are obtained, and the closed solution is derived. The results show that the fracture toughness of a penny-shaped interface crack in piezoelectric bimaterials is related to the thickness of the material, the impact time, the material characteristics, and the electric field. At the same time, it can be found that different materials have different roles in the crack propagation, so it is very important to study the crack opening displacement (COD) intensity factor of the crack for safety design.

Eringen’s model via linearization of nonlocal hyperelasticity

We consider Riesz’ fractional gradient and a truncated version of it. The equations of nonlocal nonlinear elasticity based on those gradients are known. We perform a formal linearization and arrive at the equations of linear elasticity based on those nonlocal operators. We prove the existence of solutions of the linear equations, notably, by a nonlocal version of Korn’s inequality. Finally, we show that the linearizations obtained are particular cases of Eringen’s model with singular kernels.

Lie group analysis, solitons, self-adjointness and conservation laws of the nonlinear elastic structural element equation

This study is based on the Lie group method for the nonlinear elastic structural element equation (ESE Equation). We obtain a three-dimensional Lie algebra. By utilizing this Lie algebra a four-dimensional optimal system is constructed. The governing ESE Equation is converted to nonlinear ordinary differential equations (ODEs) via symmetry reduction. We use a modified auxiliary equation (MAE) procedure to deal with nonlinear ODEs. These ODEs reveal the dynamics of the periodic and soliton solutions. We obtain soliton solutions through rational, trigonometric, and hyperbolic functions. Wolfram Mathematica simulations vividly illustrate the wave characteristics of the derived solutions, affirming their properties as singular periodic solutions, a singular solution, an optical dark soliton solution, and a singular soliton solution. We also obtain the local conservation laws by a new conservation theorem introduced by Ibragimov.

An accelerated scheme for solving inhomogeneous elasticity in phase-field method and application to Ni-based multi-component alloys

On the way towards realistic predictions of various material properties for engineering structural materials, efficient computational simulations of microstructural evolution play an indispensable role. In the present work, an in-house phase-field code with inhomogeneous elasticity has been developed aiming to deal with precipitation processes for multi-component alloys with high efficiency. Two different numerical iterative algorithms, named as Model A and B, are implemented to solve inhomogeneous elastic equations. The resulting simulated morphology and energy evolutions under different conditions show consistency for these two algorithms. In order to perform simulation efficiently, a parallel technique based on CUDA (Computer Unified Device Architecture) architecture is used, which greatly improves the computational efficiency. After estimating the computational efficiency of different acceleration techniques for different grid sizes, Model A is extended to simulate concentration and morphology evolutions in the precipitation and rafting processes of γ′ phase in a Ni-based five-component alloy of Ni−Co−Al−Mo−W. Meanwhile, the CALPHAD method and related databases are linked to provide thermodynamic and diffusivity data. Finally, the computation costs of the parallel technique for different simulation scales are explored, testifying that this technique is capable of performing larger scale simulations for multi-component multi-phase structural materials.

Coupled mechano-electro-thermal model to predict phase transition in phase change materials

An advanced coupled Mechano-Electro-Thermal Phase Field model is presented to capture the crystallization/amorphization of phase change nanofilm under external voltage. It consisted of the Allen-Cahn equation for phase transition (PT), elasticity equations to resolve stresses, heat conduction, and the Poisson equation to consider thermoelectricity. The proposed model depicts the intricate interplay between temperature, stress, pre-existing defects (surface void), and PT in model phase change material, i.e., Germanium-antimony-tellurium (GeSbTe or GST). The predicted stress generation, thermal and PT-induced, agrees with existing experimental data. Additionally, surface voids can cause remarkable temperature and PT kinetics changes, decreasing the transformation rate. Also, due to a highly heterogeneous temperature, the thermal strain shows a crucial effect on crystallization kinetics. Finally, the critical stepwise voltage (Vcr) is calculated, representing an inverse power relation with the pulse time, i.e., the maximum voltage results in complete crystallization below the melting temperature. An initial increase in the average temperature occurs at the critical voltage, followed by a constant temperature regime as crystalization evolves and a jump to higher temperatures until complete crystallization. However, the temperature rapidly jumps above the melting temperature for larger voltages (V > Vcr), leading to the reverse transformation, i.e., amorphization.

Spatial-temporal high-order rotated-staggered-grid finite-difference scheme of elastic wave equations for TTI medium

In the finite-difference (FD) solution of the elastic wave equation for TTI medium, the traditional FD scheme can improve the computational accuracy of the spatial derivatives. However, the difference discretization on the temporal partial derivatives still uses the FD scheme of second-order accuracy, which cannot guarantee the computational accuracy of the temporal partial derivatives and often causes the temporal dispersion problem. The temporal high-order FD scheme can improve the computational accuracy of the temporal partial derivatives. Nevertheless, due to the complexity of the elastic wave equation in TTI medium, using the conventional staggered-grid technique to solve the elastic wave equation, some spatial derivatives need to be approximated by spatial interpolation, which will become new error and reduce the accuracy of the solution of the equations. Therefore, the currently commonly used temporal high-order FD method for the elastic wave equation in isotropic medium cannot be directly applied to the simulation of anisotropic medium. To solve the above problems, we propose a new temporal high-order rotated-staggered-grid FD stencil and derive a new temporal high-order FD scheme under this new stencil for the TTI medium elastic wave equation. Then we estimated the FD coefficients of the new stencil using Taylor series expansions (TE) and optimized the FD coefficients based on least squares (LS). Numerical experiments of the homogeneous TTI medium model and the modified BP model demonstrate that the TE + LS-based temporal high-order rotated-staggered-grid FD scheme can well solve the temporal dispersion in the numerical simulation of TTI medium.