Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled 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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Characterizing Reactive Transport Behavior in a Three-Dimensional Discrete Fracture Network

While several studies have linked network and in-fracture scale properties to conservative transport behavior in subsurface fractured media, studies on reactive transport cases remain relatively underdeveloped. In this study, we explore the behavior of an irreversible kinetic reaction during the interaction of two solute plumes, one consisting of species A and the other species B. When the plumes converge, these species react kinetically to form a new species C via A+B→kC. This reactive system is studied using a three-dimensional discrete fracture network (DFN) model coupled with reactive Lagrangian particle tracking. We find that the interplay of network topology and chemical properties of the reactive solutes controls reactive transport processes. The network topology drives species A and B together, and the chemical properties dictate whether and how quickly a reaction occurs. Results demonstrate that reactions are most likely to occur in high-velocity fractures that make up the network backbone. The interplay between species’ chemical properties and transport is characterized by a non-dimensional Damköhler (Da) number. We show that the spatial distribution of reactions is sensitive to Da, which subsequently influences late-time tailing behavior in outlet breakthrough time distributions. The results of this study provide initial insights into how an irreversible reaction occurs during transport in a fracture network, using a methodology that can be applied to study reactive transport in a wide range of fractured media environments and contexts.

58 GEOSCIENCES↗
DOE OSTI2021

Projected Changes in Temperature and Precipitation Over the United States, Central America, and the Caribbean in CMIP6 GCMs

The Coupled Model Intercomparison Project Phase 6 (CMIP6) dataset is used to examine projected changes in temperature and precipitation over the United States (U.S.), Central America and the Caribbean. The changes are computed using an ensemble of 31 models for three future time slices (2021–2040, 2041–2060, and 2080–2099) relative to the reference period (1995–2014) under three Shared Socioeconomic Pathways (SSPs; SSP1-2.6, SSP2-4.5, and SSP5-8.5). The CMIP6 ensemble reproduces the observed annual cycle and distribution of mean annual temperature and precipitation with biases between – 0.93 and 1.27 °C and – 37.90 to 58.45%, respectively, for most of the region. However, modeled precipitation is too large over the western and Midwestern U.S. during winter and spring and over the North American monsoon region in summer, while too small over southern Central America. Temperature is projected to increase over the entire domain under all three SSPs, by as much as 6 °C under SSP5-8.5, and with more pronounced increases in the northern latitudes over the regions that receive snow in the present climate. Annual precipitation projections for the end of the twenty-first century have more uncertainty, as expected, and exhibit a meridional dipole-like pattern, with precipitation increasing by 10–30% over much of the U.S. and decreasing by 10–40% over Central America and the Caribbean, especially over the monsoon region. Seasonally, precipitation over the eastern and central subregions is projected to increase during winter and spring and decrease during summer and autumn. Over the monsoon region and Central America, precipitation is projected to decrease in all seasons except autumn. The analysis was repeated on a subset of 9 models with the best performance in the reference period; however, no significant difference was found, suggesting that model bias is not strongly influencing the projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

An ontology to represent synthetic building occupant characteristics and behavior

Since the introduction of the occupant behavior Drivers-Needs-Actions-Systems (DNAS) framework in 2013, researchers have used the framework or further developed it based on their case studies, which include efforts to collect new data on occupant behaviors. The effort is often costly for the relatively few new data points added. Problems emerge when the already collected data do not meet the modelers' interoperability requirements. Previous studies addressed this issue by developing more sophisticated ontologies that enable integration with other datasets and synthetic data methodologies that would meet unique research applications. This paper presents an extension of the DNAS framework for the representation of synthetic occupant data to support various applications and use cases across the building life cycle. An agent-based modeling application is one of our motivations that requires more elaborate characteristics of an occupant-agent or a group-of-agent. The extension, built upon a review of the literature, introduces new elements to the framework that fall into five categories, including socio-economic, geographical location, activities, subjective values, and individual and collective adaptive actions. On-going research includes identifying occupant datasets and developing data fusion methods to generate synthetic occupants, as well as to demonstrate its applications in agent-based modeling coupled with building performance simulation.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2021

Simultaneous x-ray and neutron 4D tomographic study of drying-driven hydro-mechanical behavior of cement-based materials at moderate temperatures

The drying of cement-based materials is intimately related to their durability, which has significant economic, social and environmental repercussions. The evolution of the saturation of the pore network and the associated drying shrinkage are in fact leading causes of cracking and of the ensuing penetration of aggressive chemicals. This process is highly heterogeneous, due to the thermo-hydric spatial gradients developing in the material from the exposed surfaces to its core and because of local effects, driven by the intrinsically heterogeneous micro-structure (e.g., by the distribution of pores and aggregates). It follows that macroscopic, sample-scale measurements cannot fully disclose the complexity of the underlying processes. In the last few decades, significant advances in full-field techniques have allowed an unprecedented insight into these local processes. For cement-based materials, x-ray and neutron tomography lend themselves as ideal, and highly complementary, tools for the study of their thermo-hydro-mechanical behavior. Notably, the high sensitivity to density variations of x-ray imaging gives access to the developments of fractures, in 4D (3D + time). On the other hand, neutron tomography allows the study of the evolution of the moisture field in 4D, thanks to its high hydrogen sensitivity. The combination of these two techniques provides a unique insight in thermo-hydro-mechanical couplings, e.g., the effect of cracks on the water content field. This contribution presents novel 5D datasets (3D tomographies along time, plus truly simultaneous x-ray and Neutron rapid acquisitions) in-operando of a cement paste and of a concrete sample heated at moderate temperatures (up to 140{sup ∘}C). The analysis of this 5D data-set (once aligned in time and across modalities) allows for example a quantification of the 4D moisture profiles which were found to predict an overall water loss at hydric equilibrium coherent with the corresponding analytical analysis. In the cement paste sample, the x-ray dataset captures the evolution of an extensive cracking network, opening and propagation toward the core of the sample. A novel analysis procedure is here proposed which allows the extraction of these fractures and the analysis of their interplay with local drying as captured through neutron imaging. This for example reveals the depth of penetration of drying in the vicinity of the fractures along time, which is essential for the assessment and calibration of hydro-mechanical coupled models.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Microstructure prediction for Ti-22Al-25Nb in laser powder bed fusion

This work presents a physics-informed framework for predicting solidification morphology and defect susceptibility in additively manufactured Ti–22Al–25Nb across a broad processing space. The framework integrates solidification microstructure selection (SMS) analysis with a single-track defect-based printability map to establish a unified methodology linking processing parameters to both interfacial morphology and manufacturability. Thermal gradients G and solidification rates R are first computed using the Thermo-Calc Additive Manufacturing (TC-AM) module, a finite-interface-dissipation (FID) phase-field (PF) model coupled with CALPHAD method is then employed to systematically distinguish planar and dendritic regimes as functions of $G$ and $R$. By superimposing the printability map onto the morphology projections, a comprehensive process–structure framework is obtained. Across most processing conditions, the predicted microstructure is predominantly dendritic, while planar growth emerges only under selected laser power $P$ and scan speed $v$ combinations. In addition to morphology classification, the framework quantifies the dendritic area fraction and introduces a width-based morphology descriptor to characterize the spatial extent of planar/dendritic regions within the melt pool. It provides mechanistic insight into the interplay between solidification physics and defect formation, offering practical guidance for parameter selection and microstructural control in Ti–22Al–25Nb additive manufacturing (AM).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modelling uranium dioxide corrosion under repository conditions: A pore-scale study of the chemical and thermal processes

This paper investigates the corrosion behaviors of spent uranium dioxide (UO 2 ) fuel when placed in geological repositories. We performed pore-scale reactive transport simulations of the UO 2 corrosion process in a defective fuel rod with different orientations of breach or fracture on its clad. It is found that the corrosion rate has strong dependency on the breach orientation. The highest corrosion rate of UO 2 is calculated when the angle between the groundwater flow direction and the clad fracture reaches 180°. UO 2 corrosion simulations with different flow rates are compared. The results show the higher flow rate accelerates the corrosion of UO 2 fuel. The effect of pH on the corrosion process is also determined. It is indicated that higher pH can limit the corrosion of UO 2 by reducing the reaction rate. The dependence of reaction rate and time on reactive surface area is explored. Spent fuel with lower surface area demonstrates longer lifetime under corrosion conditions. A numerical model coupling thermal conduction and chemical reactions is developed to assess the impact of temperature on the process of UO 2 corrosion. The results show that higher temperature leads to larger corrosion rates for UO 2 . Here, the predicted reaction rates are higher in comparison with the isothermal results. Hence, our results can help improve the fundamental understanding of UO 2 corrosion in geological repositories for long-term storage of spent nuclear fuels and provide guidance for the safe operations and selection of appropriate repositories.

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

Comparative verification of hydro-mechanical fracture behavior: Task G of international research project DECOVALEX–2023

Numerical simulations become a necessity when experimental approaches cannot cover the required physical and time scale of interest. One of such area is a simulation of long-term host rock behaviors for nuclear waste disposal and simulation tools involved in the assessment must go through rigorous validation tests. The DECOVALEX project (Development of COupled models and their VAlidation against EXperiments) is dedicated to this purpose by international participants. This work is part of the ongoing phase DECOVALEX–2023 (D–2023, Task G) particularly aiming to simulate fracture behaviors under various conditions. Here, in this study, we cross-verified a variety of numerical methods including continuous and discontinuous approaches against four benchmark exercises with emphasis on numerical accuracy and parameterization of the various numerical approaches. The systematic inter-comparisons of test cases highlight advantages and disadvantages of the different numerical models. Numerical details on discretization effects (e.g. mesh density and orientation) and domain size were investigated in detail for practical applications. It became evident that meticulous attention to mesh resolution and domain size is imperative for achieving accurate numerical simulations, even for static cracks. Moreover, when comparing numerical methods to closed-form solutions for static cracks, all models successfully reproduced the maximum crack opening but encountered challenges near the crack tips. Finally, the paper discusses how to convert between and therefore compare parameters of various numerical approaches. Our benchmark studies reveal that each model necessitates a distinct number of parameters, even in simple scenarios like static crack aperture benchmarks. It is generally more practical to employ fewer parameters to mitigate model over-parameterization and enhance experimental feasibility.

58 GEOSCIENCES↗
DOE OSTI2022

Predicting extreme anisotropy and shape variations in impact testing of tantalum single crystals

We report recent Taylor cylinder impact tests carried out for Ta single crystals showed strong variations in dimensional changes for different crystallographic directions aligned with the cylindrical axis. In order to capture the effect of crystallography on the deformation characteristics and final shapes of the impacted cylinders, a single crystal material subroutine is adapted and embedded in the solid mechanics/dynamics Finite Element solver Abaqus to simulate the aforementioned single crystal Ta Taylor impact experiments. Details of the coupled model implementation, and insights on the role played by single crystal anisotropic flow on the deformation behavior across a broad range of strain rates and temperatures for different single crystal orientations are presented and discussed. We demonstrate the predictive capability of the adopted crystal plasticity model to capture the significant role played by crystal orientation-induced anisotropy, as well as strain hardening and adiabatic heating, on the dynamic deformation response of crystalline materials. This re-emphasizes the need of microstructure-aware models to improve the accuracy of simulations for high-consequence engineering design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of thermal environments on fast charging Li-ion batteries

Battery thermal management systems (BTMSs) are expected to keep the battery temperature at a moderate level (~30 °C) to minimize the thermally exacerbated degradation. However, during fast charging, a strong cooling system is required to restrict the temperature rise of Li-ion batteries (LiBs), which significantly increases the cost and weight of battery packs, and induces a large temperature variation inside the battery. Here, in this work, we find that all these drawbacks could be relieved by allowing LiBs to charge at higher temperatures. Since the fast charging of a LiB only takes a tiny fraction of its lifetime, the aging rate is limited even at a charging temperature of 60 °C. Three types of thermal environments are proposed: kept constant at 30 °C, preheated to 60 °C, and adiabatic fast charging. With an experimentally validated electrochemical-thermal (ECT) coupled model, we explore the interplay between thermal management and the fast-charging performance. It is found that a gradually increasing temperature profile is the best option to balance the lithium plating and thermal management of the battery. Combining adiabatic fast charging with a preheating step, we can achieve minimal cooling need, perfect temperature uniformity within a battery, and fast-charging capability simultaneously.

25 ENERGY STORAGE↗
DOE OSTI2024

Effect of Si impurities on microstructure and tensile properties of a cast Al-Mg-Fe alloy

Al-Mg alloys are attractive for structural castings owing to their superior strength and ductility in the as-cast state. However, given the tight tolerance for impurities, Al-Mg alloys produced from secondary sources still face multiple challenges. Here, we report the effect of increased Si impurity (0.05–1.6 wt%), which is commonly found in secondary Al sources, on microstructure and tensile properties of a cast Al-4.3Mg-1.6Fe (wt%) alloy, commercially referred to as Castaduct-42 alloy. Microstructural characterization revealed that Si addition increased the volume fraction, size, and aspect ratio of primary Al13Fe4 intermetallic particles as well as the volume fraction of other binary and ternary eutectic phases. Tensile testing results demonstrated that increasing Si impurities from 0.05 to 1.6 wt% reduced ductility from 14.3 ± 1.4 % to 2.2 ± 1.0 %. A particle cracking damage accumulation model coupled with failure analysis indicated primary Al13Fe4 intermetallic particles to be the major contributing factor to the deterioration in ductility with increasing Si content. Additionally, the stabilizing effect of Si on primary Al 13 Fe 4 was inconsistent with the CALPHAD calculation results based on existing CALPHAD databases which predict a slightly decreasing primary Al13Fe4 phase fraction with increasing Si concentration. This work provides new insights into the phase stability and mechanical behavior of the lesser studied Al-Mg-Fe-Si alloy system that will contribute to the development of sustainable cast Al-Mg based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lattice misorientation evolution and grain refinement in Al-Si alloys under high-strain shear deformation

The starting microstructure of an alloy can be tailored to achieve varying degrees of microstructural refinement and enhance mechanical properties through severe plastic shear deformation during solid-phase processing. Crystal plasticity-based grain misorientation modeling, coupled with systematic pin-on-disk tribometry-based subsurface shear deformation experiments on as-cast Al-xSi alloys (x = 0, 1, 4 at %), was conducted. The post-deformation microstructural analysis, through a combined computational and experimental approach, conclusively shows that the initial volume fraction of the hard Si phase enhances the evolution of local lattice misorientation, leading to efficient grain refinement during severe plastic shear deformation. The shear–deformation-induced nanostructure resulted in more than double the nanoindentation hardness in the processed alloy.

Gwalani, Bharat↗
DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Heavy and light flavor jet quenching in heavy-ion collisions in a perturbative QCD approach

In this work, we study the flavor dependence of jet quenching in heavy-ion collisions. A next-to-leading-order perturbative QCD framework is applied to account for both quark and gluon contributions to light and heavy flavor hadron productions. A linear Boltzmann transport model coupled to hydrodynamic simulation is utilized to study the nuclear modification of heavy and light flavor jets in the dynamically evolving quark-gluon plasma. Within our jet quenching framework, we obtain a nice description of the nuclear modification factors of charged hadrons, D mesons, B mesons and B–decayed D mesons over a wide range of transverse momentum (8-300 GeV). Our result shows that perturbative QCD is sufficient to explain the color, mass and energy dependence of parton energy loss and jet quenching in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Collision-energy dependence of deuteron cumulants and proton-deuteron correlations in Au+Au collisions at RHIC

We report the first measurements of cumulants, up to 4th order, of deuteron number distributions and proton-deuteron correlations in Au+Au collisions recorded by the STAR experiment in phase-I of Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider. Deuteron cumulants, their ratios, and proton-deuteron mixed cumulants are presented for different collision centralities covering a range of center-of-mass energy per nucleon pair $\sqrt{s_{NN}} = 7.7$ to 200 GeV. It is found that the cumulant ratios at lower collision energies favor a canonical ensemble over a grand canonical ensemble in thermal models. An anti-correlation between proton and deuteron multiplicity is observed across all collision energies and centralities, consistent with the expectation from global baryon number conservation. The UrQMD model coupled with a phase-space coalescence mechanism qualitatively reproduces the collision-energy dependence of cumulant ratios and proton-deuteron correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Implication of stratospheric aerosol geoengineering on compound precipitation and temperature extremes in Africa

Three Coupled Model Intercomparison Project 5 (CMIP5) models that simulated the G4 experiment of the Geoengineering Model Intercomparison Project (GeoMIP) were used to investigate the impact of stratospheric aerosol injection (SAI) on combined temperature and precipitation extremes in Africa that can have greater negative impacts on human and the environment than individual rainfall or temperature extremes. The examined compound extremes included the dry (R warm|dry and R cold|dry ) and wet (R warm|wet and R cold|wet ) modes assessed during the injection (SAI, 2050-2069) and post-injection (postSAI, 2070-2089) periods compared with the historical period (1986-2005). We found a significant projected change in the occurrence of both wet and dry modes during SAI and postSAI related to the historical period. The magnitude and sign of this change depend on the season and the geographical location. During the SAI and postSAI, the wet (R warm|wet and R cold|wet ) modes are projected to be significantly lower while the dry modes are noted to increase in a large part of African continent depending on the season and the geographical location and may consequently leads to an increase of the droughts prone areas. Furthermore, the termination effect is noted to re duce the occurrence of dry modes, which may reduce the potential negative effects of the injection after halting. As the effect may vary from one region to another and according to the season, it suggested assessing the key sector impacts of SAI. Thus, this change in dry modes due to SAI could affect all activities which depend on water resources such as water supply, agriculture and food production, energy demand, and production with adverse effects on health, security, and sustainable development, but this needs to be assessed and quantified at regional scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Advancing the Performance of Lithium-Rich Oxides in Cooperation with Inherent Complexities: Engineering Domain Structure

The nanodomain structure of lithium- and manganese-rich, composite cathode materials has been systematically altered through synthesis conditions while keeping larger-scale morphological differences to a minimum. Clear changes in electrochemical performance across the samples studied are observed, especially with respect to the anomalous impedance at low states-of-charge. Atomic-scale modeling, coupled to electrochemical measurements and physical characterization, reveals the local consequences of the high-voltage activation charge process as a function of specific domain structures. Furthermore, the results explain the influence of different postactivation structures on the insertion of Li-ions and the associated impedance during discharge. This work adds to our series of studies on lithium- and manganese-rich oxides, demonstrating that the inherent performance of LMRs can be greatly enhanced through rational, highly controlled synthetic strategies based on an understanding of synthesis–structure–property relationships across length scales.

25 ENERGY STORAGE↗
DOE OSTI2025

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗