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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

A Tool for Low Noise Procedures Design and Community Noise Impact Assessment: The Rotorcraft Noise Model (RNM)

To improve aircraft noise impact modeling capabilities and to provide a tool to aid in the development of low noise terminal area operations for rotorcraft and tiltrotors, the Rotorcraft Noise Model (RNM) was developed by the NASA Langley Research Center and Wyle Laboratories. RNM is a simulation program that predicts how sound will propagate through the atmosphere and accumulate at receiver locations located on flat ground or varying terrain, for single and multiple vehicle flight operations. At the core of RNM are the vehicle noise sources, input as sound hemispheres. As the vehicle "flies" along its prescribed flight trajectory, the source sound propagation is simulated and accumulated at the receiver locations (single points of interest or multiple grid points) in a systematic time-based manner. These sound signals at the receiver locations may then be analyzed to obtain single event footprints, integrated noise contours, time histories, or numerous other features. RNM may also be used to generate spectral time history data over a ground mesh for the creation of single event sound animation videos. Acoustic properties of the noise source(s) are defined in terms of sound hemispheres that may be obtained from theoretical predictions, wind tunnel experimental results, flight test measurements, or a combination of the three. The sound hemispheres may contain broadband data (source levels as a function of one-third octave band) and pure-tone data (in the form of specific frequency sound pressure levels and phase). A PC executable version of RNM is publicly available and has been adopted by a number of organizations for Environmental Impact Assessment studies of rotorcraft noise. This paper provides a review of the required input data, the theoretical framework of RNM's propagation model and the output results. Code validation results are provided from a NATO helicopter noise flight test as well as a tiltrotor flight test program that used the RNM as a tool to aid in the development of low noise approach profiles.

Conner, David A.↗
NASA NTRS1988

The dynamics and control of large flexible space structures, part 11

A mathematical model is developed to predict the dynamics of the proposed Spacecraft Control Laboratory Experiment during the stationkeeping phase. The Shuttle and reflector are assumed to be rigid, while the mass connecting the Shuttle to the reflector is assumed to be flexible with elastic deformations small as compared with its length. It is seen that in the presence of gravity-gradient torques, the system assumes a new equilibrium position primarily due to the offset in the mass attachment point to the reflector from the reflector's mass center. Control is assumed to be provided through the Shuttle's three torquers and throught six actuators located by painrs at two points on the mass and at the reflector mass center. Numerical results confirm the robustness of an LQR derived control strategy during stationkeeping with maximum control efforts significantly below saturation levels. The linear regulator theory is also used to derive control laws for the linearized model of the rigidized SCOLE configuration where the mast flexibility is not included. It is seen that this same type of control strategy can be applied for the rapid single axis slewing of the SCOLE through amplitudes as large as 20 degrees. These results provide a definite trade-off between the slightly larger slewing times with the considerable reduction in over-all control effort as compared with the results of the two point boundary value problem application of Pontryagin's Maximum Principle.

Bainum, Peter M.↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

What the Polar Cap Tells Us about the Substorm Growth Phase

The polar cap region in the 30 to 60 minute period prior to the onset of the auroral substorm has been examined using global images from the Polar Ultraviolet Imager (UVI) to look for observational evidence of processes related to the substorm growth phase. In particular, the area of the polar cap has been measured to determine changes in its size in relation to the orientation of the interplanetary magnetic field (IMF). It was found that the size of the polar cap region increases during the growth phase even if the IMF has no southward component. Three phenomena have been observed to produce the increase in the size of the polar cap: (1) motion of the auroral oval to lower latitude, (2) thinning of the auroral oval, and (3) reduction of intense auroral precipitation in the polar cap region. The first phenomenon has been considered to be a result of the growth of the tail lobe magnetic field and the second is related to the thinning of the plasma sheet. Both of these have been supported by in situ observational evidence and are consistent with current models of substorm development. However, the third phenomenon appears to be unrelated to the first two and does not appear to be the result of opening of the polar cap flux tubes to the solar wind IMF. This reduction of auroral precipitation provides evidence of a growth phase process, or change in auroral precipitation processes, that is not explained by current substorm models.

Brittnacher, M. J.↗
NASA NTRS1994

Lick optical spectra of quasar HS 1946+7658 at 10 kilometers per second resolution Lyman-alpha forest and metal absorption systems

We present optical spectra of the most luminous known quasi stellar object (QSO) HS 1946+7658 (z(sub em) = 3.051). Our spectra have both full wavelength coverage, 3240-10570 A, and in selected regions, either high signal-to-noise ratio, SNR approximately equals 40-100, or unusually high approximately 10 km/sec resolution, and in parts of the Ly alpha forest and to the red of Ly alpha emission they are among the best published. We find 113 Ly alpha systems and six metal-line systems, three of which are new. The metal systems at z(sub abs) = 2.844 and 3.050 have complex velocity structure with four and three prominent components, respectively. We find that the system at z(sub abs) = 2.844 is a damped Ly alpha absorption (DLA) system, with a neutral hydrogen column density of log N(H I) = 20.2 +/- 0.4, and it is the cause of the Lyman limit break at lambda approximately equals 3520 A. We believe that most of the H I column density in this system is in z(sub abs) = 2.8443 component which shows the strongest low-ionization absorption lines. The metal abundance in the gas phase of the system is (M/H) approximately equals -2.6 +/- 0.3, with a best estimate of (M/H) = -2.8, with ionizaion parameter log gamma = -2.75, from a photoionization model. The ratios of the logarithmic abundances of C, O, Al, and Si are all within a factor of 2 of solar, which is important for two reasons. First, we believe that the gas abundances which we measure are close to the total abundances, because the ratio of aluminum to other elements is near cosmic, and Al is a refractory element which depletes very readily like chromium, in the interstellar medium. Second, we do not see the enhancement of O with respect to C of (O/C) approximately equals 0.5-0.9 reported in three partial Lyman limit systems by Reimers et al. (1992) and Vogel & Reimers (1993); we measure (O/C) = -0.06 for observed ions and (O/C) approximately equals 0.2 after ionization corrections, which is consistent with solar abundances. We see C II*(lambda 1335) offset by 15 km/sec with respect to C II(lambda 1334), presumably because the gas density varies from 2 to 8 cm(exp -3) with changing velocity in the DLA system. These densities imply that the damped component is 6-25 pc thick, which is reasonable for a single cloud in a cold spiral disk. They also imply that the cloud is relatively highly ionized with more C III than C II, more O III than O I, and log N(H I) = 20.72, which is 3 times the H I column. The system at z(abs) = 1.7382 is also believed to be damped with N(H I) approximately equals 10(exp 21) cm(exp -2), because we see Cr II, but its Ly alpha line will never be seen because it is below the Lyman limit of the other DLA system. We see a 2.6 sigma lack of Ly alpha forest lines well away from the QSO redshift, which may be a chance fluctuation. We also see a correlation between column density N(H I) and Doppler parameter b for 96 unsaturated Ly Alpha forest absorption lines, and although this correlation persists in the 36 Ly alpha lines which lie in regions where the SNR approximately equals 8-16, we agree with Rauch et al. (1993) that it is probably a bogus effect of low supernova remnant (SNR). The same applies to lines with very low b values: in regions where SNR less than or equal to 8 we see many Ly alpha lines which appear to have 10 less than or equal to b less than or equal to 20, but when 8 less than or equal to SNR less than or equal to 16 we see only one line with b less than or equal to 15 km/sec, and two others which we believe have b less than or equal to 20, with values of 20 and 16 km/sec. Traditional Ly alpha line samples which include all lines which have W/sigma(W) greater than or equal to 4 are not adequate to explore the distribution of the properties of individual clouds, because we need much higher (W/sigma(W)) and SNR to avoid the strong biases.

Fan, Xiao-Ming↗
DOE OSTI2019

Numerical Simulation and Experimental Investigation of Nitrogen Transfer Mechanism from Gas to Liquid Steel During Pressurized Electroslag Remelting Process

Since the gas–slag–metal three-phase transfer mechanism of nitrogen is not suitable to account for the nitrogen pickup during pressurized electroslag remelting process, the laboratory experiments have been carried out to study the transfer mechanism of nitrogen using different nitrogen partial pressures and electrode immersion depths. Then, two possible transfer mechanisms of nitrogen have been proposed, and a 2D transient numerical model has been developed to investigate the behavior of nitrogen based on the new proposed mechanisms. The simulated results are compared with the experiment to validate the feasibility of new proposed mechanism. The results show that the nitrogen content in ingot increases with the increasing nitrogen partial pressure and decreasing electrode immersion depth. At the current experimental condition, the electrode immersion depth seems to play a more important role in the nitrogen pickup than the nitrogen partial pressure. Due to the limits of thermodynamic and kinetic conditions, the nitrogen pickup is negligible via the transfer from gas to liquid metal through molten slag whether nitrogen chemically or physically dissolves into slag. The mechanism that nitrogen directly reacts with the partially exposed liquid metal film under the electrode tip due to the fluctuation of gas/slag interface could reasonably account for the experimental results. Furthermore, based on the mechanism, the predicted variation tendency of nitrogen content in ingot is in agreement with the experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of Heave Plate Hydrodynamic Force Parameterization on a Two-Body Wave Energy Converter

Heave plates are one approach to generating the reaction force necessary to harvest energy from ocean waves. In a Morison equation description of the hydrodynamic force, the components of drag and added mass depend primarily on the heave plate oscillation. These terms may be parameterized in three ways: (1) as a single coefficient invariant across sea state, most accurate at the reference sea state, (2) coefficients dependent on the oscillation amplitude, but invariant in phase, that are most accurate for relatively small amplitude motions, and (3) coefficients dependent on both oscillation amplitude and phase, which are accurate for all oscillation amplitudes. We validate a MATLAB model for a two-body point absorber wave energy converter against field data and a dynamical model constructed in ProteusDS. We then use the MATLAB model to evaluate the effect of these parameterizations on estimates of heave plate motion, tension between the float and heave plate, and wave energy converter electrical power output. We find that power predictions using amplitude-dependent coefficients differ by up to 30% from models using invariant coefficients for regular waves ranging in height from 0.5 to 1.9 m. Amplitude- and phase-dependent coefficients, however, yield less than a 5% change when compared with coefficients dependent on amplitude only. This suggests that amplitude-dependent coefficients can be important for accurate wave energy converter modeling, but the added complexity of phase-dependent coefficients yields little further benefit. We show similar, though less pronounced, trends in maximum tether tension, but note that heave plate motion has only a weak dependence on coefficient fidelity. Finally, we emphasize the importance of using experimentally derived added mass over that calculated from boundary element methods, which can lead to substantial under-prediction of power output and peak tether tension.

dynamical model↗
DOE OSTI2026

Effects of B 2 O 3 on the Growth, Structural, and Magneto-Optical Properties of Yttrium Iron Garnet Single-Crystal Fibers

This study explores the fabrication of yttrium iron garnet (YIG) single crystal fibers using the laser heated pedestal growth (LHPG) method with the experimental addition of B 2 O 3 . The incorporation of B 2 O 3 facilitates the fiber fabrication process by lowering the required growth temperatures and likely modifying melt viscosity behavior, consistent with the established fluxing behavior of B2O3 and the comparative viscosity trend observed in the TMA−VFT analysis, thereby improving process efficiency while maintaining fiber quality. Structural characterization using EBSD and SC-XRD reveals a transition from polycrystalline to single-crystal behavior, with improved alignment along the [111] direction without altering the garnet structure. Magnetic measurements show increases in saturation magnetization in B 2 O 3 -assisted fibers. Three-dimensional anisotropy energy modeling, based on EBSD-derived Euler angles, indicates that the enhanced crystallinity and orientation contribute to reorientation of MCA energy distribution due to improved crystallographic alignment. Faraday rotation measurements show that the B 2 O 3 -assisted sample exhibits a rotation angle closer to reported values for high-quality YIG, suggesting improved phase purity and crystallographic quality. These findings demonstrate that B 2 O 3 -assisted LHPG growth is a scalable and nontoxic approach to producing high-performance YIG fibers for integrated photonic and magnetic field sensing applications.

Crystal structure↗
DOE OSTI2026

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonparametric extensions of nuclear equations of state: Probing the breakdown scale of relativistic mean-field theory

Phenomenological calculations of the properties of dense matter, such as relativistic mean-field theories, represent a pathway to predicting the microscopic and macroscopic properties of neutron stars. However, such theories do not generically have well-controlled uncertainties and may break down within neutron stars. To faithfully represent the uncertainty in this breakdown scale, we develop a hybrid representation of the dense-matter equation of state, which assumes the form of a relativistic mean-field theory at low densities, while remaining agnostic to any nuclear theory at high densities. To achieve this, we use a nonparametric equation of state model to incorporate the correlations of the underlying relativistic mean-field theory equation of state at low pressures and transition to more flexible correlations above some chosen pressure scale. We perform astrophysical inference under various choices of the transition pressure between the theory-informed and theory-agnostic models. Here, we further study whether the chosen relativistic mean-field theory breaks down above some particular pressure and find no such evidence. Using simulated data for future astrophysical observations at about two-to-three times the precision of current constraints, we show that our method can identify the breakdown pressure associated with a potential strong phase transition.

Equations of state of nuclear matter↗
DOE OSTI2024

Three-dimensional checkerboard spin structure on a breathing pyrochlore lattice

The standard approach to realize a spin-liquid state is through magnetically frustrated states, relying on ingredients such as the lattice geometry, dimensionality, and magnetic interaction type of the spins. While Heisenberg spins on a pyrochlore lattice with only antiferromagnetic nearest-neighbor interactions are theoretically proven disordered, spins in real systems generally include longer-range interactions. The spatial correlations at longer distances typically stabilize a long-range order rather than enhancing a spin-liquid state. Both states can, however, be destroyed by short-range static correlations introduced by chemical disorder. Here, using disorder-free specimens with a clear long-range antiferromagnetic order, we refine the spin structure of the Heisenberg spinel ZnFe 2 O 4 through neutron magnetic diffraction. The unique wave vector (1, 0, $\frac{1}{2}$) leads to a spin structure that can be viewed as alternatively stacked ferromagnetic and antiferromagnetic tetrahedra in a three-dimensional checkerboard form. Stable coexistence of these opposing types of clusters is enabled by the bipartite breathing pyrochlore crystal structure, leading to a second-order phase transition at 10 K. The diffraction intensity of ZnFe 2 O 4 is an exact complement to the inelastic scattering intensity of several chromate spinel systems which are regarded as model classical spin liquids. Our results challenge this attribution, and suggest instead of the six-spin ring mode, spin excitations in chromate spinels are closely related to the (1, 0, $\frac{1}{2}$) type of spin order and the four-spin ferromagnetic cluster locally at one tetrahedron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Novel Phases and Dynamics of Multi-orbital Mott Insulators

The goal of this project is to gain theoretical insight into magnetism in the presence of spin-orbit coupling (SOC) for multi-orbital Mott insulators. Inspired by the exact solution of the Kitaev model that harbors a quantum spin liquid with novel excitations and struck by its relevance for materials with anisotropic orbital interactions, the PI will explore various new classes of 4d and 5d transition metal oxides. Starting with all electron Hamiltonians, the PI will derive minimal magnetic models. The aim will be to understand the role played by orbital frustration, even in the absence of any geometric frustration, in creating orbitally ordered as well as spin-orbital liquid phases. Orbital frustration arises primarily from the directional or anisotropic nature of d-orbitals in contrast to the isotropic nature of the spin degree of freedom. These models will be investigated by a variety of theoretical and numerical methods, including exact diagonalization, mean field theories, density matrix renormalization group and quantum Monte Carlo methods. Testable predictions for three experiments: nuclear magnetic resonance (NMR) and resonant x-ray scattering (RXS) to probe orbital ordering and pump-probe experiments to probe quantum dynamics will test the validity of these models for materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1974

Ion heating in a plasma focus

Ion acceleration and heating in a plasma focus were investigated by the numerical integration of the three-dimensional equations of motion. The electric and magnetic fields given were derived from experimental data. The results obtained show that during the collapse phase of focus formation, ions are efficiently heated to temperatures of several keV. During the phase of rapid current reduction, ions are accelerated to large velocities in the axial direction. The results obtained with the model are in general agreement with experimental results.

Hohl, F.↗
NASA NTRS1981

The 34-meter antenna-subreflector translations to maximize RF gain

The extension of the 26 meter antenna to 34 meter diameter decreased the F/D ratio. This F/D change resulted in unacceptable gain losses due to the hyperboloid's lateral deflections. A three direction translating mounting device was added to the hyperboloid. This device was controlled by a microprocessor to minimize the offsets of the phase centers in the cassegrain RF system and also compensated for boresight directions. The use of the radiation program to predict the gain losses from displacements computed by a structural computing program using an analytical model of the 34 meter reflector structure is discussed. Field test results showed accurate predictions for the Y and Z hyperboloid translations. In the X direction, the prediction value was low. However, the computed gain losses vs primary foci offsets by the radiation program were verified by field tests.

Katow, M. S.↗
NASA NTRS1983

Heating measurements on Space Shuttle Orbiter models with differentially deflected elevons

The phase change paint technique was used to make heat transfer measurements on the windward wing/elevons area of Space Shuttle Orbiter models with differentially deflected elevons. Outboard elevons were deflected windward at angles from 0 to 20 deg. The inboard elevons were deflected leeward. The models were tested in air at Mach 6 and 10 with two flow conditions and three angles of attack for both Mach numbers. Heat transfer coefficient contour maps of the wing/elevon area were obtained for all test conditions. The maximum heat transfer coefficient on the outboard elevon and on certain areas of the wing that were obtained when the elevons were deflected are compared to values obtained when the elevons were undeflected.

Wells, W. L.↗
NASA NTRS1986

Mesoscale waves as a probe of Jupiter's deep atmosphere

Images from the Voyager north/south mapping sequences were searched for waves. A remarkable class of mesoscale waves was identified, with the following features: (1) the wavetrains are usually aligned zonally, i.e., wavecrests are north-south; (2) the average wavelength is 300 km with a standard deviation of only 20%; (3) the wavetrains are long; (4) the waves occur within 25 degrees of the equator, the bulk being at the equator itself; (5) the waves are centered at the extrema (in latitude) of the zonal flow; and (6) the meridional extent of the waves is typically 1 degree of latitude. These observations are interpreted as evidence of gravity waves propagating vertically within a leaky duct. A three-level model is assumed composed of a stable duct which extends up to the base of the NH3 cloud deck near 600 mb. Above this is a thin wave-trapping region characterized by a Richardson number Ri less than 1/4 and containing a critical level, where the local value of the zonal flow velocity equals the phase speed of the wave. This in turn is overlain by a stable region, representing the tropopause region and stratosphere.

Flaser, F. M.↗
NASA NTRS1987

Application of Advanced Very High Resolution Radiometer vegetation index to study atmosphere-biosphere exchange of CO2

Normalized difference vegetation indices derived from radiances measured by the Advanced Very High Resolution Radiometer were used to prescribe the phasing of terrestrial photosynthesis. The satellite data were combined with field data on soil respiration and a global map of net primary productivity to obtain the seasonal exchange of CO2 between the atmosphere and the terrestrial biosphere. The monthly fluxes of CO2 thus obtained were employed as source/sink functions in a global three-dimensional atmospheric tracer transport model to simulate the annual oscillations of CO2 in the atmosphere. The results demonstrate that satellite data of high spatial and temporal resolution can be used to provide quantitative information about seasonal and longer-term variations of photosynthetic activity on a global scale.

Fung, I. Y.↗