Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “behavior modeling”

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

Crushing behaviors of unidirectional carbon fiber reinforced plastic composites under dynamic bending and axial crushing loading

The dynamic crushing characteristics of unidirectional carbon fiber reinforced plastic composites under two loading types, dynamic three-point bending and axial crushing, are investigated by experiment and finite element simulation. Hat section samples with two different layup orientations are tested at various impact velocities to investigate its effects. The experiment results show that delamination plays a critical role in dynamic bending deformation. Different layup orientations lead to similar crushing bending behaviors, but remarkable variations in delamination pattern evolutions. In contrast, layup orientation exerts more significant effects on failure modes and energy absorption in axial crushing. A finite element model with multiple layers of thick shell element and cohesive element is established, and simulations are performed for both dynamic three-point bending and axial crushing. The model can capture the crushing behaviors in the two loading types reasonably well. The damage evolution and energy absorption through various mechanisms at different loading conditions are discussed.

42 ENGINEERING↗
DOE OSTI2021

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 304L Stainless Steel Alloy

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The material of interest here is 304L stainless steel, and the objective of this work is to calibrate thermal-mechanical models: one for the constitutive behavior and another for failure. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of the models. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s. It also included notched tension and shear-dominated compression hat tests specifically designed to calibrate the failure model. All test specimens were extracted from a single piece of plate to maintain consistency. The constitutive model adopted was a modular $J_2$ plasticity model with isotropic hardening that included rate and temperature dependence. A criterion for failure initiation based on a critical value of equivalent plastic strain fitted the failure data appropriately and was adopted. Possible ranges of the values of the parameters of the models were determined partially on historical data from calibrations of the same alloy from other lots and are given here. The calibration of the parameters of the models were based on finite element simulations of the various material tests using relatively ne meshes and hexahedral elements. When using the model in structural finite element calculations, however, element formulations and sizes different from those in the calibration are likely to be used. A brief investigation demonstrated that the failure initiation predictions can be particularly sensitive to the element selection and provided an initial guide to compensate for the effect of element size in a specific example.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Enabling Dynamic Probabilistic Risk Assessment of Physical Security Using EMRALD and MAAP

The optimization of physical security in nuclear power plants requires sophisticated methodologies that integrate operator actions and plant behavior through advanced simulation tools. To address this, Idaho National Laboratory [JL2.1]has developed the Modeling and Analysis for Safety and Security using the Dynamic EMRALD Framework (MASS-DEF) methodology, an approach that integrates force-on-force simulations, dynamic probabilistic risk assessment, and thermal-hydraulics modeling [JL3.1]to enhance security planning while reducing costs. We developed a tool that produces reduced order models using thermal hydraulic simulations from the Modular Accident Analysis Program (MAAP) [1]. These models can quickly evaluate reactor core behavior during attack simulations, and in so doing, address two barriers of traditional methods: (1) MAAP simulations are computationally intensive, and (2) attack scenarios must be run in a secure environment, which complicates analysis and validation. By precomputing scenario outcomes for a small number of modified parameters, the reduced order model significantly decreases the computational cost and enables offsite review of the results.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Strong-coupling expansion of multi-band interacting models: Mapping onto the transverse-field J 1 - J 2 Ising model

We investigate a class of two-dimensional two-band microscopic models in which the inter-band repulsive interactions play the dominant role. We first demonstrate three different schemes of constraining the ratios between the three types of inter-band interactions – density-density, spin exchange, and pair-hopping – that render the model free of the fermionic sign-problem for any filling and, consequently, amenable to efficient Quantum Monte Carlo simulations. We then study the behavior of these sign-problem-free models in the strong-coupling regime. In the cases where spin-rotational invariance is preserved or lowered to a planar symmetry, the strong-coupling ground state is a quantum paramagnet. However, in the case where there is only a residual Ising symmetry, the strong-coupling expansion maps onto the transverse-field J 1 - J 2 Ising model, whose pseudospins are associated with local inter-band magnetic order. We show that by varying the band structure parameters within a reasonable range of values, a variety of ground states and quantum critical points can be accessed in the strong-coupling regime, some of which are not realized in the weak-coupling regime. Here, we compare these results with the case of the single-band Hubbard model, where only intra-band repulsion is present, and whose strong-coupling behavior is captured by a simple Heisenberg model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Engineering grain boundary anisotropy to elucidate grain growth behavior in alumina

Current grain growth models have evolved to account for the relationship between grain boundary energy/mobility anisotropy and the five degrees of grain boundary character. However, the role of grain boundary networks on overall growth kinetics remains poorly understood. To experimentally investigate this problem, a highly textured Al 2 O 3 was fabricated by colloidal casting in a strong magnetic field to engineer a unique spatial distribution of grain boundary character. Microstructural evolution was quantified and compared to an untextured sample. From this comparison, a prevalence of (0001)/(0001) terminated grain boundaries with anisotropic networks were identified in the textured sample. These boundaries and their networks were found to be driving grain growth at a faster rate than predicted by models. These findings will allow better modelling of grain growth in real systems by experimentally exploring the impact thereon of grain boundary plane anisotropy and relative energy/mobility differences between neighboring boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A simple three-component mixing problem for the evaluation of a new reaction rate model

A simple computational mixing problem is presented which can be utilized to assess the behavior of Reynolds-averaged reaction rate models in a problem with temporally varying mixedness. In this problem, three mixing components are homogeneously distributed but initially separated in a triply periodic domain. Further, these components are initialized within a Taylor–Green-like velocity field, which creates a mixing history evolving from the so-called “no-mix limit” to a well-mixed state. Large-eddy simulation results from this problem in configurations involving both premixed and nonpremixed reactants are then compared with zero-dimensional Reynolds-averaged Navier–Stokes results utilizing a new model for multicomponent reacting mixtures. The new model is shown to appropriately respect the no-mix limit and outperforms an earlier model (Morgan, 2022), particularly at early times when components are near the no-mix limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Water Migration and Swelling in Engineered Barrier Materials for Radioactive Waste Disposal

Deep, underground repositories are needed to isolate radioactive waste from the biosphere. Because bentonite is an integral component of many multibarrier repository systems, information on the hydraulic behavior of bentonite is crucial for modeling the long-term viability of such systems. In this paper the hydraulic behavior of bentonite samples was analyzed as a function of aggregate size, and samples were subjected to hydrothermal treatments involving contact with NaCl, KCl, and deionized water. Neutron and X-ray imaging were used to quantify water sorption into packed bentonite samples and bentonite swelling into the water column. The distance between the original clay-water interface and the wetting front was determined as a function of time. Average water uptake exhibited a square-root-of-time dependence in freshly prepared samples, but more variable rates were observed for samples previously in contact with water. The radiography was supported by small-angle neutron scattering analysis and ultra-small-angle neutron scattering analysis of aggregate size distributions and by inelastic neutron scattering to understand the physicochemical environment of the sorbed water. Results showed that hydrothermal treatment with KCl had the greatest effect of increased water transport in the bentonite, possibly as a result of the interaction of K + with smectite layers in the clay.

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

Tensor network formulation of two-dimensional gravity

We show how to formulate a lattice gauge theory whose naive continuum limit corresponds to two-dimensional (Euclidean) quantum gravity including a positive cosmological constant. More precisely the resultant continuum theory corresponds to gravity in a first-order formalism in which the local frame and spin connection are treated as independent fields. Recasting this lattice theory as a tensor network allows us to study the theory at strong coupling without encountering a sign problem. In two dimensions this tensor network is exactly soluble and we show that the system has a series of critical points that occur for pure imaginary coupling and are associated with first order phase transitions. We then augment the action with a Yang-Mills term which allows us to control the lattice spacing and show how to apply the tensor renormalization group to compute the free energy and look for critical behavior. Finally we perform an analytic continuation in the gravity coupling in this extended model and show that its critical behavior in a certain scaling limit depends only on the topology of the underlying lattice. We also show how the lattice gauge theory can be naturally generalized to generate the Polyakov or Liouville action for two dimensional quantum gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

An investigation of shock-induced phase transition in soda-lime glass

There exists a large body of evidence from experiments and molecular dynamics simulations to suggest the occurrence of phase transitions in soda-lime glass (SLG) and other silica glasses subject to shock compression to pressures above 3 GPa. In light of these findings, the current work investigated the existence of phase transition in SLG using shock and release experiments. The experiments employed symmetric SLG–SLG impact to achieve complete unloading to zero stress after shock compression to stresses in the range of 3–7 GPa. The stress–strain response and the Lagrangian release wave speed behavior of SLG obtained from these experiments are seen to reveal a mismatch between the loading and unloading paths of the pressure–strain curve for the material, which serves as compelling evidence for the occurrence of a shock-induced phase transition in the material at relatively low pressures. Furthermore, the release wave speed vs strain data obtained from experiments were used to construct a methodology for modeling the shock and release behavior of SLG. Lastly, this scheme implemented in numerical simulations was able to capture the release behavior of shock compressed SLG, for which a robust and satisfactory model was previously unavailable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Experiment and non-local crystal plasticity finite element study of nanoindentation on Al-8Ce-10Mg alloy

Cerium and magnesium strengthened aluminum alloys, or Al-Ce-Mg, is a recently developed alloy family that exhibits good mechanical properties at elevated temperatures (~300 °C). To examine the single-crystal properties of Al-Ce-Mg alloys, nanoindentation experiments are conducted in this study. A crystal plasticity finite element model (CPFEM) with the evolution of geometrically necessary dislocations (GNDs) is applied to simulate the indentation in individual grains. A parametric study is carried out to investigate the sensitivity of each crystal plasticity model parameter to the indentation behavior. The highly sensitive parameters are calibrated by matching the indentation load-depth curves, while the rest parameters are obtained from bulk polycrystal uniaxial tension tests. Overall, satisfactory matching between experiment and simulation is obtained for each individual grain. The calculated hardness, as determined from the experiment, shows the dependence on indent depth, which is captured by the GND model. Furthermore, the effect of grain orientation and neighboring grains to nanoindentation behavior have been discussed with the comparison between the simulation and experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High‐pressure oxidation of n ‐butane

Abstract The oxidation of n ‐butane at elevated pressures has been investigated by experiments in a laminar flow reactor at 100 bar and temperatures of 450–900 K. The onset temperature for reaction increased from 550 K under oxidizing conditions (Φ = 0.02) to 625 K under reducing conditions (Φ = 13). NTC behavior was observed at 600–650 K (Φ = 0.02) and 625–675 K (Φ = 1.0). A detailed chemical kinetic model for the oxidation of n ‐butane was established. The present model and those suggested in literature were evaluated against the present experimental results and literature data at elevated pressures. None of the tested models could accurately reproduce the NTC behavior of n ‐butane under stoichiometric conditions of the present study, but all evaluated models could reproduce experimental data from literature with different levels of accuracy.

Hashemi, Hamid↗
DOE OSTI2021

Fracture Mechanics—Theory, Modeling and Applications

The field of fracture mechanics was developed during the throes of World War II, and since then, it has been a very active area of research. This is a very energetic and vibrant field because fracture processes are relevant for a vast range of engineering applications spanning different temporal and length scales. Depending on the application, material fracture and fragmentation may be either a desirable feature (for example, for the purpose of hydraulic fracturing operations), or something that is to be avoided at all costs (for example, in key components of mechanical systems). In some other instances, fractures are already present and must be accommodated through novel design (for example, in the case of earthquake ruptures occurring in the Earth’s crust). Because of this wide spectra of motivation, there is a continuous need for the fracture mechanics research community to gain more insight on how fractures are created, evolve, and interact with themselves and with material microstructures. In addition, there are still many remaining questions about how fracture propagation, growth, and interaction impact overall material behavior, and how to develop accurate models to predict this behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Probabilistic data fusion and physics-informed machine learning: A new paradigm for modeling under uncertainty, and its application to accelerating the discovery of new materials

In this report we summarize the work conducted by PI Perdikaris and his group under this Early Career project DE–SC0019116 during the period of 09/01/2018 – 08/31/2023. The central aim of the work was to introduce a new paradigm for scientific data analysis that can seamlessly synthesize rigorous mathematical modeling with data of variable fidelity (e.g., measurements at multiple scales/resolutions or predictions of variable fidelity models) and multiple modalities (e.g., images, time–series, or scattered measurements). The setting we are interested in involves complex systems that are partially observed and whose dynamical behavior could be hard to model or totally unknown. The inherent uncertainty associated with this setting necessitates a departure from the classical deterministic realm of modeling and scientific computation, and, consequently, our main building blocks can no longer be crisp deterministic numbers and governing laws, but instead we must operate with probabilistic models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Characterize traction–separation relation and interfacial imperfections by data-driven machine learning models

Abstract Interfacial mechanical properties are important in composite materials and their applications, including vehicle structures, soft robotics, and aerospace. Determination of traction–separation (T–S) relations at interfaces in composites can lead to evaluations of structural reliability, mechanical robustness, and failures criteria. Accurate measurements on T–S relations remain challenging, since the interface interaction generally happens at microscale. With the emergence of machine learning (ML), data-driven model becomes an efficient method to predict the interfacial behaviors of composite materials and establish their mechanical models. Here, we combine ML, finite element analysis (FEA), and empirical experiments to develop data-driven models that characterize interfacial mechanical properties precisely. Specifically, eXtreme Gradient Boosting (XGBoost) multi-output regressions and classifier models are harnessed to investigate T–S relations and identify the imperfection locations at interface, respectively. The ML models are trained by macroscale force–displacement curves, which can be obtained from FEA and standard mechanical tests. The results show accurate predictions of T–S relations ( R 2 = 0.988) and identification of imperfection locations with 81% accuracy. Our models are experimentally validated by 3D printed double cantilever beam specimens from different materials. Furthermore, we provide a code package containing trained ML models, allowing other researchers to establish T–S relations for different material interfaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Predicting Atomistic Transitions with Transformers

Accurate knowledge of the atomistic transition pathways in materials and material surfaces is crucial for many material science problems. However, conventional simulation techniques used to find these transitions are extremely computationally intensive. Even with large-scale, accelerated material simulations, the computational cost constrains the applicable domain in practice. Machine learning models, with the potential to learn the complex emergent behaviors governing atomistic transitions as a fast surrogate model, have great promise to predict transitions with a vastly reduced computational cost. Here, we demonstrate how transformers can be trained to predict atomistic transitions in nano-clusters. We show how we evaluate physical validity of the predictions and how a multitude of additional, different microstates can be generated by slightly varying the data provided to the model.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Integrative Circuit-Host Modeling of a Genetic Switch in Varying Environments

Synthetic biology is advancing into a new phase where real-world applications are emphasized. There is hence an urgent need for mathematical modeling that can quantitatively describe the behaviors of genetic devices in natural, fluctuating environments. We utilize an integrative circuit-host modeling framework to examine the dynamics of a genetic switch and its host cell in varying environments. For both steady-state and transient cases, we find increasing nutrient reduces the bistability region of the phase space and eventually drives the switch from bistability to monostability. In response, cellular growth and proteome partitioning experience the same transition. Antibiotic perturbations cause the similar circuit and host responses as nutrient variations. However, one difference is the trend of growth rate, which augments with nutrient but declines with antibiotic levels. The framework provides a mechanistic scheme to account for both the dynamic and static characteristics of the circuit-host system upon environmental perturbations, underscoring the intimacy of gene circuits and their hosts and elucidating the complexity of circuit behaviors arising from environmental variations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Orthotropic Viscoelastic Modeling of Polymeric Battery Separator

The separator is one of the most important components for battery safety since its failure can cause short-circuit that may lead to thermal runaway. As such, understanding the mechanical behaviors of separators is critical to improve battery safety. Nevertheless, modeling the mechanical behaviors of separators is challenging since they usually depend on multiple factors such as orientation, temperature and strain rate due to the microstructure of separators. Here we develop an orthotropic linear viscoelastic model for separators, which takes both anisotropy and strain rate effects into account. This model is the first step of a future model that considers more factors such as temperature and large deformation. The model has been implemented as a user material model in a finite element package, where a discretization algorithm is developed to evaluate the stiffness-based hereditary integral with a kernel of Prony series. The user model has been applied to a polypropylene (PP) separator Celgard ® 2400. The simulation results agree well with experimental data, which include the rate dependent uniaxial tensile behavior of the separator in the two in-plane material directions and in shear, and the strain contour patterns of every strain component during biaxial tension.

25 ENERGY STORAGE↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗