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

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Effects of plasma squareness on RMP response in MAST-U

The effect of the squareness of the plasma boundary shape on the plasma response to the applied n = 2 (n is the toroidal mode number) resonant magnetic perturbation, and consequently on controlling the edge-localized mode (ELM) in MAST-U in the double-null configuration, is numerically systematically investigated. The modeling is performed utilizing three different plasma response models: the single-fluid model, the two-fluid model and the magnetohydrodynamic-kinetic hybrid model. The effect of the plasma squareness on ELM control is found to depend on the equilbrium kinetic profiles. Fixing the equilibrium radial profiles to that from a recent MAST-U discharge, the experimental plasma shape is found to be close to optimal for ELM control in terms of squareness. The optimal coil phasing does not significantly vary with varying plasma squareness. These findings robustly hold with (i) two figures of merit (FoMs) defined for predicting the ELM control feasibility, and (ii) three plasma response models employed for evaluating these FoMs.

MAST-U↗
DOE OSTI2026

Predicting fracture behavior in single crystal nickel using a coupled crystal plasticity phase field damage approach

Understanding the fracture behavior of single crystal metals is critical for predicting material performance under mechanical loading. Here, in this study, we investigate the fracture characteristics of single crystal nickel tensile bars using a crystal plasticity coupled phase field damage (CP-PFD) model. Experimental tensile tests were conducted on 15 specimens spanning five crystallographic orientations and three thickness variants per orientation. The results revealed two distinct fracture modes: brittle fractures with 45-degree angled surfaces and ductile fractures characterized by significant necking. The CP-PFD model successfully replicated these fracture behaviors, demonstrating strong agreement with experimental observations. The model effectively predicted the strain at which necking and fracture occurred, as well as the orientation-dependent fracture mechanisms. By comparing experimental and simulated fracture surfaces, we establish the CP-PFD model as a robust tool for predicting single crystal behavior and damage evolution. This work provides insight into the microstructural dependence of fracture behavior and establishes a predictive framework for modeling orientation-dependent damage evolution in single-crystal nickel.

Crystal plasticity↗
DOE OSTI2022

Control of a Three-Phase Grid-Connected Voltage-Sourced Converter Using Long Short-Term Memory Networks

With the rise of inverter-based resources (IBRs) within the power system, the control of grid-connected converters (GCCs) has become pertinent due to the fact they interface IBRs to the grid. The conventional method of control for a GCC such as the voltage-sourced converter (VSC) is through a decoupled control loop in the synchronous reference frame. However, this model-based control method is sensitive to parameter changes causing deterioration in controller performance. Data-driven approaches such as machine learning can be utilized to design controllers that are capable of operating GCCs in various system conditions. This work explores a deep learning-based control method for a three-phase grid-connected VSC, specifically utilizing a long short-term memory (LSTM) network for robust control. Simulations of a conventional controlled VSC are conducted using Simulink to collect data for training the LSTM-based controller. The LSTM model is built and trained using the Keras and TensorFlow libraries in Python and tested in Simulink. The performance of the LSTM-based controller is evaluated under different case studies and compared to the conventional method of control. Simulation results demonstrate the effectiveness of this approach by outperforming the conventional controller and maintaining stability under different system parameter changes.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Developing a Consistent Travel-Time Framework for Comparing Three-Dimensional Velocity Models for Seismic Location Accuracy

Abstract Location algorithms have historically relied on simple, one-dimensional (1D) velocity models for fast seismic event locations. 1D models are generally used as travel-time lookup tables, one for each seismic phase, with travel-times pre-calculated for event distance and depth. These travel-time lookup tables are extremely fast to use and this fast computational speed makes them the preferred type of velocity model for operational needs. Higher-dimensional (i.e., three-dimensional—3D) seismic velocity models are becoming readily available and provide more accurate event locations over 1D models. The computational requirements of these 3D models tend to make their operational use prohibitive. Additionally, comparing location accuracy for 3D seismic velocity models tends to be problematic, as each model is determined using different ray-tracing algorithms. Attempting to use a different algorithm than the one used to develop a model usually results in poor travel-time prediction. We demonstrate and test a framework to create first-P and first-S 3D travel-time correction surfaces using an open-source framework ( PCalc + GeoTess , https://www.sandia.gov/salsa3d/software/geotess ) that easily stores 3D travel-time and uncertainty data. This framework produces fast travel-time and uncertainty predictions and overcomes the ray-tracing algorithm hurdle because the lookup tables can be generated using the exact ray-tracing algorithm that is preferred for a model.

3D velocity models↗
DOE OSTI2021

Microstructure-dependent rate theory model of defect segregation and phase stability in irradiated polycrystalline LiAlO 2

We report gamma lithium aluminate (LiAlO 2 ) is a breeder material for tritium and is one of key components in a tritium-producing burnable absorber rod (TPBAR). Dissolution and precipitation of second phases such as LiAl 5 O 8 and voids are observed in irradiated LiAlO 2 . Such microstructure changes cause the degradation of thermomechanical properties of LiAlO 2 and affect tritium retention and release kinetics, and hence, the TPBAR performance. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed for investigating the accumulation of species and phase stability in polycrystalline LiAlO 2 structures under irradiation. Three sublattices (i.e. [Li, Al, V] I [O, V o ] II [Li i , Al i , O i , V i ] III ), and concentrations of six diffusive species (i.e. Li; vacancy of Li or Al at [Li, Al, V] I sublattice, O vacancy at [O, V o ] II sublattice, and Li, Al and O interstitials at [Li i , Al i , O i , V i ] III interstitial sublattices; are used to describe spatial and temporal distributions of defects and chemistry. Microstructure-dependent thermodynamic and kinetic properties including the generation, reaction, and chemical potentials of defects and defect mobility are taken into account in the model. The parametric studies demonstrated the capability of the developed RIS model to assess the effect of thermodynamic and kinetic properties of defects on the segregation and depletion of species in polycrystalline structures and to explain the phase stability observed in irradiated LiAlO 2 samples. The developed RIS model will be extended to study the precipitation of LiAl 5 O 8 and voids and tritium retention by integrating the phase-field method.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetic states of the quasi-one-dimensional iron chalcogenide Ba 2 FeS 3

Quasi-one-dimensional iron-based ladders and chains, with the 3d iron electronic density n=6, are attracting considerable attention. Recently, a new iron chain system Ba 2 FeS 3 , also with n=6, was prepared under high-pressure and high-temperature conditions. Here the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered ↑-↓-↑-↓ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard U coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large U/W (W=bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba 2 FeS 3 provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material.

1-dimensional systems↗
DOE OSTI2026

Trisodium Phosphate Phases and Solubility in Alkaline Solutions Relevant to Radioactive Waste Processing

The U.S. Department of Energy’s Hanford Site faces significant challenges in managing millions of gallons of legacy radioactive waste, where phosphate precipitation can obstruct pipelines during retrieval and processing. To resolve long-standing inconsistencies in reported solubility and clarify factors governing trisodium phosphate hydrates, we examined the Na3PO4:NaOH:H2O system. Powder and single-crystal X-ray diffraction revealed that commercial precursors undergo transformations that produce multiple hydrates, including three previously unreported phases comprising an ordered polymorph of Na3(PO4)·12H2O·1/6NaOH, Na3PO4·5H2O, and Na3PO4·9H2O. Computational modeling indicated that the ordered dodecahydrate is more stable than its disordered counterpart, suggesting kinetic persistence of structural disorder. Solubility measurements were conducted with solutions prepared either from anhydrous Na3PO4 or from Na3(PO4)·12H2O·1/6NaOH and revealed that release of interstitial NaOH from the hydrate precursor elevated solution alkalinity, thereby significantly reducing phosphate solubility relative to solutions prepared with anhydrous Na3PO4 across 20–44 °C. These results begin to reconcile inconsistencies in prior solubility data and clarify how phase composition dictates phosphate precipitation under alkaline conditions.

Graham, Trenton R. (ORCID:0000000189078004)↗
DOE OSTI2021

Improving the MJO Forecast of S2S Operation Models by Correcting Their Biases in Linear Dynamics

The operational dynamic subseasonal to seasonal (S2S) models for Madden-Julian oscillation (MJO) forecasting mostly still suffer from systematic errors in capturing the MJO's key dynamic features, such as its growth rate and propagation speed. By deriving the linear dynamic operators using the linear inverse modeling (LIM) approach, we propose a method to partly correct the errors in MJO linear dynamic operators to improve the MJO predictions of three operational dynamic S2S models. Correcting the deficiencies of the too-fast decay rates and the unrealistic propagating phase speeds lead to MJO prediction skills being extended by approximately 2–4 days. The improvements are more significant for the models with larger biases in MJO amplitude and propagation. This approach in principle may be extendable to predictions of other types of climate variability such as ENSO on one hand, and possible inclusions of nonlinear dynamics effects on the other hand.

58 GEOSCIENCES↗
DOE OSTI2026

Phase-field approach to cellular blebbing

Bulges in the plasma membrane of cells known as blebs can form spontaneously in a wide range of biological processes, but what controls their shape and stability remains incompletely understood. Here, to address this we introduce a dual phase-field model with coupled order parameters representing the cell cortex and plasma membrane that can quantitatively model blebbing in three dimensions. Simulations and sharp-interface analyses reveal that, depending on whether blebbing occurs by detachment of the plasma membrane or rupture of the actin cortex, blebs can form discontinuously through a saddle-node bifurcation or continuously with increasing cortical tension. The model predictions are in good quantitative agreement with existing experimental data for laser-induced cortex rupture.

Biological and medical sciences↗
DOE OSTI2021

Smart Inverter Twin Model for Anomaly Detection

Smart inverters connected to a communication network are vulnerable to various anomalies in the form of cyberattacks. In this paper, a self-security approach is implemented using the digital twin concept for smart inverters. The digital twin is formed using the inverter’s dynamic model. Then, the incoming setpoints are autonomously examined using the digital twin, and only the safe setpoints are engaged to the inverter’s local controller. This paper demonstrates the details of the self-security algorithm and how the inverter’s digital twin is formed. In particular, the stable and unstable operation region is experimentally verified by changing the power setpoints engaged to the local controller, using a laboratory setup including a three-phase 1.5-kVA SiC-MOSFET inverter and a 12-kW NHR 9410 regenerative power grid emulator. The results demonstrate that the digital twin model can potentially protect inverters from abnormal operation by examining the incoming commands (new setpoints) using the inverter’s digital twin before engaging the setpoints to the local controller.

Hossen, Tareq↗
DOE OSTI2025

DC-Link Capacitor Design for a Neutral-Point-Less Three-Level Dual-Phase Inverter for Traction Application

Conventional multilevel inverters, such as neutral￾point-clamped and T-type inverters, have gained popularity in electric vehicle applications due to their advantages, including high voltage range, high power capability, low switching losses, low total harmonic distortion, and low electromagnetic interference. However, these traditional multilevel inverters require a neutral point connection to generate a zero-voltage vector. The neutral point current oscillates at three times the fundamental frequency, leading to voltage imbalance and overvoltage stress on the power modules. Additionally, the use of two stacked DC-link capacitors increases the volume required for the same overall capacitance and complicates packaging due to ripple current and heat dissipation from separate components. This is a significant concern for traction drive units, where space is limited. In this paper, a neutral-point-less multilevel dual three￾phase inverter topology is investigated for traction inverter applications. Simulation results demonstrate that the proposed topology effectively retains the benefits of multilevel operation while utilizing a single DC-link capacitor. The inverter model was simulated in conjunction with an industry-standard battery model to evaluate the potential for capacitor size reduction compared to a conventional three-level inverter.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Erosion and cathodic arc plasma of Nb–Al cathodes: composite versus intermetallic

Many properties of cathodic arcs from single-element cathodes show a correlation to the cohesive energy of the cathode material. For example, the burning voltage, the erosion rate, or, to a lesser extent, plasma properties like electron temperatures or average ion energy and charge states. For multi-element cathodes, various phases with different cohesive energies can initially be present in the cathode, or form due to arc exposure, complicating the evaluation of such correlations. To test the influence of morphology and phase composition of multi-element cathodes on cathodic arc properties, a Nb–Al cathode model system was used that includes: pure Nb and Al cathodes; intermetallic Nb 3 Al, Nb 2 Al and NbAl 3 cathodes; and three composite Nb–Al cathodes with atomic ratios corresponding to the stoichiometric ratios of the intermetallic phases. Pulsed cathodic arc plasmas from these cathodes were examined using a mass-per-charge and energy-per-charge analyzer, showing that charge-state-resolved ion energy distributions of plasmas from the intermetallic and corresponding composite cathodes are nearly identical. An examination of converted layers of eroded cathodes using x-ray diffraction and scanning electron microscopy indicates the formation of a surface layer with similar phase composition for intermetallic and their corresponding composite cathode types. The average arc voltages do not follow the trend of cohesive energies of Nb, Al and intermetallic Nb–Al phases, which have been calculated using density functional theory. Possible reasons for this effect are discussed based on the current knowledge of multi-element arc cathodes and their arc plasma available in literature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nudging allows direct evaluation of coupled climate models with in situ observations: a case study from the MOSAiC expedition

Comparing the output of general circulation models to observations is essential for assessing and improving the quality of models. While numerical weather prediction models are routinely assessed against a large array of observations, comparing climate models and observations usually requires long time series to build robust statistics. Here, we show that by nudging the large-scale atmospheric circulation in coupled climate models, model output can be compared to local observations for individual days. We illustrate this for three climate models during a period in April 2020 when a warm air intrusion reached the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition in the central Arctic. Radiosondes, cloud remote sensing and surface flux observations from the MOSAiC expedition serve as reference observations. The climate models AWI-CM1/ECHAM and AWI-CM3/IFS miss the diurnal cycle of surface temperature in spring, likely because both models assume the snowpack on ice to have a uniform temperature. CAM6, a model that uses three layers to represent snow temperature, represents the diurnal cycle more realistically. During a cold and dry period with pervasive thin mixed-phase clouds, AWI-CM1/ECHAM only produces partial cloud cover and overestimates downwelling shortwave radiation at the surface. AWI-CM3/IFS produces a closed cloud cover but misses cloud liquid water. Our results show that nudging the large-scale circulation to the observed state allows a meaningful comparison of climate model output even to short-term observational campaigns. We suggest that nudging can simplify and accelerate the pathway from observations to climate model improvements and substantially extends the range of observations suitable for model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Global quantum phase diagram and non-Abelian chiral spin liquid in a spin- 3 2 square-lattice antiferromagnet

Since strong quantum fluctuations are essential for the emergence of quantum spin liquids, there have been extensive exploration and identification of spin liquid candidates in spin-$\frac{1}{2}$ systems, while such activities are rare in higher spin systems. Here we report an example of non-Abelian chiral spin liquid emerging in a spin-$\frac{3}{2}$ Heisenberg model on a square lattice. By tuning Heisenberg exchange interaction and scalar chirality interaction, we map out a quantum phase diagram enclosing three conventional magnetic orders and a chiral spin liquid based on density-matrix renormalization group studies. The nature of the spin liquid is identified as a long-sought bosonic version of the Read-Rezayi state that supports non-Abelian Fibonacci anyonic statistics, identified by the ground state entanglement spectrum. Significantly, we establish that the non-Abelian chiral spin liquid emerges through the enlarged local degrees of freedom and enhanced quantum fluctuations near the classical phase boundaries of competing magnetic orders. Finally, our numerical discovery of an exotic quantum spin liquid in a spin-$\frac{3}{2}$ system suggests a route for discovering fractionalized quantum phases in frustrated higher spin magnetic compounds.

2-dimensional systems↗
DOE OSTI2023

Phase coexistence implications of violating Newton’s third law

Newton’s third law, action = reaction, is a foundational statement of classical mechanics. However, in natural and living systems, this law appears to be routinely violated for constituents interacting in a nonequilibrium environment. Here, we use computer simulations to explore the macroscopic phase behavior implications of breaking microscopic interaction reciprocity for a simple model system. We consider a binary mixture of attractive particles and introduce a parameter that is a continuous measure of the degree to which interaction reciprocity is broken. In the reciprocal limit, the species are indistinguishable, and the system phase separates into domains with distinct densities and identical compositions. Furthermore, increasing nonreciprocity is found to drive the system to explore a rich assortment of phases, including phases with strong composition asymmetries and three-phase coexistence. Many of the states induced by these forces, including traveling crystals and liquids, have no equilibrium analogs. By mapping the complete phase diagram for this model system and characterizing these unique phases, our findings offer a concrete path forward toward understanding how nonreciprocity shapes the structures found in living systems and how this might be leveraged in the design of synthetic materials.

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

Ultrahigh-pressure crystallographic passage towards metallic hydrogen

The structural evolution of molecular hydrogen H 2 under multi-megabar compression and its relation to atomic metallic hydrogen is a key unsolved problem in condensed-matter physics. Although dozens of crystal structures have been proposed by theory, only one, the simple hexagonal-close-packed (hcp) structure of only spherical disordered H 2 , has been previously confirmed in experiments. Through advancing nano-focused synchrotron X-ray probes, here we report the observation of the transition from hcp H 2 to a post-hcp structure with a six-fold larger supercell at pressures above 212 GPa, indicating the change of spherical H 2 to various ordered configurations. Theoretical calculations based on our XRD results found a time-averaged structure model in the space group $P\bar{6}2c$ with alternating layers of spherically disordered H 2 and new graphene-like layers consisting of H 2 trimers (H 6 ) formed by the association of three H 2 molecules. Here, this supercell has not been reported by any previous theoretical study for the post-hcp phase, but is close to a number of theoretical models with mixed-layer structures. The evidence of a structural transition beyond hcp establishes the trend of H 2 molecular association towards polymerization at extreme pressures, giving clues about the nature of the molecular-to-atomic transition of metallic hydrogen. Considering the spectroscopic behaviours that show strong vibrational and bending peaks of H 2 up to 400 GPa, it would be prudent to speculate the continuation of hydrogen molecular polymerization up to its metallization.

Ji, Cheng [Center for High Pressure Science and Te↗