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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hydro-Code Implementation and Testing of a Kinetic Phase Transition Framework

In this report we describe the Kinetic Phase Transition (KPT) framework that has been worked out over the last 10 years (from around 2014) and the implementation of it into three different codes, the one-dimensional hydro- LASLO and the three-dimensional magneto-hydro- ALEGRA, Sandia codes, via subroutines in the LAMBDA Equations of State and constitutive models package, and Flag, an arbitrary Lagrangian-Eulerian multiphysics code developed within the Lagrangian Applications project (LAP) at LANL. We discuss the introduction of phase mass (and/or volume) fractions that are needed in a code for it to be ‘phase aware’, that is, not only the thermodynamic state is known in each point but also the mixture of the materials’ phases in that point. Further we point to the need of a full Equations of State for each phase in a material to achieve phase awareness and we review the equilibrium phase model, where a phase mixture is at its lowest Gibbs free energy state, to make this point clear. Contrasting the kinetic phase transition to this equilibrium model seamlessly introduce us to the KPT framework that is subsequently thoroughly discussed. While the determination of the total state and the states and mass fractions of phases in each point is a problem that can borrow many of its numerical details from Eulerian codes and mixture of materials (not phases), the update of mass fractions with time in a KPT framework needs a new set of considerations. General for any update model is that we need to prevent mass fractions from becoming unphysical (negative or their sum to be larger than one). We have solved this problem by implementing a subdivision of the hydro time step that prevents the phase from being fully present to not present at all in one subdivided time step by limiting the size of the subdivided time step. This scheme also corrects numerical problems from abrupt changes in parameter values, the so called Gibbs phenomena, that gives rise to slushing between phases in the KPT framework. Interspersed throughout the report are discussions on different thermodynamics considerations. EOS validity windows, limitations on the EOS phase space, are needed for the KPT framework and are discussed separately and exemplified. The KPT framework described in this report has been verified by code comparison, but validation is still an active area of research. There is room for improvement in the update model, both in the model for determination of rates and in how to prevent the mass fractions from becoming unphysical. In addition, the parameters in the KPT update model and the placement of the phase boundary in the EOS phase space, and interactions with other constitutive models, are closely related and interfering with each other. One possible way forward is to simultaneously develop KPT parameters, EOS, and constitutive models for each material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Three-dimensional microstructural characterization of FBR MOX fuel and the contribution of microstructural features to the thermal conductivity of the fuel

Combination of microstructural characterization, property measurements, and phase field modeling is used to investigate fast breeder reactor (FBR) mixed oxide (MOX) fuel irradiated to burnup of 13.7% fissions per initial metal atom (FIMA). Here, the fuel was characterized at different radial locations, which revealed that grey phase can be present in the central region if it nucleates on five metal precipitates (FMPs). In addition, in the mid-radial region FMPs do not diffuse out of the region once formed and the size of Pd–Te precipitates is dictated by the porosity present in the region. Thermal conductivity measurements were conducted as a function of radial location and the microstructure of the fuel was correlated with the observed trend. Reconstructions of the 3D solid and gaseous fission product structures in different regions of the fuel were used to simulate the effective thermal conductivity (ETC) of the respective regions and determine which microstructural feature has the strongest impact on thermal conductivity. Based on conducted assessment, FMPs and Pd-Te precipitates improve local conductivity of central and mid-radial regions even in the presence of grey phase, but defects are primary contributor to the degradation of thermal conductivity on the periphery of the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Block Mott insulating state induced by next-nearest-neighbor hopping in the S = 3 2 zigzag chain BaCoTe 2 O 7

Quasi-one-dimensional correlated electronic multiorbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-field splitting, the electronic correlations (Hubbard repulsion U and Hund coupling J H ), and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe 2 ⁢ O 7 , with electronic density n=7, was prepared experimentally. Here, we systematically study the electronic and magnetic properties of this quasi-one-dimensional compound from the theoretical perspective. Based on first-principles density functional theory calculations, strongly anisotropic one-dimensional electronic Co 3⁢d bands were found near the Fermi level. By evaluating the relevant hopping amplitudes, we provide the magnitude and origin of the nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping matrices in BaCoTe 2 ⁢ O 7 . With this information, we constructed a three-orbital electronic Hubbard model for this zigzag chain system, and studied two cases: with only a NN hopping matrix, and with NN plus NNN hopping matrices. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group method, we constructed the ground-state phase diagram. A robust staggered ↑ - ↓ - ↑ - ↓ antiferromagnetic (AFM) region was found when only the NN hopping matrix in the chain direction was employed. However, for the realistic case where the NNN hopping matrix is also included, the dominant state becomes instead a block AFM ↑ - ↑ - ↓ - ↓ order, in agreement with experiments. The system displays Mott insulator characteristics with three half-filled orbitals, when the block AFM order is stable. In conclusion, our results for BaCoTe 2 ⁢O 7 provide guidance to experimentalists and theorists working on this zigzag one-dimensional chain and related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

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

OC6 Phase III Definition Document

The objective of Phase III of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project was to evaluate the accuracy of aerodynamic load predictions by offshore wind modeling tools for a floating offshore wind turbine as it experiences large surge-translational and pitch-rotational motion, as would occur during normal operation. A variety of models were examined in the project, with the goal of performing a three-way validation between engineering-level tools, higher-fidelity tools, and measurement data from two wind tunnel experimental campaigns. This document provides the necessary information to model and simulate the system studied in the OC6 Phase III project. This information is not only useful for project participants, but also future modelers who want to replicate this work for verification and validation purposes.

17 WIND ENERGY↗
DOE OSTI2023

Experimental and modeling studies of IPDI-based polyurea elastomers – The role of hard segment fraction

Segmented polyureas (PUa) are industrially important class of polymers widely used in coatings, sealant, and adhesive applications. Here, in this paper, we report synthesis, characterization, and modeling of Isophorone Diisocyanate-Diethyl-Toluene-Diamine-Polyether amine (IPDI-DETDA-PO PUa) with varied hard segment contents of 20, 30, and 40 weight percent. For each of the three materials, we study its structure and phase behavior using FTIR, DSC, and TEM, and clearly show the presence of microphase separation between the hard and soft nanodomains. We then measure the linear viscoelastic response of the PUa-s using DMA (frequency sweeps at multiple temperatures). The DMA data are shown to obey the time-temperature superposition. Finally, we develop a new micromechanical model describing the DMA results; the model describes a phase-separated PUa as two “Fractional-order Maxwell gels” branches, connected in parallel, with the first FMG branch representing the “percolated hard phase" and the second one modeling the “filled soft phase". In agreement with the earlier thermodynamic theories, the volume-fraction of the percolated hard phase is related to the hard segment weight-fraction (HSWF), defined as the combined mass of IPDI and DETDA normalized to the total mass of the polymer. The data and model are found to be in a good qualitative and quantitative agreement.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An Optimal Power Control Strategy for Grid-Following Inverters in a Synchronous Frame

This work proposes a power control strategy based on the linear quadratic regulator with optimal reference tracking (LQR-ORT) for a three-phase inverter-based generator (IBG) using an LCL filter. The use of an LQR-ORT controller increases robustness margins and reduces the quadratic value of the power error and control inputs during transient response. A model in a synchronous reference frame that integrates power sharing and voltage–current (V–I) dynamics is also proposed. This model allows for analyzing closed-loop eigenvalue location and robustness margins. The proposed controller was compared against a classical droop approach using proportional-resonant controllers for the inner loops. Mathematical analysis and hardware-in-the-loop (HIL) experiments under variations in the LCL filter components demonstrate fulfillment of robustness and performance bounds of the LQR-ORT controller. Experimental results demonstrate accuracy of the proposed model and the effectiveness of the LQR-ORT controller in improving transient response, robustness, and power decoupling.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Prediction of exotic magnetic states in the alkali-metal quasi-one-dimensional iron selenide compound Na 2 FeSe 2

The magnetic and electronic phase diagram of a model for the quasi-one-dimensional alkali-metal iron selenide compound Na 2 FeSe 2 is presented. The novelty of this material is that the valence of iron is Fe 2+ , contrary to most other iron-chain compounds with valence Fe 3+ . Using first-principles techniques, we developed a three-orbital tight-binding model that reproduces the ab initio band structure near the Fermi level. Including Hubbard and Hund couplings and studying the model via the density-matrix renormalization group and Lanczos methods, we constructed the ground-state phase diagram. A robust region where the block state ↑↑↓↓↑↑↓↓ is stabilized was unveiled. The analog state in iron ladders, employing 2×2 ferromagnetic blocks, is by now well established, but in chains a block magnetic order has not been observed yet in real materials. The phase diagram also contains a large region of canonical staggered spin order ↑↓↑↓↑↓↑ at very large Hubbard repulsion. At the block-to-staggered transition region, an exotic phase is stabilized with a mixture of both states: an inhomogeneous orbital-selective charge density wave with the exotic spin configuration ↑↑↓↑↓↓↑↓. Our predictions for Na 2 FeSe 2 may guide crystal growers and neutron-scattering experimentalists towards the realization of block states in one-dimensional iron selenide chain materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Emergent quasi-one-dimensional antiferromagnetism in the distorted kagome magnet CePtPb

CePtPb hosts a distorted kagome lattice of Ce 3+ ions, providing a clean platform to investigate how reduced local symmetry and strong spin-orbit coupling reshape frustrated magnetism. Magnetization, specific heat, and magnetocaloric effect measurements, combined with a symmetry analysis of the single-ion anisotropy, demonstrate that the local 𝑚⁢2⁢𝑚 site symmetry selects a nearly Ising-like Kramers doublet with easy axes lying within the 𝑎⁢𝑏 plane. This results in three distinct in-plane Ising directions and an overall easy-plane anisotropy. The low-energy magnetic response is well captured by a three-sublattice Ising model, which quantitatively reproduces the saturation magnetization for arbitrary in-plane field orientations, including $[110]$ and $[1\bar{⁢1}⁢0]$, as well as the ratio of the field-induced critical fields. For 𝐵∥$[110]$, the phase diagram exhibits two quantum critical points at 𝐵 c⁢1 = 0.25T and 𝐵 c⁢2 = 0.55T, arising from the sequential polarization of the three Ising sublattices. In conclusion, these results reveal that the system develops quasi-one-dimensional spin chains along the 𝑐 axis, emerging from the nominally three-dimensional crystal structure composed of stacked kagome layers, and illustrate how reduced local symmetry can drive effective dimensional reduction in rare-earth Ising magnets.

Li, Fangli [Southern University of Science and Tec↗
DOE OSTI2025

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Discovery and Timing of Three Millisecond Pulsars in Radio and Gamma-Rays with the Giant Metrewave Radio Telescope and Fermi Large Area Telescope

In this work, we performed deep observations to search for radio pulsations in the directions of 375 unassociated Fermi Large Area Telescope γ-ray sources using the Giant Metrewave Radio Telescope (GMRT) at 322 and 607 MHz. In this paper we report the discovery of three millisecond pulsars (MSPs), PSR J0248+4230, PSR J1207–5050, and PSR J1536–4948. We conducted follow-up timing observations for ~5 yr with the GMRT and derived phase-coherent timing models for these MSPs. PSR J0248+4230 and J1207–5050 are isolated MSPs having periodicities of 2.60 ms and 4.84 ms. PSR J1536–4948 is a 3.07 ms pulsar in a binary system with an orbital period of ~62 days about a companion of a minimum mass of 0.32 M⊙. We also present multifrequency pulse profiles of these MSPs from the GMRT observations. PSR J1536–4948 is an MSP with an extremely wide pulse profile having multiple components. Using the radio timing ephemeris we subsequently detected γ-ray pulsations from these three MSPs, confirming them as the sources powering the γ-ray emission. For PSR J1536–4948 we performed combined radio–γ-ray timing using ~11.6 yr of γ-ray pulse times of arrival (TOAs) along with the radio TOAs. PSR J1536–4948 also shows evidence for pulsed γ-ray emission out to above 25 GeV, confirming earlier associations of this MSP with a ≥10 GeV point source. The multiwavelength pulse profiles of all three MSPs offer challenges to models of radio and γ-ray emission in pulsar magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Probing CP and flavor violation in neutral kaon decays with ALPs

We analyze the three-body decays of the long-lived neutral kaon K L → ππa, where a is an axion-like particle (ALP), and compare them to the two-body decay K L → π 0 a. While the latter requires both flavor violation (FV) and CP violation (CPV), the former can proceed via FV alone, allowing the ratio of decay rates to serve as a probe of CPV of the underlying UV theory. We emphasize the importance of weak-interaction-induced contributions, often neglected in recent calculations. We explore both minimal and non-minimal flavor-violating scenarios, and identify classes of models where ALP production from neutral three-body decays is comparable to — or even dominates over — the two-body decay, despite its reduced phase space. Finally, we discuss the phenomenological implications of our results and show how these decays can provide complementary probes of ALP couplings beyond those accessible via charged kaon channels.

Axions and ALPs↗
DOE OSTI2021

Robust hierarchical dispatch for residential distribution network management considering home thermal flexibility and model predictive control

In the transactive energy (TE) paradigm, the devices of participative consumers, or prosumers, may be aggregated and employed to drive operational objectives at the network level. Home heating, ventilation and air-conditioning (HVAC) systems in particular are well-suited to modulate their behaviours based on both home thermal flexibility and requests from the utility grid. This paper develops a robust, hierarchical power dispatch scheme in the context of a residential distribution network. The formulation couples a unique, multiphase linear distribution optimal power flow (OPF) at the upper level with model predictive control (MPC)-based HVAC fleet controllers at the lower level. The proposed approach is tested on nearly 2000 homes with a three-phase distribution network in an intraday market setting, where two major applications are explored and analysed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Properties of first-order hadron-quark phase transition from inverting neutron star observables

By inverting the observational data of several neutron star observables in the three-dimensional parameter space of the constant speed of sound (CSS) model while fixing all hadronic equation-of-state parameters at their currently known most probable values, we constrain the three parameters of the CSS model and their correlations. Using two lower radius limits of $R$ $2.01$ = $11.41$ km and $R$ $2.01$ = $12.2$ km for PSR J0740+6620 obtained from two independent analyses using different approaches by the Neutron Star Interior Composition Explorer (NICER) Collaboration, the speed of sound squared $c$ $^2_{QM}$ in quark matter is found to have a lower limit of 0.35 or 0.43 in units of $c$ $2$ , above its conformal limit of $c$ $^2_{QM}$ < $1/3$. Furthermore, an approximately linear correlation between the first-order hadron-quark transition density $ρ$ $t$ and its strength $Δε$ is found. Moreover, the presence of twin stars is deemed improbable by the present work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Dynamic Restoration Strategy for Distribution System Resilience Enhancement

In electric power distribution systems, distributed energy resources (DERs) can act as controllable power sources and support utility operators to minimize power outages after extreme weather events (e.g., hurricane, earthquake, wildfire) and thus help enhance the grid's resilience. Meanwhile, the influences of extreme events and the capabilities of DERs are dynamic and difficult to predict. Hence, the desired distribution system restoration strategy should be able to evolve according to real-time fault/dis-turbance information and the availability of DERs. In this paper, we propose a new dynamic restoration strategy for distribution systems to enhance system resilience against potential hazards. An efficient reconfiguration algorithm is developed to eliminate the use of integer variables to relieve the computational burden. Model predictive control is implemented to adjust the system topology and DER operation set points based on the updated fault information and DER forecasts. The effectiveness of the proposed restoration model in enhancing distribution system resilience is validated through an IEEE 123-bus test system. Simulation results also validate that the proposed restoration model can mitigate the occurrence of unexpected events and the fluctuations of DERs.

distributed energy resources (DERs)↗
DOE OSTI2024

Convolution Neural Network for Fault Identification in Distribution Feeder with High Penetration Solar PV

Identification and zonal classification of the faults is a decisive factor in the relay’s decision to trip or not. Different types of fault like three-phase, line-to-line-to-ground and single-line-to-ground can occur at various locations in the feeder. These faults are seen as the variation in the instantaneous values of three-phase voltages and currents, i.e., waveforms, that are measured at the relay location. The objective of this work is to develop a machine learning model that can identify a fault and classify it to various protection zones based on measured waveforms. In this work, a data-driven relay based on Convolutional Neural Network (CNN) is proposed for fault identification in distribution feeders with high penetration solar PV. The proposed CNN model takes local current and voltage waveforms as input and classify it into fault, no-fault or a capacitor switching. Further, the CNN also attempts to identify fault zones based on the images of waveforms. The overall testing accuracy of the trained model exceeds 95%.

Ramesh, Meghana↗
DOE OSTI2025

Development of Innovative Non-Concentrated Alloy (NCA) Fuel Cladding for Advanced Nuclear Reactors

As the nuclear energy sector advances toward next-generation reactors, the need for high-performance fuel cladding materials has become increasingly urgent. Traditional alloys like zirconium and stainless steel are reaching their performance limits under higher temperatures, more corrosive coolants, and extended irradiation. This report presents the development of a new class of fuel cladding materials based on Non-Concentrated Alloys (NCAs) ? multi-element systems designed to deliver enhanced mechanical strength, corrosion resistance, and radiation tolerance. Through a combination of computational modeling (Computer Coupling of Phase Diagrams and Thermochemistry: CALPHAD), simulation-guided alloy selection, and experimental fabrication via arc melting and spark plasma sintering (SPS), three strategic alloy design paths were explored: (1) FeCrAl-based NCAs, (2) refractory-lean neutron-efficient alloys, and (3) equimolar high-entropy compositions. Microstructural analysis confirmed the formation of stable body-centered cubic BCC_A2 phases, while mechanical testing demonstrated hardness values significantly exceeding those of conventional cladding materials. The results highlight the tunability of NCA systems and their potential for balancing strength and ductility ? a critical consideration for in-reactor performance. Looking forward, future work will focus on thermomechanical optimization, CALPHAD refinement, and benchmarking against industry standards to enable scalable deployment. This work not only advances the science of nuclear materials but also supports broader goals in nuclear safety, performance, and nuclear energy innovation.

36 - MATERIALS SCIENCE↗