Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Three-Dimensional Unsteady Simulation of Aerodynamics and Heat Transfer in a Modern High Pressure Turbine Stage

Unsteady 3-D RANS simulations have been performed on a highly loaded transonic turbine stage and results are compared to steady calculations as well as to experiment. A low Reynolds number k-epsilon turbulence model is employed to provide closure for the RANS system. A phase-lag boundary condition is used in the tangential direction. This allows the unsteady simulation to be performed by using only one blade from each of the two rows. The objective of this work is to study the effect of unsteadiness on rotor heat transfer and to glean any insight into unsteady flow physics. The role of the stator wake passing on the pressure distribution at the leading edge is also studied. The simulated heat transfer and pressure results agreed favorably with experiment. The time-averaged heat transfer predicted by the unsteady simulation is higher than the heat transfer predicted by the steady simulation everywhere except at the leading edge. The shock structure formed due to stator-rotor interaction was analyzed. Heat transfer and pressure at the hub and casing were also studied. Thermal segregation was observed that leads to the heat transfer patterns predicted by steady and unsteady simulations to be different.

Shyam, Vikram↗
NASA NTRS2003

Influence of Ice-phase of Hydrometeors on Moist-Convection

Climate models often ignore the influence of ice-phase physics (IPP) of hydrometeors as a second order effect. This has also been true for McRAS (Microphysics of clouds with Relaxed Arakawa Schubert Scheme) developed by the authors. Recognizing that the temperature sounding is critical for moist-convection, and, that IPP would modify it, we investigated the influence of introducing IPP into McRAS coupled to FvGCM (finite volume General Circulation Model with NCAR physics). We analyzed three 3-yr long simulations; the first called Control Case, CC and had no IPP; the other two called Experiments El and E2 had IPP introduced with two different in-cloud freezing assumptions. Simulation El assumed that all hydrometeors remain liquid in the updraft and freeze upon detrainment. Simulation E2 invoked the in-cloud freezing of new condensate generated at subfreezing temperatures in the updraft while old cloud water continued to ascend as liquid. Upon detrainment, this cloud water also froze like in E1. With these assumptions, about 50% of hydrometeors froze in the tower and the rest froze in the anvil. However, in both El and E2, the frozen hydrometeors melted during fall at the first encounter of above freezing ambient temperature. Comparative analysis revealed that El simulated far more mid-level and far less deep clouds while E2 had modified deep and more mid-level clouds as compared to CC along with some major changes around the melt-level. We infer that IPP produced a more realistic response in E2. At the basic level, the results show that ice-phase processes influence convective detrainment at mid- and deep levels in accord with TOGAGOARE observations. The results suggest that IPP can help to mitigate less-than-observed mid-level and over-abundance of deep convective clouds in McRAS.

Sud, Y. C.↗
NASA NTRS2003

The Growth of Central Black Hole and the Ionization Instability of Quasar Disk

A possible accretion model associated with the ionization instability of quasar disks is proposed to address the growth of the central black hole harbored in the host galaxy. The evolution of quasars in cosmic time is assumed to change from a highly active state to a quiescent state triggered by the S-shaped ionization instability of the quasar accretion disk. For a given external mass transfer rate supplied by the quasar host galaxy, ionization instability can modify accretion rate in the disk and separates the accretion flows of the disk into three different phases, like a S-shape. We suggest that the bright quasars observed today are those quasars with disks in the upper branch of S-shaped instability, and the faint or 'dormant' quasars are simply the system in the lower branch. The middle branch is the transition state which is unstable. We assume the quasar disk evolves according to the advection-dominated inflow-outflow solutions (ADIOS) configuration in the stable lower branch of S-shaped instability, and Eddington accretion rate is used to constrain the accretion rate in each phase. The mass ratio between black hole and its host galactic bulge is a nature consequence of ADIOS. Our model also demonstrates that a seed black hole (BH) similar to those found in spiral galaxies today is needed to produce a BH with a final mass 2 x 10(exp 8) solar mases.

Lu, Ye↗
DOE OSTI2021

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Thermodynamic Modeling of Aqueous LiCl, LiBr, LiI, and LiNO 3 Solutions

Thermodynamic models are essential to facilitate the advancement of process design, optimization, and operation of electrolyte systems. In this work, a comprehensive thermodynamic framework based on the Electrolyte Nonrandom Two-Liquid model is developed to calculate phase equilibria behavior and salt solubility of aqueous LiCl, LiBr, LiI, and LiNO 3 solutions. The model describes the non-ideality of the electrolyte solutions by using two binary interaction parameters for each electrolyte-molecule pair in the system. Each binary interaction parameter is further expressed with up to three temperature coefficients which are regressed from experimental data. To take into account the hydration of lithium ion, two separate chemistries for the dissociation of lithium salts are investigated. In the first case, the lithium ion is considered as a bare ion, Li + , while in the second case hydration of the lithium ion from Li + to Li(H 2 O) + is considered. Here, the calculated thermodynamic properties compare adequately with the experimental data for both sets of chemistries for concentrations up to saturation and temperatures from 273.15 K up to 623.15 K. Moderate to significant improvements are observed with the incorporation of the hydration chemistry for aqueous LiCl, LiBr, and LiI solutions when compared to the non-hydrated lithium ion model results.

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

Analysis of ultrafast magnetization switching dynamics in exchange-coupled ferromagnet–ferrimagnet heterostructures

Magnetization switching in ferromagnets has so far been limited to the current-induced spin–orbit–torque effects. Recent observation of helicity-independent all-optical magnetization switching (HI-AOS) in an exchange-coupled ferromagnet–ferrimagnet (FM-FEM) heterostructures expanded the range and applicability of such ultrafast heat-driven magnetization switching. Here we report the element-resolved HI-AOS dynamics of such an exchange-coupled system, using a modified microscopic three-temperature model. We have studied the effect of (i) the Curie temperature of the FM, (ii) FEM composition, (iii) the long-range Ruderman–Kittel–Kasuya–Yosida (RKKY) exchange-coupling strength, and (iv) the absorbed optical energy on the element-specific time-resolved magnetization dynamics. The phase-space of magnetization illustrates how the RKKY coupling strength and the absorbed optical energy influence the switching time. Our analysis demonstrates that the threshold switching energy depends on the composition of the FEM and the switching time depends on the Curie temperature of the FM as well as RKKY coupling strength. Importantly, this simulation anticipates new insights into developing faster and more energy-efficient spintronics devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamics of Liquid–Liquid Phase Separation in Submicrometer Aerosol

Nanoscale materials, when compared to their bulk components, possess unique properties. In particular, shifts in phase transitions can occur for submicron particles. For instance, small particles do not undergo the process of liquid-liquid phase separation (LLPS). LLPS has applications in emulsions such as Janus particles, controllable morphology to create drug-rich phases during drug delivery, and is often observed in atmospheric aqueous aerosol particles. In atmospheric particles, LLPS is tracked as a function of particle water activity, which is equivalent to the relative humidity (RH) at equilibrium. We probed three organic/inorganic aerosol systems in the range of RH over which phase separation occurs (SRH). Our findings indicate that SRH for submicron aerosol particles is lower than for micron-sized droplets. These findings show that it may be necessary to update the representation of phase transitions in aerosol particles in climate models. The vast majority of organic/inorganic aerosol particles have submicron diameters, and a decrease in SRH for submicron particles indicates that the current estimation of phase separated aerosols may be overestimated. Furthermore, understanding the properties of LLPS at the nanoscale can provide key parameters to describe these systems and may lead to better control of phase separation in submicron particles.

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

Theory of criticality for quantum ferroelectric metals

A variety of compounds, for example, doped paraelectrics and polar metals, exhibit both ferroelectricity and correlated electronic phenomena such as low-density superconductivity and anomalous transport. Characterizing such properties is tied to understanding the quantum dynamics of inversion symmetry breaking in the presence of itinerant electrons. Here, we present a comprehensive analysis of the properties of a metal near a quantum critical transition to a ferroelectric state, in both two and three dimensions. In this study, starting from a minimal model of electrons coupled to a transverse polar phonon via a Rashba-type spin-orbit interaction, we compute the dynamical response of both electrons and phonons. We find that the system can evince both Fermi and non-Fermi liquid phases, as well as enhanced pairing in both singlet and triplet channels. Furthermore, we systematically compute corrections to one-loop theory and find a tendency to quantum order-by-disorder, leading to a phase diagram that can include second-order, first-order, and finite-momentum phase transitions. Finally, we show that the entire phase diagram can be controlled via application of external strain, either compressive or volume-preserving. Our results provide a map of the dynamical and thermodynamical phase space of quantum ferroelectic metals, which can serve in characterizing existing materials and in seeking applications for quantum technologies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A Process Intensification Approach for CO 2 Absorption Using Amino Acid Solutions and a Guanidine Compound

Environmentally friendly amino-acid salt solutions are used for the absorption of carbon dioxide from concentrated flue-gas streams via chemical absorption. Process intensification reduces operating and capital costs by combining chemical reactions and separation operations. Here, we present a new process-intensification approach that combines the CO 2 capture and the amino-acid regeneration steps into a single process carried out in a slurry three-phase reactor. The absorbed CO 2 precipitates as a solid carbonated guanidine compound. The cycle is completed by separation of the solid precipitate to strip the CO 2 and regenerate the guanidine compound, while the liquid solution is recycled to the slurry reactor. The process was studied by modifying a model developed by the authors for a gas-liquid bubble column without the presence of the guanidine compound. The guanidine precipitation reaction was accounted for using kinetic parameters calculated by the authors in another study. The proposed model was implemented by modifying an existing computer code used for the simulation of gas-liquid bubble columns. The calculated results showed that the proposed cycle can significantly reduce energy, equipment, and operating costs and can make an important contribution to developing a competitive cost-effective large-scale process for CO 2 capture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1974

Current system in an auroral substorm

An instrumented sounding rocket was launched at 0040 LT into the expansive phase of an auroral substorm on March 16, 1973, from Poker Flat, Alaska. This report describes coordinated measurements of electric field, magnetic field, and energetic electrons. By comparing the observed magnetic field with that predicted by model currents, a current system, which includes field-aligned currents, is deduced. The model indicates three currents: a southward current, a westward current, and two field-aligned currents. A southward electric field of 60 mV/m identifies the southward current as a Pedersen current and the westward current as a Hall current. The field-aligned current measured from particle detectors shows a structure similar to that in the magnetic model. The total field-aligned electron current measured by particle detectors agrees within 10% of the current model.

Kintner, P. M.↗
NASA NTRS1989

A distributed finite-element modeling and control approach for large flexible structures

An unconventional framework is described for the design of decentralized controllers for large flexible structures. In contrast to conventional control system design practice which begins with a model of the open loop plant, the controlled plant is assembled from controlled components in which the modeling phase and the control design phase are integrated at the component level. The developed framework is called controlled component synthesis (CCS) to reflect that it is motivated by the well developed Component Mode Synthesis (CMS) methods which were demonstrated to be effective for solving large complex structural analysis problems for almost three decades. The design philosophy behind CCS is also closely related to that of the subsystem decomposition approach in decentralized control.

Young, K. D.↗
DOE OSTI2024

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗
NASA NTRS1990

Runway Exit Designs for Capacity Improvement Demonstrations. Phase 1: Algorithm Development

A description and results are presented of a study to locate and design rapid runway exits under realistic airport conditions. The study developed a PC-based computer simulation-optimization program called REDIM (runway exit design interactive model) to help future airport designers and planners to locate optimal exits under various airport conditions. The model addresses three sets of problems typically arising during runway exit design evaluations. These are the evaluations of existing runway configurations, addition of new rapid runway turnoffs, and the design of new runway facilities. The model is highly interactive and allows a quick estimation of the expected value of runway occupancy time. Aircraft populations and airport environmental conditions are among the multiple inputs to the model to execute a viable runway location and geometric design solution. The results presented suggest that possible reductions on runway occupancy time (ROT) can be achieved with the use of optimally tailored rapid runway designs for a given aircraft population. Reductions of up to 9 to 6 seconds are possible with the implementation of 30 m/sec variable geometry exits.

Trani, A. A.↗
DOE OSTI2024

Synthesis of results for Brine Availability Test in Salt (BATS) DECOVALEX-2023 Task E

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1996

Diamond Research Overview and a Model for Lab Experiments Using Oxyacetylene Torch

High pressure synthetic diamonds have now been a commercial product for over 40 years. Russian scientists invented the chemical vapor deposition (CVD) process 30 years ago, while the Japanese followed 10 years later and the U.S. was introduced to it 10 years after that. The new syntheses focus is on liquid and solid phase approaches. There are three CVD processes: microwave plasma, hot filament, and oxy-acetylene torch. The oxy-acetylene torch is an excellent materials synthesis lab experiment, emphasizing the simplicity of the science.

Roy, Rustum↗
DOE OSTI2020

Atmospheric ice nuclei concentration measurements over a high altitude-station in the Western Ghats, India

Mixed-phase and ice clouds play a significant role in modulating the Indian summer monsoon rainfall, and modeling studies suggest that inaccurate representation of ice nucleating particles (INP) is one of the factors, which have resulted in a deficiency in the prediction of cloud properties and hydrological cycle. We do not understand the ice nucleation mechanisms to represent the ice formation in the model, and therefore more INP measurements at various supercooled temperatures are needed. This study reports for the first time real-time INP concentration at a high altitude station in the Western Ghats region of India, and these results will be used to constrain the ice nucleation parameterizations in the model. The variation of INP for any particular cloud condition can have strong microphysical responses of mixed-phase precipitation formation processes. Also, to improve the representation of heterogeneous ice nucleation, INP measurements of various types of aerosols commonly observed over the Indian region are needed. This study presents, the INP concentrations and daily concentration of non-refractory chemical composition (organics, sulfate, nitrate, ammonium, and chloride), cloud condensation nuclei (CCN), and aerosol size distribution. Spectrometer for Ice Nuclei (SPIN) was operated to measure INP concentration at the mountain site from August to December 2018. INP measurements were performed at three different temperatures (-25, -30, and - 34 °C) in the immersion freezing mode relevant for mixed-phase cloud conditions. The average INP concentrations approximately varied from 0.18 to 12.4 L-1, 0.39 to 24 L-1, 1.1 to 40.2 L-1 at -25, -30, and - 34 °C, respectively. These concentrations are within the concentrations reported in the literature from different regions across the globe. The air mass back trajectory analysis showed that continental air mass contains more INP compared to maritime air mass. The information on INP concentration and chemical characteristics of aerosol will help to improve heterogeneous ice nucleation parameterization in numerical models.

Kumar, Anil V.↗
DOE OSTI2022

Numerical analysis of soot emissions from gasoline-ethanol and gasoline-butanol blends under gasoline compression ignition conditions

In the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine were performed to investigate the impact of blending two biofuels, ethanol and n-butanol, with gasoline on the trade-off between combustion phasing and soot emissions under low load conditions. Here, in order to represent market gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF) + ethanol (ETPRF) surrogate (with 20% ethanol by mole; E20) was formulated using a neural network based octane predictor such that the surrogate had the same ethanol content, Research Octane Number (RON) and Octane Sensitivity (S). In addition, a novel skeletal kinetic mechanism for ETPRF and TPRF + n-butanol (BTPRF) blends, incorporating polycyclic aromatic hydrocarbon (PAH) chemistry, was developed. A three-dimensional (3D) engine CFD formulation employing the skeletal mechanism, adaptive mesh refinement (AMR), finite-rate chemistry approach, and hybrid method of moments (HMOM) was adopted to capture the in-cylinder combustion phenomena and soot emissions. The engine CFD model was validated against RD5-87 experimental data for a broad range of start-of-injection (SOI) timings (-21/-27/-36/-45 crank angle degrees (CAD) after top-dead center (aTDC)), with respect to in-cylinder pressure, heat release rate, combustion phasing, and soot emissions. The closed-cycle simulation results were analyzed to elucidate the non-monotonic trend of soot emissions versus SOI timing: SOI-36 > SOI-45 > SOI-21 > SOI-27. Thereafter, the validated CFD model was employed to simulate the combustion of a gasoline-ethanol blend with 45% (by mole) ethanol (E45) and a gasoline-butanol blend with 45% (by mole) n-butanol (B45) under the same operating conditions to study the effects of fuel composition and SOI timing on combustion phasing and soot emissions. The sooting propensity followed the trend: B45 > E20 > E45 at all SOI timings. Overall, it was observed that the autoignition propensity was primarily related to fuel chemistry. On the other hand, sooting propensity showed strong coupling with both fuel chemistry and physical properties, with greater impact of fuel physical properties at advanced SOI timings.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

On CCN Effects upon Convective Cold Pool Timing and Features

Abstract Cold pools produced by deep convection can initiate new convection, and their representation in larger-scale weather and climate models could improve prediction of the extent and timing of upscale growth. Cold pools originate from latent cooling from precipitation changing phase, but little attention has been paid to microphysical influences on cold pool characteristics, particularly CCN effects. Datasets obtained from the CACTI and RELAMPAGO field campaigns, along with idealized numerical modeling, are utilized to investigate the hypothesis that convective storms forming in higher-CCN environments generate their first surface rainfall later, delaying cold pool initiation. Aircraft observations of CCN and shallow convection on 9 days do suggest a CCN effect. Those ingesting more CCN contained fewer drizzle drops, although a decreased cloud depth with increasing CCN was also likely a limiting factor. In three of those cases that later developed into deep convection, the timing of cold pool onset was not ubiquitously delayed in environments with more CCN. Idealized numerical simulations suggest that an ordinary thunderstorm can experience small delays in cold pool onset with increasing CCN due to changes in graupel production, but CCN effects on the cold pool from a supercell thunderstorm can be easily overpowered by its unique dynamics. A strong inverse relationship between cold pool strength, expansion rate, and depth with increasing CCN is suggested by the results of the ordinary thunderstorm simulation. Further consideration of CCN appears warranted for future cold pool parameterization development, but other environmental factors affecting storm morphology and precipitation cannot be ignored.

Meteorology & Atmospheric Sciences↗