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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

The Postshock Chemical Lifetimes of Outflow Tracers and a Possible New Mechanism to Produce Water Ice Mantles

We have used a coupled time-dependent chemical and dynamical model to investigate the lifetime of the chemical legacy in the wake of C-type shocks. We concentrate this study on the chemistry of H2O and O2, two molecules which are predicted to have abundances that are significantly affected in shock-heated gas. Two models are presented: (1) a three-stage model of preshock, shocked, and postshock gas; and (2) a Monte Carlo cloud simulation where we explore the effects of stochastic shock activity on molecular gas over a cloud lifetime. For both models we separately examine the pure gas-phase chemistry as well as the chemistry including the interactions of molecules with grain surfaces. In agreement with previous studies, we find that shock velocities in excess of 10 km/s are required to convert all of the oxygen not locked in CO into H2O before the gas has an opportunity to cool. For pure gas phase models the lifetime of the high water abundances, or "H2O legacy," in the postshock gas is approximately (4-7) x 10(exp 5) yr, independent of the gas density. A density dependence for the lifetime of H2O is found in gas-grain models as the water molecules deplete onto grains at the depletion timescale. Through the Monte Carlo cloud simulation we demonstrate that the time-average abundance of H2O, the weighted average of the amount of time gas spends in preshock, shock, and postshock stages, is a sensitive function of the frequency of shocks. Thus we predict that the abundance of H2O, and to a lesser extent O2, can be used to trace the history of shock activity in molecular gas. We use previous large-scale surveys of molecular outflows to constrain the frequency of 10 km/s shocks in regions with varying star formation properties and discuss the observations required to test these results. We discuss the postshock lifetimes for other possible outflow tracers (e.g., SiO and CH3OH) and show that the differences between the lifetimes for various tracers can produce potentially observable chemical variations between younger and older outflows. For gas-grain models we find that the abundance of water-ice on grain surfaces can be quite large and is comparable to that observed in molecular clouds. This offers a possible alternative method to create water mantles without resorting to grain surface chemistry: gas heating and chemical modification due to a C-type shock and subsequent depletion of the gas-phase species onto grain mantles.

Bergin, Edwin A.↗
NASA NTRS2006

Implications of Cometary Water: Deep Impact, Stardust and Hayabusa

Three recent in situ spacecraft missions have explored comets or asteroids, producing data in conflict with the standard comet paradigm, the Whipple Dirty Snowball Model (DSM). We have developed an alternative Wet Comet Model (WCM) which proposes that comets undergo an irreversible phase change to a wet comet when they enter within Mars orbit. The WCM may explain some of the observational discrepancies seen by Deep Impact, Stardust and Hayabusa. In particular, it accurately predicted Deep Impact observation of organics, biominerals, and meltwater temperatures. Predictions concerning Stardust s returned cometary dust particles have yet to be falsified, but if comets are largely composed of the silicates seen by Stardust, there may be a cometary explanation for Itokawa s low density rubble-pile observed by Hayabusa.

Sheldon, Robert B.↗
NASA NTRS2004

Theoretical Investigation of OCN(-) Charge Transfer Complexes in Condensed Phase Media: Spectroscopic Properties in Amorphous Ice

Density functional theory (DFT) calculations of cyanate (OCN(-)) charge-transfer complexes were performed to model the "XCN" feature observed in interstellar icy grain mantles. OCN(-) charge-transfer complexes were formed from precursor combinations of HNCO or HOCN with either NH3 or H2O. Three different solvation strategies for realistically modeling the ice matrix environment were explored, including (1) continuum solvation, (2) pure DFT cluster calculations, and (3) an ONIOM DFT/PM3 cluster calculation. The model complexes were evaluated by their ability to reproduce seven spectroscopic measurements associated with XCN: the band origin of the OCN(-) asymmetric stretching mode, shifts in that frequency due to isotopic substitutions of C, N, O, and H, plus two weak features. The continuum solvent field method produced results consistent with some of the experimental data but failed to account for other behavior due to its limited capacity to describe molecular interactions with solvent. DFT cluster calculations successfully reproduced the available spectroscopic measurements very well. In particular, the deuterium shift showed excellent agreement in complexes where OCN(-) was fully solvated. Detailed studies of representative complexes including from two to twelve water molecules allowed the exploration of various possible solvation structures and provided insights into solvation trends. Moreover, complexes arising from cyanic or isocyanic acid in pure water suggested an alternative mechanism for the formation of OCN(-) charge-transfer complexes without the need for a strong base such as NH3 to be present. An extended ONIOM (B3LYP/PM3) cluster calculation was also performed to assess the impact of a more realistic environment on HNCO dissociation in pure water.

Park, Jin-Young↗
NASA NTRS2015

Validating and Verifying Biomathematical Models of Human Fatigue

Airline pilots experience acute and chronic sleep deprivation, sleep inertia, and circadian desynchrony due to the need to schedule flight operations around the clock. This sleep loss and circadian desynchrony gives rise to cognitive impairments, reduced vigilance and inconsistent performance. Several biomathematical models, based principally on patterns observed in circadian rhythms and homeostatic drive, have been developed to predict a pilots levels of fatigue or alertness. These models allow for the Federal Aviation Administration (FAA) and commercial airlines to make decisions about pilot capabilities and flight schedules. Although these models have been validated in a laboratory setting, they have not been thoroughly tested in operational environments where uncontrolled factors, such as environmental sleep disrupters, caffeine use and napping, may impact actual pilot alertness and performance. We will compare the predictions of three prominent biomathematical fatigue models (McCauley Model, Harvard Model, and the privately-sold SAFTE-FAST Model) to actual measures of alertness and performance. We collected sleep logs, movement and light recordings, psychomotor vigilance task (PVT), and urinary melatonin (a marker of circadian phase) from 44 pilots in a short-haul commercial airline over one month. We will statistically compare with the model predictions to lapses on the PVT and circadian phase. We will calculate the sensitivity and specificity of each model prediction under different scheduling conditions. Our findings will aid operational decision-makers in determining the reliability of each model under real-world scheduling situations.

fatigue↗
DOE OSTI2021

Mixed-phase clouds and precipitation in Southern Ocean cyclones and cloud systems observed poleward of 64°S by ship-based cloud radar and lidar

Mixed-phase clouds (MPCs), composed of both liquid and ice, are prevalent in Southern Ocean cyclones. A characterisation of these clouds on fine vertical scales is required in order to understand the microphysical processes within these clouds, and for model and satellite evaluations over this region. We investigated three examples of cloud systems collected by ship-mounted remote sensing instruments adjacent to East Antarctica at latitudes between 64°S and 69°S. These cases allow us to examine the properties of mid-level MPCs, with cloud tops between 2 and 6 km. Mid-level MPCs contain multiple layers of supercooled liquid water (SLW) embedded within ice during the passage of cyclones. SLW layers are capped by strong temperature inversions and are observed at temperatures as low as - 31°C. Convective generating cells (GCs) are present inside supercooled liquid-topped mid-level MPCs. The horizontal extent, vertical extent and maximum upward Doppler velocity of these GCs were 0.6 – 3.6 km, 0.7 – 1.0 km and 0.5 – 1.0 m s -1 , respectively, and are consistent with observations from previous lower-latitude studies. Ice precipitation is nearly ubiquitous, except in the thinnest clouds at the trailing end of the observed systems. Seeding of lower SLW layers from above leads to periods with either larger ice particles or greater ice precipitation rates. Additionally, periods of supercooled drizzle lasting up to two hours were observed toward the end of two of the three cyclone systems. This supercooled drizzle turns into predominantly ice precipitation as the result of seeding by ice clouds located above the precipitating SLW layer.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Automation of reverse engineering process in aircraft modeling and related optimization problems

During the year of 1994, the engineering problems in aircraft modeling were studied. The initial concern was to obtain a surface model with desirable geometric characteristics. Much of the effort during the first half of the year was to find an efficient way of solving a computationally difficult optimization model. Since the smoothing technique in the proposal 'Surface Modeling and Optimization Studies of Aerodynamic Configurations' requires solutions of a sequence of large-scale quadratic programming problems, it is important to design algorithms that can solve each quadratic program in a few interactions. This research led to three papers by Dr. W. Li, which were submitted to SIAM Journal on Optimization and Mathematical Programming. Two of these papers have been accepted for publication. Even though significant progress has been made during this phase of research and computation times was reduced from 30 min. to 2 min. for a sample problem, it was not good enough for on-line processing of digitized data points. After discussion with Dr. Robert E. Smith Jr., it was decided not to enforce shape constraints in order in order to simplify the model. As a consequence, P. Dierckx's nonparametric spline fitting approach was adopted, where one has only one control parameter for the fitting process - the error tolerance. At the same time the surface modeling software developed by Imageware was tested. Research indicated a substantially improved fitting of digitalized data points can be achieved if a proper parameterization of the spline surface is chosen. A winning strategy is to incorporate Dierckx's surface fitting with a natural parameterization for aircraft parts. The report consists of 4 chapters. Chapter 1 provides an overview of reverse engineering related to aircraft modeling and some preliminary findings of the effort in the second half of the year. Chapters 2-4 are the research results by Dr. W. Li on penalty functions and conjugate gradient methods for quadratic programming problems.

Li, W.↗
DOE OSTI2020

Multi-scale dynamics of atmospheric-pressure discharges ignited over liquid electrodes

A multiphase computational model of atmospheric-pressure gas discharges ignited over liquid electrodes is developed. The model takes into account both the liquid electrode motion under applied electric fields and plasma generation in the gas phase. The influence of the applied voltage and liquid properties (density, viscosity, and surface tension) on the liquid surface deformation is analyzed. It is shown that the shape of dynamic cones formed on the liquid surface differs from the static Taylor's cones. The influence of the liquid surface protrusions on gas breakdown dynamics is demonstrated. It is shown that the breakdown develops in two stages: first, a fast ionization wave propagates from the cathode to the anode, and then the cathode sheath collapses due to secondary electron emission from the liquid surface. It is shown that dynamics of the multiphase system containing gas, plasma, and liquid states is characterized by three disparate time scales: the fast electron time scale, the intermediate ion time scale, and the slow liquid dynamics time scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Extending a scaling equation of state to QCD

Whether quantum chromodynamics (QCD) exhibits a phase transition at finite temperature and density is an open question. It is important for hydrodynamic modeling of heavy ion collisions and neutron-star mergers. Lattice QCD simulations have definitively shown that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Here, we combine the parametric scaling equation of state, usually associated with the three-dimensional Ising model, with a background equation of state based on a smooth crossover from hadrons to quarks and gluons. Comparison to experimental data from the Beam Energy Scan II at the Relativistic Heavy Ion Collider or in heavy ion experiments at other accelerators may allow the critical exponents and amplitudes in the scaling equation of state to be determined for QCD if a critical point exists.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Modelling the limiter ramp-up of WEST for addressing the future challenges of ITER

This paper presents a joint experimental and numerical investigation into the physics of long limited plasma ramp-up in tokamaks with tungsten (W) first walls, a critical phase for ITER operation. The comparison between the average plasma quantities simulated using the SolEdge-HDG code and the measurements taken during three successive WEST discharges after boronisation shows how challenging it is to predict this phase. While simulations reproduce the general trends at the midplane, with reasonable match in density profiles, they consistently underestimate core temperatures possibly due to too large perpendicular heat conductivity. On the contrary, at the high field side (HFS) limiter, simulations overestimate the measured quantities, and highlights the limitation of using Bohm boundary conditions at grazing magnetic angles. Experimental measurements reveal that the boron layer is rapidly eroded, on a timescale comparable to a single ITER discharge. The subsequent transition from a boron-coated to a tungsten wall increases recycling and significantly degrades the core plasma, reducing the electron temperature by nearly half due to W contamination, despite wall parameters remaining stable. Furthermore, comparisons with Langmuir probes, bolometry, reflectometry, and spectroscopy indicate that the experimental far scrape-Off layer (SOL) is significantly wider than simulated. This wide SOL implies that boron erosion extends along the entire HFS limiter rather than being confined to the contact point. This work highlights some characteristics of the plasma during this phase of the discharge and emphasizes the current modelling issues that need to be resolved in order to obtain reliable predictions concerning the ITER ramp-up.

ITER↗
DOE OSTI2022

Snowmass Neutrino Frontier: DUNE Physics Summary

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long-baseline neutrino oscillation experiment with a primary physics goal of observing neutrino and antineutrino oscillation patterns to precisely measure the parameters governing long-baseline neutrino oscillation in a single experiment, and to test the three-flavor paradigm. DUNE's design has been developed by a large, international collaboration of scientists and engineers to have unique capability to measure neutrino oscillation as a function of energy in a broadband beam, to resolve degeneracy among oscillation parameters, and to control systematic uncertainty using the exquisite imaging capability of massive LArTPC far detector modules and an argon-based near detector. DUNE's neutrino oscillation measurements will unambiguously resolve the neutrino mass ordering and provide the sensitivity to discover CP violation in neutrinos for a wide range of possible values of $\delta_{CP}$. DUNE is also uniquely sensitive to electron neutrinos from a galactic supernova burst, and to a broad range of physics beyond the Standard Model (BSM), including nucleon decays. DUNE is anticipated to begin collecting physics data with Phase I, an initial experiment configuration consisting of two far detector modules and a minimal suite of near detector components, with a 1.2 MW proton beam. To realize its extensive, world-leading physics potential requires the full scope of DUNE be completed in Phase II. The three Phase II upgrades are all necessary to achieve DUNE's physics goals: (1) addition of far detector modules three and four for a total FD fiducial mass of at least 40 kt, (2) upgrade of the proton beam power from 1.2 MW to 2.4 MW, and (3) replacement of the near detector's temporary muon spectrometer with a magnetized, high-pressure gaseous argon TPC and calorimeter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2000

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗
NASA NTRS1997

Surface Modeling and Grid Generation of Orbital Sciences X34 Vehicle

The surface modeling and grid generation requirements, motivations, and methods used to develop Computational Fluid Dynamic volume grids for the X34-Phase 1 are presented. The requirements set forth by the Aerothermodynamics Branch at the NASA Langley Research Center serve as the basis for the final techniques used in the construction of all volume grids, including grids for parametric studies of the X34. The Integrated Computer Engineering and Manufacturing code for Computational Fluid Dynamics (ICEM/CFD), the Grid Generation code (GRIDGEN), the Three-Dimensional Multi-block Advanced Grid Generation System (3DMAGGS) code, and Volume Grid Manipulator (VGM) code are used to enable the necessary surface modeling, surface grid generation, volume grid generation, and grid alterations, respectively. All volume grids generated for the X34, as outlined in this paper, were used for CFD simulations within the Aerothermodynamics Branch.

Alter, Stephen J.↗
DOE OSTI2025

Enhancing nanoscale charged colloid crystallization near a metastable liquid binodal

Achieving predictive control over crystallization using non-classical nucleation while avoiding kinetic traps would be a step towards designing materials with new functionalities. We address these challenges by inducing the bottom-up assembly of nanocrystals into ordered arrays, or superlattices. Using electrostatics—rather than density—to tune the interactions between particles, we watch self-assembly proceed through a metastable liquid phase. Here, we systematically investigate the phase behaviour as a function of quench conditions in situ and in real time using small-angle X-ray scattering. By fitting to colloid, liquid and superlattice models, we extract the time evolution of each phase and the system phase diagram, which we find to be consistent with short-range attractive interactions. Using the predictive power of the phase diagram, we establish control of the self-assembly rate over three orders of magnitude, and we identify one- and two-step self-assembly regimes, with only the latter implicating the metastable liquid as an intermediate. The presence of the metastable liquid increases the superlattice formation rate relative to the equivalent one-step pathway, and the superlattice order increases with the rate, revealing a generalizable kinetic strategy for promoting and enhancing ordered assembly.

Tanner, Christian P. N. [University of California,↗
DOE OSTI2022

Comparison of Power System Current Sensors via Playback of Electrical Disturbances

The need to accurately measure high-frequency content in power system voltage and current phenomena is increasingly becoming more of a priority. As the amount of distributed energy resources (DER) and nonlinear loads penetrating the grid increases, so do challenges associated with traditional measurement and metering applications. In this paper, three commercially-available medium-voltage-level current sensors are characterized in terms of their harmonic amplitude and phase performance against reference signals that are “played back” through the sensors through the use of an arbitrary waveform generator. It is shown that none of the three sensors studied are able to faithfully replicate all of the input signals completely, though there are advantages and disadvantages to each in terms of noise, resonance, and induced phase drift. Additionally, the Goodness-of-Fit metric, typically used for PMU model validation, is used to generate side-by-side comparisons of sensor accuracy over a small window around the events under study.

Wilson, Aaron↗
DOE OSTI2025

CyRRL (Cyber Resilient Reinforcement Learning for grid voltage control) [SWR-24-115]

This codebase contains a multi-agent, actor-critic reinforcement learning implementation for cyber-resilient grid voltage control. It uses a 123-bus OpenDSS system as the environment, with three-phase power flow translating nodal power injections into solved nodal voltages. The reward function penalizes deviations from nominal voltage as well as reactive power dispatch, while encouraging agents to take actions that result in fast convergence to nominal conditions. The codebase models false data injection attacks and includes functionality for training, testing, hyper-parameter tuning, and visualization.

Murphy, Sinnott [National Renewable Energy Laborat↗
DOE OSTI2021

Post-explosion Evolution of Core-collapse Supernovae

Here we investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 explosion models during several seconds. We present a broad study based on three progenitors (11.2, 15, and 27 M ⊙ ), different neutrino heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind (NDW) that survives for several seconds. This indicates that NDWs are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 M ⊙ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Time dependences of the 0.51 and 2.2 MeV lines in solar flares

The time dependence of the 0.51 MeV line in solar flares is determined by the following factors: variation with time of the rate of nuclear reactions which produce the positron emitters; mean lives of the positron emitters; slowing-down of the positrons which depends on the ambient density and magnetic field; and annihilation or positronium formation times of the positrons with free or bound electrons. The discussion is limited to the observed intensity-time profile of a 0.35-8 MeV gamma-ray emission and of a 0.51-MeV line calculated for three different ambient densities in the annihilation region. Major conclusions are that (1) MeV electrons and protons of at least 30 MeV/nuc have similar time dependences and hence a common origin in the flare region, which support a two-phase acceleration model in solar flares; and (2) the delayed character of the 0.51-MeV emission is clearly defined, where at a time when the rate of nuclear reactions is zero, there still is considerable 0.51-MeV emission.

Ramaty, R.↗
NASA NTRS1988

Utility of Emulation and Simulation Computer Modeling of Space Station Environmental Control and Life Support Systems

Over the years, computer modeling has been used extensively in many disciplines to solve engineering problems. A set of computer program tools is proposed to assist the engineer in the various phases of the Space Station program from technology selection through flight operations. The development and application of emulation and simulation transient performance modeling tools for life support systems are examined. The results of the development and the demonstration of the utility of three computer models are presented. The first model is a detailed computer model (emulation) of a solid amine water desorbed (SAWD) CO2 removal subsystem combined with much less detailed models (simulations) of a cabin, crew, and heat exchangers. This model was used in parallel with the hardware design and test of this CO2 removal subsystem. The second model is a simulation of an air revitalization system combined with a wastewater processing system to demonstrate the capabilities to study subsystem integration. The third model is that of a Space Station total air revitalization system. The station configuration consists of a habitat module, a lab module, two crews, and four connecting nodes.

Yanosy, James L.↗