Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical modelling”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Advances in parameterization of cumulus convection in numerical models

Observational studies over the last decade have revealed new sources of difficulty in the already complex problem of parameterization of precipitating cumulus convection in large scale numerical weather prediction models. The structure of mesoscale updrafts and downdrafts and the difficulties they present for cumulus parameterization is reviewed. An approach is proposed which follows Kuo (1974) but includes cumulus and mesoscale downdrafts in the vertical distributions of heating and moistening. The development of a realistic life cycle behavior in the numerical prediction of a mesoscale convective complex will be shown.

Molinari, J.↗
NASA NTRS1992

A numerical model including PID control of a multizone crystal growth furnace

This paper presents a 2D axisymmetric combined conduction and radiation model of a multizone crystal growth furnace. The model is based on a programmable multizone furnace (PMZF) designed and built at NASA Lewis Research Center for growing high quality semiconductor crystals. A novel feature of this model is a control algorithm which automatically adjusts the power in any number of independently controlled heaters to establish the desired crystal temperatures in the furnace model. The control algorithm eliminates the need for numerous trial and error runs previously required to obtain the same results. The finite element code, FIDAP, used to develop the furnace model, was modified to directly incorporate the control algorithm. This algorithm, which presently uses PID control, and the associated heat transfer model are briefly discussed. Together, they have been used to predict the heater power distributions for a variety of furnace configurations and desired temperature profiles. Examples are included to demonstrate the effectiveness of the PID controlled model in establishing isothermal, Bridgman, and other complicated temperature profies in the sample. Finally, an example is given to show how the algorithm can be used to change the desired profile with time according to a prescribed temperature-time evolution.

Panzarella, Charles H.↗
NASA NTRS2019

Air Pollution Inputs to the Mojave Desert By Fusing Surface Mobile and Airborne in Situ and Airborne and Satellite Remote Sensing: A Case Study of Interbasin Transport With Numerical Model Validation

Deserts are fragile and highly sensitive ecosystems that increasingly are affected by upwind urban areas and industrial activities. The Los Angeles Basin (LAB) contributes to poor air quality in downwind deserts including the Mojave Desert. Few studies have investigated potential air pollution inputs to the Mojave, whose fragile ecosystem includes endangered plant and animal species. Data were collected on 19 August 2015 by a mobile air quality laboratory, AMOG (AutoMObile trace Gas) Surveyor, that observed inputs can arise from the LAB as well as the San Joaquin Valley (SJV), California. The campaign used a strong methane (CH4) plume as a tracer for the downwind fate of emissions from Bakersfield area petroleum production and also measured ozone (O3). Additional in situ concurrent airborne GHG and O3 data were collected by AJAX - Alpha Jet Atmospheric eXperiment. Both AMOG and AJAX measure winds. Mojave Desert air quality was very poor (visibility ~4 km). Based on the winds, an additional source was inferred beyond the LAB and SJV Basins. Numerical transport modeling and analysis of aerosol lidar data collected the same day by the Cloud Profiling LiDAR onboard the Earth Research-2 stratospheric airplane demonstrated that fires in Northern California were responsible, with prevailing winds transporting air southwards along the eastern Sierra Nevada Range (Bishop Valley) to the Mojave. Whereas the southern and eastern Mojave are impacted by SJV and LAB outflow, the north Mojave generally avoids these inputs. This study shows it can be affected by even distant wildfires, which likely will increase in occurrence and intensity from climate change. Thus, regulatory efforts to reduce air quality impacts on the endangered Mojave ecosystem must include wildfires and also account for the significant differences between different regions of the Mojave. Currently, there is a paucity of studies, highlighting the critical need for field research.

Ira Leifer↗
NASA NTRS1992

Numerical modeling of fluid and electromagnetic phenomena in an arcjet

An explicit numerical technique is used to solve the axisymmetric reduced electromagnetic field equation. The effect of an electrical arc on a viscous, axisymmetric flow is approximated using an implicit thin layer Navier-Stokes solver with additional electromagnetic source terms in conjunction with the explicit finite difference code.

Flowe, Anita C.↗
NASA NTRS1991

Two-dimensional numerical modeling and solution of convection heat transfer in turbulent He II

Numerical schemes are employed to investigate heat transfer in the turbulent flow of He II. FEM is used to solve a set of equations governing the heat transfer and hydrodynamics of He II in the turbulent regime. Numerical results are compared with available experimental data and interpreted in terms of conventional heat transfer parameters such as the Prandtl number, the Peclet number, and the Nusselt number. Within the prescribed Reynolds number domain, the Gorter-Mellink thermal counterflow mechanism becomes less significant, and He II acts like an ordinary fluid. The convection heat transfer characteristics of He II in the highly turbulent regime can be successfully described by using the conventional turbulence and heat transfer theories.

Zhang, Burt X.↗
NASA NTRS2017

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗
NASA NTRS1988

Numerical models of thermally and mechanically coupled two-layer convection of highly viscous fluids

Thermal convection in the earth mantle is investigated by means of numerical simulations. The mantle models comprise two horizontal layers of viscous incompressible fluid (with identical or differing properties) separated by a fixed horizontal interface; an isothermal, fixed-heat-flux, or insulating lower boundary; and heat supplied either from below or internally (in the lower layer only, equally in both layers, or primarily in the upper layer). The mathematical formulation of the models is explained, and results for linear stability and finite-amplitude convection are presented in extensive tables and graphs and discussed in detail. Particular attention is given to the presence of thermal coupling without interface distortion in many cases, the predominance of long-wavelength cells in the finite-amplitude models, and the large (50 percent) temperature difference across the interface in all cases.

Ellsworth, Kirk↗
NASA NTRS2011

Numerical Modeling of Cometary Meteoroid Streams Encountering Mars and Venus

We have simulated numerically the existence of meteoroid streams that encounter the orbits of Mars and Venus, potentially producing meteor showers at those planets. We find that 17 known comets can produce such showers, the intensity of which can be determined through observations. Six of these streams contain dense dust trails capable of producing meteor outbursts.

Christou, A. A.↗
NASA NTRS1981

Regional eddy vorticity transport and the equilibrium vorticity budgets of a numerical model ocean circulation

A mean vorticity budget analysis is presented of Holland's (1978) numerical ocean general circulation experiment. The stable budgets are compared with classical circulation theory to emphasize the ways in which the mesoscale motions of the model alter (or leave unaltered) classical vorticity balances. The basinwide meridional transports of vorticity by the mean flow and by the mesoscale flow in the mean are evaluated to establish the role(s) of the mesoscale in the larger scale equilibrium vorticity transports. The vorticity equation for this model fluid system is presented and the budget analysis method is described. Vorticity budgets over the selected regions and on a larger scale are given, and a summary of budget results is provided along with remarks about the utility of this type of analysis.

Harrison, D. E.↗
NASA NTRS1977

Comparisons among second-order, fourth-order, and pseudospectral advection techniques in a mesoscale numerical model

The reported investigation is concerned with mesoscale squall-line simulation utilizing numerical techniques. Particular attention is given to fourth-order results and their implications for the dynamical evolution of mesoscale squall-line systems which contain severe local storms. An adiabatic inviscid set of prognostic equations is utilized in a z coordinate system for the fundamental experiment. The fourth-order advection is employed for all time-dependent equations. The Euler-backward time marching scheme is utilized with a 60 second time step. The horizontal mesh length is 42 km. Horizontal diffusion is accomplished by utilizing a smoother-desmoother. Lateral boundary conditions are designed to reduce the development of strong gradients in dependent variables near the boundaries by bringing interior values from the grid to the boundaries.

Kaplan, M.↗
NASA NTRS1988

Numerical modeling of chemical vapor deposition (CVD) in a horizontal reactor

In the present numerical prediction of the deposition rate of silicon from silane in a CVD process, the conservation equations for mass, momentum, energy, and chemical species are solved on a staggered grid using the SIMPLE algorithm, while the rate of chemical reactions in the gas phase and on the susceptor surface is obtained from an Arrhenius rate equation. Predicted deposition rates as a function of position along the susceptor with and without the gas phase chemical reaction are compared with the available experimental and numerical data; agreement is excellent except at the leading edge of the susceptor, where the deposition rate is overpredicted.

Sheikholeslami, M. Z.↗
NASA NTRS1999

Numerical Modeling of Three-Dimensional Fluid Flow with Phase Change

We present a numerical method to compute phase change dynamics of three-dimensional deformable bubbles. The full Navier-Stokes and energy equations are solved for both phases by a front tracking/finite difference technique. The fluid boundary is explicitly tracked by discrete points that are connected by triangular elements to form a front that is used to keep the stratification of material properties sharp and to calculate the interfacial source terms. Two simulations are presented to show robustness of the method in handling complex phase boundaries. In the first case, growth of a vapor bubble in zero gravity is studied where large volume increase of the bubble is managed by adaptively increasing the front resolution. In the second case, growth of a bubble under high gravity is studied where indentation at the rear of the bubble results in a region of large curvature which challenges the front tracking in three dimensions.

Esmaeeli, Asghar↗
NASA NTRS1994

Numerical modeling of coronal mass ejections based on various pre-event model atmospheres

We examine how the initial state (pre-event corona) affects the numerical MHD simulation for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homogeneous density and temperature distribution at lower boundary (i.e. solar surface) have been used to analyze the role of streamer properties in determining the characteristics of loop-like transients. The present paper extends these studies to show how a broader class of global coronal properties leads not only to different types of CME's, but also modifies the adjacent quiet corona and/or coronal holes. We consider four pre-event coronal cases: (1) Constant boundary conditions and a polytropic gas with gamma = 1.05; (2) Non-constant (latitude dependent) boundary conditions and a polytropic gas with gamma = 1.05; (3) Constant boundary conditions with a volumetric energy source and gamma = 1.67; (4) Non-constant (latitude dependent) boundary conditions with a volumetric energy source and gamma = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions (streamers) and open field regions. The CME's initiation is simulated by introducing at the base of the corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values are determined using MHD characteristic theory. The simulations show how different CME's, including loop-like transients, clouds, and bright rays, might occur. There are significant new features in comparison to published results. We conclude that the pre-event corona is a crucial factor in dictating CME's properties.

Suess, S. T.↗
NASA NTRS1995

Numerical Modeling of Coronal Mass Ejections Based on Various Pre-event Model Atmospheres

We examine how the initial state (pre-event corona) affects the numerical MHD simulation for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homogeneous density and temperature distribution, at the lower boundary (i.e., solar surface) have been used to analyze the role of streamer properties in determining the characteristics of loop-like transients. The present paper extends these studies to show how a broader class of global coronal properties leads not only to different types of CME's, but also modifies the adjacent quiet corona and/or coronal holes. We consider four pre-event coronal cases: (1) constant boundary conditions and a polytropic gas with gamma = 1.05; (2) non-constant (latitude dependent) boundary conditions and a polytropic gas with gamma = 1.05; (3) constant boundary conditions with a volumetric energy source and gamma = 1.67; (4) non-constant (latitude dependent) boundary conditions with a volumetric energy source and gamma = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions (streamers) and open field regions. The CME's initiation is simulated by introducing at the base of the corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values are determined using magnetohydrodynamic (MHD) characteristic theory. The simulations show how different CME's, including loop-like transients, clouds and bright rays, might occur. There are significant new features in comparison to published results. We conclude that the pre-event corona is a crucial factor in dictating CME's properties.

Wang, A. H.↗
NASA NTRS2002

Simultaneous Laser Raman-rayleigh-lif Measurements and Numerical Modeling Results of a Lifted Turbulent H2/N2 Jet Flame in a Vitiated Coflow

An experiment and numerical investigation is presented of a lifted turbulent H2/N2 jet flame in a coflow of hot, vitiated gases. The vitiated coflow burner emulates the coupling of turbulent mixing and chemical kinetics exemplary of the reacting flow in the recirculation region of advanced combustors. It also simplifies numerical investigation of this coupled problem by removing the complexity of recirculating flow. Scalar measurements are reported for a lifted turbulent jet flame of H2/N2 (Re = 23,600, H/d = 10) in a coflow of hot combustion products from a lean H2/Air flame ((empty set) = 0.25, T = 1,045 K). The combination of Rayleigh scattering, Raman scattering, and laser-induced fluorescence is used to obtain simultaneous measurements of temperature and concentrations of the major species, OH, and NO. The data attest to the success of the experimental design in providing a uniform vitiated coflow throughout the entire test region. Two combustion models (PDF: joint scalar Probability Density Function and EDC: Eddy Dissipation Concept) are used in conjunction with various turbulence models to predict the lift-off height (H(sub PDF)/d = 7,H(sub EDC)/d = 8.5). Kalghatgi's classic phenomenological theory, which is based on scaling arguments, yields a reasonably accurate prediction (H(sub K)/d = 11.4) of the lift-off height for the present flame. The vitiated coflow admits the possibility of auto-ignition of mixed fluid, and the success of the present parabolic implementation of the PDF model in predicting a stable lifted flame is attributable to such ignition. The measurements indicate a thickened turbulent reaction zone at the flame base. Experimental results and numerical investigations support the plausibility of turbulent premixed flame propagation by small scale (on the order of the flame thickness) recirculation and mixing of hot products into reactants and subsequent rapid ignition of the mixture.

Cabra, R.↗
NASA NTRS1986

Numerical modeling of water-vapor transport during pre-storm and COHMEX

Initial conditions are designed for numerical simulation of mesocale processes in the atmosphere using the Limited Area Mesoscale Prediction System (LAMPS) model. These initial conditions represent an idealized baroclinic wave in which the transport of water vapor can be simulated. The constructed atmosphere has two homogeneous air masses, polar front, polar jet stream and a stratosphere. All these simulate the basic structure of the earth's atmosphere. The hydrostatic and geostrophic balances make it possible to evaluate mutually consistent fields of wind and of the height of isobaric surfaces.

Djuric, Dusan↗
NASA NTRS1987

A comparison of analytical and numerical models for steadily driven magnetic reconnection

The effect of boundary conditions in both analytical and numerical solutions of steadily driven reconnection on the reconnection of antiparallel magnetic fields is considered, with special attention given to the mathematical problem of properly specifying boundary conditions for the MHD equation. A unified formulation developed by Priest and Forbes (1986) for steady state reconnection, which includes the Petschek solution and the Sonnerup (1970) solution as special cases, is used to reinterpret the previous numerical experiments of steadily driven reconnection. It is shown that many contradictory features of these experiments were caused by the use of boundary conditions which are different from those required by Petschek's (1964) theory.

Forbes, T. G.↗
NASA NTRS2002

Simultaneous Raman-Rayleigh-LIF Measurements and Numerical Modeling Results of a Lifted H2/N2 Turbulent Jet Flame in a Vitiated Coflow

An experimental and numerical investigation is presented of a H2/N2 turbulent jet flame burner that has a novel vitiated coflow. The vitiated coflow emulates the recirculation region of most combustors, such as gas turbines or furnaces. Additionally, since the vitiated gases are coflowing, the burner allows for exploration of recirculation chemistry without the corresponding fluid mechanics of recirculation. Thus the vitiated coflow burner design facilitates the development of chemical kinetic combustion models without the added complexity of recirculation fluid mechanics. Scalar measurements are reported for a turbulent jet flame of H2/N2 in a coflow of combustion products from a lean ((empty set) = 0.25) H2/Air flame. The combination of laser-induced fluorescence, Rayleigh scattering, and Raman scattering is used to obtain simultaneous measurements of the temperature, major species, as well as OH and NO. Laminar flame calculation with equal diffusivity do agree when the premixing and preheating that occurs prior to flame stabilization is accounted for in the boundary conditions. Also presented is an exploratory pdf model that predicts the flame's axial profiles fairly well, but does not accurately predict the lift-off height.

Cabra, R.↗