SEARCH · Engineering Papers

Results for “elastic fields”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Mechanical, electronic, optical, piezoelectric and ferroic properties of strained graphene and other strained monolayers and multilayers: an update

Abstract This is an update of a previous review (Naumis et al 2017 Rep. Prog. Phys. 80 096501). Experimental and theoretical advances for straining graphene and other metallic, insulating, ferroelectric, ferroelastic, ferromagnetic and multiferroic 2D materials were considered. We surveyed (i) methods to induce valley and sublattice polarisation ( P ) in graphene, (ii) time-dependent strain and its impact on graphene’s electronic properties, (iii) the role of local and global strain on superconductivity and other highly correlated and/or topological phases of graphene, (iv) inducing polarisation P on hexagonal boron nitride monolayers via strain, (v) modifying the optoelectronic properties of transition metal dichalcogenide monolayers through strain, (vi) ferroic 2D materials with intrinsic elastic ( σ ), electric ( P ) and magnetic ( M ) polarisation under strain, as well as incipient 2D multiferroics and (vii) moiré bilayers exhibiting flat electronic bands and exotic quantum phase diagrams, and other bilayer or few-layer systems exhibiting ferroic orders tunable by rotations and shear strain. The update features the experimental realisations of a tunable two-dimensional Quantum Spin Hall effect in germanene, of elemental 2D ferroelectric bismuth, and 2D multiferroic NiI 2 . The document was structured for a discussion of effects taking place in monolayers first, followed by discussions concerning bilayers and few-layers, and it represents an up-to-date overview of exciting and newest developments on the fast-paced field of 2D materials.

Physics↗

DOE OSTI2020

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Mercury Crustal and Lithospheric Properties from High-resolution Gravity Field Models

We have analyzed all the available tracking data from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) mission to determine high-resolution gravity field models. We have applied a technique based on line-of-sight accelerations that results in models with improved correlations with topography in areas where the spacecraft altitude was low. We have used these models in an analysis of admittance to determine properties of the crust and lithosphere. We focus on four areas where correlations between gravity and topography are high. These locations represent different areas on Mercury: the high-Mg region, the Strindberg crater plus some lobate scarps, heavily cratered terrain, and smooth plains (the northern rise). We find crustal densities around 2600 kg m-3 for the first three areas, and a higher density for the northern rise, indicative of volcanic materials. Elastic thickness is generally low and varies between 11 and 30 km.

Sander Johannes Goossens↗

DOE OSTI2021

Inelastic peridynamic model for molecular crystal particles

The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. Finally, the material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

A Validated Finite Element Modeling Tool for Hydrodynamic Loading and Structural Analysis of Ocean Deployed Macroalgae Farms

The objective of this project was to develop and validate a fine-tuned, 3D computational modeling tool for high-fidelity simulations of macroalgae cultivation and harvest systems. To develop this tool, a full understanding of the geometric, material and hydrodynamic properties associated with farming kelp in exposed environments was necessary. Experiments and numerical modeling were done at increasing scales from 1) tank tests with a 1 m model, 2) to field tests with a 61 m culture line at a sheltered site, 3) to field tests with a 122 m culture line at an exposed site and to 4) a 5- line array at the exposed site with 186 m of culture length. The results of the experimental tank tests yielded a set of hydrodynamic drag-area coefficients for a 1 m aggregate based on densely grown kelp at 16 kg/m with a length up to 3 m as described in Fredriksson et al. This work provided the basis to model kelp as a 1 m aggregate based on in-situ geometric and material properties to include measurements of kelp blade length and width, number of blades, yield per m, material mass density and cantilever tests to obtain modulus of elasticity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗

DOE OSTI2021

First principles modeling of polymer encapsulant degradation in Si photovoltaic modules

An outstanding issue in the longevity of photovoltaic (PV) modules is the accelerated degradation caused by the presence of moisture. Moisture leads to interfacial instability, de-adhesion, encapsulant decomposition, and contact corrosion. However, experimental characterization of moisture in PV modules is not trivial and its impacts can take years or decades to establish in the field, presenting a major obstacle to designing high-reliability modules. First principles calculations provide an alternative way to study the ingress of water and its detrimental effect on the structure and decomposition of the polymer encapsulant and interfaces between the encapsulant and the semiconductor, the metal contacts, or the dielectric layer. Here, we use density functional theory (DFT) computations to model single chain, crystalline and cross-linked structures, infrared (IR) signatures, and degradation mechanisms of ethylene vinyl acetate (EVA), the most common polymer encapsulant used in Si PV modules. IR-active modes computed for low energy EVA structures and possible decomposition products match well with reported experiments. The EVA decomposition energy barriers computed using the Nudged Elastic Band (NEB) method show a preference for acetic acid formation as compared to acetaldehyde, are lowered in the presence of a water solvent or hydroxyl ion catalyst, and match well with reported experimental activation energies. This systematic study leads to a clear picture of the hydrolysis-driven decomposition of EVA in terms of energetically favorable mechanisms, possible intermediate structures, and IR signatures of reactants and products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2017

A Review of the Proposed K (sub Isi) Offset-Secant Method for Size-Independent Linear-Elastic Toughness Evaluation

The proposed size-independent linear-elastic fracture toughness, K (sub Isi), for potential inclusion in ASTM E399 targets a consistent 0.5 millimeters crack extension for all specimen sizes through an offset secant that is a function of the specimen ligament length. The K (sub Isi) method also includes an increase in allowable deformation, and the removal of the P (sub max)/P (sub Q) criterion. A finite element study of the K (sub Isi) test method confirms the viability of the increased deformation limit, but has also revealed a few areas of concern. Findings: 1. The deformation limit, b (sub o) greater than or equal to 1.1 times (K (sub I) divided by delta (sub ys) squared) maintains a K-dominant crack tip field with limited plastic contribution to the fracture energy; 2. The three dimensional effects on compliance and the shape of the force versus CMOD (Crack-Mouth Opening Displacement) trace are significant compared to a plane strain assumption; 3. The non-linearity in the force versus CMOD trace at deformations higher than the current limit of 2.5 times (K (sub I) divided by delta (sub ys) squared) is sufficient to introduce error or even "false calls" regarding crack extension when using a constant offset secant line. This issue is more significant for specimens with W (width) greater than or equal to 2 inches; 4. A non-linear plasticity correction factor in the offset secant may improve the viability of the method at deformations between 2.5 times (K (sub I) divided by delta (sub ys) squared) and 1.1 times (K (sub I) divided by delta (sub ys) squared).

James, Mark↗

DOE OSTI2022

Bootstrapping 2d ϕ 4 theory with Hamiltonian truncation data

We combine the methods of Hamiltonian Truncation and the recently proposed generalisation of the S-matrix bootstrap that includes local operators to determine the two-particle scattering amplitude and the two-particle form factor of the stress tensor at s > 0 in the 2d Φ 4 theory. We use the form factor of the stress tensor at s ≤ 0 and its spectral density computed using Lightcone Conformal Truncation (LCT), and inject them into the generalized S-matrix bootstrap set-up. The obtained results for the scattering amplitude and the form factor are fully reliable only in the elastic regime. We independently construct the “pure” S-matrix bootstrap bounds (bootstrap without including matrix elements of local operators), and find that the sinh-Gordon model and its analytic continuation the “staircase model” saturate these bounds. Surprisingly, the Φ 4 two-particle scattering amplitude also very nearly saturates these bounds, and moreover is extremely close to that of the sinh-Gordon/staircase model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗

NASA NTRS2004

Nonlinear Analysis of Bonded Composite Single-LAP Joints

This study presents a semi-analytical solution method to analyze the geometrically nonlinear response of bonded composite single-lap joints with tapered adherend edges under uniaxial tension. The solution method provides the transverse shear and normal stresses in the adhesive and in-plane stress resultants and bending moments in the adherends. The method utilizes the principle of virtual work in conjunction with von Karman s nonlinear plate theory to model the adherends and the shear lag model to represent the kinematics of the thin adhesive layer between the adherends. Furthermore, the method accounts for the bilinear elastic material behavior of the adhesive while maintaining a linear stress-strain relationship in the adherends. In order to account for the stiffness changes due to thickness variation of the adherends along the tapered edges, their in-plane and bending stiffness matrices are varied as a function of thickness along the tapered region. The combination of these complexities results in a system of nonlinear governing equilibrium equations. This approach represents a computationally efficient alternative to finite element method. Comparisons are made with corresponding results obtained from finite-element analysis. The results confirm the validity of the solution method. The numerical results present the effects of taper angle, adherend overlap length, and the bilinear adhesive material on the stress fields in the adherends, as well as the adhesive, of a single-lap joint

Oterkus, E.↗

NASA NTRS2011

Validation of Methods to Predict Vibration of a Panel in the Near Field of a Hot Supersonic Rocket Plume

This paper describes the measurement and analysis of surface fluctuating pressure level (FPL) data and vibration data from a plume impingement aero-acoustic and vibration (PIAAV) test to validate NASA s physics-based modeling methods for prediction of panel vibration in the near field of a hot supersonic rocket plume. For this test - reported more fully in a companion paper by Osterholt & Knox at 26th Aerospace Testing Seminar, 2011 - the flexible panel was located 2.4 nozzle diameters from the plume centerline and 4.3 nozzle diameters downstream from the nozzle exit. The FPL loading is analyzed in terms of its auto spectrum, its cross spectrum, its spatial correlation parameters and its statistical properties. The panel vibration data is used to estimate the in-situ damping under plume FPL loading conditions and to validate both finite element analysis (FEA) and statistical energy analysis (SEA) methods for prediction of panel response. An assessment is also made of the effects of non-linearity in the panel elasticity.

Bremner, P. G.↗

NASA NTRS1990

Experimental flaw detection by scattering of plate waves

The surface displacements of elastic waves scattered from a surface-breaking crack are measured with a capacitive transducer on a 2D plate, and comparisons with theoretical predictions are reported. The noncontacting detector has a wide and flat frequency response which measures the response of a glass plate to the propagation of Rayleigh-Lamb waves. A theoretical simulation of the surface displacements is generated by multiplying experimental traces with a cosine window, and Green's function is used to compute the source function. Power spectra are also recorded from the experiment to compare the results with a finite-element model of the crack, and good agreement is reported for results in both the time and frequency domains. The crack can be detected in the far field when the spacing of the receiver locations is small, and the transducer/probe appears to be a practical flaw-detection instrument.

Spetzler, Hartmut↗

NASA NTRS1994

Interlaminar shear test method development for long term durability testing of composites

The high speed civil transport is a commercial aircraft that is expected to carry 300 passengers at Mach 2.4 over a range of more than 6000 nautical miles. With the existing commercial structural material technology (i.e., aluminum) the performance characteristics of the high speed civil transport would not be realized. Therefore there has been a concerted effort in the development of light weight materials capable of withstanding elevated temperatures for long duration. Thermoplastic composite materials are such candidate materials and the understanding of how these materials perform over the long term under harsh environments is essential to safe and effective design. The matrix dominated properties of thermoplastic composites are most affected by both time and temperature. There is currently an effort to perform short term testing to predict long term behavior of in-plane mechanical properties E22 (transverse modulus of elasticity) and G12 (shear modulus). Out-of-plane properties such as E33, G13, and G23 are inherently more difficult to characterize. This is especially true for the out-of-plane shear modulus G23 and hence there is no existing acceptable standard test method. Since G23 is the most matrix dominated property, it is essential that a test method be developed. A shear test methodology is developed to do just that. The test method, called the double notched specimen, along with the previously developed shear gage was tested at room temperature. Mechanical testing confirmed the attributes of the methodology. A finite element parametric study was conducted for specimen optimization. Moire interferometry, a high sensitivity laser optical method, was used for full-field analysis of the specimen. From this work, material parameters will be determined and thus enable the prediction of long term material behavior of laminates subjected to general loading states.

Ifju, Peter G.↗

NASA NTRS2024

Out of Plane Stresses Using a Beam-Based Elastic Foundation Bonded Joint Model

Efficient, semi-analytical models of bonded joints are utilized for initial joint design and sizing due to their convenience and speed. Many of these models are based on beam and plate theories, which are most applicable when in-plane stresses/strains dominate. While this can be a sensible assumption for obtaining the deformations of the joint, out of plane (transverse) stresses are often of critical importance for failure in bonded joints composed of fiber reinforced composite materials. Since through-thickness strains are neglected in these semi-analytical theories, the present study shows a method to back-calculate the out of plane stresses based on stress balance equations, which use the in-plane stresses and displacements that have already been obtained. The method was implemented within the Joint Element Designer (JED) software, and verification and validation examples have been shown. As is typical for these types of models, corner stress singularities are not captured by the JED model, but the general trends of the stress and deformation fields agreed with a higher fidelity finite element model. As such, this method for obtaining out of plane stresses, when coupled with the JED code, is very useful for the complete stress analysis of bonded composite joints.

composites↗

DOE OSTI2021

Enhanced relaxed physical factorization preconditioner for coupled poromechanics

The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled poromechanics. For its application, however, it is necessary to invert blocks with the algebraic form C^ = (C + βFF T ), where C is a symmetric positive definite matrix, FF T a rank-deficient term, and β a real non-negative coefficient. The inversion of C^, performed in an inexact way, can become unstable for large values of β, as it usually occurs at some stages of a full poromechanical simulation. In this work, we propose a family of algebraic techniques to stabilize the inexact solve with C^. This strategy can prove useful in other problems as well where such an issue might arise, such as augmented Lagrangian preconditioning techniques for Navier-Stokes or incompressible elasticity. First, we introduce an iterative scheme obtained by a natural splitting of matrix C^. Second, we develop a technique based on the use of a proper projection operator annihilating the near-kernel modes of C^. Both approaches give rise to a novel class of preconditioners denoted as Enhanced RPF (ERPF). Furthermore, effectiveness and robustness of the proposed algorithms are demonstrated in both theoretical benchmarks and real-world large-size applications, outperforming the native RPF preconditioner.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Disentangling Coexisting Structural Order Through Phase Lock-In Analysis of Atomic-Resolution STEM Data

As a real-space technique, atomic-resolution STEM imaging contains both amplitude and geometric phase information about structural order in materials, with the latter encoding important information about local variations and heterogeneities present in crystalline lattices. Such phase information can be extracted using geometric phase analysis (GPA), a method which has generally focused on spatially mapping elastic strain. Here we demonstrate an alternative phase demodulation technique and its application to reveal complex structural phenomena in correlated quantum materials. As with other methods of image phase analysis, the phase lock-in approach can be implemented to extract detailed information about structural order and disorder, including dislocations and compound defects in crystals. Extending the application of this phase analysis to Fourier components that encode periodic modulations of the crystalline lattice, such as superlattice or secondary frequency peaks, we extract the behavior of multiple distinct order parameters within the same image, yielding insights into not only the crystalline heterogeneity but also subtle emergent order parameters such as antipolar displacements. When applied to atomic-resolution images spanning large (~0.5 × 0.5 μm 2 ) fields of view, this approach enables vivid visualizations of the spatial interplay between various structural orders in novel materials.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗