Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase 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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Vacancy and interstitial interactions with crystal/amorphous, metal/covalent interfaces

In this work, we use atomistic simulations to investigate the interaction of vacancies and interstitials with interfaces between a crystalline metal and an amorphous, covalently-bonded solid. We select the gold (Au)/silicon (Si) binary system as a model material and construct interface models along two different facets of crystalline Au and with amorphous Si (a-Si) created at three different quench rates. We compute formation energies of vacancies, self-interstitials, and interstitial impurities as a function of position relative to the interface and find that they have markedly lower values near the interface than in the interior of the adjoining phases. We conclude that crystal/amorphous, metal/covalent interfaces may be as effective at removing radiation-induced point defects as interfaces in polycrystalline metals composites. Moreover, irrespective of interface character, the average formation energies of all point defects at all the Au/a-Si interfaces we investigated are comparable. Thus, unlike in polycrystalline metals, where an interface’s crystallographic character has a marked effect on its interactions with point defects, all interface types in crystal/amorphous, metal/covalent composites may be equally effective at absorbing all radiation-induced defects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiterminal Inverse AC Josephson Effect

When a Josephson junction is exposed to microwave radiation, it undergoes the inverse AC Josephson effect–the phase of the junction locks to the drive frequency. As a result, the I–V curves of the junction acquire “Shapiro steps” of quantized voltage. If the junction has three or more superconducting contacts, coupling between different pairs of terminals must be taken into account and the state of the junction evolves in a phase space of higher dimensionality. Here, we study the multiterminal inverse AC Josephson effect in a graphene sample with three superconducting terminals. We observe robust fractional Shapiro steps and correlated switching events, which can only be explained by considering the device as a completely connected Josephson network. We successfully simulate the observed behaviors using a modified two-dimensional RCSJ model. Our results suggest that multiterminal Josephson junctions are a playground to study highly connected nonlinear networks with novel topologies.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Magnetized laser–plasma interactions in high-energy-density systems: Parallel propagation

In this study, we investigate parametric processes in magnetized plasmas, driven by a large-amplitude pump light wave. Our focus is on laser–plasma interactions relevant to high-energy-density (HED) systems, such as the National Ignition Facility and the Sandia MagLIF concept. We present a self-contained derivation of a “parametric” dispersion relation for magnetized three-wave interactions, meaning the pump wave is included in the equilibrium, similar to the unmagnetized work of Drake et al., Phys. Fluids 17, 778 (1974). For this, we use a multi-species plasma fluid model and Maxwell's equations. The application of an external B field causes right- and left-polarized light waves to propagate with differing phase velocities. This leads to Faraday rotation of the polarization, which can be significant in HED conditions. Phase-matching and linear wave dispersion relations show that Raman and Brillouin scattering have modified spectra due to the background B field, though this effect is usually small in systems of current practical interest. We study a scattering process we call stimulated whistler scattering, where a light wave decays to an electromagnetic whistler wave (ω≲ω ce ) and a Langmuir wave. This only occurs in the presence of an external B field, which is required for the whistler wave to exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Petrographic and spectral study of hydrothermal mineralization in drill core from Hawaii: A potential analog to alteration in the martian subsurface

Abstract Continuous rock core was collected for 1764 m (5786’) on the Pohakuloa Army Training base near the center of the big island of Hawaii. The core traverses basaltic lava flows from the volcano's shield-building phase, and perched aquifers and higher temperature groundwaters were encountered. The collected samples record water-rock interactions of basaltic materials in a setting that may be a model for groundwater interactions on Mars. We collected visible and infrared point spectra of materials in the lowest portion of the core, where alteration was noted to become more prominent. We identified three types of phyllosilicate spectral signatures and three types of zeolites. The phyllosilicates show similarity to those identified on Mars using data from the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM). Based on the field survey, 25 depths were selected for sampling and petrographic analysis of thin sections. The spectral data and thin section work have a strong agreement in the types of materials identified by the two different techniques. Both the spectral and petrographic data indicate low to moderate temperature geothermal alteration occurred in the lower half of the core. The field spectra are a useful reconnaissance tool for selecting mineralogic diversity for subsequent higher resolution and more time-consuming laboratory analysis.

Geochemistry & Geophysics↗
DOE OSTI2025

Aerosol Thermodynamics (Aerosol Liquid Water Concentration and Aerosol pH) at S2 during CoURAGE

These data report aerosol liquid water content (ALWC) and aerosol pH for the S2 (Mt. Airy) site during CoURAGE for April - June 2025. The ISORROPIA-II aerosol thermodynamic equilibrium model was used to compute ALWC and pH. The model was run in forward mode, with solids formation disabled (metastable mode), according to Pye et al. (Atmospheric Chemistry and Physics, 2020). Information about the model can be found at https://www.epfl.ch/labs/lapi/models-and-software/isorropia/, and in Fountoukis and Nenes (Atmospheric Chemistry and Physics, 2007). Inputs to the model were: (1) measured Temperature and relative humidity, both from the DoE ARM “metwxt” data products at S2, (2) aerosol chemical composition from the ARM ACSM, and (3) gas-phase ammonia concentrations. For ISORROPIA-II, ambient RH values above 0.995 were input to the model as 0.995; therefore, caution should be exercised interpreting ALWC when RH was above 0.995. Aerosol composition or NH3(g) concentrations below the LOD were input into the model as 0.5*LOD. Model runs were only conducted for data points at 30-min resolution with all three inputs (meteorology, ACSM, and NH3). If any of the three model inputs were missing, aerosol thermodynamic outputs are flagged as -999. Similarly, aerosol thermodynamic outputs are flagged as -999 if measured NH3(g) and aerosol ammonium (NH4+) were simultaneously below the respective LODs. The ALWC includes liquid water from inorganic aerosol components, calculated directly by ISORROPIA-II, as well as liquid water from organics, estimated from the ACSM organics assuming a kappa value of 0.12 according to Guo et al. (Atmospheric Chemistry and Physics, 2015).

Aerosol Liquid Water Content↗
NASA NTRS1979

Siple station experiments on wave-particle interactions in the magnetosphere

Natural and controlled whistler-mode signals have been used to study nonlinear mechanisms of wave growth and wave-wave interactions (WWI) in the magnetosphere; three general classes of WWI (triggering, suppression, and entrainment) are identified and interpreted in terms of a model based on cyclotron resonance interaction. This model is also used to estimate the wave field intensity associated with different types of WWI. A new type of triggered emission, the band-limited impulse (BLI) is interpreted in terms of the switching of phase-bunched currents. In addition, an experiment to find a threshold for the excitation of the coherent wave instability is discussed, and observed VLF wave-induced transient bursts of X-rays, light and E-region ionization enhancements are considered with regard to a step function wave interacting with all resonant electrons in a given energy range.

Helliwell, R. A.↗
NASA NTRS1988

Space Station Active Thermal Control System modeling

The Space Station Active Thermal Control System (ATCS) has been modeled using modified SINDA/SINFLO programs to solve two-phase Thermo-fluid problems. The modifications include changes in several subroutines to incorporate implicit solution which allows larger time step as compared to that for explicit solutions. Larger time step saves computer time but involves larger computational error. Several runs were made using various time steps for the ATCS model. It has been found that for a reasonable approach, three times larger time step as compared to that used in explicit method is a good value which will reduce the computer time by approximately 50 percent and still maintain the accuracy of the output data to within 90 percent of the explicit values.

Hye, Abdul↗
NASA NTRS1994

Ground Based Studies of Gas-Liquid Flows in Microgravity Using Learjet Trajectories

A 1.27 cm diameter two phase gas-liquid flow experiment has been developed with the NASA Lewis Research Center to study two-phase flows in microgravity. The experiment allows for the measurement of void fraction, pressure drop, film thickness and bubble and wave velocities as well as for high speed photography. Three liquids were used to study the effects of liquid viscosity and surface tension, and flow pattern maps are presented for each. The experimental results are used to develop mechanistically based models to predict void fraction, bubble velocity, pressure drop and flow pattern transitions in microgravity.

Bousman, W. S.↗
NASA NTRS2004

The Role of Rhenium on the Temporal Evolution of the Nanostructure of a Model Ni-Al-Cr-Re Superalloy

Rhenium (2 at.%) additions to a model Ni-8.5 at.% Cr-10 at.% Al alloy are studied with respect to its effects on the temporal evolution of the nanostructure and the partitioning behavior of the four elements between the gamma (fcc) and gamma' (L1(sub 2)) phases. Chemical evolution of this quaternary alloy aged at 1073 K from 0.25 to 264 h, is investigated by three-dimensional atom-probe (3DAP) microscopy. The morphology of gamma'-precipitates remains spheroidal, even at an aging time of 264 h. The results demonstrate that Re slows the coarsening of gamma'-precipitates, in comparison to the ternary Ni-10 at.% A1-8.5 at.% Cr alloy at 1073 K.

Yoon, Kevin E.↗
DOE OSTI2025

Experimental observation of nonlinear relation between pressure and water flux is consistent with the solution-diffusion model

In several recent studies, it has been proposed that the fundamental understanding of penetrant transport in dense polymer membranes occurring via the solution-diffusion model, which has been the generally accepted theoretical framework for describing penetrant transport in such materials for the past several decades, is flawed. An alternate mechanistic framework based on the idea of two-phase flow in a porous medium (i.e., pore-flow) has been broadly advanced instead, with proponents of this approach claiming that the pore-flow theoretical framework provides the necessary mechanistic insight to design novel polymeric membrane materials for emerging applications. In this study, we show experimental results for hydraulic permeation of water that are entirely consistent with the solution-diffusion theory, without modification, for three dense polymeric membranes: crosslinked poly(ethylene glycol diacrylate) (XLPEGDA), Nafion 117 ionomer in the sodium counterion form (Nafion 117-Na), and cellulose acetate (CA). By measuring water flux at transmembrane pressures up to 240 bar, we observe a nonlinear relationship between the transmembrane pressure (TMP) and water flux, J w , for XLPEGDA and Nafion 117-Na, while this relationship is linear for CA. We demonstrate that the behavior of these three materials is described via the solution-diffusion model. According to the solution-diffusion model, flux is, to a good approximation, proportional to the transmembrane concentration difference induced by the pressure difference across the membrane, rather than to TMP itself. Water sorption isotherms are reported for all three materials. They further justify the nonlinear relationship between TMP and J w observed in XLPEGDA and Nafion 117-Na, emphasizing that the nonlinearity in the flux/TMP relationship stems from nonlinearities in the sorption isotherm with pressure. Additionally, the relationship between water flux and TMP can be predicted, a priori, with no adjustable parameters when a predictive model for the diffusion coefficient of water is employed in conjunction with the experimental water sorption isotherms in the solution-diffusion model. Furthermore, our results demonstrate the validity of the solution-diffusion model to describe transport of penetrants in dense polymer membranes, while highlighting the sensitivity of the solution-diffusion model to the many physical and mathematical simplifications commonly applied to the theory in literature.

materials↗
DOE OSTI2026

A careful examination of closure models in Euler–Lagrange Simulations of compressible multiphase flow in a planar shock particle curtain problem

In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.

Compressible flow↗
DOE OSTI2020

Effect of heat treatment on the tensile behavior of selective laser melted Ti-6Al-4V by in situ X-ray characterization

Selective Laser Melted Ti-6Al-4V (as-SLMed) exhibits decreased yield strength, increased work hardening, and increased ductility after heat treatment at 730 degrees C (HT-730) or 900 degrees C (HT-900) for 2 h. To understand the change of mechanical properties, in situ high energy X-ray diffraction (HEXRD) is used to examine the phase composition, load partitioning, slip system activity, and dislocation density evolution in all three specimens. The as-SLMed specimen is dominated by martensitic α'. After heat treatment, α' partly or fully decomposes into α+β, reducing the yield strength. In HT-730, beta precipitates with confined size show much higher lattice strain than the α'/α matrix during deformation; in HT-900, the lattice strain difference is mostly eliminated. This is a key reason for the increased ductility in HT-900. From the anisotropic lattice strain development, basal slip is identified as the easiest slip system in α'/α. Using an elasto-plastic self-consistent (EPSC) model, the critical resolved shear stress ratio between prismatic slip and basal slip (CRSS prism /CRSS basal ) is estimated to be 1.31 and 1.16 in the as -SLMed and the HT-900 specimens, respectively. Further, the alpha phase in HT-900 is able to activate multiple slip systems and accumulate more dislocations during plastic deformation. This explains why HT-900 has better ductility and higher work hardening rate than the other two specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Artificial-Intelligence Aided Design and Synthesis of Novel Layered 2D Multi-Principal Element Materials for Energy Storage (Final Report)

This DOE-EPSCoR project aimed to predict, synthesize, and characterize novel layered two-dimensional (2D) high-entropy materials (HEMs). These 2D-HEMs, composed of multiple principal elements in nearly equal concentrations, are distinct from traditional 2D materials (typically containing two or three elements) and conventional alloys (dominated by a single primary element with minor secondary additions). Their unique structural and compositional features enable significant lattice strain accommodation, resulting in enhanced electrode performance and potential applications in catalysis, hydrogen storage, sensing, quantum information technologies, and flexible electronics. The research focused on addressing four fundamental questions: (i) What combinations of elements can form stable and synthesizable 2D-HEMs? (ii) What mechanisms drive the stability and synthesizability of crystalline single-phase 2D-HEMs? (iii) How do local chemical disorder and defects influence the macroscopic electronic and mechanical properties? and (iv) What charge storage mechanisms are active in selectively synthesized 2D-HEMs for battery and supercapacitor electrode applications? To achieve these goals, the project employed an integrated theory-experiment approach, incorporating high-throughput first-principles calculations, theoretical modeling, data mining, experimental synthesis, and advanced characterization techniques. The advanced computing resources and state-of-the-art experimental characterization facilities at Oak Ridge National Laboratory (ORNL) were leveraged through collaboration. Beyond scientific advancements, the project contributed to workforce development. Two postdoctoral researchers and three graduate students at the University of Maine were trained through co-advising by ORNL scientists and collaborative interactions, strengthening their expertise in cutting-edge materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensedphase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi’s golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid–porphyrin–C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor–acceptor, donor–ground, and acceptor–ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC 60 /THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC 60 /THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Integrated control of edge localized modes and divertor flux using mixed toroidal harmonic resonant magnetic perturbations in EAST

Mixed harmonic resonant magnetic perturbations (RMPs) for integrated edge localized modes (ELMs) and divertor flux control are demonstrated in EAST target plasmas of low input torque and normalized beta β N ~1.7–1.9, which are close to the equivalent value in ITER high Q operation. The applied RMPs are designed to combine a static harmonic of the toroidal mode number n = 3 with a static or rotating harmonic of n = 2. ELM suppression is achieved without a drop of plasma energy confinement, and tungsten concentration is effectively reduced during the application of RMPs. With mixed harmonics, the toroidal varying steady state heat and particle fluxes on the divertor target can be modified with the rotating n = 2 harmonic, which agrees with the numerical modeling of three-dimensional magnetic topology, with plasma responses being taken into account. ELM suppression correlates with the times of larger n = 3 response with mixed n = 2 and n = 3 RMPs. The mixture of harmonics and the rotating n = 2 harmonic does not require additional coil current because the variation is only in the upper-lower coil current phase space. Furthermore, these results further affirm the effectiveness of integrated ELM and divertor flux control using RMPs with mixed harmonics and improve the understanding of the role of plasma responses in ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Maximum CO 2 diffusion inside leaves is limited by the scaling of cell size and genome size

Maintaining high rates of photosynthesis in leaves requires efficient movement of CO 2 from the atmosphere to the mesophyll cells inside the leaf where CO 2 is converted into sugar. CO 2 diffusion inside the leaf depends directly on the structure of the mesophyll cells and their surrounding airspace, which have been difficult to characterize because of their inherently three-dimensional organization. Yet faster CO 2 diffusion inside the leaf was probably critical in elevating rates of photosynthesis that occurred among angiosperm lineages. Here we characterize the three-dimensional surface area of the leaf mesophyll across vascular plants. We show that genome size determines the sizes and packing densities of cells in all leaf tissues and that smaller cells enable more mesophyll surface area to be packed into the leaf volume, facilitating higher CO 2 diffusion. Measurements and modelling revealed that the spongy mesophyll layer better facilitates gaseous phase diffusion while the palisade mesophyll layer better facilitates liquid-phase diffusion. Our results demonstrate that genome downsizing among the angiosperms was critical to restructuring the entire pathway of CO 2 diffusion into and through the leaf, maintaining high rates of CO 2 supply to the leaf mesophyll despite declining atmospheric CO 2 levels during the Cretaceous.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Interaction Effects on the Dynamical Anderson Metal-Insulator Transition Using Kicked Quantum Gases

Understanding the interplay of interaction and disorder in quantum transport poses long-standing scientific challenges for theory and experiment. While highly controlled ultracold atomic platforms combining atomic interactions with spatially disordered lattices have led to remarkable advances, the extension of such controlled studies to phenomena in high-dimensional disordered systems, such as the three-dimensional Anderson metal-insulator transition has been limited. Kicked quantum gases provide an alternate experimental platform that captures the Anderson model in momentum space and features dynamical localization as the analog of Anderson localization. Here, we utilize a momentum space lattice platform using quasiperiodically kicked ultracold atomic gases to experimentally investigate interaction effects on the three-dimensional dynamical Anderson metal-insulator transition. Here, we observe interaction-driven subdiffusion and a divergence of delocalization onset time on approaching the phase boundary. Mean-field numerical simulations show qualitative agreement with experimental observations, but with significant quantitative deviations.

74 ATOMIC AND MOLECULAR PHYSICS↗