Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field model”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗
NASA NTRS1990

Gravity field determination for Mars Observer

Mars Observer will be the first near circular, low altitude and short-periodic orbiter of Mars. From previous Mars orbiters, such as Mariner 9 and Vikings 1 and 2, a large gravitational oblateness and in general a gravity field ten times stronger than the earth's field have been determined. Because of these, the gravity field will dominate the evolution of the spacecraft's orbital motion and is the principal error source for reconstructing the spacecraft's motion. Thus a gravity calibration (GC) period has been established prior to the mapping phase in order to refine the gravitational field model. This paper describes the plan being made by Navigation to implement the GC strategy. By analyzing simulated Doppler data acquired during the seven day GC period, we give (1) the expected improvement in the gravity field model and (2) the resulting improvement in reconstructed spacecraft orbital motion. In addition, we have developed a strategy for converging the GC gravity model.

Esposito, Pasquale B.↗
NASA NTRS1998

Current Status on the use of Parallel Computing in Turbulent Reacting Flow Computations Involving Sprays, Monte Carlo PDF and Unstructured Grids

The state of the art in multidimensional combustor modeling as evidenced by the level of sophistication employed in terms of modeling and numerical accuracy considerations, is also dictated by the available computer memory and turnaround times afforded by present-day computers. With the aim of advancing the current multi-dimensional computational tools used in the design of advanced technology combustors, a solution procedure is developed that combines the novelty of the coupled CFD/spray/scalar Monte Carlo PDF (Probability Density Function) computations on unstructured grids with the ability to run on parallel architectures. In this approach, the mean gas-phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The gas-turbine combustor flows are often characterized by a complex interaction between various physical processes associated with the interaction between the liquid and gas phases, droplet vaporization, turbulent mixing, heat release associated with chemical kinetics, radiative heat transfer associated with highly absorbing and radiating species, among others. The rate controlling processes often interact with each other at various disparate time 1 and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and liquid phase evaporation in many practical combustion devices.

Raju, M. S.↗
NASA NTRS1982

Superalloy composition modeling

Development of a predictive method for determination of the gamma/gamma prime phase fields, i.e., gamma prime volume fraction as a function of the multicomponent composition, is described. The cluster variation method used for binary alloys in which the precipitated phase is coherent with the matrix phase is extended for application to the multicomponent coherent gamma/gamma prime nickel-based superalloys. It is shown that the cluster variation method can accurately describe the equilibrium (incoherent) gamma/gamma prime phase fields in the binary Ni-Al phase diagram. The gamma/gamma prime phase field for the Ni-Cr-Al ternary phase diagram is computed as a function of temperature. A reasonable fit results between the calculated and the experimental diagrams. The modeling of the six-component Ni-Cr-Al-Co-Mo-Ti base superalloy and the effect of Ni substitution of Co are discussed.

Barefoot, J.↗
NASA NTRS1994

Hybird state of the tail mangetic configuration during steady convection events

Previous observations have shown that during periods of steady magnetospheric convection (SMC) a large amount of magnetic flux crosses the plasma sheet (corresponding to approximately 10 deg wide auroral oval at the nightside) and that the magnetic configuration in the midtail is relaxed (the curent sheet is thick and contains enhanced B(sub Z). These signatures are typical for the substorm recovery phase. Using near-geostationary magnetic field data, magnetic field modeling and a noval diagostic technique (isotropic boundary algorithm), we show that in the near-Earth tail the magnetic confirguration is very stretched during the SMC events. This stretching is caused by an intense, thin westward current. Because of the srongly depressed B(sub Z), there is a large radial gradient in the near-tail magetic field. These signatures have been peviously associated only with the substorm growth phase. Our results indicate that during the SMC periods the magnetic configuration is very peculiar, with co-existing thin near-Earth current sheet and thick midtail plasma sheet. The deep local minimum of the equatorial B(sub Z) that devleops at R approximately 12 R(sub E) is consistent with steady, adiabatic, Earthward convection in the midtail. These results impose contraints on the existing substorm theories, and call for an explanation of how such a stressed configuration can persist for such a long time without tail current disruptions that occur at the end of a substorm growth phase.

Sergeev, V. A.↗
NASA NTRS1985

A study of backscattering and emission from closely packed inhomogeneous media

The effects of close spacing between small scattering spheres were examined by keeping the distance-dependent terms in the expressions for the transverse scattered fields. The phase matrix was then derived from these fields and was used in the radiative transfer formulation to model scattering and emission from a densely populated, inhomogeneous layer. Computed results were compared with those obtained when the phase matrix was specialized to the far-field condition. It was found that the use of the far-zone condition tended to underestimate both the level of the copolarized backscattering and the cross-polarized backscattering. In emission computations, the use of the far-zone condition overestimated the level of the brightness temperature. These effects decreased with a decrease in the volume fraction or an increase in the exploring frequency, as expected. An improvement on the snow parameter (density and crystal size) estimation was shown to be possible when this new phase matrix was used.

Fung, A. K.↗
NASA NTRS1994

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
NASA NTRS1995

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
NASA NTRS1996

Photon Number-Phase Uncertainty Relation in the Evolution of the Field in a Kerr-Like Medium

A model of a single-mode field, initially prepared in a coherent state, coupled to a two-level atom surrounded by a nonlinear Kerr-like medium contained inside a very good quality cavity is considered. We derive the photon number-phase uncertainty relation in the evolution of the field for a weak and strong nonlinear coupling respectively, within the Hermitian phase operator formalism of Pegg and Barnett, and discuss the effects of nonlinear coupling of the Kerr-like medium on photon number-phase uncertainty relation of the field.

Fan, An-Fu↗
NASA NTRS1985

Extension of Ostwald Ripening Theory

The objective is to develop models based on the mean field approximation of Ostwald ripening to describe the growth of second phase droplets or crystallites. The models will include time variations in nucleation rate, control of saturation through addition of solute, precipitating agents, changes in temperature, and various surface kinetic effects. Numerical integration schemes have been developed and tested against the asymptotic solution of Liftshitz, Slyozov and Wagner (LSW). A second attractor (in addition to the LSW distribution) has been found and, contrary to the LSW theory, the final distribution is dependent on the initial distribution. A series of microgravity experiments is being planned to test this and other results from this work.

Baird, J.↗
NASA NTRS1983

The structure of the heliospheric current sheet - 1978-1982

Continuing the study near solar minimum in 1976-1977 (Hoeksema et al., 1982), the configuration of the heliospheric magnetic field for the period 1976-1982 is calculated using a potential field model. Particular attention is given to the structure during the rising phase, maximum, and early decline of sunspot cycle 21, from 1978 to 1982. Four warps in the current sheet (the boundary between interplanetary magnetic field toward and away from the sun) are seen early in this interval; these give rise to a four-sector structure in the interplanetary magnetic field observed at earth. It is noted that the location of the current sheet changes slowly and extends to a heliographic latitude of approximately 50 deg. The strength of the polar field correction throughout this period is determined and included in the model calculations. The lower latitude magnetic fields become much stronger as the polar fields weaken and reverse polarity near maximum, decreasing the effect of the polar field correction.

Hoeksema, J. T.↗
NASA NTRS2014

Multi-Scale Computational Modeling of Two-Phased Metal Using GMC Method

A multi-scale computational model for determining plastic behavior in two-phased CMSX-4 Ni-based superalloys is developed on a finite element analysis (FEA) framework employing crystal plasticity constitutive model that can capture the microstructural scale stress field. The generalized method of cells (GMC) micromechanics model is used for homogenizing the local field quantities. At first, GMC as stand-alone is validated by analyzing a repeating unit cell (RUC) as a two-phased sample with 72.9% volume fraction of gamma'-precipitate in the gamma-matrix phase and comparing the results with those predicted by finite element analysis (FEA) models incorporating the same crystal plasticity constitutive model. The global stress-strain behavior and the local field quantity distributions predicted by GMC demonstrated good agreement with FEA. High computational saving, at the expense of some accuracy in the components of local tensor field quantities, was obtained with GMC. Finally, the capability of the developed multi-scale model linking FEA and GMC to solve real life sized structures is demonstrated by analyzing an engine disc component and determining the microstructural scale details of the field quantities.

nickel alloys↗
NASA NTRS2002

A model of hollow cathode plasma chemistry

We have developed a new model of hollow cathode plasma chemistry based on the observation that xenon ion mobility is diffusion limited due to resonant charge exchange reactions. The model shows that vapor phase barium atoms are ionized almost immediately and electric fields accelerate the ions upstream from the emission zone. We have also applied the model to the orifice region, where the resultant ion generation profile correlates with previously reported orifice erosion.

hollow cathode↗
NASA NTRS1980

Magnetic variations during substorms

Extensive observations of magnetic field variations during auroral substorms have been used to construct phenomenological models of the substorm process. In one model, three substorm phases, growth, expansion and recovery correspond respectively to energy storage, energy release and return to equilibrium. New data and more detailed studies show the simple model required modification. Consideration is given to a number of recent studies of substorm magnetic field variations which illustrate the limitations of the simple model. These include the existence of several substorm expansions following a single growth phase; the frequent occurrence of multiple onsets of a single substorm; plasma sheet thinning prior to each onset; correlations between northward fluctuations of the IMF and expansion onset. These observations are incorporated into a revised phenomenological model.

Mcpherron, R. L.↗
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 NTRS1977

A design handbook for phase change thermal control and energy storage devices

Comprehensive survey is given of the thermal aspects of phase change material devices. Fundamental mechanisms of heat transfer within the phase change device are discussed. Performance in zero-g and one-g fields are examined as it relates to such a device. Computer models for phase change materials, with metal fillers, undergoing conductive and convective processes are detailed. Using these models, extensive parametric data are presented for a hypothetical configuration with a rectangular phase change housing, using straight fins as the filler, and paraffin as the phase change material. These data are generated over a range of realistic sizes, material properties, and thermal boundary conditions. A number of illustrative examples are given to demonstrate use of the parametric data. Also, a complete listing of phase change material property data are reproduced herein as an aid to the reader.

Humphries, W. R.↗
NASA NTRS1983

The structure of the heliospheric current sheet, 1978 - 1982

The structure of the heliospheric magnetic field changes substantially during the 11 year sunspot cycle. Its configuration for the period 1976 through 1982 using a potential field model was calculated. The structure during the rising phase, maximum, and early decline of sunspot cycle 21, from 1978 to 1982 is considered.

Hoeksema, J. T.↗