Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solid 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Chromium isotope fractionation during reduction of Chromium(VI) by Iron(II/III)-bearing clay minerals

Chromium stable isotope ratios are used to trace the reduction of Cr(VI) to Cr(III) in both ancient and modern systems. However, quantitative interpretation of Cr isotopic signatures has been stymied by the large variability in isotopic fractionation factors for Cr(VI) reduction by different reductants. Here we determine Cr isotope fractionation factors during Cr(VI) reduction by Fe(II/III)-bearing clay minerals, which are abundant in subsurface environments. Several variables were tested: pH, total Fe content of the clay, and the fraction of reduced Fe within the clay (Fe(II)/Fe(total)). The latter controls the standard reduction potential of the clay. Our results demonstrate that neither pH nor total Fe content of the clay have major effects on isotopic fractionation. In contrast, as the effective standard reduction potential of the clay and thus the standard free energy of Cr(VI) reduction become more negative, Cr isotope fractionation factors decrease in magnitude from - 4.9 to -1.3‰ according to a linear free energy relationship. This linear free energy relationship can be predicted from Marcus electron transfer theory and allows first-order predictions of Cr isotope fractionation factors to be made from the standard reduction potential or Fe(II)/Fe(total) of a clay, potentially improving our ability to model Cr isotope signatures in geochemical systems. Chromium is the first isotope system to show such a linear free energy relationship over a diverse range of reductants, including both aqueous and solid-phase reductants, and may provide a model for determining other redox-driven kinetic isotope effects in environmentally important isotope systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗
DOE OSTI2025

Developing Stable Critical Materials and Microstructure for High-Flux and Efficient Hydrogen Production through Reversible Solid Oxide Cells

Reversible Solid Oxide Cells (RSOCs), which operate as either Solid Oxide Fuel Cells (SOFCs) or Solid Oxide Electrolytic Cells (SOECs), hold great promise for clean, high-efficiency energy conversion and hydrogen production. However, their commercial potential is hindered by stability issues arising from temperature-induced materials degradation. While substantial progress has been made in developing durable, reduced-temperature SOFCs - bringing them closer to commercial deployment, SOECs still exhibit substantially higher degradation rates at both the cell and stack levels under practical operating conditions.

08 HYDROGEN↗
DOE OSTI2020

Foams in Hohlraums

Foam materials are starting to find application in laser-heated hohlraums used to drive inertial confinement fusion (ICF) implosions. The disparate scales of the solid foam elements (nm) compared to the typical problem scale (mm) make it expensive to resolve the individual foam elements in hydrodynamic calculations. In this project we developed two new models for treating foam materials – a 3D hydrodynamic model that treats the foam as bits of pre-expanded solid material separated by voids, and a subgrid model that captures the essential physics of laser-heated foams without requiring resolution of the individual solid elements. We also explored the microphysics of the transition from interpenetrating kinetic behavior to fully-collisional fluid-like behavior during the early time heating and expansion of foam elements through a combination of single-fluid and multi-fluid hydrodynamics calculations and particle-in-cell calculations. The model was benchmarked against these calculations and against experimental data on laser-heated chemical and additive-manufactured (AM) foams. Since there was no existing data on laser-heated AM foams, we performed the first experiments on such foams at the Jupiter Laser Facility. Samples of three different types of printed AM foams were heated using a single 0.5 micron laser. The foam densities ranged from 10 to 100 mg/cc (supercritical for 0.5 micron light). We measured the backscattered light (power and spectrum), the transmitted light, and the time-integrated temperature of the foam samples. This benchmarked subgrid foam model is a new capability that should allow routine calculation of ICF designs using foam elements, potentially opening up a new design space.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Double Paddle‐Wheel Enhanced Sodium Ion Conduction in an Antiperovskite Solid Electrolyte

Abstract Antiperovskite structure compounds (X 3 AB, where X is an alkali cation and A and B are anions) have the potential for highly correlated motion between the cation and a cluster anion on the A or B site. This so‐called “paddle‐wheel” mechanism may be the basis for enhanced cation mobility in solid electrolytes. Through combined experiments and modeling, the first instance of a double paddle‐wheel mechanism, leading to fast sodium ion conduction in the antiperovskite Na 3− x O 1− x (NH 2 ) x (BH 4 ), is shown. As the concentration of amide (NH 2 − ) cluster anions is increased, large positive deviations in ionic conductivity above that predicted from a vacancy diffusion model are observed. Using electrochemical impedance spectroscopy, powder X‐ray diffraction, synchrotron X‐ray diffraction, neutron diffraction, ab initio molecular dynamics simulations, and NMR, the cluster anion rotational dynamics are characterized and it is found that cation mobility is influenced by the rotation of both NH 2 − and BH 4 − species, resulting in sodium ion conductivity a factor of 10 2 higher at x = 1 than expected for the vacancy mechanism alone. Generalization of this phenomenon to other compounds could accelerate fast ion conductor exploration and design.

25 ENERGY STORAGE↗
DOE OSTI2019

Surface Engineering by Simultaneous Action of Multiple External Fields (Final Technical Report)

This research project aims at enabling surface engineering strategies based on the fundamental understanding and predictive modeling of the surface morphological response of solid materials subjected to the combined action of multiple external forces. Surface morphological response to mismatch strain in strained-layer heterostructures also is analyzed. Our research focuses on the surface morphological evolution and stability of electrically conducting (metallic, such as Cu) and semiconducting (such as Si) solids under the simultaneous application of mechanical stresses, electric fields, and temperature gradients. We emphasize on identifying the conditions under which the multiply driven surface morphology is stable and on exploring the complexity of the corresponding various morphologically stable surface patterns. We also emphasize on the driven evolution of nanoscale features aiming at stabilizing and controlling nanoscale patterns on surfaces by their manipulation through simultaneously applied multiple external forces. Specifically, we address systematically the morphological response to the combined action of multiple external forces of bulk solid surfaces, surfaces of thin films grown epitaxially on thick or thin substrates, as well as nanoscale surface features such as coherently strained islands grown epitaxially on substrate surfaces. The research is based on a modeling approach that combines theoretical analyses of surface morphological stability with self-consistent dynamical simulations of surface evolution based on properly parameterized continuum and multiscale surface transport models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Morphology and Transport of Multivalent Cation-Exchanged Ionomer Membranes Using Perfluorosulfonic Acid–Ce Z+ as a Model System

Perfluorosulfonic acids (PFSAs) are commonly used as solid polymer electrolyte membranes (PEMs) in electrochemical energy devices, where they are vulnerable to attack by radical species during operation. A popular strategy to combat this problem is to introduce radical scavengers like cerium (Ce) ions that neutralize these species before they attack the PFSA. Such cation doping creates a multi-ion system, in which understanding the mechanisms of cation solvation and transport becomes important for the effective design and utilization of PFSA–cation systems. Ce ions also provide a representative model system for multication-exchanged ionomers in electrochemical systems. In this study, hydration and conductivity measurements, along with X-ray fluorescence and scattering, are employed to elucidate how Ce ion exchange alters PFSA’s ionic solvation, as well as nano- and mesoscale morphologies, which ultimately control its ion transport properties. A molecular transport model is used to deconvolute the impact of Ce ions on the local solvation structure of water in the membrane from mesoscale changes of the transport pathways. The combined experimental and theoretical analysis reveals a nonlinear decrease in conductivity driven by cation solvation at the molecular level and morphological changes at longer length scales. Migration–diffusion coupling, its nonlinear dependence on ion exchange and hydration, and its overall implications for ionomer performance are also discussed. Finally, these findings have the potential to be translated into other mixed cation–ionomer systems for a wide range of energy and environmental devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mesoscale modeling and semi-analytical approach for the microstructure-aware effective thermal conductivity of porous polygranular materials

Here we established a comprehensive modeling approach for investigating the microstructure-aware effective thermal conductivity ($κ_{eff}$) for porous microstructures containing solid particles and gaseous pores. Our approach combines the mesoscale computational modeling framework and the semi-analytical method, allowing for efficient prediction of $κ_{eff}$ for realistic porous microstructures, while considering complicated microstructural thermal conduction pathways effectively in the prediction. We used the diffuse-interface mesoscale computational model to generate extensive simulated $κ_{eff}$ data for realistic digital representations of microstructures with wide ranges of porosity ($f_p$), thermal conductivity of the gas phase ($κ_g$), and thermal conductivity of the solid phase ($κ_s$). From the simulated data, we identified two property variation regimes for $κ_{eff}$: (1) a slow $κ_{eff}$ increase for $κ_s ~ κ_g$; and (2) a faster $κ_{eff}$ increase for $κ_s \gg κ_g$. To capture the key features of the relationship between the microstructure and $κ_{eff}$, we derived a semi-analytical model by introducing structure and intensification factors. The two new factors incorporate the calibrated effective contribution of the solid volume with $κ_s$ and additional interfacial effects into the prediction of $κ_{eff}$, respectively, allowing for consideration of parallel, serial, and interfacial conduction mechanisms effectively. Using the selected simulation data, we quantified key model parameters within the semi-analytical model and verified that the parameterized model exhibits excellent agreement with simulated $κ_{eff}$ for the entire range of the parameter space.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Planet Formation by Gas-assisted Accretion of Small Solids

We compute the accretion efficiency of small solids, with radii 1 cm ≤ R s ≤ 10 m, on planets embedded in gaseous disks. Planets have masses 3 ≤ M p ≤ 20 Earth masses (M ⊕ ) and orbit within 10 au of a solar mass star. Disk thermodynamics is modeled via 3D radiation-hydrodynamics calculations that typically resolve the planetary envelopes. Both icy and rocky solids are considered, explicitly modeling their thermodynamic evolution. The maximum efficiencies of 1 ≤ R s ≤ 100 cm particles are generally ≲10%, whereas 10 m solids tend to accrete efficiently or be segregated beyond the planet’s orbit. A simplified approach is applied to compute the accretion efficiency of small cores, with masses M p ≤ 1 M ⊕ and without envelopes, for which efficiencies are approximately proportional to $M^{2/3}_{p}$. The mass flux of solids, estimated from unperturbed drag-induced drift velocities, provides typical accretion rates dM p /dt ≲ 10 -5 M ⊕ yr -1 . In representative disk models with an initial gas-to-dust mass ratio of 70–100 and total mass of 0.05–0.06 M ⊙ , the solids’ accretion falls below 10 -6 M ⊕ yr -1 after 1–1.5 Myr. The derived accretion rates, as functions of time and planet mass, are applied to formation calculations that compute dust opacity self-consistently with the delivery of solids to the envelope. Assuming dust-to-solid coagulation times of ≈0.3 Myr and disk lifetimes of ≈3.5 Myr, heavy-element inventories in the range 3–7 M ⊕ require that ≈90–150 M ⊕ of solids cross the planet’s orbit. The formation calculations encompass a variety of outcomes, from planets a few times M ⊕ , predominantly composed of heavy elements, to giant planets. The peak luminosities during the epoch of the solids’ accretion range from ≈10 -7 to ≈10 -6 L ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Performance assessment of active insulation systems in residential buildings for energy savings and peak demand reduction

Active insulation systems (AISs) in buildings are envelopes that integrate thermal insulation, thermal energy storage, and controls. Although different designs for AISs have been proposed in the literature, a comprehensive analysis of feasible AISs is lacking. This paper discusses the energy performance, peak demand reduction potential, and performance characteristics of an AIS that uses a concrete wall as thermal mass sandwiched between two solid-state thermal switches (STSs). These STSs change their thermal conductivity using an on/off metal switch to create or break a thermal bridge across the STS. This paper first describes the experimental setup, used to determine the ratio of thermal resistance during R-high (low thermal conductivity) and R-low (high thermal conductivity) states of the STSs. This ratio was then used in whole-building energy simulations to evaluate the performance of AIS walls across different climate zones with/without a freeze timer of 60 min. The timer was added to reduce the number of switches of STSs from one state to another, and hence the energy needed for these switches. Analysis of the switching frequency and interval of STSs, thermal conductivity of walls, impact of wall orientation, and heat transfer through the wall from the use of AIS at different climate zones/locations were performed. The simulation results show that the AIS can achieve energy savings ranging from ~980 to 2,290 kWh in a single-family home with a floor area of ~220 m 2 compared with an IECC 2018 baseline. We found the energy savings was higher in dry climate zones which represent 17% of residential buildings in the United States, compared to humid or marine climate.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗
DOE OSTI2023

Multiphase Modeling in a Parallel Plate Fluidized Bed Receiver for Concentrating Solar Power

A novel high temperature particle solar receiver is developed by using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a packed bed at the lower end, leaving the system from the bottom at a constant flow rate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built and a near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2K.

CSP↗
DOE OSTI2020

IC Project: w19_hossfault “Modelling of stick-slip behavior in sheared granular fault gouge & nonlinear elasticity behavior in cracked solid”

Low frequency earthquakes, non-volcanic tremor, and acoustic emissions are examples of weak seismic signals that may help detect major seismic events, i.e., earthquakes. In a laboratory setting, in which stick-slip events are simulated, acoustic emissions are detected far from the stick-slip events. In both, field and laboratory scale cases, the acoustic emission signals are sourced or detected in a volume remote from the volume that spawns the earthquake. Therefore, it is relevant to establish the causal relationship between signals detected on passive, remote monitors and the dynamics of the elastic structures that launch important seismic events. The work conducted under this IC allocation allowed us to utilize a numerical model that let us follow this causality, i.e., examine and connect the dynamics in a granular system (fault gouge), to signals detected on passive remote monitors. It was demonstrated that stress chains are key in the dynamics of granular systems and that their evolution is the source of the acoustic emission.

58 GEOSCIENCES↗
DOE OSTI2026

Development and Application of High-Fidelity Models for Heterogeneous CO2 Frost Formation

Carbon America has developed a cryogenic carbon capture technology ("FrostCC") that separates CO2 from point source emissions by solidifying it at cold temperatures through preferential desublimation. Cooling is achieved through a series of interlinked compression, heat exchange, and expansion operations. In the current system, frosting of CO2 happens in heat exchangers, followed by CO2 recovery in a separate extraction step. In this work, multiphysics computational fluid dynamics (CFD) models are developed and validated for compressible and low Mach flows to simulate the formation of solid CO2 in flue gas flowing in a heat exchanger geometry. The models track the mass transfer rate of CO2 from gas phase to solid phase, heat released from desublimation, and the evolution of the solid CO2 layer. Simulations are used to answer scientific questions related to the angle of heat exchanger pipes, where buoyancy effects from flow velocity and pipe orientation influence CO2 frosting. Results show that upwardly angled pipes produce notably different flow structures compared to horizontal or vertical configurations, and that carbon capture efficiency correlates with buoyancy effects for pipe angles within plus or minus 23 degrees of horizontal.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo -7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation

The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.

02 PETROLEUM↗
DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗