Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase Prediction”

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 271 records · Page 15

Metastable sound speed in gas-liquid mixtures

A new method of calculating speed of sound for two-phase flow is presented. The new equation assumes no phase change during the propagation of an acoustic disturbance and assumes that only the total entropy of the mixture remains constant during the process. The new equation predicts single-phase values for the speed of sound in the limit of all gas or all liquid and agrees with available two-phase, air-water sound speed data. Other expressions used in the two-phase flow literature for calculating two-phase, metastable sound speed are reviewed and discussed. Comparisons are made between the new expression and several of the previous expressions -- most notably a triply isentropic equation as used, a triply isentropic equation as used, among others, by Karplus and by Wallis. Appropriate differences are pointed out and a thermodynamic criterion is derived which must be satisfied in order for the triply isentropic expression to be thermodynamically consistent. This criterion is not satisfied for the cases examined, which included two-phase nitrogen, air-water, two-phase parahydrogen, and steam-water. Consequently, the new equation derived is found to be superior to the other equations reviewed.

Bursik, J. W.↗

Phase Stability of a Powder Metallurgy Disk Superalloy

Advanced powder metallurgy superalloy disks in aerospace turbine engines now entering service can be exposed to temperatures approaching 700 C, higher than those previously encountered. They also have higher levels of refractory elements, which can increase mechanical properties at these temperatures but can also encourage phase instabilities during service. Microstructural changes including precipitation of topological close pack phase precipitation and coarsening of existing gamma' precipitates can be slow at these temperatures, yet potentially significant for anticipated disk service times exceeding 1,000 h. The ability to quantify and predict such potential phase instabilities and degradation of capabilities is needed to insure structural integrity and air worthiness of propulsion systems over the full life cycle. A prototypical advanced disk superalloy was subjected to high temperature exposures, and then evaluated. Microstructural changes and corresponding changes in mechanical properties were quantified. The results will be compared to predictions of microstructure modeling software.

Gabb, Timothy P.↗

Experimental Demonstration and System Analysis for Plasmonic Force Propulsion

This project developed a higher-fidelity model of a recently envisioned small spacecraft propulsion system for precision pointing and proximity control. Plasmonic force propulsion harnesses solar light focused onto plasmon reactive subwavelength nanostructures to accelerate and expel nanoparticle propellant via strong optical forces. The goal of the project was to show that plasmonic space propulsion can provide the level of proximity and attitude control envisioned for future NASA nano/picosatellite missions, a level that is better than state-of-the-art approaches. We achieved this goal by showing that plasmonic force thrusters are feasible for a range of advanced mission concepts requiring swarm formations in a deep space environment. We performed three case studies that evaluated the performance of the plasmonic force propulsion thruster in a deep space, microsatellite swarm formation. These case studies assumed the propulsion system could generate thrust at the level predicted from our Phase 1 study (1.6 μN). Through these cases we were able to analyze the concept within a mission specific context through detailed orbital dynamics calculations. Results indicate that, with the Phase 1 estimated thrust level, the approach is promising for providing attitude control to swarm formation spacecraft. Further, we achieved goals related to technology development. Specifically, we experimentally demonstrated nanoparticle acceleration due to plasmonic forces with asymmetric nanostructures excited by focused laser light. Additionally, we investigated the thrust sensitivity and nanoparticle propellant injection dependencies upon thermal effects. As a result of our study, plasmonic force propulsion is at an early TRL 3. Active research and design has been conducted analytically and in the laboratory. Furthermore, practical applications such as the three case studies have been identified for the scientific basic principles that were observed. Future efforts related to fundamental understanding of these techniques should focus on 1) developing a standalone array of asymmetric nanostructures that can effectively interact with a stream or reservoir of particles or 2) experimentally evaluate a dielectrophoretic injector for nanoparticle propellant. The main limitation discovered about plasmonic propulsion regards performance estimates significantly below the Phase 1 estimations. Specifically, original assumptions in the Phase 1 project (notably, a linear array of asymmetric nanostructures) is not a viable approach to achieving significant acceleration, high exhaust velocity, of nanoparticles. More specifically, we assumed in Phase 1 that nanoparticles would be accelerated in series by a long linear array of asymmetric nanostructures. That is, the acceleration of the nanoparticle would build and increase with the kick received by each subsequent nanostructure. This is fundamentally flawed. The potential profile of a single nanostructure is such that it prohibits this phenomenon. The potential energy associated with the plasmon-generated dielectrophoretic force is a potential well, which is good for trapping nanoparticles, but cannot provide significant acceleration of particles to expel them out and away from the nanostructure. Further, a nanoparticle expelled from the first nanostructure would need to overcome the potential barrier for entry into the next nanostructure accelerating stage. Fundamentally, this effect means that a linear array of nanostructures is not a viable accelerating structure. Correspondingly then, acceleration can, or should, only be provided by one nanostructure, and the net acceleration and thrust force of a single nanostructure is small (~cm/s exhaust velocities, sub-nN level thrust vs. the 100’s m/s, μN originally envisioned). While our experiments demonstrated acceleration and manipulation of a nanoparticle using laser light in aqueous environment, the achievable energy and momentum addition to the nanoparticle from a single nanostructure stage is too low for useful propulsion. In terms of thrust prediction, the estimated thrust of 1.6 μN in Phase 1 is reduced to a few nN of thrust with this new insight and understanding of the concept. This thrust level is too small to achieve attitude control of swarms as originally envisioned.

Plasmonic Force↗

A method of mapping compact structure in radio sources using VLBI observations

A two-part technique is described for determining the angular structure of the compact components of radio sources from VLBI observations. With this technique, the source structure is first approximated, using both amplitudes and closure phases, by a model consisting of circularly symmetric Gaussian components located on a grid of positions on the sky. The second part begins by employing this model to predict the visibility phase corresponding to each observed visibility amplitude; these estimated visibility phases are then adjusted to agree with the observed closure phases. The resulting estimates of the visibility phases and the observed visibility amplitudes are then combined in a direct Fourier transform to produce a 'dirty' source map that is deconvolved via the CLEAN procedure on the basis of the point-source response. Some examples based on data generated from test models are provided.

Cotton, W. D.↗

Parallel Grand Canonical Monte Carlo (ParaGrandMC) Simulation Code

This report provides an overview of the Parallel Grand Canonical Monte Carlo (ParaGrandMC) simulation code. This is a highly scalable parallel FORTRAN code for simulating the thermodynamic evolution of metal alloy systems at the atomic level, and predicting the thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces.

Vesselin I Yamakov↗

Semiannual variation of geomagnetic activity.

It is proposed that the semiannual variation in geomagnetic activity is caused by a semiannual variation in the effective southward component of the interplanetary field. The southward field arises because the interplanetary field is ordered in the solar equatorial coordinate system, whereas the interaction with the magnetosphere is controlled by a magnetospheric system. Several simple models utilizing this effective modulation of the southward component of the interplanetary field are examined. One of these closely predicts the observed phase and amplitude of the semiannual variation. This model assumes that northward interplanetary fields are noninteracting and that the interaction with southward fields is ordered in solar magnetospheric coordinates. The prediction of the diurnal variation of the strength of the interaction at the magnetopause by this model, does not, however, match the diurnal variation of geomagnetic activity as derived from ground-based data.

Russell, C. T.↗

Global Positioning System Navigation Above 76,000 km for NASA's Magnetospheric Multiscale Mission

NASA's Magnetospheric Multiscale (MMS) mission, launched in March of 2015, consists of a controlled formation of four spin-stabilized spacecraft in similar highly elliptic orbits reaching apogee at radial distances of 12 and 25 Earth radii (RE) in the first and second phases of the mission. Navigation for MMS is achieved independently on-board each spacecraft by processing Global Positioning System (GPS) observables using NASA Goddard Space Flight Center (GSFC)'s Navigator GPS receiver and the Goddard Enhanced Onboard Navigation System (GEONS) extended Kalman filter software. To our knowledge, MMS constitutes, by far, the highest-altitude operational use of GPS to date and represents a high point of over a decade of high-altitude GPS navigation research and development at GSFC. In this paper we will briefly describe past and ongoing high-altitude GPS research efforts at NASA GSFC and elsewhere, provide details on the design of the MMS GPS navigation system, and present on-orbit performance data from the first phase. We extrapolate these results to predict performance in the second phase orbit, and conclude with a discussion of the implications of the MMS results for future high-altitude GPS navigation, which we believe to be broad and far-reaching.

GPS↗

Global Positioning System Navigation Above 76,000 km for NASA's Magnetospheric Multiscale Mission

NASA's Magnetospheric Multiscale (MMS) mission, launched in March of 2015, consists of a controlled formation of four spin-stabilized spacecraft in similar highly elliptic orbits reaching apogee at radial distances of 12 and 25 Earth radii (RE) in the first and second phases of the mission. Navigation for MMSis achieved independently on-board each spacecraft by processing Global Positioning System (GPS) observables using NASA Goddard Space Flight Center (GSFC)'s Navigator GPS receiver and the Goddard Enhanced Onboard Navigation System (GEONS) extended Kalman filter software. To our knowledge, MMS constitutes, by far, the highest-altitude operational use of GPS to date and represents a high point of over a decade of high-altitude GPS navigation research and development at GSFC. In this paper we will briefly describe past and ongoing high-altitude GPS research efforts at NASA GSFC and elsewhere, provide details on the design of the MMS GPS navigation system, and present on-orbit performance data from the first phase. We extrapolate these results to predict performance in the second phase orbit, and conclude with a discussion of the implications of the MMS results for future high-altitude GPS navigation, which we believe to be broad and far-reaching.

navigation↗

A liquidus phase diagram for a primitive shergottite

To see if there is any relationship between primitive shergottites such as Eg and evolved shergottites such as Shergotty and Zagami, we performed one-bar experiments on the Eg composition. Broadly, our experimental results compare favorably with prediction. Our inferred phase diagram and comparison to Shergotty and Zagami melting experiments of Stolper and McSween are given. It does not appear possible to derive bulk Shergotty or Zagami by either equilibrium or fractional crystallization of Eg. However, if Shergotty and Zagami are cumulates, it may be possible to derive the inferred interstitial liquid from a composition such as Eg.

Jurewicz, A. J. G.↗

Calculation of Hub Loads at Low Airspeeds with Active Control

The effect of individual blade control (IBC) on the full-scale, low airspeed, level flight UH-60A oscillatory fixed system 4P hub loads and the rotating system blade bending moments is studied. The effect of a single 3P IBC input has been considered in this analytical study. At the low speed under consideration, it has been found that convergence of the comprehensive analysis is important for obtaining good predictions. Good correlation has been obtained with the measured full-scale wind tunnel data for the shapes of the fixed system 4P hub loads variations with the 3P IBC input phase, and also for the "best" phase of the 3P input (for minimum hub loads). The blade bending moment comparison shows mixed results. The 3P lead lag and the 4P flap bending moment trends with the 3P IBC input phase are reasonably predicted, whereas the 5P lead lag bending moment trend is not predicted well. Finally, the prediction of the baseline (no IBC) bending moments needs further study.

Kottapalli, Sesi↗

Parallel Grand-Canonical Monte Carlo (ParaGrandMC) User’s Manual Version 2.0

This manual describes the commands and command line options for the Parallel Grand Canonical Monte Carlo version 2.0 (ParaGrandMC.2.0) simulation code. This is a highly scalable parallel FORTRAN 2003 code for simulating the thermodynamic evolution of materials at the atomic level, and predicting their thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is specifically designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces. In addition to performing Monte Carlo (MC) simulations, the code can also perform Molecular Dynamics (MD) and Langevin Dynamics (LD) simulations, which can be combined and interchanged with MC for faster and more efficient system evolution. A detailed description of the MC part of the code is provided in the NASA ParaGrandMC report: NASA/CR–2016-219202; http://www.sti.nasa.gov.

High performance computing↗

CFD Modeling of Tank Pressurization and Axial Jet Mixing Experiments with and without Non-Condensable Gas in the Ullage

A two-phase CFD model for tank pressurization and following jet mixing of a cryogenic storage tank is presented using VOF approach for representing the phase boundary and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. The CFD model was validated against pressurization and liquid jet mixing data for a 110-inch diameter tank provided by Bullard1. Cases with like-gas and non-condensable gas pressurization are studied. The results of the cases with like-gas pressurization and following mixing are presented first, focusing on the effects of turbulence at the vapor-liquid interface on the tank pressure and phase change rates predictions. The second part of this paper is devoted to testing and validating the CFD model against non-condensable gas pressurization and following mixing. Tank pressures predicted in both cases are compared with each other and with the experimental data.

Computational Fluid Dynamics↗

CFD Modeling of Tank Pressurization and Axial Jet Mixing Experiments with and without Non-Condensable Gas in the Ullage

A two-phase CFD model for tank pressurization and following jet mixing of a cryogenic storage tank is presented using VOF approach for representing the phase boundary and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. The CFD model was validated against pressurization and liquid jet mixing data for a 110-inch diameter tank provided by Bullard. Cases with like-gas and non-condensable gas pressurization are studied. The results of the cases with like-gas pressurization and following mixing are presented first, focusing on the effects of turbulence at the vapor-liquid interface on the tank pressure and phase change rates predictions. The second part of this paper is devoted to testing and validating the CFD model against non-condensable gas pressurization and following mixing. Tank pressures predicted in both cases are compared with each other and with the experimental data.

Computational Fluid Dynamics↗

Numerical Predictions of Wind Turbine Power and Aerodynamic Loads for the NREL Phase II and IV Combined Experiment Rotor

Accurate, reliable and robust numerical predictions of wind turbine rotor power remain a challenge to the wind energy industry. The literature reports various methods that compare predictions to experiments. The methods vary from Blade Element Momentum Theory (BEM), Vortex Lattice (VL), to variants of Reynolds-averaged Navier-Stokes (RaNS). The BEM and VL methods consistently show discrepancies in predicting rotor power at higher wind speeds mainly due to inadequacies with inboard stall and stall delay models. The RaNS methodologies show promise in predicting blade stall. However, inaccurate rotor vortex wake convection, boundary layer turbulence modeling and grid resolution has limited their accuracy. In addition, the inherently unsteady stalled flow conditions become computationally expensive for even the best endowed research labs. Although numerical power predictions have been compared to experiment. The availability of good wind turbine data sufficient for code validation experimental data that has been extracted from the IEA Annex XIV download site for the NREL Combined Experiment phase II and phase IV rotor. In addition, the comparisons will show data that has been further reduced into steady wind and zero yaw conditions suitable for comparisons to "steady wind" rotor power predictions. In summary, the paper will present and discuss the capabilities and limitations of the three numerical methods and make available a database of experimental data suitable to help other numerical methods practitioners validate their own work.

Duque, Earl P. N.↗

A Computer Model for Analyzing Volatile Removal Assembly

A computer model simulates reactional gas/liquid two-phase flow processes in porous media. A typical process is the oxygen/wastewater flow in the Volatile Removal Assembly (VRA) in the Closed Environment Life Support System (CELSS) installed in the International Space Station (ISS). The volatile organics in the wastewater are combusted by oxygen gas to form clean water and carbon dioxide, which is solved in the water phase. The model predicts the oxygen gas concentration profile in the reactor, which is an indicator of reactor performance. In this innovation, a mathematical model is included in the computer model for calculating the mass transfer from the gas phase to the liquid phase. The amount of mass transfer depends on several factors, including gas-phase concentration, distribution, and reaction rate. For a given reactor dimension, these factors depend on pressure and temperature in the reactor and composition and flow rate of the influent.

Guo, Boyun↗

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗

LINAC Longitudinal Simulation and Measurement of Output Energy

The Fermilab Linac, a pivotal and historic accelerator at Fermilab, is crucial to the laboratory's operations, supplying a 400 MeV beam to various acceleration facilities. Due to daily variations in the Linac's output energy, precise monitoring and machine tuning are important to ensure the exiting energy meets the Booster's acceptance criteria. To address this, Beam Position Monitors (BPMs) are employed to assess and adjust the beam's energy, providing essential data on both the transverse position and longitudinal phase of the beam. We also aim to regulate the longitudinal phase profile to match our simulation predictions. By utilizing a Python-based simulation and beam data from three BPMs, we can determine the optimal phasing correction required for the final RF stage to achieve the desired energy. Currently, the calculations of longitudinal phase profile are based on simulations, so additional research is needed to verify if the predicted longitudinal phase distribution aligns with actual real-world data.

Safaryan, Milena↗

Numerical Formulation for the Prediction of Solid/Liquid Change of a Binary Alloy

A computational model is presented for the prediction of solid/liquid phase change energy transport including the influence of free convection fluid flow in the liquid phase region. The computational model considers the velocity components of all non-liquid phase change material control volumes to be zero but fully solves the coupled mass-momentum problem within the liquid region. The thermal energy model includes the entire domain and uses an enthalpy like model and a recently developed method for handling the phase change interface nonlinearity. Convergence studies are performed and comparisons made with experimental data for two different problem specifications. The convergence studies indicate that grid independence was achieved and the comparison with experimental data indicates excellent quantitative prediction of the melt fraction evolution. Qualitative data is also provided in the form of velocity vector diagrams and isotherm plots for selected times in the evolution of both problems. The computational costs incurred are quite low by comparison with previous efforts on solving these problems.

Schneider, G. E.↗