Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Thermal loading effects on chalk hydromechanical behavior for nuclear waste disposal

Safe disposal of heat-generating nuclear waste depends on host rock stability under thermal, hydrological, and mechanical stresses. This study investigates the effect of thermal loading on mechanical behavior of the shallowly buried Ghareb formation chalk through triaxial and hydrostatic constant strain rate and creep tests at temperatures up to 100 ˚C and effective pressures up to 20.7 MPa. Experimental results show that thermal loading reduces the elastic moduli of chalk by 50–75%, and a transition occurs above 60 ˚C where creep rates increase rapidly. Water saturation nearly doubles the thermally induced strain compared to dry conditions and strongly decreases material rigidity. Thermal loading also leads to significant pore pressure increases under undrained conditions and leads to reductions in the apparent permeability during drained conditions. Laboratory experimental data were used to parameterize and develop a preliminary constitutive model for predicting future deformation during repository operations in the Ghareb. The strongly coupled effects – mechanical weakening, fluid pressure fluctuations, and permeability modification – demonstrate that elevated repository temperatures will have a pronounced effect on the near field Ghareb behavior during waste disposal operations. The findings indicate that the coupled interactions must be considered in predictive models and repository design to ensure long-term nuclear waste isolation and safety.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Shock compression behavior of stainless steel 316L octet-truss lattice structures

Lattice structures offer desirable mechanical properties for applications of energy absorption and impact mitigation but limited research has been carried out on their shock compression behavior. In this work, the shock compression behavior of stainless steel 316L (SS316L) octet-truss lattice structures was investigated through experimental techniques and numerical simulations. Plate impact experiments with high-speed imaging were conducted at impact velocities of 270 – 390 m/s on lattice specimens with 5x5x10 unit cell geometries additively manufactured (AM) using direct metal laser sintering. High-speed imaging together with digital image correlation was used to extract full-field measurements and define a two-wave structure consisting of an elastic wave and planar compaction (shock) wave which propagated along the impact direction. A linear shock velocity versus particle velocity relation was found to approximate the measurements with a unit slope and a linear fit constant equal to the crushing speed. Furthermore, the shock velocity versus particle velocity relation, full-field measurements, and elastic limit together with the Eulerian form of the Rankine-Hugoniot jump conditions were used to find relations for the stress and internal energy behind the shock. Stress behind the shock increased with relative density and particle velocity, and specific internal energy converged to a single curve similar to that of bulk AM SS316L. Explicit finite element analysis using the Johnson-Cook constitutive model demonstrated similar shock behavior observed in experiments and a linear shock velocity versus particle velocity relation and corresponding Hugoniot calculations were found to be in agreement with experimental results. Numerical simulations confirmed negligible effects of exterior versus interior measurements and further validated the application of one-dimensional shock theory.

42 ENGINEERING↗
DOE OSTI2020

Localization of deformation in fluid-saturated sandstone

Herein we discuss a laboratory study of inelastic deformation and localization was conducted to evaluate microseismicity and fracture in fluid-saturated rock. Plane strain compression experiments were performed on oil- and water-saturated Berea sandstone under dry, drained, undrained, and quasi-unjacketed conditions. Parameters associated with the elastic and inelastic material response were accurately measured and used in a constitutive model to predict localization of deformation under the drained and undrained conditions. Acoustic emission (AE) activity was recorded and the onset of inelastic response in fluid-saturated rock, similar to dry rock, coincided with an increase in AE rate. Fracture initiation and propagation were successfully monitored by AE locations. Clustering of AE events and the start of inhomogeneous deformation inferred from displacement measurements are shown to be well correlated with the model prediction of the onset of localized deformation under various rock-fluid interaction conditions.

42 ENGINEERING↗
DOE OSTI2022

Data-driven analysis of neutron diffraction line profiles: application to plastically deformed Ta

Abstract Non-destructive evaluation of plastically deformed metals, particularly diffraction line profile analysis (DLPA), is valuable both to estimate dislocation densities and arrangements and to validate microstructure-aware constitutive models. To date, the interpretation of whole line diffraction profiles relies on the use of semi-analytical models such as the extended convolutional multiple whole profile (eCMWP) method. This study introduces and validates two data-driven DLPA models to extract dislocation densities from experimentally gathered whole line diffraction profiles. Using two distinct virtual diffraction models accounting for both strain and instrument induced broadening, a database of virtual diffraction whole line profiles of Ta single crystals is generated using discrete dislocation dynamics. The databases are mined to create Gaussian process regression-based surrogate models, allowing dislocation densities to be extracted from experimental profiles. The method is validated against 11 experimentally gathered whole line diffraction profiles from plastically deformed Ta polycrystals. The newly proposed model predicts dislocation densities consistent with estimates from eCMWP. Advantageously, this data driven LPA model can distinguish broadening originating from the instrument and from the dislocation content even at low dislocation densities. Finally, the data-driven model is used to explore the effect of heterogeneous dislocation densities in microstructures containing grains, which may lead to more accurate data-driven predictions of dislocation density in plastically deformed polycrystals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tracking nanoparticle growth in pulsed carbon arc discharge

The dynamics of nanoparticle growth in pulsed anodic arc discharge has been studied in time-resolved mode. To this end, a fast moving probe was employed to extract the material generated in a pulsed arc plasma held between two graphite electrodes. The probe motion was synchronized with the pulse phase and the exposure time to the plasma was set to 10 ms. The graphite anode was eroded in a helium atmosphere (300 Torr) by an arc plasma pulsed at 1 Hz with a 10% duty cycle and showing 250 A of peak current. The structure and morphology of the probe depositions were characterized by Raman spectroscopy, scanning electron microscopy, and focused ion beam. A maximal deposition rate of 260 μm/s was measured 5 mm away from the arc core during the active 0.1 s of the pulse. Such a rate yields a growth flux of 1.3 × 10 21 cm –2 s –1 , rich in carbon nanostructures (graphene platelets, nanotubes) with a characteristic aggregate size within 1–10 μm. The deposition during the inactive 0.9 s of the pulse was several orders of magnitude slower and consisted of amorphous carbon traces. Furthermore, the nanoparticle distribution along the collecting probe is correlated with the pulse phase, thereby providing information on particle transport. Pulsed nanosynthesis can be modeled as a periodical growth process, where the volume and propagation velocity of the growth region can be adjusted through modulation of the pulse signal waveform. The proposed model constitutes a suitable framework to investigate the pulsed arc synthesis of nanomaterials with tailored physical and chemical properties.

42 ENGINEERING↗
DOE OSTI2020

Experimental evaluation of shear modulus scaling of dynamic strength at extreme pressures

Recent progress in the development of dynamic strength experimental platforms is allowing for unprecedented insight into the assumptions used to construct constitutive models operating in extreme conditions. In this work, we make a quantitative assessment of how tantalum strength scales with its shear modulus to pressures of hundreds of gigapascals through a cross-platform examination of three dynamic strength experiments. Specifically, we make use of Split–Hopkinson pressure bar and Richtmyer–Meshkov instability experiments to assess the low-pressure strain and strain rate dependence. Concurrent examination of magnetically driven ramp-release experiments up to pressures of 350GPa allows us to examine the pressure dependence. In this work, using a modern description of the shear modulus, validated against both ab initio theory and experimental measurements, we then assess how the experimentally measured pressure dependence scales with shear modulus. Furthermore, we find that the common assumption of scaling strength linearly with the shear modulus is too soft at high pressures and offer discussion as to how descriptions of slip mediated plasticity could result in an alternative scaling that is consistent with the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. I. Inverse Lagrangian analysis

Absolute measurements of solid-material compressibility by magnetically driven shockless dynamic compression experiments to multi-megabar pressures have the potential to greatly improve the accuracy and precision of pressure calibration standards for use in diamond anvil cell experiments. Here, to this end, we apply characteristics-based inverse Lagrangian analysis (ILA) to 11 sets of ramp-compression data on pure platinum (Pt) metal and then reduce the resulting weighted-mean stress–strain curve to the principal isentrope and room-temperature isotherm using simple models for yield stress and Grüneisen parameter. We introduce several improvements to methods for ILA and quasi-isentrope reduction, the latter including calculation of corrections in wave speed instead of stress and pressure to render results largely independent of initial yield stress while enforcing thermodynamic consistency near zero pressure. More importantly, we quantify in detail the propagation of experimental uncertainty through ILA and model uncertainty through quasi-isentrope reduction, considering all potential sources of error except the electrode and window material models used in ILA. Compared to previous approaches, we find larger uncertainty in longitudinal stress. Monte Carlo analysis demonstrates that uncertainty in the yield-stress model constitutes by far the largest contribution to uncertainty in quasi-isentrope reduction corrections. We present a new room-temperature isotherm for Pt up to 444 GPa, with 1-sigma uncertainty at that pressure of just under ±1.2%; the latter is about a factor of three smaller than uncertainty previously reported for multi-megabar ramp-compression experiments on Pt. The result is well represented by a Vinet-form compression curve with (isothermal) bulk modulus K 0 = 270.3 ± 3.8 GPa, pressure derivative K$^{'}_{0}$= 5.66 ± 0.10, and correlation coefficient $R_{K_{0},K^{'}_{0}}$= –0.843.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Implementation of compound refractive lenses for large field-of-view x-ray phase-contrast imaging during hypervelocity impact experiments

Synchrotron x-ray phase-contrast imaging (XPCI) offers time-resolved visualization of dynamic compression phenomena, but its intrinsically small field-of-view (FOV) limits the time that key features remain in frame. A novel approach to enlarge the FOV is achieved by positioning a two-dimensional parabolic compound refractive lens (CRL) upstream of the sample to deliberately defocus the white beam. Ray-tracing simulations and XPCI measurements show that this CRL configuration can expand the beam by ∼50% vertically and ∼15% horizontally based on the full width at half-maximum of the beam. Implementing the CRL, however, attenuates the photon flux and lowers signal-to-noise ratio (SNR). Task-based analysis using a calibration grid (30 μm dots) showed that both setups fail to consistently meet the Rose criterion (SNR ≥ 5) for features of this size in single-bunch imaging. Extrapolating the measured SNR Rose values suggests that the minimum consistently detectable feature lies closer to 30–40 μm for the standard XPCI setup and above 40 μm for CRL-XPCI. Despite this limitation, the CRL configuration nearly doubles the illuminated area, enabling simultaneous tracking of front and rear observations of boron carbide targets subjected to rod and sphere impacts at 1.0–2.6 km/s. Image tracking algorithms and photonic Doppler velocimetry were used to measure penetration and rear-surface velocity histories. Together, these measurements capture crack fronts, penetration, and material breakout, offering new benchmark data for validating high-strain-rate constitutive models of ceramic materials.

Ceramic materials↗
DOE OSTI2021

Time–connectivity superposition and the gel/glass duality of weak colloidal gels

Colloidal gels result from the aggregation of Brownian particles suspended in a solvent. Gelation is induced by attractive interactions between individual particles that drive the formation of clusters, which in turn aggregate to form a space-spanning structure. We study this process in aluminosilicate colloidal gels through time-resolved structural and mechanical spectroscopy. Using the time–connectivity superposition principle a series of rapidly acquired linear viscoelastic spectra, measured throughout the gelation process by applying an exponential chirp protocol, are rescaled onto a universal master curve that spans over eight orders of magnitude in reduced frequency. This analysis reveals that the underlying relaxation time spectrum of the colloidal gel is symmetric in time with power-law tails characterized by a single exponent that is set at the gel point. The microstructural mechanical network has a dual character; at short length scales and fast times it appears glassy, whereas at longer times and larger scales it is gel-like. These results can be captured by a simple three-parameter constitutive model and demonstrate that the microstructure of a mature colloidal gel bears the residual skeleton of the original sample-spanning network that is created at the gel point. Our conclusions are confirmed by applying the same technique to another well-known colloidal gel system composed of attractive silica nanoparticles. The results illustrate the power of the time–connectivity superposition principle for this class of soft glassy materials and provide a compact description for the dichotomous viscoelastic nature of weak colloidal gels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiaxial Plastic Deformation of Zircaloy-4 Nuclear Fuel Cladding Tubes

Here, this work is motivated by the desire to devise an internal pressure test that can mimic a displacement-controlled loading scenario and demonstrate how to apply the multiaxial stress and strain data from the test to develop an elastic/plastic constitutive model for a thin-walled tubular component. This is achieved by conducting simultaneous measurements of tangential and axial strain during the pressure test and integrating these strain measures into a feedback loop with the pressure controller. It is shown how data from such a test can be used to develop a large mechanical property data set relevant to biaxial loading conditions. The data obtained have high confidence evidenced by their low variability and alignment with other literature studies. Additionally, data from these internal pressure tests combined with full-tube axial tensile tests allow for the derivation of the Hill anisotropic yield function. The developed Hill yield function is validated by comparing the plastic strain ratios from the full-tube tension tests and by comparing the predicted yield stress in the tangential direction with measured values from ring tension tests in a previous study.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

The Role of Unit-Cell Topology in Modulating the Compaction Response of Additively Manufactured Cellular Materials using Simulations and Validation Experiments

Additive manufacturing has enabled a transformational ability to create cellular structures (or foams) with tailored topology. Compared to their monolithic polymer counterparts, cellular structures are potentially suitable for systems requiring materials with high specific energy-absorbing capability to provide enhanced damping. In this work, we demonstrate the utility of controlling unit-cell topology with the intent of obtaining a desired stress–strain response and energy density. Using mesoscale simulations that resolve the unit-cell sub-structures, we validate the role of unit-cell topology in selectively activating a buckling mode and thereby modulating the characteristic stress–strain response. Simulations incorporate a linear viscoelastic constitutive model and a hyperelastic model for simulating large deformation of the polymer under both tension and compression. Simulated results for nine different cellular structures are compared with experimental data to gain insights into three different modes of buckling and the corresponding stress–strain response.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Neutrino-assisted early dark energy: theory and cosmology

Abstract The tension between measurements of the Hubble constant obtained at different redshifts may provide a hint of new physics active in the relatively early universe, around the epoch of matter-radiation equality. A leading paradigm to resolve the tension is a period of early dark energy, in which a scalar field contributes a subdominant part of the energy budget of the universe at this time. This scenario faces significant fine-tuning problems which can be ameliorated by a non-trivial coupling of the scalar to the standard model neutrinos. These become non-relativistic close to the time of matter-radiation equality, resulting in an energy injection into the scalar that kick-starts the early dark energy phase, explaining its coincidence with this seemingly unrelated epoch. We present a minimal version of this neutrino-assisted early dark energy model, and perform a detailed analysis of its predictions and theoretical constraints. We consider both particle physics constraints — that the model constitute a well-behaved effective field theory for which the quantum corrections are under control, so that the relevant predictions are within its regime of validity — and the constraints provided by requiring a consistent cosmological evolution from early through to late times. Our work paves the way for testing this scenario using cosmological data sets.

Astronomy & Astrophysics↗
DOE OSTI2021

Impact of confined geometries on hopping and trapping of motile bacteria in porous media

We use a random walk particle-tracking (RWPT) approach to elucidate the impact of porous media confinement and cell-cell interactions on bacterial transport. The model employs stochastic alternating motility states consisting of hopping movement and trapping reorientation. The stochastic motility patterns are defined based on direct visualization of individual trajectory data. We validate our model against experimental data, at single-cell resolution, of bacterial E. coli motion in three-dimensional confined porous media. Results show that the model is able to efficiently simulate the spreading dynamics of motile bacteria as it captures the impact of cell-cell interaction and pore confinement, which marks the transition to a late-time subdiffusive regime. Furthermore, the model is able to qualitatively reproduce the observed directional persistence. Lastly, our RWPT model constitutes a meshless simple method which is easy to implement and does not invoke ad hoc assumptions but represents the basis for a multiscale approach to the study of bacterial dispersal in porous systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Analyses of internal structures and defects in materials using physics-informed neural networks

Characterizing internal structures and defects in materials is a challenging task, often requiring solutions to inverse problems with unknown topology, geometry, material properties, and nonlinear deformation. Here, we present a general framework based on physics-informed neural networks for identifying unknown geometric and material parameters. By using a mesh-free method, we parameterize the geometry of the material using a differentiable and trainable method that can identify multiple structural features. We validate this approach for materials with internal voids/inclusions using constitutive models that encompass the spectrum of linear elasticity, hyperelasticity, and plasticity. We predict the size, shape, and location of the internal void/inclusion as well as the elastic modulus of the inclusion. Our general framework can be applied to other inverse problems in different applications that involve unknown material properties and highly deformable geometries, targeting material characterization, quality assurance, and structural design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale analysis in solids with unseparated scales: fine-scale recovery, error estimation, and coarse-scale adaptivity

There are several engineering applications in which the assumptions of homogenization and scale separation may be violated, in particular, for metallic structures constructed through additive manufacturing. Instead of resorting to direct numerical simulation of the macroscale system with an embedded fine scale, an alternative approach is to use an approximate macroscale constitutive model, but then estimate the model-form error using a posteriori error estimation techniques and subsequently adapt the macroscale model to reduce the error for a given boundary value problem and quantity of interest. Here, we investigate this approach to multiscale analysis in solids with unseparated scales using the example of an additively manufactured metallic structure consisting of a polycrystalline microstructure that is neither periodic nor statistically homogeneous. As a first step to the general nonlinear case, we focus here on linear elasticity in which each grain within the polycrystal is linear elastic but anisotropic.

42 ENGINEERING↗
DOE OSTI2021

Accelerated Creep Testing of Inconel 718 Using the Stepped Isostress Method (SSM)

This paper demonstrates the stepped isostress method (SSM); an accelerated creep test (ACT) for the rapid assessment of metallic materials. The SSM test is based on the time-temperature-stress superposition principle (TTSSP) where temperature and/or stress are step increased to accelerate the time-to-rupture. The SSM test has proven successful for the ACT of polymers and polymeric composites but has yet to be proven for metallic materials. In this study, new test matrix design rules for the SSM of metals are established based on deformation mechanism, time-temperature transformation, and time-temperature precipitation maps. A test matrix of SSM and conventional creep tests (CCTs) are executed for alloy Inconel 718 at 750°C (1382°F) with stress levels ranging from 100 to 350 MPa. Validation CCT data is gathered from the Japan National Institute of Material Science (NIMS). Material constants for the Sine-hyperbolic (Sinh) constitutive model are calibrated using the SSM data and employed to predict the conventional creep response. Here, when blindly compared to the CCT data, the SSM calibrated Sinh model can accurately predict the conventional creep response across logarithmic decades and thus accelerate the capture of conventional creep data. Fractography indicates ductile fracture by transgranular microvoid coalescence (TMVC) with the same fracture mode observed in both SSM and CCT specimens. Creep cavitation is indicated by the population, smoothness, and size of microvoids.

36 MATERIALS SCIENCE↗