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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2025

Modeling of melting and erosion in high velocity microparticle impacts

High velocity collisions can cause melting, a situation involving elasticity, plasticity, and phase transformations in tandem, which is very challenging for mechanical modeling. Here, this paper advances the use of smoothed particle hydrodynamics (SPH) simulations for this complex situation, investigating stainless steel microspheres impacting a tin substrate. After extensive validation against a range of mechanical and geometrical measurements of experimental impacts, the model provides a detailed, quantitative view of melt volume production over a broad range of impact velocities. It also enables direct decomposition of the impact energy partitioning, including analysis of the energy absorbed by plasticity in the impactor, the system kinetic energy, and the adiabatic heating of the substrate in addition to the melting itself. Since these quantities are not available experimentally, SPH modeling substantially enhances the quantitative understanding of such impacts. Finally, the model speaks to the phase characteristics of the ejecta formed during impact, confirming the experimental inference that the ejected material likely contains a significant fraction of solid material.

Constitutive modeling↗
DOE OSTI2021

USAMP Low-Cost Mg Sheet Component Development and Demonstration Project

The overarching objective of this USAMP project was to develop and demonstrate door panels made from magnesium (Mg) sheet with a cost penalty over conventional steel stampings of no more than $\$5.50$/kg saved. The technical approach integrated experiments with advanced computational tools based on Integrated Computational Materials Engineering (ICME) methods to develop new alloy chemistries and their thermomechanical processing that promise improved formability and lower forming temperatures. A penultimate task before finally forming the stampings was to incorporate actual microstructure into models that would enable formability simulations. This approach would, for the first time, account for individual magnesium grains moving in an anisotropic fashion unlike that for aluminum or steel that have isotropic properties upon which the current simulation tools are based. In separate activities, new coatings and lubricants to facilitate forming and improved corrosion protection and joining strategies, were developed to ensure that the door could be produced with stated product requirements. A technical cost model, which included parts production, assembly, and paint for a door specifically designed for Mg sheet, showed the cost penalty to be between $\$4.26$ to $\$6.60$/kg saved, which enveloped the project’s cost targets. The cost of the coated Mg sheet was identified as the key driver for the cost penalty. The mass of the Mg-intensive door was 7.9 kg, which was 54% less than the baseline steel door.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Low-Cost Magnesium Alloy Sheet Component Development and Demonstration Project

The overarching objective of this USAMP project was to develop and demonstrate door panels made from magnesium (Mg) sheet with a cost penalty over conventional steel stampings of no more than $5.50/kg saved. The technical approach integrated experiments with advanced computational tools based on Integrated Computational Materials Engineering (ICME) methods to develop new alloy chemistries and their thermomechanical processing that promise improved formability and lower forming temperatures. A penultimate task before finally forming the stampings was to incorporate actual microstructure into models that would enable formability simulations. This approach would, for the first time, account for individual magnesium grains moving in an anisotropic fashion unlike that for aluminum or steel that have isotropic properties upon which the current simulation tools are based. In separate activities, new coatings and lubricants to facilitate forming and improved corrosion protection and joining strategies, were developed to ensure that the door could be produced with stated product requirements. A technical cost model, which included parts production, assembly, and paint for a door specifically designed for Mg sheet, showed the cost penalty to be between 4.26 USD to 6.60 USD/kg saved, which enveloped the project’s cost targets. The cost of the coated Mg sheet was identified as the key driver for the cost penalty. The mass of the Mg-intensive door was 7.9 kg, which was 54% less than the baseline steel door.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Model for the Softening Factor within Stages of Work Hardening

A formulation for a dimensionless coefficient c b is derived that represents a scale of microstructural softening for alloys which follow Kocks-Mecking (K-M) work hardening behavior. The variation of the true plastic strain ε p between the proportional limit σ y and the strength σ u at the instability is determined using the Considère criterion. Parameterization of the model is limited to variables expressly measured within tensile experiments. Further development is now made for the softening factor c bi through the individual and sequential stages 3 and 4 of plastic deformation during K-M work hardening Q behavior. Application is shown for tensile test results of Ti-6Al-4V made by different additively manufactured (AM) processes. It is found that the variation in plastic strain as a function of c bi produces a continuous curve representative of the alloy system. The results of data analysis indicate that Θ o3 and c b3 increase, while Θ o4 and c b4 decrease, as the total plastic strain ε p increases. Furthermore, formulations derived for K-M stages 3 and 4 enable the evaluation of other material parameters such as the activation volume ν* for the onset of plastic deformation. This activation volume is found to be near constant, at a ν*-value of 0.353±0.036 nm 3 as computed using a strain-rate sensitivity of strength exponent m of 0.014, irrespective of the AM method used to produce the Ti-6Al-4V alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Non-linear elastic-plastic behavior of the invert glass lithium phosphorous oxynitride (LiPON)

Micro and nano-scale mechanical behavior of network/ionic glasses is dictated by their composition and the number of constraints per structural unit. From this perspective, the glass LiPON presents the opportunity for enhancement of the microscale ductility due to its rather unconstrained orthophosphate structure. Here, we use instrumented nanoindentation with different tip geometries to investigate the mechanical response of LiPON glass. The results reveal that the elastic modulus of LiPON is not constant and depends on pressure. With the method utilizing spherical nanoindentation and continuous stiffness measurement (CSM) we determine the yield point of LiPON. We propose the Drucker-Prager type of yield criterion for LiPON and estimate its yield stress in compression as 2.4 GPa. There exists another stress threshold, however, around 490 MPa, at which the elastic deformation becomes non-linear, and this can be mistaken for the yield stress.

Constitutive model↗
DOE OSTI2020

Discrepancy in the Critical State Void Ratio of Poorly Graded Sand due to Shear Strain Localization

The critical state (CS) concept is a theoretical framework that models the constitutive behavior of soils, including sand and other granular materials. It supports the notion of a unique postfailure state, where the soil ultimately experiences continuous shearing with no change in the plastic volumetric strain. However, the published literature has frequently noted the nonconvergence of sand specimens with different initial densities to a unique CS in the compression plane due to many factors such as specimen fabric, particle morphology, breakage, and grain size distribution. This paper examines the CS for poorly graded (uniform) glass beads and 3 different types of silica sands using 50 conventional triaxial compression (CTC) experiments, 12 oedometer tests, and in situ synchrotron microcomputed tomography (SMT) scans for 10 CTC experiments. The results of the 50 CTC experiments revealed a diffused CS zone in the compression plane, which was further examined using the in situ SMT scans. A thorough three-dimensional image analysis of the SMT scans accurately quantified the evolution of the local void ratio (elocal) versus axial compression within zones of intensive shearing toward the center of the specimen. The evolution of the void ratio was also measured using the entire volume of the specimen (eglobal). At the CS, the elocal/eglobal ratio was assessed to be ~1.25 when a single shear band developed within the scanned specimens and ~1.1-1.15 for specimens that failed via external bulging that was internally manifested by the development of multiple shear bands. This finding suggests that the CS zone in the compression plane can be attributed to the common wrong consideration of eglobal evolution in lieu of elocal within the developing shear bands. Furthermore, the lack of shear band development in uniaxial compression has made the results of the oedometer test reliable in quantifying the CS parameters in the compression plane.

58 GEOSCIENCES↗
DOE OSTI2021

Genomic Signatures of a Major Adaptive Event in the Pathogenic Fungus Melampsora larici-populina

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host–pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype–phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

pnnl/MCRASTA

McRasta (Markov Chain Rate and State Analysis) was developed to estimate parameter uncertainty in constitutive friction models via Bayesian inverse and Markov Chain Monte Carlo (MCMC) methods.

Fichera, Marissa [Pacific Northwest National Labor↗
DOE OSTI2021

A Constitutive Structural Parameter for the Work Hardening Behavior of Additively Manufactured Ti-6Al-4V

The mechanical behavior of Ti-6Al-4V produced by additive manufacturing processes is assessed based on a formulation developed from the Kocks-Mecking relationship. A constitutive parameter derived for the microstructure is characteristic of the work hardening behavior determined by the plastic strain between the proportional limit and the strength at the instability point. The varied plasticity behavior associated with surface and build direction effects can be evaluated with this approach as presented for the case of Ti-6Al-4V under uniaxial tension.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On the peridynamic effective force state and multiphase constitutive correspondence principle

This article concerns modeling unsaturated deformable porous media as an equivalent single-phase and single-force state peridynamic material through the effective force state. The balance equations of linear momentum and mass of unsaturated porous media are presented by defining relevant peridynamic states. The energy balance of unsaturated porous media is utilized to derive the effective force state for the solid skeleton that is an energy conjugate to the nonlocal deformation state of the solid, and the suction force state. Through an energy equivalence, a multiphase constitutive correspondence principle is built between classical unsaturated poromechanics and peridynamic unsaturated poromechanics. The multiphase correspondence principle provides a means to incorporate advanced constitutive models in classical unsaturated porous theory directly into unsaturated peridynamic poromechanics. Finally, numerical simulations of localized failure in unsaturated porous media under different matric suctions are presented to demonstrate the feasibility of modeling the mechanical behavior of such three-phase materials as an equivalent single-phase peridynamic material through the effective force state concept.

42 ENGINEERING↗
DOE OSTI2021

Assessment Approach to Advanced Fuel Models

This document constitutes completion of the NEAMS milestone M3MS-21IN0201016, which is titled: Assessment of advanced fuel improvements. In this report we present: (1) simulations of separate effects creep tests for uranium alloy and HT-9; (2) examples of standardizing and updating metallic fuel assessments to include the latest models and key output figures of merit, which is coordinated with LANL to better evaluate advanced models; (3) examples of how the FIPD database at ANL is integrated into metallic fuel assessment cases; and finally (4) documentation of fuel-specific figures of merit and the BISON review process for assessments with best practices for preparing and reviewing finite element simulations in the appendix.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

On the inference of viscoelastic constants from stress relaxation experiments

Several constitutive theories have been proposed in the literature to model the viscoelastic response of materials, including widely used rheological constitutive models. These models are characterized by certain parameters ("time constants") that define the time scales over which the material relaxes. These parameters are primarily obtained from stress relaxation experiments using curve-fitting techniques. However, the question of how best to estimate these time constants remains open. As a step towards answering this question, we propose an optimal experimental design approach based on ideas from information geometry, namely Fisher information and Kullback-Leibler divergence. The material is modeled as a spring element in parallel with multiple Maxwell elements and described using a one- or two-term Prony series. Treating the time constants as unknowns, we develop expressions for the Fisher information and Kullback-Leibler divergence that allow us to maximize information gain from experimental data. Based on the results of this study, we propose that the largest time constant estimated from a stress relaxation experiment for a linear viscoelastic material should be at most one-fifth of the total time of the experiment in order to maximize information gain. Our results also provide confirmation that the equilibrium modulus of the material cannot be reliably determined from curve-fitting to data from a stress relaxation experiment.

Fisher information↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A reactive electrochemomechanical theory for growth and remodeling of polyelectrolyte hydrogels and application to dynamic polymerization of $\text{DNA}$ hydrogels

Here this study develops a framework for growth and remodeling of active polyelectrolyte hydrogels that accounts for effects of compositional changes on the mechanical response. By developing a reactive electrochemomechanical theory, thermodynamical constraints upon reactive and remodeling processes are elucidated within a general framework that allows any number of chemical reactions to evolve the response of the gel and transfer mass and charge between constituents. Fully coupled, nonlinear constitutive relations are adopted for molar fluxes, allowing exploration of effects including cross-diffusion, electrophoresis, and electro-osmosis. A robust finite element implementation is developed in the open source FEBio software (febio.org) by exploiting an equivalence between electrochemomechanics and mixture theory. The implementation is verified against analytical solutions for free swelling, and a proper reduction to a prior chemomechanical theory is demonstrated for neutral gels swollen only by a solvent with no solutes. The theory and implementation are then applied to model the tunable large swelling achieved through dynamic polymerization of DNA crosslinkers seen in our recently developed experimental hydrogel system (Cangialosi et al., 2017). A novel constitutive model for reaction-driven evolution of the locking stretch λ L in a non-Gaussian mechanical free energy was developed, where the increasing concentration of DNA crosslinkers makes further swelling energetically favorable. With a single free parameter, excellent agreement was found between measured and predicted equilibrium swelling ratios. This study demonstrated the ability to extend the electrochemomechanical framework to include chemical reactions and composition-aware constitutive models, and showed that development of reactive models allows simulation of complex dynamic polymerization phenomena not treated before. The theoretical frame here can be further expanded in scope to incorporate additional non-ideal and nonlinear phenomena.

42 ENGINEERING↗
DOE OSTI2021

Applicability of discrete element method with spherical and clumped particles for constitutive study of granular materials

Discrete element method (DEM) has been intensively used to study the constitutive behaviour of granular materials. However, to what extent a real granular material can be reproduced by virtual DEM simulations remains unclear. This study attempts to answer this question by comparing DEM simulations with typical features of experimental granular materials. Three groups of models with spherical and clumped particles are investigated from four perspectives: (i) deviatoric stress and volumetric behaviour; (ii) critical state behaviour; (iii) stress-dilatancy relationship; and (iv) the evolution of principal stress ratio against axial strain. The results demonstrate that DEM with spherical or clumped particles is capable of qualitatively describing macroscopic deviatoric stress responses, volumetric behaviour, and critical state behaviour observed in experiments for granular materials. On the other hand, some qualitative deviations between experiments and the investigated DEM simulations are also observed, in terms of the stress-dilatancy behaviour and principal stress ratio against axial strain, which are proven to be critical for constitutive modelling. The results demonstrate that DEM with spherical or clumped particles may not necessarily fully capture experimental features of granular materials even from a qualitative perspective. It is thus encouraged to thoroughly validate DEM with experiments when developing constitutive models based on DEM observations.

58 GEOSCIENCES↗
DOE OSTI2024

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗