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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Comparison of excess free energy at an interface according to the applied interpolation scheme for elasticity: A phase-field method

Phase-field modeling is an effective simulation technique for modeling microstructure evolution of elastically anisotropic systems. To introduce the elastic energy contribution in a phase field model, an interpolation scheme is used to define the mechanical properties within the phases and across the continuous interface. Several existing interpolation schemes introduce a potential excess elastic energy at the interface, which undesirable effect on microstructure evolution needs to be evaluated. In this study, we focused on three interpolation schemes including Khachaturyan’ scheme (KHS), Voigt–Taylor’s scheme (VTS), and Steinbach–Apel’s scheme (SAS). Comparisons of these schemes’ performances were performed in three configuration types using the MOOSE (Multiphysics Object-Oriented Simulation Environment) framework: bi-crystal, isotropic particle-matrix and anisotropic particle-matrix. The contribution of excess elastic energy on the interface energy as a function of interface width and the computational time to steady-state were evaluated in these three configurations. SAS introduces the lowest excess elastic energy contribution and the VTS has the biggest contribution amongst the considered schemes. Moreover, when modeling precipitation in an anisotropic elastic material, the SAS approach seems to predict more physical convex shapes during growth, making it preferable to KHS and VTS. Finally, as currently implemented, SAS requires the largest computational time and KHS requires the smallest time to reach steady-state amongst the considered schemes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Physics-Informed Gaussian Process Regression for States Estimation and Forecasting in Power Grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method. We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power. For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Kinetic Pathways Templated by Low-Temperature Intermediates during Solid-State Synthesis of Layered Oxides

Layered oxides have been the dominant cathodes in lithium-ion batteries, and among them, high-nickel (Ni) systems are attractive because of their high capacity. For practical use, synthetic control of stoichiometry and structural ordering is crucial but has been nontrivial due to the complexity inherent to synthesis reactions, which often proceed via nonequilibrium pathways. In this paper, we report here a combined in situ synchrotron X-ray diffraction and ab initio study of solid-state synthesis of layered oxides starting from acetate precursors for LiCoO 2 and LiNiO 2 and their solid solution LiNi 0.8 Co 0.2 O 2 . While all three systems ultimately evolve into the same thermodynamically stable layered phase ($R\overline{3}m$), each chemistry involves distinct metastable intermediates. We explain the phase progressions using a structural template model, demonstrating that during the synthesis of LiCoO 2 , the formed metastable spinel polymorph (Li 2 Co 2 O 4 ; $Fd\overline{3}m$) is a kinetically facile lithiation product of spinel Co 3 O 4 —the low-temperature (LT) intermediate from the decomposed Co-acetate. Similarly, in the Ni-based systems, the acetate decomposition products, rocksalts (Ni,Co)O, topotactically template the kinetic pathways of forming disordered rocksalts (Li x (Ni,Co) 2–x O 2 ; $Fm\overline{3}m$), consequently leading to off-stoichiometric Li x (Ni,Co)O 2 with undesired high Li/Ni mixing. These findings highlight new opportunities for engineering precursors to form LT intermediates that template the synthesis of target phases and structural properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Signatures of QCD conductivities in heavy-ion collisions

Dissipative processes are pivotal for understanding the hydrodynamic evolution of hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. The interplay of multiple conserved charges—net baryon, strangeness, and electric charge—is of particular interest. Here, we simulate the longitudinal hydrodynamic evolution with the three diffusion currents in a hydrodynamic model with a lattice-QCD-based equation of state, NEOS -4 D , and estimate rapidity distributions including diffusive corrections to the phase-space distribution in the presence of multiple charges, which ensure charge conservation at particlization. We determine the response of particle yields at midrapidity to changes in the diagonal and off-diagonal conductivities. Inversely, we find that most components of the conductivity matrix can be constrained experimentally using identified particle multiplicities at different collision energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Deep Learning-enhanced Block-Diagram Modeling of Solar Power Systems

Data-driven models of power system inverter-based resources are desired to run simulations faster than with detailed electromagnetic transient models, to hide proprietary design details, to support control system design applications, and to aggregate the effects of distributed energy resources. This paper applies a customized Hammerstein Wiener framework to train block diagram models from thousands of electromagnetic transient simulations or experimental test records. The block diagram models integrate with larger grid simulations as voltagecontrolled current sources or current-controlled voltage sources for several simulators. Guidelines for block architecture and training are presented. Three-phase balanced, three-phase unbalanced, and single-phase examples all achieve an acceptable root mean square error of no more than 0.05 per-unit.

Mcdermott, Thomas E. [Private consulting company]↗
DOE OSTI2022

Novel Power-Hardware-in-the-Loop Interface Method for Grid-Forming Inverter Systems

Power-hardware-in-the-loop (PHIL) simulations of grid-forming (GFM) inverter systems facilitate the testing of drastic scenarios, such as on-grid to off-grid transitions and islanded microgrid operations without a stiff grid. To the authors’ best knowledge, most studies in the literature focus on PHIL simulations for grid-following inverter systems. Only a few studies focus on GFM inverters, and those are challenging and problematic, especially for high-power applications. This article proposes a novel PHIL simulation platform that enables interfacing high-power GFM inverter systems. The paper proposes the concept of a virtual GFM inverter as a part of the proposed PHIL interface. This addition of a virtual GFM inverter in the PHIL interface expands the conventional ideal transformer model (ITM) method and enables it to overcome the issues of instability of existing ITM methods. In the validation stage, a PHIL experiment is conducted on a three-phase, 480-V, 125-kVA GFM inverter system with the proposed interfacing method. The results corroborate that the proposed PHIL simulation method performs well and is stable for GFM inverter systems.

droop control↗
DOE OSTI2023

Virtual Flow Solver - Geophysics: A 3D Incompressible Navier-Stokes Solver

Virtual Flow Solver - Geophysics (VFS-Geophysics) is a three-dimensional (3D) incompressible Navier-Stokes solver based on the Curvilinear Immersed Boundary (CURVIB) method. The CURVIB is a sharp interface type of immersed boundary (IB) method that enables the simulation of fluid flows in the presence of geometrically complex moving bodies. The CURVIB method can be applied to wind/MHK turbine simulations and energy applications. VFS-Geophysics is the result of many years of research work by several graduate students, post-docs, and research associates that have been involved in the Computational Hydrodynamics and Biofluids Laboratory directed by Professor Fotis Sotiropoulos. The preparation of the present manual has been supported by the U.S. Department of Energy (DE-EE 0005482).

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Predictive Modeling and Uncertainty Quantification in Condition Monitoring of Active Components: A Reactor Coolant Pump Use Case

This work develops data-driven models for onset of thermal barrier leakage in reactor coolant pumps. It incorporates uncertainty quantification to enhance the reliability and robustness of pre- dictions. Using synthetic data generated by the Generic Pressurized Water Reactor simulator, realistic degradation scenarios were simulated across lifecycle stages—beginning, middle, and end of life. Key variables, including differential pressure, flow rate, vibration, and temperatures, were analyzed using machine learning framework. The fully connected neural network models demonstrated exceptional performance, achieving R2 scores exceeding 0.99 and root mean square errors as low as around 8.23 × 10-2 gallon per minute (gpm) for the three stages of the lifecy- cle. UQ analysis further validated the model’s robustness, with narrow uncertainty bounds during steady-state operations and appropriately wider bounds during transitional phases, reflecting the physical behavior of the system. This work addresses important gaps in real-time condition moni- toring and regulatory compliance by integrating advanced condition monitoring technologies with UQ into IST programs. The ability to detect thermal barrier leakage early and quantify prediction reliability supports optimizing maintenance strategies while ensuring nuclear power plants’ safe and reliable operation.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Modeling and Harmonic Stability of MMC-HVDC With Passive Circulating Current Filters

A modular multilevel converter (MMC) is an emerging converter technology widely used in wind power applications using high-voltage direct current (HVDC) transmission, and it has been largely investigated in medium-voltage solar harvesting as well as electric motor drives. Different from studies dealing with active circulating current control loops, this paper focuses on impedance modeling and harmonic stability studies for MMCs with passive circulating current filters (PCCFs). A power stage circuit including a PCCF is transformed and remodeled to obtain small-signal formulations. This transformation leads to additional high-order matrix computing complexity due to the added impedance subnetwork matrix from the PCCF as well as relevant crossed-frequency impacts. To simplify the computational burden, the proposed model aims to solve one arm equation instead of all six MMC arm equations. To achieve this result, additional challenges occur when the MMC is connected to a renewable energy source instead of the grid voltage, where low-level zero-sequence components might exist in the neutral point of three phases, especially in the case of no integrated advanced modular voltage balancing control, which prevents computational simplification. To overcome this challenge, further impedance matrix adjustments are conducted in this paper to theoretically suppress the impact of zero-sequence harmonics when obtaining small-signal impedance, taking into account the frequency coupling effect. Finally, simulations are carried out to validate the developed impedance model under different operating scenarios. Harmonic stability case studies of MMCs with PCCFs connected to renewable energy current sources are also presented, where frequency-domain stability analyses based on Bode diagrams are compared to time-domain simulation waveforms and FFTs, which validate the effectiveness of the proposed impedance model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Crack-Free 30% Chromium-Nickel Alloy Welding Products for Nuclear Service

We report prior research in the development of 30% chromium-nickel alloy nuclear welding wires has resulted in the resolution of primary water stress corrosion cracking (PWSCC) and ductility dip cracking (DDC) as well as improvement in solidification cracking (SC) resistance. The resolution of DDC exhibits some Laves phase, which has a negative effect on SC resistance. In this study, the use of an alternate carbide former, tantalum (Ta), in combination with niobium (Nb) was researched. Three heats of recently designed Filler Metal 52MSS-Ta (i.e., HV1648, HV1673A, and VX131WXW) were melted, fabricated, and systematically studied. DDC and SC were evaluated with thermodynamic modeling using the Scheil solidification simulation model, two types of varestraint tests, and strain-to-fracture (STF) testing. The varestraint and STF test results showed an improved SC resistance with reduced Laves phase and concurrent excellent DDC resistance. Optimized compositions with low Laves phase also exhibited high threshold strain values (TSVs) in the STF test. VX131WXW — which contains 2.81 wt-% Ta, 0.6 wt-% Nb, and 6 wt-% iron (Fe) - exhibited a TSV of 24%. Thermo-Calc computed the Laves phase to be 0.24% for VX131WXW compared to 0.06% in HV1673A. This difference in Laves phase resulted in the lower SC resistance of VX131WXW compared to HV1673A when measured with longitudinal varestraint testing. The maximum crack distance for HV1673A was about 0.6 mm while that of heat VX131WXW was about 1.0 mm. The typical diluted weld deposit made with VX131WXW was also resistant to PWSCC due to the chromium content exceeding 24%. These simultaneous results mark progress toward crack-free welds and provide direction for further optimization of Ta-containing filler metals.

30% chromium-nickel alloy↗
DOE OSTI2025

Novel Results Visualization for Dynamic PSA and New Modeling Features in EMRALD

The Event Modeling Risk Assessment Linked Diagram (EMRALD) tool, developed at the Idaho National Laboratory (INL), was designed to simplify the creation of dynamic models and support various research projects. One of the primary goals of EMRALD was to provide visual methods for modeling. EMRALD consists of two main components: a web-based user interface for model development and a solve engine for running model simulations. Over time, it has evolved to meet the diverse needs of its users. Initially, EMRALD's results were simple text outputs with final key state percentages and uncertainty bounds. However, because EMRALD utilizes a three-phase discrete event simulation and tracks the paths of each simulation run leading to a key state, there is significant potential to analyze large sets of path results data, including state paths, events, and timing. Visualizing this data meaningfully posed a challenge. To address this, a novel time-based Sankey diagram was developed. EMRALD exports results data in a format that can be opened in this Sankey viewer, allowing users to visualize paths, occurrences, events, and probability data for the entire simulation run in a single diagram. Moreover, when EMRALD was first created, there were limited tools capable of meeting its graphical requirements, many of which are no longer supported. In 2024, a new web-based interface was developed using modern graphing tools, enabling additional modeling features. This paper discusses the new dynamic PSA results visualization capability and the enhanced modeling tools available in EMRALD.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

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

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

First Images of the Protoplanetary Disk around PDS 201

Scattered light imaging has revealed nearly a dozen circumstellar disks around young Herbig Ae/Be stars—enabling studies of structures in the upper disk layers as potential signs of ongoing planet formation. We present the first images of the disk around the variable Herbig Ae star PDS 201 (V* V351 Ori) and an analysis of the images and spectral energy distribution through 3D Monte Carlo radiative transfer simulations and forward modeling. The disk is detected in three data sets with the Large Binocular Telescope Interferometer/Large Binocular Telescope mid-infrared camera at the Large Binocular Telescope, including direct observations in the Ks and L′ filters, and an L′ observation with the 360° vector apodizing phase plate coronagraph. The scattered light disk extends to a very large radius of ∼250 au, which places it among the largest of such disks. Exterior to the disk, we establish detection limits on substellar companions down to ∼5 M {sub Jup} at ≳1.″5 (≳500 au), assuming the Baraffe et al. models. The images show a radial gap extending to ∼0.″4 (∼140 au at a distance of 340 pc) that is also evident in the spectral energy distribution. The large gap is a possible signpost of multiple high-mass giant planets at orbital distances (∼60–100 au) that are unusually massive and widely separated compared to those of planet populations previously inferred from protoplanetary disk substructures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Improved pore network models to simulate single-phase flow in porous media by coupling with lattice Boltzmann method

In this paper, different pore network models to simulate single-phase flow in porous media are built and their accuracy is evaluated. In addition to the conventional pore network model (CPNM) which consists of regular pore bodies and throat bonds, three improved pore network models (IPNMs) are developed allowing to better describing the real pore and throat geometry. The first improved pore network model (IPNM1) replaces the regular throat bond with a throat bond showing the real throat cross section. The second improvement (IPNM2) uses a series of sub-throat bonds with varying cross sections to better describe the real throat geometry, which is firstly proposed in this paper. The third model (IPNM3) extracts the real pore-throat-pore geometry without simplification. The conductance of fluid flow through these more realistic throat bonds is calculated by the lattice Boltzmann method (LBM). The accuracy and computational efficiency of the different pore network models are evaluated taking the LBM simulation over the whole porous medium as reference solution. The global permeability and detailed pressure distributions in the pores for the different pore network models are validated. The results show that the accuracy of the pore network model increases from CPNM to IPNM3, but at the expense of increasing computational cost. This study suggests that IPNM3 can replace a whole-domain LBM simulation with similar accuracy but much lower computational cost. As a first-order approximation the newly proposed IPNM2 is suggested as good compromise between accuracy and computational cost.

42 ENGINEERING↗
DOE OSTI2022

Modeling of high intensity orbital angular momentum beams for laser–plasma interactions

In this work, we explore the field of high intensity orbital angular momentum (OAM) beams, their generation with spiral phase mirrors, and the theory behind modeling both ideal and realistic beams. We explore OAM beam asymmetries introduced by aberrations in the beam, manufacturing defects, and bandwidth. A full three-dimensional description of the paraxial Laguerre–Gaussian (LG) modes suitable for modeling beams down to f/2 focal geometries is derived. A perturbative approach to modeling asymmetric OAM beams is introduced showing that only three LG modes are sufficient to model a wide variety of OAM asymmetries. The models are compared with experimental results followed by a discussion on the future of high intensity OAM beams in plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗
DOE OSTI2020

Elliptic flow coalescence to identify the $f_0$(980) content

Here we use a simple coalescence model to generate $f_0$(980) particles for three configurations: a $s\bar{s}$ meson, a $u\bar{u}s\bar{s}$ tetraquark, and a $K^+K^–$ molecule. The phase-space information of the coalescing constituents is taken from a multiphase transport (AMPT) simulation of heavy-ion collisions. It is shown that the number of constituent quarks scaling of the elliptic flow anisotropy can be used to discern $s\bar{s}$ from $u\bar{u}s\bar{s}$ and $K^+K^–$ configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗