Search NASA⌕ Search

SEARCH · Search NASA

Results for “Elasticity equations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

10-th order of accuracy for numerical solution of 3-D elasticity equations for heterogeneous materials on unfitted Cartesian meshes

We have developed the Optimal Local Truncation Error Method (OLTEM) with 10-th order of accuracy on unfitted Cartesian meshes for a system of 3-D elasticity equations with smooth irregular interfaces. 5 x 5 x 5 = 125-point stencils (similar to those for quadratic finite elements) for elastic heterogeneous materials are used for OLTEM. There are no unknowns at the interface points between different materials; the structure of the global discrete equations is the same for homogeneous and heterogeneous materials. The calculation of unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal 10-th order of accuracy for OLTEM on unfitted Cartesian meshes, i.e., the increase by 7 orders in accuracy compared to quadratic finite elements on conformal meshes. A new post-processing procedure provides the 9-th order of accuracy for stresses in the 3-D case. Similar to basic computations it uses OLTEM with the 125-point stencils, the interface conditions and the elasticity equations. It was shown that the use of the elasticity equations for post-processing improves the accuracy of 0.1% stresses by 6 orders compared to post-processing without the use of PDEs. At an accuracy of for stresses, OLTEM with the new post-processing procedure reduces the number of degrees of freedom by 360 - 8000 times compared to quadratic finite elements with similar stencils. OLTEM with the 125-point stencils yields even more accurate results than high-order finite elements with much wider stencils. OLTEM provides accurate numerical results for compressible and nearly incompressible materials.

elasticity equations↗
DOE OSTI2018

Convergence of the multigrid reduction in time algorithm for the linear elasticity equations

This paper presents some recent advances for parallel-in-time methods applied to linear elasticity. With recent computer architecture changes leading to stagnant clock speeds, but ever increasing numbers of cores, future speedups will be available through increased concurrency. Thus, sequential algorithms, such as time stepping, will suffer a bottleneck. This paper explores multigrid reduction in time (MGRIT) for an important application area, linear elasticity. Previously, efforts at parallel-in-time for elasticity have experienced difficulties, for example, the beating phenomenon. As a result, practical parallel-in-time algorithms for this application area currently do not exist. This paper proposes some solutions made possible by MGRIT (e.g., slow temporal coarsening and FCF-relaxation) and, more importantly, a different formulation of the problem that is more amenable to parallel-in-time methods. Using a recently developed convergence theory for MGRIT and Parareal, we show that the changed formulation of the problem avoids the instability issues and allows the reduction of the error using two temporal grids. We then extend our approach to the multilevel case, where we demonstrate how slow temporal coarsening improves convergence. Finally, the paper ends with supporting numerical results showing a practical algorithm enjoying speedup benefits over the sequential algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Super-Grid Modeling of the Elastic Wave Equation in Semi-Bounded Domains

Abstract We develop a super-grid modeling technique for solving the elastic wave equation in semi-bounded two- and three-dimensional spatial domains. In this method, waves are slowed down and dissipated in sponge layers near the far-field boundaries. Mathematically, this is equivalent to a coordinate mapping that transforms a very large physical domain to a significantly smaller computational domain, where the elastic wave equation is solved numerically on a regular grid. To damp out waves that become poorly resolved because of the coordinate mapping, a high order artificial dissipation operator is added in layers near the boundaries of the computational domain. We prove by energy estimates that the super-grid modeling leads to a stable numerical method with decreasing energy, which is valid for heterogeneous material properties and a free surface boundary condition on one side of the domain. Our spatial discretization is based on a fourth order accurate finite difference method, which satisfies the principle of summation by parts. We show that the discrete energy estimate holds also when a centered finite difference stencil is combined with homogeneous Dirichlet conditions at several ghost points outside of the far-field boundaries. Therefore, the coefficients in the finite difference stencils need only be boundary modified near the free surface. This allows for improved computational efficiency and significant simplifications of the implementation of the proposed method in multi-dimensional domains. Numerical experiments in three space dimensions show that the modeling error from truncating the domain can be made very small by choosing a sufficiently wide super-grid damping layer. The numerical accuracy is first evaluated against analytical solutions of Lamb’s problem, where fourth order accuracy is observed with a sixth order artificial dissipation. We then use successive grid refinements to study the numerical accuracy in the more complicated motion due to a point moment tensor source in a regularized layered material.

58 GEOSCIENCES↗
DOE OSTI2018

Comparison of pseudo-static limit equilibrium and elastic wave equation analyses of dynamic earth pressures on retaining structures

The seismic earth pressure increment is typically computed using either pseudo-static limit equilibrium methods or elastic wave equation analyses of the interaction between a retaining structure and backfill material, yet current interpretations of the two methods provide conflicting recommendations. The focus of this study is to compare the seismic earth pressure increment computed using the two methods. This approach is demonstrated by subjecting an initially uniform prototype site selected from standard site classifications to harmonic excitation in one-dimensional equivalent linear analyses. Then, the seismic earth pressure resultant for a rigid wall is computed using the two methods. The limit equilibrium approach utilizes the acceleration records from the equivalent linear analysis to compute a seismic coefficient, whereas the elastic solution incorporates the reduced modulus and damping from the final iteration of the analysis, as well as the relative displacement records. The results presented herein corroborate the findings of recent centrifuge experiments and associated analyses.

58 GEOSCIENCES↗
DOE OSTI2025

High-pressure elasticity and equation of state of the fluoroelastomer Viton® A-500

Viton® A is a semi-crystalline copolymer of polyvinylidene fluoride and hexafluoropropylene used in various engineering applications due to its mechanical properties and chemical inertness. In situ ultrasonic spectroscopy and x-ray radiography measurements were performed in a Paris–Edinburgh press to measure the pressure dependence of the transverse and longitudinal acoustic velocity of the fluoroelastomer A-500 from 2.7 to 5.7 GPa at 296 K. In addition, we performed high-pressure Brillouin scattering measurements to obtain acoustic velocities from ambient pressure to 5.7 GPa to supplement the ultrasonic measurements, especially at low pressures. The acoustic velocities were then used to calculate a pressure–volume (P–V) equation of state, the bulk and shear moduli, and the Poisson's ratio. These quantities are compared with the reported pressure-dependent properties of related polymers over this range of pressures.

Acoustical properties↗
DOE OSTI2023

Acoustic properties, elasticity, and equation of state of glycerol under pressure

Herein we employed high-pressure Brillouin scattering to study the pressure dependencies of acoustic modes of glycerol up to 14 GPa at 300 K. We observed longitudinal acoustic velocities and transverse acoustic velocities for the first time from 5 to 14 GPa. The results allow the determination of a complete set of elastic properties and an accurate determination of the pressure–volume (P–V) equation of state (EOS). EOS parameters, K 0 = 14.9 ± 1.8 GPa and K' 0 = 5.6 ± 0.5, were determined from fits to the data from ambient pressure to 14 GPa. Direct volume measurements of the P–V EOS are consistent with those determined by Brillouin scattering. A deviation from a Cauchy-like relationship for elastic properties was observed, and the pressure dependencies of the photoelastic constants and relaxation times were documented from 5 to 14 GPa. These results have broad implications for glass-forming liquids, viscoelastic theory, and mode coupling theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Optimal local truncation error method for 3-D elasticity interface problems

The paper deals with a new effective numerical technique on unfitted Cartesian meshes for simulations of heterogeneous elastic materials. Here, we develop the optimal local truncation error method (OLTEM) with 27- point stencils (similar to those for linear finite elements) for the 3-D time-independent elasticity equations with irregular interfaces. Only displacement unknowns at each internal Cartesian grid point are used. The interface conditions are added to the expression for the local truncation error and do not change the width of the stencils. The unknown stencil coefficients are calculated by the minimization of the local truncation error of the stencil equations and yield the optimal second order of accuracy for OLTEM with the 27-point stencils on unfitted Cartesian meshes. A new post-processing procedure for accurate stress calculations has been developed. Similar to basic computations it uses OLTEM with the 27-point stencils and the elasticity equations. The post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing numerical techniques (e.g., with finite elements). Numerical experiments show that at an accuracy of 0.1% for stresses, OLTEM with the new post-processing procedure significantly (by 10 5 -10 9 times) reduces the number of degrees of freedom compared to linear finite elements. OLTEM with the 27-point stencils yields even more accurate results than high-order finite elements with wider stencils.

42 ENGINEERING↗
DOE OSTI2025

Finite domain solution of a hydraulic fracture in a permeable rock

In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media.

Domain-based method↗
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

The Enskog Equation for Confined Elastic Hard Spheres

A kinetic equation for a system of elastic hard spheres or disks confined by a hard wall of arbitrary shape is derived. It is a generalization of the modified Enskog equation in which the effects of the confinement are taken into account and it is supposed to be valid up to moderate densities. From the equation, balance equations for the hydrodynamic fields are derived, identifying the collisional transfer contributions to the pressure tensor and heat flux. A Lyapunov functional, H[f], is identified. For any solution of the kinetic equation, H decays monotonically in time until the system reaches the inhomogeneous equilibrium distribution, that is a Maxwellian distribution with a density field consistent with equilibrium statistical mechanics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

AMPE

AMPE is a phase-field code. It evolves in time the Allen-Cahn equations associated with a phase­field model, coupled with the time-evolution of a composition field (in an alloy) and a temperature field which are described by diffusion equations. It can also solve the elasticity equations assuming mechanical equilibrium at every time-step. Grain orientations are described by normalized quaternions which also evolve in time according to Allen-Cahn type equations.

Wickett, Mike↗
DOE OSTI2019

Simulation of unsteady blood flows in a patient-specific compliant pulmonary artery with a highly parallel monolithically coupled fluid-structure interaction algorithm

Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient-specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems, but need improvements when the problems become harder in the case when (1) the geometry becomes complex (from a few branches to a full pulmonary artery), (2) the model becomes more complex (from fluid-only calculation to coupled fluid-structure interaction calculation), (3) both the fluid and wall models become highly nonlinear, and (4) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid-structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient-specific, full size pulmonary artery obtained from CT (Computed Tomography) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier-Stokes equations and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. We report the proposed numerical algorithm and software scale well beyond 10,000 processor cores on a supercomputer for solving the fluid-structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Study on 3D Internal Magnetic Field Distribution and Dynamic Mechanics of a Giant Magnetostrictive Actuator

The three-dimensional (3D) coupled magneto-elastic dynamic model of a giant magnetostrictive actuator (GMA) was built based on Maxwell’s equations, elasticity, and vibration theories. The general analytical expressions for the 3D internal magnetic field distribution of GMA were derived, and the corresponding analytical solutions were obtained. This analytical model is able to predict the internal magnetic field distribution and the output displacement of GMA system with different physical parameters under external excitations. The dynamic behaviors, such as the amplitude vs. frequency responses and the time domain waveforms of the GMA, were also studied. Analytical solutions showed that the internal magnetic field has a prominent end effect along the axial direction of a giant magnetostrictive material (GMM) rod, whereas a minimal skin effect along its radial direction, which are matched very well with the FEM results. The amplitude-frequency relationship exhibited prominent resonance and anti-resonance characteristic, which is a typical hysteresis feature of the GMA system. In order to testify the validity of the model prediction, a GMA testing system was also established. Excellent agreement was found between the experiment measurements and model predictions in the range of our research (frequency is lower than 300 Hz). The analytical results of the present study have potential applications in the fields of active vibration control and other precision machining of the GMA system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Effects of internal stresses and intermediate phases on the coarsening of coherent precipitates: A phase-field study

Phase stability, topology and size evolution of precipitates are important factors in determining the mechanical properties of crystalline materials. In this article, the Cahn-Hilliard type of phase-field model was coupled to elasticity equations within a mixed-order Galerkin finite element framework to study the coarsening morphology of coherent precipitates. The effects of capillarity, particle size and fraction, compositional strain, and inhomogeneous elasticity on the kinetics and kinematics of coherent precipitates in a binary dual phase crystal admitting a third intermediate stable/meta-stable phase were investigated. The results demonstrated the ability of the model to simulate coarsening under the concomitant action of Ostwald ripening and mismatch elastic strain mechanisms. Using a phenomenological coarsening power law, coarsening rates were determined to depend on precipitate size and volume fraction, compositional strain, and strain mismatch between precipitates and the matrix. Results also showed that the necking incubation time between two neighboring precipitates depends inversely on the precipitate's initial sizes; however, under fixed volume fraction of precipitates, any increase in the initial sizes of the precipitates mitigates the coarsening. Meanwhile, the compositional strain and the growth of the intermediate stable/meta-stable phase leads to substantial enhancements of precipitate coarsening.

36 MATERIALS SCIENCE↗