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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Calculating n-point charge correlations in evolving systems

Background: In dynamic systems such as heavy-ion collisions, charge susceptibilities and local charge correlations change with time. These changes are accompanied by nonlocal correlations, which spread diffusively with time and are constrained by local charge conservation. Such correlations have been measured at the Relativistic Heavy-Ion Collider (RHIC) and have been linked to the chemical evolution of the matter and possible phase separation or critical phenomena. The evolution of such correlations has been modeled for two-point correlators superimposed onto hydrodynamics, but not for three- or four-point correlators. These higher-order correlators represent an essential basis for calculating the kurtosis and skewness of charge distributions measured at RHIC. For short-range correlations, such as those between charges on the same particle, and the associated charge-balance correlations a theoretical formalism was lacking for handling higher-order correlations. In particular, one needed to understand how higher-order charge fluctuations was split onto lower numbers of particles. For example, correlations of order three have contributions with charge all on one particle, split onto two particle, or split onto three. By assuming local chemical equilibrium of short-range correlations, the sources of the various correlators were found to be uniquely determined. The relevant theoretical foundation is presented here, including a diagrammatic technique to correctly account for all the possible terms to higher-order correlations. Here, we found a consistent theoretical treatment and its applicability, viability, and tractability are discussed. The formalism derived here enables realistic and quantitative modeling of the three- and four-point charge correlations necessary for understanding measurements of the kurtosis and skewness of charge distributions at RHIC. These observables have been promoted as signals of phase transitions or measures of chemical evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Three-state Potts nematic order in stacked frustrated spin models with SO(3) symmetry

Here, we propose stacked two-dimensional lattice designs of frustrated and SO(3) symmetric spin models consisting of antiferromagnetic (AF) triangular and ferromagnetic (FM) sixfold symmetric sublattices that realize emergent $\mathbb{Z}_{3}$ Potts nematic order. Considering bilinear-biquadratic spin interactions, our models describe an SO(3)-symmetric triangular lattice AF subject to a fluctuating magnetization arising from the FM coupled sublattice. We focus on the classical AFM-FM windmill model and map out the zero- and finite-temperature phase diagram using Monte Carlo simulations and analytical calculations. We discover a state with composite Potts nematic order above the ferrimagnetic three-sublattice up-up-down ground state and relate it to Potts phases in SO(3)-broken Heisenberg and Ising AFMs in external magnetic fields. Finally, we show that the biquadratic exchange in our model is automatically induced by thermal and quantum fluctuations in the purely bilinear Heisenberg model, easing the requirements for realizing these lattice designs experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2006

Scalar Dissipation Modeling for Passive and Active Scalars: a priori Study Using Direct Numerical Simulation

Transitional databases from Direct Numerical Simulation (DNS) of three-dimensional mixing layers for single-phase flows and two-phase flows with evaporation are analyzed and used to examine the typical hypothesis that the scalar dissipation Probability Distribution Function (PDF) may be modeled as a Gaussian. The databases encompass a single-component fuel and four multicomponent fuels, two initial Reynolds numbers (Re), two mass loadings for two-phase flows and two free-stream gas temperatures. Using the DNS calculated moments of the scalar-dissipation PDF, it is shown, consistent with existing experimental information on single-phase flows, that the Gaussian is a modest approximation of the DNS-extracted PDF, particularly poor in the range of the high scalar-dissipation values, which are significant for turbulent reaction rate modeling in non-premixed flows using flamelet models. With the same DNS calculated moments of the scalar-dissipation PDF and making a change of variables, a model of this PDF is proposed in the form of the (beta)-PDF which is shown to approximate much better the DNS-extracted PDF, particularly in the regime of the high scalar-dissipation values. Several types of statistical measures are calculated over the ensemble of the fourteen databases. For each statistical measure, the proposed (beta)-PDF model is shown to be much superior to the Gaussian in approximating the DNS-extracted PDF. Additionally, the agreement between the DNS-extracted PDF and the (beta)-PDF even improves when the comparison is performed for higher initial Re layers, whereas the comparison with the Gaussian is independent of the initial Re values. For two-phase flows, the comparison between the DNS-extracted PDF and the (beta)-PDF also improves with increasing free-stream gas temperature and mass loading. The higher fidelity approximation of the DNS-extracted PDF by the (beta)-PDF with increasing Re, gas temperature and mass loading bodes well for turbulent reaction rate modeling.

scalar dissipation modeling↗

DOE OSTI2021

Enhancement of Distribution System State Estimation Using Pruned Physics-Aware Neural Networks

Realizing complete observability in the three-phase distribution system remains a challenge that hinders the implementation of classic state estimation algorithms. In this paper, a new method, called the pruned physics-aware neural network (P2N2), is developed to improve the voltage estimation accuracy in the distribution system. The method relies on the physical grid topology, which is used to design the connections between different hidden layers of a neural network model. To verify the proposed method, a numerical simulation based on one- year smart meter data of load consumptions for three-phase power flow is developed to generate the measurement and voltage state data. The IEEE 123-node system is selected as the test network to benchmark the proposed algorithm against the classic weighted least squares (WLS). Numerical results show that P2N2 outperforms WLS in terms of data redundancy and estimation accuracy.

distribution system state estimation↗

DOE OSTI2022

Development and Validation of Resistance-Capacitance Model (RCM) for Phase Change Material (PCM) Embedded in 3D Periodic Structures

The low thermal conductivity of Phase Change Materials (PCM) can be improved with extended surfaces such as additively manufactured 3D periodic lattice structures. Three different aluminum alloy-based lattices (base sizes 10, 20, and 40 mm) with average porosity of 0.95 filled with paraffin wax, with a nominal phase change temperature of 55°C, were experimentally investigated. In this work, a computationally efficient 2D Resistance Capacitance-based model (RCM) was developed for predicting the thermal characteristics of these geometries. Non-uniform porosity in the PCM-metal domain was estimated using image processing and served as model input. The solver does not solve for higher-order physics as in CFD but still can provide a good prediction of thermal resistance and energy storage at a very low computational cost. The simulation-to-real-time factor for this geometry is of the order of 10-4, while CFD simulations typically have a real-time factor greater than 1. The model was validated against the experimental data for melting under three different heat fluxes (6250 W/m2, 12500 W/m2, and 18750 W/m2). The mean deviation of the predicted average PCM temperature was between 1.34 K-2.81 K for different cases. The maximum average temperature deviation of 5.45 K was observed for the 20 mm geometry at the highest heat flux test condition. The effects of natural convection were neglected in the model, but the predicted PCM temperature and energy storage still showed good agreement with the experimental data.

25 ENERGY STORAGE↗

DOE OSTI2021

Evaluating the grain-scale deformation behavior of a single-phase FCC high entropy alloy using synchrotron high energy diffraction microscopy

We report although the deformation behavior of high-entropy alloys (HEAs) has been extensively studied at the macroscale, many important properties have yet to be explored for these alloys at the microscale, thus hampering accurate prediction of damage and failure processes. Synchrotron high-energy diffraction microscopy (HEDM) and fast-Fourier transform-based crystal plasticity modeling was conducted to investigate the three-dimensional (3D) grain-resolved micromechanical response for approximately 1,900 constituent grains within a single-phase FCC HEA up to 1% applied strain. The evolution of grain-resolved elastic strains, lattice reorientations, and maximum resolved shear stresses (mRSS) were evaluated to quantify elastic, yield, and fully plastic behavior. Overall, the initial critical resolved shear stress (CRSS), determined via in situ HEDM and companion modeling, was found to be > 20% higher than estimated using the classical polycrystalline Taylor factor (M = 3.06). However, a descriptive parameter based on the average grain-resolved Taylor factor ($\overline{M}$) was found to show excellent agreement with plastic yielding behavior observed within HEDM datasets. Noticeable deviations in HEDM lattice reorientations compared to both EVP-FFT simulations and classical predictions for FCC polycrystals were discovered, highlighting the complexity in correlating local lattice reorientations, Taylor, and Schmid factors with plastic response for this material at the grain-scale. Therefore, it is anticipated that the overall trends and parameter identification of 3D grain-resolved properties in this study can serve as an important foundation for continued mesoscale investigation on both well-established and newly developed Cantor-like HEAs.

36 MATERIALS SCIENCE↗

NASA NTRS2003

A Micro-Ultrastable Oscillator (micro-US0) for Micro/Nano Sciencecraft

Syntonics LLC developed a prototype micro-Ultra Stable Oscillator (micro-USO) under a Space Base Technology Grant (NAGS-10395). Syntonics conducted the micro-USO Program in two phases. In Phase I, we developed a set of verified analytical models (including thermal, electrical, and control models) for a baseline USO, conducted a series of six technology studies, and built three approx.9OOg prototype units. These prototypes provided a tool for evaluating competing design topologies. In Phase II we prepared the conceptual design of a approx.100-15Og micro-USO.

Source record↗

NASA NTRS1993

Fuel injector: Air swirl characterization aerothermal modeling, phase 2, volume 1

A well integrated experimental/analytical investigation was conducted to provide benchmark quality relevant to a prefilming type airblast fuel nozzle and its interaction with the combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM), and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

NASA NTRS1993

Fuel Injector: Air swirl characterization aerothermal modeling, phase 2, volume 2

A well integrated experimental/analytical investigation was conducted to provide benchmark quality data relevant to prefilming type airblast fuel nozzle and its interaction with combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) equipment was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM) and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

DOE OSTI2022

Prediction of orbital-selective Mott phases and block magnetic states in the quasi-one-dimensional iron chain Ce 2 O 2 FeSe 2 under hole and electron doping

The recent detailed study of quasi-one-dimensional iron-based ladders, with the 3d iron electronic density n=6, has unveiled surprises, such as orbital-selective phases. However, similar studies for n=6 iron chains are still rare. Here a three-orbital electronic Hubbard model was constructed to study the magnetic and electronic properties of the quasi-one-dimensional n=6 iron chain Ce 2 O 2 FeSe 2 , with focus on the effect of doping. Specifically, introducing the Hubbard U and Hund J H couplings and studying the model via the density matrix renormalization group, we report the ground-state phase diagram varying the electronic density away from n=6. For the realistic Hund coupling J H /U=1/4, several electronic phases were obtained, including a metal, orbital-selective Mott, and Mott insulating phases. Doping away from the parent phase, the competition of many tendencies leads to a variety of magnetic states, such as ferromagnetism, as well as several antiferromagnetic and magnetic “block” phases. In the hole-doping region, two different interesting orbital-selective Mott phases were found: OSMP1 (with one localized orbital and two itinerant orbitals) and OSMP2 (with two localized orbitals and one itinerant orbital). Moreover, charge disproportionation phenomena were found in special doping regions. We argue that our predictions can be tested by simple modifications in the original chemical formula of Ce 2 O 2 FeSe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Coupling solidification kinetics with phase-behavior computations in hydrodynamic simulations of high-pressure, dynamic-compression processes

In this study, we report a numerical scheme to integrate models for the kinetics of solidification processes together with phase-behavior computations in the context of continuum-scale hydrodynamic simulations. The objective of the phase-behavior computations is to determine the pressure and temperature, given the following three sets of inputs: (1) an appropriate equation of state to describe our system, (2) the phase fraction(s) produced by the kinetic models, (3) and the volume and internal energy obtained by solving the conservation equations that govern the hydrodynamic behavior. The kinetics are assumed to be governed by the Kolmogorov–Johnson–Mehl–Avrami equation, and the nucleation and growth rates that enter into that equation are functions of the pressure and temperature produced by the phase-behavior computations. Our formulation allows for the fluid and solid phases to be at different temperatures (thermal nonequilibrium) and pressures (arising from surface-tension-induced Laplace contributions). The formulation is presented in a fairly general setting that is independent of any particular material, although we demonstrate it in some examples related to high-energy-density science applications where materials are rapidly compressed to pressures exceeding several gigapascals in less than a microsecond. We conclude with a critical evaluation of our approach and provide suggestions for future work to improve the predictive capabilities and generality of the models.

42 ENGINEERING↗

DOE OSTI2021

The Fundamental Creep Behavior Model of Gr.91 Alloy by ICME Approach (Final Report)

The current project mainly focused on the fundamental creep cracking mechanism of the Grade 91 system during the operation conditions in the advanced technologies FE power plants and build the link among Composition-Processing Parameters-Phase Stability-Microstructure-Creep Resistance. A model based on computational thermodynamics and diffusion kinetics will be developed to provide guidance on how to improve the creep resistance of the alloy system. The PI mainly investigate the Gr.91 base alloy and weldment with the Integrated Computational Materials Engineering (ICME) approach. The long-term goal of the proposed program is to develop a model for creep and fatigue resistant alloys with different elemental systems, compositions, and processing parameters of weldment and heat-treatment. In addition, new creep and fatigue resistant alloys are to be designed based on the predictions from the model. The computational as well as experimental results include three sections: 1) The compositional optimization and secondary phases evaluation regarding the creep resistance in Grade 91 steel through the CALPHAD approach. In this section, the formation of the critical secondary phases, M23C6, MX and Z-phase were well predicted based on many conditions like the temperature as well as compositions. Meanwhile, the role of different alloying elements was also considered as they will directly affect the formation of the secondary phases mentioned above based on the compositional evolution of the alloying elements. 2) An investigation of Creep Resistance in Grade 91 Steel through Computational Thermodynamics. We systematically considered the evolution of the four critical temperatures, Ac3, Ac1, the threshold of M23C6 and Z-phase, on basis of the computational approach. Meanwhile, the equilibrium cooling as well as Scheil simulations were also considered to perform further predictions based on the different cooling rates. 3) Creep lifetime and microstructural evolution of the Grade 91 steel. In this part, mainly experimental approach was performed on the real creep test from the beginning of the commercial Gr.91 alloys. The alloys were normalized, tempered, welded (Gleeble) and PWHT and finally subjected to the creep lifetime analyses. Finally, all the samples after each process were further characterized based on the OM, XRD and SEM technique to draw final conclusions.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Design Principles to Govern Electrode Fabrication for the Lithium Trivanadate Cathode

A full depth of discharge mathematical model for the lithium trivanadate cathode, considering lithiation of the layered α-phase, phase change, and lithiation of the rock-salt like β-phase at lower potentials, is developed. The coupled electrode-scale and crystal-scale model is fit to electrochemical data, and additionally validated with operando EDXRD. There is good agreement between the simulated and measured spatial variation of the volume fraction of the β-phase. This mathematical model is used to guide electrode fabrication, accounting for both ionic and electronic transport effects. Values of three design parameters—electrode thickness, porosity, and volume fraction of conductor—are identified, and the sensitivity of the energy density to these design parameters is quantified. The model is also used to investigate electrode design to create electrodes that deliver the maximum achievable energy density under the constraint that the α to β-phase transition is avoided, since phase change has been demonstrated to reduce cycle life. The energy density sacrificed to avoid phase change decreases at higher discharge rates, but the target values for electrode fabrication remain the same as those when optimizing the electrode for the full depth of discharge.

25 ENERGY STORAGE↗

DOE OSTI2020

Cooperative Model Predictive Control Scheme for Dispersed Smart Inverters at the Grid Edge

Growing number of distributed energy resources (DERs) requires cutting-edge control schemes to ensure uniform but not identical involvement of DERs to address rapid fluctuations on voltage of the grid clusters. Designing such control schemes for a grid with high penetration of DERs connected various clusters is complex and depends on physical configuration of the clusters themselves. This paper proposes a two-layered cooperative model predictive control scheme for fleet of dispersed smart inverters connected to a three-phase grid cluster located at the edge of a grid cluster. The upper layer of the controller ensures cooperative operation of fleet of smart inverters to perform voltage regulation and address voltage fluctuations in the grid cluster; while the lower level controller guarantees each individual inverter obeys the dictated references for active/reactive powers by the upper layer controller. In addition, the proposed twolayered controller has low voltage ride through (LVRT) capability to support the grid clusters during voltage sags. By employing the proposed control approach, the smart inverters will inject desired active/reactive powers based on circumstance of the entire cluster while guaranteeing proper power curtailment among the generation units. The control scheme is verified by simulation.

Fard, Amin Y.↗

DOE OSTI2020

Phase Shift Control of a Three-Phase Inverter for Balanced Secondary Currents in Misaligned Three-Phase Inductive Power Transfer Systems

This paper presents a strategy for generating fundamental reference voltage commands for a double-sided LCC tuned three-phase wireless power transfer (WPT) system to achieve balanced output currents when the receiver and transmitter are not perfectly aligned. Formulas for translating the normalized fundamental current commands into phase-shift commands for two of the phase legs in a two-level three-phase inverter are derived. The methodology is validated through circuit simulations performed using coupler models derived from finite-element simulations. Simulation results for a 1.5 k W system design are included, which indicate near perfect correction of imbalance in the secondary currents for a typical misalignment case.

Pries, Jason↗

DOE OSTI2022

Predicting beam transmission using 2-dimensional phase space projections of hadron accelerators

We present a method to compress the 2D transverse phase space projections from a hadron accelerator and use that information to predict the beam transmission. This method assumes that obtaining at least three projections of the 4D transverse phase space is possible and that an accurate simulation model is available for the beamline. Using a simulated model, we show that—a computer can train a convolutional autoencoder to reduce phase-space information which can later be used to predict the beam transmission. Finally, we argue that although using projections from a realistic nonlinear distribution produces less accurate results, the method still generalizes well.

43 PARTICLE ACCELERATORS↗

NASA NTRS1984

Structure of nonevaporating sprays - Measurements and predictions

Structure measurements were completed within the dilute portion of axisymmetric nonevaporating sprays (SMD of 30 and 87 microns) injected into a still air environment, including: mean and fluctuating gas velocities and Reynolds stress using laser-Doppler anemometry; mean liquid fluxes using isokinetic sampling; drop sizes using slide impaction; and drop sizes and velocities using multiflash photography. The new measurements were used to evaluate three representative models of sprays: (1) a locally homogeneous flow (LHF) model, where slip between the phases was neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop interaction with turbulent fluctuations were ignored; and (3) a stochastic separated flow (SSF) model, where effects of both interphase slip and turbulent fluctuations were considered using random sampling for turbulence properties in conjunction with random-walk computations for drop motion. The LHF and DSF models were unsatisfactory for present test conditions-both underestimating flow widths and the rate of spread of drops. In contrast, the SSF model provided reasonably accurate predictions, including effects of enhanced spreading rates of sprays due to drop dispersion by turbulence, with all empirical parameters fixed from earlier work.

Solomon, A. S. P.↗