Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field simulations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Adaptive-Grid Methods for Phase Field Models of Microstructure Development

In this work the authors show how the phase field model can be solved in a computationally efficient manner that opens a new large-scale simulational window on solidification physics. Our method uses a finite element, adaptive-grid formulation, and exploits the fact that the phase and temperature fields vary significantly only near the interface. We illustrate how our method allows efficient simulation of phase-field models in very large systems, and verify the predictions of solvability theory at intermediate undercooling. We then present new results at low undercoolings that suggest that solvability theory may not give the correct tip speed in that regime. We model solidification using the phase-field model used by Karma and Rappel.

Provatas, Nikolas↗
NASA NTRS2009

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗
NASA NTRS2012

Method for utilizing properties of the sinc(x) function for phase retrieval on nyquist-under-sampled data

Disclosed herein are systems, methods, and non-transitory computer-readable storage media for simulating propagation of an electromagnetic field, performing phase retrieval, or sampling a band-limited function. A system practicing the method generates transformed data using a discrete Fourier transform which samples a band-limited function f(x) without interpolating or modifying received data associated with the function f(x), wherein an interval between repeated copies in a periodic extension of the function f(x) obtained from the discrete Fourier transform is associated with a sampling ratio Q, defined as a ratio of a sampling frequency to a band-limited frequency, and wherein Q is assigned a value between 1 and 2 such that substantially no aliasing occurs in the transformed data, and retrieves a phase in the received data based on the transformed data, wherein the phase is used as feedback to an optical system.

Dean, Bruce H.↗
NASA NTRS1976

Summary of nozzle-exhaust plume flowfield analyses performed in support of space shuttle applications

Results are summarized of various analyses of the exhaust plumes of the space shuttle. Specific topics discussed included: the development of the two-phase plume flow field model, including finite rate chemistry and free molecular effects; analyses of exhaust plume simulation studies; verification of the analytical two-phase plume flow field model; and complete documentation of the two-phase plume computer code.

Smith, S. D.↗
NASA NTRS2003

Convection Effects in Three-dimensional Dendritic Growth

A phase-field model is developed to simulate free dendritic growth coupled with fluid flow for a pure material in three dimensions. The preliminary results presented here illustrate the strong influence of convection on the three-dimensional (3D) dendrite growth morphology. The detailed knowledge of the flow and temperature fields in the melt around the dendrite from the simulations allows for a detailed understanding of the convection effects on dendritic growth.

Lu, Yili↗
NASA NTRS2018

(GO)2-SIM: a GCM-Oriented Ground-Observation Forward-Simulator Framework for Objective Evaluation of Cloud and Precipitation Phase

General circulation model (GCM) evaluation using ground-based observations is complicated by inconsistencies in hydrometeor and phase definitions. Here we describe (GO)2-SIM, a forward simulator designed for objective hydrometeor-phase evaluation, and assess its performance over the North Slope of Alaska using a 1-year GCM simulation. For uncertainty assessment, 18 empirical relationships are used to convert model grid-average hydrometeor (liquid and ice, cloud, and precipitation) water contents to zenith polarimetric micropulse lidar and Ka-band Doppler radar measurements, producing an ensemble of 576 forward-simulation realizations. Sensor limitations are represented in forward space to objectively remove from consideration model grid cells with undetectable hydrometeor mixing ratios, some of which may correspond to numerical noise.Phase classification in forward space is complicated by the inability of sensors to measure ice and liquid signals distinctly. However, signatures exist in lidar–radar space such that thresholds on observables can be objectively estimated and related to hydrometeor phase. The proposed phase-classification technique leads to misclassification in fewer than 8% of hydrometeor-containing grid cells. Such misclassifications arise because, while the radar is capable of detecting mixed-phase conditions, it can mistake water- for ice-dominated layers. However, applying the same classification algorithm to forward-simulated and observed fields should generate hydrometeor-phase statistics with similar uncertainty. Alternatively, choosing to disregard how sensors define hydrometeor phase leads to frequency of occurrence discrepancies of up to 40%. So, while hydrometeor-phase maps determined in forward space are very different from model "reality" they capture the information sensors can provide and thereby enable objective model evaluation.

Lamer, K.↗
NASA NTRS2008

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗
NASA NTRS1971

Evaluation of acoustic testing techniques for spacecraft systems

External acoustic environments, structural responses, noise reductions, and the internal acoustic environments have been predicted for a typical shroud/spacecraft system during lift-off and various critical stages of flight. Spacecraft responses caused by energy transmission from the shroud via mechanical and acoustic paths have been compared and the importance of the mechanical path has been evaluated. Theoretical predictions have been compared extensively with available laboratory and in-flight measurements. Equivalent laboratory acoustic fields for simulation of shroud response during the various phases of flight have been derived and compared in detail. Techniques for varying the time-space correlations of laboratory acoustic fields have been examined, together with methods for varying the time and spatial distribution of acoustic amplitudes. Possible acoustic testing configurations for shroud/spacecraft systems have been suggested and trade-off considerations have been reviewed. The problem of simulating the acoustic environments versus simulating the structural responses has been considered and techniques for testing without the shroud installed have been discussed.

Cockburn, J. A.↗
NASA NTRS2023

Breaking Wave Field Statistics With A Multi-Layer Model

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The Λ(c) distribution, defined as the length of breaking fronts (per unit area) moving with speed c to c + dc following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the Λ(c) ∝ c−6 scaling without wind forcing for sufficiently steep wave fields. A scaling of Λ(c) based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Jiarong Wu↗
NASA NTRS1993

Models of Neptune's smooth recurrent radio emission

The quantitative response of the Planetary Radio Astronomy (PRA) instrument to a wave with polarization ellipse of arbitrary shape and orientation, arriving at the antennas from any direction, can be determined. This capability is used to model the time variation of intensity and circular polarization over a range of radio frequencies for proposed radio-source locations and emission characteristics at Neptune. At frequencies below 400 kHz the observed variation of intensity, polarization, and phase are closely simulated in an offset tilted dipole magnetic field by conjugate sources at midlatitude with filled emission cones. The phase of emission at higher frequencies is reproduced by sources at lower latitude. Modeled wide-cone emission does not reach the spacecraft at the observed phase nor have the polarization sense observed before closest approach. Source-surface maps of apparent polarization for the period before closest approach when instrumental response is especially sensitive to source location is presented. The method is capable of extension to more realistic models of the magnetic field.

Sawyer, Constance↗
NASA NTRS1997

Mars Global Surveyor: Aerobraking and Observations Support Using a Mars Global Circulation Model

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. Using a global atmospheric circulation model for Mars, the focus of this JRI has been to provide support for the Mars Global Surveyor (MGS) spacecraft aerobraking activities and interpretation guidance of preliminary observations. The primary atmospheric model applied in this investigation has been a high-top version of the NASA Ames Mars general circulation model (MGCM). Comparisons with an atmospheric model designed primarily for engineering purposes (Mars-GRAM) has also been carried out. From a suite of MGCM simulations, we have assessed plausible spatial and temporal variability in atmospheric density at high altitudes (e.g., 70-110 km) for seasonal dates and locations during Phase I aerobraking. Diagnostic tools have been developed to analyze circulation fields from the MGCM simulations, and these tools have been applied in the creation of a Mars climate catalogue database. Throughout Phase I aerobraking activities, analysis products have been provided to the MGS aerobraking atmospheric advisory group (AAG). Analyses of circulation variability at the coupling level between the MGCM and a Mars thermospheric global circulation model (MTGCM) has also been assessed. Finally, using a quasi-geostrophic dynamical formulation with the MGCM simulations, diagnosis of breaking planetary (Rossby) waves in Mars' middle atmosphere has been carried out. Titles of papers presented at scientific workshops and seminars, and a publication in the scientific literature are provided.

Hollingsworth, Jeffery L.↗
NASA NTRS1997

Mars Global Surveyor: Aerobraking and Observations Support Using a Mars Global Circulation Model

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. Using a global atmospheric circulation model for Mars, the focus of this JRI has been to provide support for the Mars Global Surveyor (MGS) spacecraft aerobraking activities and interpretation guidance of preliminary observations. ne primary atmospheric model applied in this investigation has been a high-top version of the NASA Ames Mars general circulation model (MGCM). Comparisons with an atmospheric model designed primarily for engineering purposes (Mars-GRAM) has also been carried out. From a suite of MGCM simulations, we have assessed plausible spatial and temporal variability in atmospheric density at high altitudes (e.g., 70-110 km) for seasonal dates and locations during Phase 1 aerobraking. Diagnostic tools have been developed to analyze circulation fields from the MGCM simulations, and these tools have been applied in the creation of a Mars climate catalogue database. Throughout Phase 1 aerobraking activities, analysis products have been provided to the MGS aerobraking atmospheric advisory group (AAG). Analyses of circulation variability at the coupling level between the MGCM and a Mars thermospheric global circulation model (MTGCM) has also been assessed. Finally, using a quasi-geostrophic dynamical formulation with the MGCM simulations, diagnosis of breaking planetary (Rossby) waves in Mars' middle atmosphere has been carried out. Titles of papers presented at scientific workshops and seminars, and a publication in the scientific literature are provided.

Hollingsworth, Jeffery L.↗
NASA NTRS1997

Mars Global Surveyor: Aerobraking and Observations Support Using a Mars Global Circulation Model

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. Using a global atmospheric circulation model for Mars, the focus of this JRI has been to provide support for the Mars Global Surveyor (MGS) spacecraft aerobraking activities and interpretation guidance of preliminary observations. The primary atmospheric model applied in this investigation has been a high-top version of the NASA Ames Mars general circulation model (MGCM). Comparisons with an atmospheric model designed primarily for engineering purposes (Mars-GRAM) has also been carried out. From a suite of MGCM simulations, we have assessed plausible spatial and temporal variability in atmospheric density at high altitudes (e.g., 70-1 10 km) for seasonal dates and locations during Phase 1 aerobraking. Diagnostic tools have been developed to analyze circulation fields from the MGCM simulations, and these tools have been applied in the creation of a Mars climate catalogue database. Throughout Phase I aerobraking activities, analysis products have been provided to the MGS aerobraking atmospheric advisory group (AAG). Analyses of circulation variability at the coupling level between the MGCM and a Mars thermospheric global circulation model (MTGCM) has also been assessed. Finally, using a quasi-geostrophic dynamical formulation with the MGCM simulations, diagnosis of breaking planetary (Rossby) waves in Mars middle atmosphere has been carried out. Titles of papers presented at scientific workshops and seminars, and a publication in the scientific literature are provided.

Hollingsworth, Jeffery L.↗
NASA NTRS2025

Simulations of Large-Scale Zero Boiloff, Densification, and Solidification of Hydrogen

An Integrated Refrigeration and Storage (IRAS) experimental system called the Ground Operations and Demonstration Unit for Liquid Hydrogen (GODU-LH2) demonstrated the ability store cryogens in a zero boiloff (ZBO) process, and to densify liquid hydrogen by reducing the temperature and pressure down to the triple point resulting in solidification. The incompressible Navier-Stokes was utilized to simulate ZBO, densification, and solidification of liquid hydrogen (LH2) in the IRAS tank. The simulations were performed using a commercially available Computational Fluid Dynamics (CFD) pressure-based mass and momentum flow model and an enthalpy-porosity energy model. Results demonstrated the simulation’s ability to predict time-dependent flow and temperature fields and solid-liquid phase locations for hydrogen during ZBO, densification, and solidification. The simulations showed good agreement with experimental data, with errors within acceptable ranges for temperature and pressure predictions, provided detailed insights into natural convection and solid-liquid phase change dynamics. These findings are critical for the design of future cryogenic fluid management systems.

Computational Fluid Dynamics↗
NASA NTRS2021

Preconditioning of overcast-to-broken cloud transitions by riming in marine cold air outbreaks

Marine cold air outbreaks (CAOs) commonly form overcast cloud decks that transition into broken cloud fields downwind, dramatically altering the local radiation budget. In this study, we investigate the impact of frozen hydrometeors on these transitions. We focus on a CAO case in the NW Atlantic, the location of the multi-year flight campaign ACTIVATE (Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment). We use MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, version 2) reanalysis fields to drive large eddy simulations with mixed-phase two-moment microphysics in a Lagrangian framework. We find that transitions are triggered by substantial rain (rainwater paths >25 g m−2), and only simulations that allow for aerosol depletion result in sustained breakups, as observed. Using a range of diagnostic ice nucleating particle concentrations, Ninp, we find that increasing ice progressively accelerates transitions, thus abbreviating the overcast state. Ice particles affect the cloud-topped boundary layer evolution, primarily through riming-related processes prior to substantial rain, leading to (1) a reduction in cloud liquid water, (2) early consumption of cloud condensation nuclei, and (3) early and light precipitation cooling and moistening below cloud. We refer to these three effects collectively as “preconditioning by riming”. Greater boundary layer aerosol concentrations available as cloud condensation nuclei (CCN) delay the onset of substantial rain. However, cloud breakup and low CCN concentration final stages are found to be inevitable in this case, due, primarily, to liquid water path buildup. An ice-modulated cloud transition speed suggests the possibility of a negative cloud–climate feedback. To address prevailing uncertainties in the model representation of mixed-phase processes, the magnitude of ice formation and riming impacts and, thereby, the strength of an associated negative cloud–climate feedback process, requires further observational evaluation by targeting riming hot spots with in situ imaging probes that allow for both the characterization of ice particles and abundance of supercooled droplets.

Marine cold air outbreaks↗
NASA NTRS1993

Three dimensional modeling of cirrus during the 1991 FIRE IFO 2: Detailed process study

A three-dimensional model of cirrus cloud formation and evolution, including microphysical, dynamical, and radiative processes, was used to simulate cirrus observed in the FIRE Phase 2 Cirrus field program (13 Nov. - 7 Dec. 1991). Sulfate aerosols, solution drops, ice crystals, and water vapor are all treated as interactive elements in the model. Ice crystal size distributions are fully resolved based on calculations of homogeneous freezing of solution drops, growth by water vapor deposition, evaporation, aggregation, and vertical transport. Visible and infrared radiative fluxes, and radiative heating rates are calculated using the two-stream algorithm described by Toon et al. Wind velocities, diffusion coefficients, and temperatures were taken from the MAPS analyses and the MM4 mesoscale model simulations. Within the model, moisture is transported and converted to liquid or vapor by the microphysical processes. The simulated cloud bulk and microphysical properties are shown in detail for the Nov. 26 and Dec. 5 case studies. Comparisons with lidar, radar, and in situ data are used to determine how well the simulations reproduced the observed cirrus. The roles played by various processes in the model are described in detail. The potential modes of nucleation are evaluated, and the importance of small-scale variations in temperature and humidity are discussed. The importance of competing ice crystal growth mechanisms (water vapor deposition and aggregation) are evaluated based on model simulations. Finally, the importance of ice crystal shape for crystal growth and vertical transport of ice are discussed.

Jensen, Eric J.↗