Search NASA⌕ Search

SEARCH · Search NASA

Results for “Displacement”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

The North-South Asymmetry of the Heliospheric Current Sheet: Results of an MHD Simulation

A displacement of the heliospheric current sheet (HCS) south of the helioequator by approx.10deg was proposed by Simpson et al. (1996) as a possible explanation of the north-south asymmetry in the galactic cosmic rays observed by Ulysses during its first fast transit in 1994-1995. The idea was not supported by magnetic field measurements on Ulysses and, on this ground, was dismissed by Simpson et al. (1996). In addition, Erdos & Balogh (1998) argued that any north-south symmetry was unlikely as there should be flux balance between the magnetic sectors of opposite polarity. Nonetheless, many in the scientific community have accepted the original suggestion of Simpson et al. (1996) that a displacement of the HCS was responsible for the cosmic ray asymmetry. In this paper, using a magnetohydrodynamic model of the solar corona and solar wind that includes both dipole and quadrupole magnetic source terms, we show that a north-south asymmetry of the magnetic field on the Sun does not give rise to a displacement of the HCS. The lack of displacement of the HCS results from a latitudinal redistribution of magnetic flux near the Sun where the plasma beta much < 1. The latitudinal redistribution is a direct consequence of the magnetic field gradient between pole and equator. Near the Sun, the latitudinal gradient in magnetic field generates meridional flows directed equatorward that tend to relax the gradient in the magnetic field (to make it more latitude-independent) as heliocentric distance increases. If there is an asymmetry between north and south magnetic field strength then the meridional flows are also asymmetric (i.e., stronger in the hemisphere of stronger magnetic field). Because the magnetic fluxes (positive and negative) in the hemispheres must be equal, the redistribution shifts the HCS into balance by approx. 16 R(solar mass). At larger distances, where the magnetic field is relatively weak (beta much > 1), the HCS can be displaced if there is a difference in total pressure between the hemispheres.

solar wind↗
NASA NTRS2014

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post-­‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi-­‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two-­‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x-­‐acceleration peak response correlated with the human response at 8 and 10-­‐G 100 ms but not 10-­‐G 70 ms. The phase lagged the human response. Head z-­‐acceleration was not correlated. Chest x-­‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x-­‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x-­‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x-­‐acceleration was not well correlated. Head and chest z-­‐acceleration was in phase but had a higher peak response. Chest z-­‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z-­‐displacement but the THOR was in phase and was comparable to the mean peak response. Head x-­‐ and z-­‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, Nate↗
NASA NTRS2018

Evacuated Airship for Mars Missions: NIAC Phase I, 2017

An evacuated or vacuum airship relies on the same principle of buoyancy used by standard balloons. However, unlike a balloon which uses a lighter than air gas to displace air and provide lift, the vacuum airship leverages a rigid structure to maintain a vacuum and displace air, thereby providing buoyancy. This method is similar to how a ship uses a rigid structure to displace water and fill the space with air; an evacuated airship uses the same mechanism, except air is displaced and the space remains vacant. Using this method, the evacuated airship is capable of utilizing the full potential of the displaced mass of air, which has interesting implications in the Martian atmosphere. Unlike other aerial vehicles, which are at a disadvantage in Martian atmospheric conditions, the evacuated airship benefits from the Martian atmosphere by virtue of the temperature and molecular composition. As a result, the evacuated airship offers an unprecedented payload capacity and, if implemented, may be used to transport current and future scientific instruments, other vehicles, rovers, and possibly even human habitations. A standard dirigible or balloon for Mars would have a severely limited span of operation and a very narrow field of study, nearly exclusively the atmosphere, but a vacuum airship can be used as a long term tool for many different missions: transportation, ground study, communications, atmospheric study, etcetera, thereby making it a far more economically sensible choice

Clarke, John-Paul↗
NASA NTRS2016

Phase Retrieval Implementation for the WFIRST Coronagraph Development Testbed

The Wide-Field Infrared Survey Telescope (WFIRST) is a 2.4m diameter space telescope NASA program. The payload will include two instruments, including a coronagraph instrument (CGI). The CGI will be the first active, high-contrast stellar coronagraph in space intended for imaging, discovery, and spectral characterization of exoplanets and debris discs. The CGI designs under development use deformable mirrors (DM) to create a point spread function (PSF) with a dark region around the obscured star object. Each DM consists of an array of actuators, and each actuator causes a local displacement to a thin, mirror-coated, glass face sheet. The actuator displacements are designed to create the PSF. Typically, phase retrieval is used to set the DM actuator displacements to their designed pattern, and then the actuators are further modified by an iterative control loop to correct for system non-idealities. To observe an exoplanet near the star object, the contrast in the annular region must be better than 10^8 (10^9 after post processing) relative to the unobstructed image of the star [1,2]. The coronagraph testbed at Jet Propulsion Laboratory is used to test coronagraph design and control algorithms. The optical system (Fig. 1) is maintained in a vacuum tank so that the PSF is stable for long integration times. Phase retrieval is an important tool for the testbed and serves several purposes, in addition to setting initial DM actuator displacements. Other uses include evaluation of systematic wavefront error introduced by system components, and alignment of the DMs and pupil masks. For example, the location of specific DM actuators relative to the Lyot mask can be determined by displacing the actuator, then measuring the location of the phase disturbance in the pupil plane with phase retrieval. In the next section we further describe the optical train of the testbed and the collection of source and pupil images available as inputs for phase retrieval, as well as one phase retrieval algorithm employed. Section 3 shows example phase retrieval results.

Marx, David↗
NASA NTRS2018

Adaptive Load Control of Flexible Aircraft Wings Using Fiber Optic Sensing

Over the past century aircraft wing design has transformed from the morphing wing used on the Wright Flyer to rigid wings with little to no shape-tailoring abilities. Modern day wings are designed to fly at a single trim condition and optimized to have a maximum aerodynamic efficiency at only this condition. Shape morphing wings on the other hand have the potential to undergo geometric changes allowing them to adapt to their mission profiles. Several flight demonstrations have been conducted over the decades using morphing-wing technologies. Active wing-twist control was demonstrated on the X-53 Active Aeroelastic Wing (AAW) research project by utilizing multiple leading- and trailing-edge control surfaces. Passive morphing technology has been demonstrated on the Rockwell RPRV-870 Highly Maneuverable Aircraft Technology (HiMAT) aircraft. In the current study, the wings of a small unmanned aerial system (sUAS) were modified to have segmented control surfaces (SCS). The modifications include segmenting the original wing control surfaces (one flap and one aileron per wing) into 44 individual sections, each section having its own independent servo control motor. The wings were also instrumented with a network of over 1800 fiber-optic strain sensors (on four sensing fibers distributed over the top and bottom surfaces of the wing) monitoring the strain response of the wing to aerodynamic loading. The SCS positions were manipulated in real time to modify the spanwise lift distribution of the wings on the sUAS. The change in the structural response of the wings caused by load redistribution was quantified by measuring the bending strains on the upper and lower wing surfaces using an on-board compact fiber-optic strain sensing (cFOSS) system. A feedback controller was developed to control the SCS positions using strain-based shape estimations from the Displacement Transfer Function (DTF). Post-processing of the strain data allowed for the transverse displacement distributions and load distributions to be compared for the conventional and segmented control surface cases using displacements and loads algorithms developed by Richards and Ko at AFRC (refs. 5-10). While the current study focused on the shifting of the spanwise aerodynamic loads as quantified by displacements, future applications for loads or displacement control might include active gust alleviation and flutter suppression.

Pena, Francisco↗
NASA NTRSDate not supplied

Method for Real-Time State Estimation of Structural Modes for an Aeroelastic Wind Tunnel Model

A method for estimating displacements, velocities, and accelerations of structural modes in generalized coordinates from measured sensor data in real time is developed and demonstrated. Data from conventional strain gauges and fiber optic strain sensors (FOSS) were combined with strain mode shapes to produce least-squares estimates of the structural mode displacements. Similarly, accelerometer data were combined with displacement mode shapes to estimate structural mode accelerations. Estimates were then combined using a Kalman filter to refine the displacement estimates and produce structural mode velocity estimates. The approach was demonstrated using simulation data for the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project in both a stable open-loop condition and an unstable condition where estimated displacement and velocity states were used for feedback control. Results supported the feasibility of using this approach for feedback control and system identification applications for wind tunnel tests.

Aeroservoelasticity↗
NASA NTRS2015

Spectral Eclipse Timing

We utilize multi-dimensional simulations of varying equatorial jet strength to predict wavelength-dependent variations in the eclipse times of gas-giant planets. A displaced hot spot introduces an asymmetry in the secondary eclipse light curve that manifests itself as a measured offset in the timing of the center of eclipse. A multiwavelength observation of secondary eclipse, one probing the timing of barycentric eclipse at short wavelengths and another probing at longer wavelengths, will reveal the longitudinal displacement of the hot spot and break the degeneracy between this effect and that associated with the asymmetry due to an eccentric orbit. The effect of time offsets was first explored in the IRAC wavebands by Williams et al. Here we improve upon their methodology, extend to a broad range of wavelengths, and demonstrate our technique on a series of multi-dimensional radiative-hydrodynamical simulations of HD 209458b with varying equatorial jet strength and hot-spot displacement. Simulations with the largest hot-spot displacement result in timing offsets of up to 100 s in the infrared. Though we utilize a particular radiative hydrodynamical model to demonstrate this effect, the technique is model independent. This technique should allow a much larger survey of hot-spot displacements with the James Webb Space Telescope than currently accessible with time-intensive phase curves, hopefully shedding light on the physical mechanisms associated with thermal energy advection in irradiated gas giants.

Hydrodynamics↗
DOE OSTI2023

In situ mechanical characterization of a single cell-cell adhesion interface under large strain

A method of measuring a stress-strain curve in a cell-cell adhesion interface, the method including: providing a structure including a first movable island supported by a first beam, a second movable island supported by a second beam, and a gap therebetween connected by a pair of cells forming a junction, the pair of cells comprising a cell-cell adhesion interface having an initial length defined by a distance between nuclei of the pair of cells; moving the second movable island with a defined displacement; determining a displacement of the first movable island based on moving the second movable island; calculating a difference between the displacement of the first movable island and the defined displacement of the second movable island based on moving the second movable island; determining an applied strain in the cell-cell adhesion interface between the pair of cells based on the difference divided by the initial length of the cell-cell adhesion interface; calculating a force between the cell-cell adhesion interface of the pair of cells based on the displacement of the first movable island; calculating a stress in the cell-cell adhesion interface between the pair of cells based on the force; and determining the stress-strain curve of the cell-cell adhesion interface between the pair of cells by plotting the calculated stress against the applied strain.

Yang, Ruiguo↗
DOE OSTI2026

HDG-1 Experiment Irradiation Monitoring Data Qualification Final Report

SUMMARY The U.S. Department of Energy (DOE) Advanced Reactor Technologies (ART) Graphite Research and Development (GRD) Program is conducting a series of six experiments to quantify the effects of irradiation on nuclear-grade graphite. This report documents the qualification of irradiation monitoring data for the fifth experiment, High Dose Graphite-1 (HDG-1). Qualified monitoring data are required by the ART program to support the design and licensing of the first high-temperature reactor (HTR) nuclear plant. Data are classified as Qualified if they meet the usage requirements described in the experiment planning and quality assurance (QA) documents, Failed if they do not meet those requirements and provide no usable information, or Trend if they do not fully meet all requirements but still provide useful information subject to an assessment of how any deficiencies may affect a particular use of the data. HDG-1 irradiation began with Advanced Test Reactor (ATR) Cycle 168B on August 24, 2020, and concluded after Cycle 173C on January 27, 2025. The HDG-1 capsule was removed from the reactor core twice—during core internal change (CIC) Cycle 170A and powered axial locator mechanism (PALM) Cycle 172A—to prevent overheating of the graphite specimens during high-power PALM cycles. The capsule was therefore irradiated during a total of seven normal ATR cycles: 168B, 169A, 171A, 171B, 173A, 173B, and 173C. Irradiation monitoring data evaluated in this report include thermocouple (TC) temperature, gas flow rate, gas moisture, gas pressure, specimen load, and graphite stack displacement. Temperature. A total of 14,508,065 TC temperature records were captured. Of these, 13,901,785 (95.8%) are Qualified and 606,280 (4.2%) are Failed. The principal source of failed temperature data was the instrument failure of TC-9 (Zone 2) on June 24, 2024, and TC-10 (Zone 1) on July 5, 2024, near the end of Cycle 173A, which resulted in 595,554 Failed readings. An additional 379 missing values and 10,347 slightly negative values from TC-13 during ATR outages are also Failed. Neither TC-9 nor TC-10 was used as a temperature-control TC, and their failures did not compromise capsule condition monitoring. Correlation analysis of all 13 TCs found no evidence of virtual junction formation. Control chart analysis revealed clear downward drift of approximately 80°C for TC-6 (Zone 3) relative to other stable TCs, and possible downward drift of approximately 60°C for TC-13 relative to the Zone 5 control TC (TC-1), though TC-13 remained consistent with the Zone 2 control TC (TC-12). Gas flow. A total of 20,088,090 gas flow rate records were captured. Of these, 19,941,463 (99.3%) are Qualified and 146,627 (0.7%) are Failed due to missing values. All argon, helium, and total gas flow data were within expected ranges throughout the irradiation. Gas moisture. A total of 1,116,005 outlet gas moisture values were captured. Of these, 1,101,421 (98.7%) are Qualified and 14,584 (1.3%) are Failed, comprising 14,556 out-of-range values and 28 missing values. The out-of-range moisture values exceeded 22,000 ppmv for approximately 1 week at the beginning of Cycle 173A, when accumulated moisture evaporated after the capsule was retrieved from water storage during PALM Cycle 172A and reinserted into the east flux trap. Moisture levels returned to below 25 ppmv for the remaining three cycles, and the transient high-moisture event did not affect the integrity of specimen irradiation. Gas pressure. A total of 7,812,035 gas pressure values were captured. Of these, 6,642,048 (85.0%) are Qualified and 1,169,987 (15.0%) outlet pressure values are Failed, comprising 718,537 zero outlet pressure values due to sensor failure from Cycle 168B through Cycle 171B, 54,550 missing values, and 396,900 too-low outlet pressure values, ranging from 1.1 to 1.6 psia after sensor replacement during Cycle 173A. Load. A total of 6,696,030 load values were captured. Of these, 6,694,580 (99.98%) are Qualified and 1,450 (0.02%) are Failed due to missing values. Applied loads to the six specimen stacks were stable throughout the irradiation. Stack displacement. A total of 6,696,030 displacement values were captured. Of these, 5,713,297 (85.32%) are Qualified and 3,781 (0.06%) are Failed due to missing values. Stack displacement increased consistently throughout the irradiation, reaching approximately 3.08 in. for Channels 5 and 6 by the end of irradiation. 978,952 (14.62%) substantially elevated displacements observed for Channel 6 beginning in Cycle 171A and for Channel 5 beginning in Cycle 173A are assigned Trend status. Raising pressure. A total of 1,115,999 raising pressure values were captured. Of these, 1,115,430 (99.95%) are Qualified and 569 (0.05%) are Failed due to missing values. Ram pressure. A total of 6,696,030 ram pressure values were captured. Of these, 6,692,249 (99.95%) are Qualified and 3,484 (0.05%) are Failed due to missing values. Stack raising was perf

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Reduced-order atomistic cascade method for simulating radiation damage in metals

Atomistic modeling of radiation damage through displacement cascades is deceptively non-trivial. Because of the high energy and stochastic nature of atomic collisions, individual primary knock-on atom (PKA) cascade simulations are computationally expensive and ill-suited for length and dose upscaling. Here, we propose a reduced-order atomistic cascade model capable of predicting and replicating radiation events in metals across a wide range of recoil energies. Our methodology approximates cascade and displacement damage production by modeling the cascade as a core-shell atomic structure composed of two damage production estimators, namely an athermal recombination corrected displacements per atom (arc-dpa) in the shell and a replacements per atom (rpa) representing atomic mixing in the core. These estimators are calibrated from explicit PKA simulations and a standard displacement damage model that incorporates cascade defect production efficiency and mixing effects. We highlight the predictability and accuracy of our reduced-order atomistic cascade method for the cases of copper and niobium by comparing its results with those from full PKA simulations in terms of defect production as well as the resulting cascade evolution and structure. We offer examples for simulating high energy cascade fragmentation and large dose ion-bombardment to demonstrate its possible applicability. Finally, we discuss the various practical considerations and challenges associated with this methodology especially when simulating subcascade formation and dose effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Utah FORGE: Fault Shear Reactivation Experimental Data for Fluid Injection-Rate Controls on Seismic Moment

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Effect of Stoichiometry on the Structure and Polarization of BaTiO 3

Barium titanate (BaTiO 3 ) is a material of interest for photonic device applications due to its strong optical non-linearity. However, BaTiO 3 -based devices have not found widespread adoption, in part due to the challenges associated with synthesizing high quality thin-films. Here, high-resolution scanning transmission electron microscope (STEM) imaging is used to investigate the atomic structure of both on- and off-stoichiometric BaTiO 3 synthesized by molecular beam epitaxy (MBE). Here, this investigation reveals an asymmetry in the way the BaTiO 3 atomic lattice accommodates off-stoichiometry growth and unveils features beyond what is expected from diffraction or surface characterization techniques. Excess titanium incorporates into the BaTiO 3 lattice to form pervasive defects despite titanium-rich films having a low surface roughness and high-quality appearance in diffraction. Excess barium forms a rough, water-soluble surface layer but does not significantly impact the quality of the BaTiO 3 lattice. STEM is used to map titanium atom displacement in real-space. The average displacement distance is 30–60 pm in the strained thin-films, higher than the <20 pm displacement in bulk BaTiO 3 . Additionally, the titanium atom displacement direction deviates from the c-axis of the unit cell, which may have implications for the material's electro-optic tensor and thus for electro-optic device design.

Cavanagh, Ashley E. [Harvard Univ., Cambridge, MA ↗
DOE OSTI2020

Interlaboratory Study of Digital Volume Correlation Error Due to X-Ray Computed Tomography Equipment and Scan Parameters: an Update from the DVC Challenge

Background: The quality of Digital Volume Correlation (DVC) full-field displacement measurements depends directly on the characteristics of the X-ray Computed Tomography (XCT) equipment, and scan procedures used to acquire the tomographic images. Objective: In this work, we seek to experimentally study the effects of XCT equipment and tomographic scan procedures on the quality of these images for DVC analysis, and to survey the level of DVC error that may be achieved using standard XCT operating procedures. Methods: Six participants in an interlaboratory study acquired high-quality XCT scans of a syntactic foam before and after rigid body motion. The resulting images were correlated using commercial DVC software to quantify error sources due to random image noise, reconstruction artifacts, as well as systematic spatial or temporal distortion. Results: In the absence of rigid body motion, the standard deviation of the displacement measurements ranged from 0.012 to 0.043 voxels using a moderate subvolume size, indicating that subvoxel measurement resolution could readily be achieved with a variety of XCT equipment and scan recipes. Comparison of consecutive scans without rigid body motion showed transient dilatational displacement gradients due to self-heating of the X-ray source and/or thermal expansion of the foam. Evaluation of the scans after rigid body motion showed significant, machine-specific spatial distortion in the displacement fields of up to 0.5 voxels; new approaches to remove this error need to be developed. Conclusions: Analysis of the scan protocols used in the interlaboratory study, as well as a complementary parametric sensitivity study, showed that the DVC error was strongly influenced by the XCT equipment, but could be mitigated by adjusting the total scan duration.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Propagation of Noise Uncertainty Through Virtual Strain Gauge Formulations for 2D Digital Image Correlation

The effect of displacement uncertainty is examined on 2-dimensional strain, calculated using linear surfaces fitted to the displacement field. A classical engineering error propagation method is used to calculate uncertainty in Green-Lagrangian strain calculations. The derived uncertainty is compared to a Monte Carlo simulation and discrepancies under 2% are seen between these two methods. The effect of virtual strain gauge size, displacement uncertainty, and boundaries on the region of interest on the strain uncertainty are considered. Here, an exponential decay relationship is observed between strain uncertainty and virtual strain gauge size, while a linear relationship is seen between strain and displacement uncertainty. For boundaries in the region of interest, strain uncertainty is affected by the reduced number of points available to perform the regression.

42 ENGINEERING↗
DOE OSTI2019

The Coupled Adjoint-State Equation in forward and inverse linear elasticity: Incompressible plane stress

A persistent challenge present in inverse or parameter estimation problems with interior data is how to deal with uncertainty in the boundary conditions employed in the forward or state model. In this work we focus on a linear plane stress inverse elasticity problem with measured displacement data where one component of the measured displacement field is known with considerably greater precision than the other. This situation is commonly encountered when the displacement field is measured using ultrasound or optical coherence tomography. We present a novel computational formulation in which no displacement or traction boundary conditions are assumed. The formulation results in coupling the state and adjoint equations, that are typically uncoupled when a well-posed state model is available. Two variants of residual-based stabilization are added. Our approach is applied to a simulated data set and experimental data from an ultrasound phantom.

42 ENGINEERING↗
DOE OSTI2020

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗
DOE OSTI2024

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗
DOE OSTI2019

Features of mechanical behavior of EV battery modules under high deformation rate

Here, we observe changes in material behavior of electric vehicle (EV) Li-ion battery cell modules when the loading speed increases from quasi-static to high speed (close to 1 m/s). A series of out-of-plane indentation experiments were performed on stacks of pouch cells placed in the specially designed enclosure to represent conditions in battery modules. Over one hundred of large format automotive pouch cells were tested under displacement rates differing by orders of magnitude. The details of internal failure were studied by X-ray tomography (XCT). We determine a shift in the force–displacement response from a parabolic behavior characteristic of particulate materials under compression to the response normal of fully dense materials with yield point. This shift is associated with dynamic effects in active materials which result in deep propagation of damage into the battery module, despite the fact that the internal short circuit is triggered when the displacement of the indenter is one half of that under the slow displacement rate.

25 ENERGY STORAGE↗