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↗
NASA NTRS2011

An In-Depth Tutorial on Constitutive Equations for Elastic Anisotropic Materials

An in-depth tutorial on the constitutive equations for elastic, anisotropic materials is presented. Basic concepts are introduced that are used to characterize materials, and notions about how anisotropic material deform are presented. Hooke s law and the Duhamel-Neuman law for isotropic materials are presented and discussed. Then, the most general form of Hooke s law for elastic anisotropic materials is presented and symmetry requirements are given. A similar presentation is also given for the generalized Duhamel-Neuman law for elastic, anisotropic materials that includes thermal effects. Transformation equations for stress and strains are presented and the most general form of the transformation equations for the constitutive matrices are given. Then, specialized transformation equations are presented for dextral rotations about the coordinate axes. Next, concepts of material symmetry are introduced and criteria for material symmetries are presented. Additionally, engineering constants of fully anisotropic, elastic materials are derived from first principles and the specialized to several cases of practical importance.

Nemeth, Michael P.↗
NASA NTRS1983

Thermal elastohydrodynamic lubrication of line contacts

A numerical solution to the problem of thermal elastohydrodynamic lubrication of line contacts was obtained by using a finite difference formulation. The solution procedure consists of simultaneous solution of the thermal Reynolds equation, the elasticity equation, and the energy equation subject to appropriate boundary conditions. Pressure distribution, film shape, and temperature distribution were obtained for fully flooded conjunctions, a paraffinic lubricant, and various dimensionless speed parameters while the dimensionless load and materials parameters were held constant. Reduction in the minimum film thickness due to thermal effects (as a ratio of thermal to isothermal minimum film thickness) is given by a simple formula as a function of the thermal loading parameter Q: H(min)/H(min,I) = 10/10+ Q(0.4). Plots of pressure distribution, film shape, temperature distribution, and flow are shown for some representative cases.

Ghosh, M. K.↗
NASA NTRS1985

Thermal elastohydrodynamic lubrication of line contacts

A numerical solution to the problem of thermal elastohydrodynamic lubrication of line contacts was obtained by using a finite difference formulation. The solution procedure consists of simultaneous solution of the thermal Reynolds equation, the elasticity equation, and the energy equation subject to appropriate boundary conditions. Pressure distribution, film shape, and temperature distribution were obtained for fully flooded conjunctions, a paraffinic lubricant, and various dimensionless speed parameters while the dimensionless load and materials parameters were held constant. Reduction in the minimum film thickness due to thermal effects (as a ratio of thermal to isothermal minimum film thickness) is given by a simple formula as a function of the thermal loading parameter Q: H(min)/H(min, I) = 10/10 + Q(0.4). Plots of pressure distribution, film shape, temperature distribution, and flow are shown for some representative cases.

Ghosh, M. K.↗
NASA NTRS1992

Acoustic intensity calculations for axisymmetrically modeled fluid regions

An algorithm for calculating acoustic intensities from a time harmonic pressure field in an axisymmetric fluid region is presented. Acoustic pressures are computed in a mesh of NASTRAN triangular finite elements of revolution (TRIAAX) using an analogy between the scalar wave equation and elasticity equations. Acoustic intensities are then calculated from pressures and pressure derivatives taken over the mesh of TRIAAX elements. Intensities are displayed as vectors indicating the directions and magnitudes of energy flow at all mesh points in the acoustic field. A prolate spheroidal shell is modeled with axisymmetric shell elements (CONEAX) and submerged in a fluid region of TRIAAX elements. The model is analyzed to illustrate the acoustic intensity method and the usefulness of energy flow paths in the understanding of the response of fluid-structure interaction problems. The structural-acoustic analogy used is summarized for completeness. This study uncovered a NASTRAN limitation involving numerical precision issues in the CONEAX stiffness calculation causing large errors in the system matrices for nearly cylindrical cones.

Hambric, Stephen A.↗
NASA NTRS1986

Multigrid method for the equilibrium equations of elasticity using a compact scheme

A compact difference scheme is derived for treating the equilibrium equations of elasticity. The scheme is inconsistent and unstable. A multigrid method which takes into account these properties is described. The solution of the discrete equations, up to the level of discretization errors, is obtained by this method in just two multigrid cycles.

Taasan, S.↗
NASA NTRS1987

Multigrid method for the equilibrium equations of elasticity using a compact scheme

A compact difference scheme is derived for treating the equilibrium equations of elasticity. The scheme is inconsistent and unstable. A multigrid method which takes into account these properties is described. The solution of the discrete equations, up to the level of discretization errors, is obtained by this method in just two multigrid cycles.

Taasan, Shlomo↗
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↗
NASA NTRS1973

Effect of surface asperity on elastohydrodynamic lubrication

The important aspects of elastohydrodynamic lubrication, with a single, one-dimensional asperity, have been found by solving numerically the coupled transient Reynolds equation and the elasticity equation. Even though the assumption of a single asperity is highly ideal, this study sheds some light on the effect of surface roughness on elastohydrodynamic lubrication. The results show that the film pressure tends to increase more than the steady state pressure, and in particular, the increase in pressure reaches a maximum as the asperity approaches the inlet of the contact region. The asperity height and the pressure increase above the steady state pressure are closely related to each other; the higher the asperity height, the larger the pressure increase. In the pure rolling case, it has been found that a local pressure peak is not developed. However, in the cases of sliding and rolling, a small, local pressure peak is developed on the pressure profile when the asperity moves into the contact region. In general, the overall film thickness profile increases with increasing asperity height, but is not significantly affected by the asperity width. Moreover, the slope of the overall film thickness profile for the transient cases is much greater than the steady state profile, which is approximately constant across the contact width. The increase in the center film thickness also depends upon the width and height of the asperity.

Lee, K.↗
NASA NTRS1985

Elastohydrodynamic lubrication calculations used as a tool to study scuffing

A new Reynolds equation was developed that takes into account the nonlinear viscous behavior of the fluid. The new Reynolds equation considers the nonlinear viscous fluid model of Eyring, the equilibrium equation, the constant mass flow, and the kinematic boundary condition. The new Reynolds equation and the elasticity equation are solved simultaneously by using a system approach and a Newton-Raphson technique. Comparisons are made with results obtained from the classical Reynolds equation. The effects of sliding speed and introducing a bump or a groove within the conjunction are studied. Results are shown for both moderate and heavy loads.

Houpert, L. G.↗
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↗