SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1992

A new numerical model of the middle atmosphere. I - Dynamics and transport of tropospheric source gases

Attention is given to a new model of the middle atmosphere which includes, in addition to the equations governing the zonal mean state, a potential vorticity equation for a single planetary-scale Rossby wave, and an IR radiative transfer code for the stratosphere and lower mesosphere, which replaces the Newtonian cooling parameterization used previously. It is shown that explicit computation of the planetary-scale wave field yields a more realistic representation of the zonal mean dynamics and the distribution of trace chemical species. Wave breaking produces a well-mixed 'surf zone' equatorward of the polar night vortex and drives a meridional circulation with downwelling on the poleward side of the vortex. This combination of mixing and downwelling produces shallow meridional gradients of trace gases in the subtropics and middle latitudes, and very steep gradients at the edge of the polar vortex. Computed distributions of methane and nitrous oxide are shown to agree well with observations.

Garcia, Rolando R.↗

NASA NTRS1993

Nonhydrostatic effects in numerical modeling of mesoscale convective systems and baroclinic waves

The present investigation is concerned with the role of convection upon mesoscale modeling results, particularly when the grid resolution becomes small enough that there is not a clear scale separation between the explicitly resolved circulations and the parameterized clouds. In those situations, the vertical accelerations in explicitly resolved circulations become strong enough that the hydrostatic assumption may no longer be valid. These concerns arise from interests in improving mesoscale modeling per se and in improving the subgrid-scale parameterizations in global models. The hydrostatic and the nonhydrostatic options of the Colorado State University Regional Atmospheric Modeling System were used to simulate dry gravity currents in two dimensions, using several different horizontal grid sizes. With horizontal grid intervals of 10 km or less, nonhydrostatic simulations produce wider and colder heads and weaker but wider forced updrafts than do the hydrostatic simulations. Comparing the hydrostatic and nonhydrostatic models show that the difference between the vertical mass fluxes is much less than the difference between the vertical velocities. When the grid is fine enough to resolve the head of the gravity current, horizontal convergence at the gust front extends upwards almost to the head of the cold air. Vertical mass flux in the forced updraft at the front varies with horizontal grid size mainly as a function of the height of the simulated head. For coarser grids, which do not resolve the head, vertical mass flux at all heights decreases with increasing horizontal grid size. A comparison on nonhydrostatic simulations with horizontal grid intervals of 1 km and 2 km illustrates how decreasing the grid size does not necessarily increase the intensity of the resolved circulation. The smaller grid enables the simulated gravity current to entrain a bubble of warm air behind the head, which results in a weaker circulation with a shorter head and weaker updraft.

Cohen, Charles↗

NASA NTRS1993

The role of global cloud climatologies in validating numerical models

The purpose of this work is to estimate sampling errors of area-time averaged rain rate due to temporal samplings by satellites. In particular, the sampling errors of the proposed low inclination orbit satellite of the Tropical Rainfall Measuring Mission (TRMM) (35 deg inclination and 350 km altitude), one of the sun synchronous polar orbiting satellites of NOAA series (98.89 deg inclination and 833 km altitude), and two simultaneous sun synchronous polar orbiting satellites--assumed to carry a perfect passive microwave sensor for direct rainfall measurements--will be estimated. This estimate is done by performing a study of the satellite orbits and the autocovariance function of the area-averaged rain rate time series. A model based on an exponential fit of the autocovariance function is used for actual calculations. Varying visiting intervals and total coverage of averaging area on each visit by the satellites are taken into account in the model. The data are generated by a General Circulation Model (GCM). The model has a diurnal cycle and parameterized convective processes. A special run of the GCM was made at NASA/GSFC in which the rainfall and precipitable water fields were retained globally for every hour of the run for the whole year.

HARSHVARDHAN↗

NASA NTRS1992

A numerical model of the large-scale solar wind in the outer heliosphere

We present a solar wind model developed primarily for the study of cosmic ray modulation in the heliosphere. The model is steady-state in the corotating frame, the wind speed is assumed to vary with heliomagnetic latitude, slower wind is taken near the current sheet. A wavy sheet will lead to dynamical effects: corotating interaction regions (CIR) and merged interaction regions (MIR) arise in a natural way. As a simplifying approximation, the wind speed is kept radial, so the calculation is basically 1-dimensional including plasma density, pressure, bulk speed, and the azimuthal component of the magnetic field. Calculations are carried out for various tilt angles of the neutral sheet illustrating various phases of the solar cycle. Possible cosmic-ray effects are discussed.

Kota, J.↗

NASA NTRS1991

Radar and microphysical characteristics of convective storms simulated from a numerical model using a new microphysical parameterization

The basic features of a new and improved bulk-microphysical parameterization capable of simulating the hydrometeor structure of convective systems in all types of large-scale environments (with minimal adjustment of coefficients) are studied. Reflectivities simulated from the model are compared with radar observations of an intense midlatitude convective system. Simulated reflectivities using the novel four-class ice scheme with a microphysical parameterization rain distribution at 105 min are illustrated. Preliminary results indicate that this new ice scheme works efficiently in simulating midlatitude continental storms.

Ferrier, Brad S.↗

NASA NTRS1993

Numerical modeling of crystal growth on a centrifuge for unstable natural convection configurations

The fluid mechanics associated with crystal growth processes on centrifuges is modeled using 2D and 3D models. Two-dimensional calculations show that flow bifurcations exist in such crystal growth configurations where the ampoule is oriented in the same direction as the resultant gravity vector and a temperature gradient is imposed on the melt. A scaling analysis is formulated to predict the flow transition point from the natural convection dominated regime to the Coriolis force dominated regime. Results of 3D calculations are presented for two thermal configurations of the crystal growth cell: top heated and bottom heated with respect to the centrifugal acceleration. In the top heated configuration, a substantial reduction in the convection intensity within the melt can be attained by centrifuge operations, and close to steady diffusion-limited thermal conditions can be achieved over a narrow range of the imposed microgravity level. In the bottom heated configuration the Coriolis force has a stabilizing effect on fluid motion by delaying the onset of unsteady convection.

Ramachandran, N.↗

NASA NTRS1993

Numerical modeling and passive thermal control of external lighting systems for Space Station Freedom

Consideration is given to three generic families of luminaries with lamp power ranging from 11 to 150 watts. A concept of an equivalent radiation node boundary temperature was used to impose worst hot and cold environments, and transient finite difference models were developed to study the effects of geometry and optical properties of thermal control coatings. Minimum and maximum transient temperatures were computed at the critical location during 90 minute orbit and were compared with allowable limits. Results show that with the proper choice of optical properties, the luminaries can be passively controlled to within acceptable limits.

Buck, Gregory A.↗

NASA NTRS1993

Numerical modeling of physical vapor transport under microgravity conditions: Effect of thermal creep and stress

One of the most promising applications of microgravity (micro-g) environments is the manufacture of exotic and high-quality crystals in closed cylindrical ampoules using physical vapor transport (PVT) processes. The quality enhancements are believed to be due to the absence of buoyant convection in the weightless environment - resulting in diffusion-limited transport of the vapor. In a typical experiment, solid-phase sample material is initially contained at one end of the ampoule. The sample is made to sublime into the vapor phase and deposit onto the opposite end by maintaining the source at an elevated temperature with respect to the deposit. Identification of the physical factors governing both the rates and uniformity of crystal growth, and the optimization of the micro-g technology, will require an accurate modeling of the vapor transport within the ampoule. Previous micro-g modeling efforts have approached the problem from a 'classical' convective/diffusion formulation, in which convection is driven by the action of buoyancy on thermal and solutal density differences. The general conclusion of these works have been that in low gravity environments the effect of buoyancy on vapor transport is negligible, and vapor transport occurs in a diffusion-limited mode. However, it has been recently recognized than in the non-isothermal (and often low total pressure) conditions encountered in ampoules, the commonly-assumed no-slip boundary condition to the differential equations governing fluid motion can be grossly unrepresentative of the actual situation. Specifically, the temperature gradients can give rise to thermal creep flows at the ampoule side walls. In addition, temperature gradients in the vapor itself can, through the action of thermal stress, lead to bulk fluid convection.

Mackowski, Daniel W.↗

NASA NTRS1994

Numerical modeling studies of wake vortex transport and evolution within the planetary boundary layer

The proposed research involves four tasks. The first of these is to simulate accurately the turbulent processes in the atmospheric boundary layer. TASS was originally developed to study meso-gamma scale phenomena, such as tornadic storms, microbursts and windshear effects in terminal areas. Simulation of wake vortex evolution, however, will rely on appropriate representation of the physical processes in the surface layer and mixed layer. This involves two parts. First, a specified heat flux boundary condition must be implemented at the surface. Using this boundary condition, simulation results will be compared to experimental data and to other model results for validation. At this point, any necessary changes to the model will be implemented. Next, a surface energy budget parameterization will be added to the model. This will enable calculation of the surface fluxes by accounting for the radiative heat transfer to and from the ground and heat loss to the soil rather than simple specification of the fluxes. The second task involves running TASS with prescribed wake vortices in the initial condition. The vortex models will be supplied by NASA Langley Research Center. Sensitivity tests will be performed on different meteorological environments in the atmospheric boundary layer, which include stable, neutral, and unstable stratifications, calm and severe wind conditions, and dry and wet conditions. Vortex strength may be varied as well. Relevant non-dimensional parameters will include the following: Richardson number or Froude number, Bowen ratio, and height to length scale ratios. The model output will be analyzed and visualized to better understand the transport, decay, and growth rates of the wake vortices. The third task involves running simulations using observed data. MIT Lincoln Labs is currently planning field experiments at the Memphis airport to measure both meteorological conditions and wake vortex characteristics. Once this data becomes available, it can be used to validate the model for vortex behavior under different atmospheric conditions. The fourth task will be to simulate the wake in a more realistic environment covering a wider area. This will involve grid nesting, since high resolution will be required in the wake region but a larger total domain will be used. During the first allocation year, most of the first task will be accomplished.

Lin, Yuh-Lang↗

NASA NTRS1992

Three-dimensional numerical modeling of indium phosphide Point-Contact Solar Cells

The Point-Contact Solar Cell (PCSC) geometry has proven very effective for silicon cells. To date, it has not been implemented in III-V materials. In addition, modeling such a geometry is very difficult because of its three-dimensional nature. We have developed a three-dimensional finite element modeling code (FIESTA ROC). In this paper, we present results from a three-dimensional modeling study of InP point-contact solar cells.

Clark, Ralph O.↗

NASA NTRS1994

Numerical model for the weakly nonlinear propagation of sound through turbulence

When finite amplitude (or intense) sound, such as a sonic boom, propagates through a turbulent atmosphere, the propagation is strongly affected by the turbulence. The interaction between sound and turbulence has mostly been studied as a linear phenomenon, i.e., the nonlinear behavior of the intense sound has been neglected. It has been shown that turbulence has an effect on the perceived loudness of sonic booms, mainly by changing its peak pressure and rise time. Peak pressure and rise time are important factors that determine the loudness of the sonic boom when heard outdoors. However, the interaction between turbulence and nonlinear effects has mostly not been included in propagation studies of sonic booms. It is therefore important to investigate the influence of acoustical nonlinearity on the interaction of intense sound with turbulence.

Lipkens, Bart↗

NASA NTRS1994

Quantitative numerical modeling of a coronal mass ejection

The results of an investigation on how the initial state (pre-event corona) quantitatively affects the magnetohydrodynamic model output of a coronal mass ejection (CME) are presented. The extension of the study on the identification of pre-event streamer properties which are crucial in determining the observed characteristics of CMEs, to more realistic configurations, also taking into account the coronal hole regions, adjacent to quiet coronal streamers, whose physical conditions are not adequately reproduced by global coronal models, is addressed. Results from the present model show three distinct types of CME's (loop, plasmoid and spike) which correspond to the classes of CME's identified in observations made with white light coronographs. The aim is a better understanding of the intercorrelations between streamers and adjacent coronal hole properties.

Wang, A. H.↗

NASA NTRS1994

A new numerical model of the middle atmosphere. 2: Ozone and related species

A new two-dimensional model with detailed photochemistry is presented. The model includes descriptions of planetary wave and gravity wave propagation and dissipation to characterize the wave forcing and associated mixing in the stratosphere and mesosphere. Such a representation allows for explicit calculation of the regions of strong mixing in the middle atmosphere required for accurate simulation of trace gas transport. The new model also includes a detailed description of photochemical processes in the stratosphere and mesosphere. The downward transport of H2, H2O, and NO(y) from the mesosphere to the stratosphere is examined, and it is shown that mesospheric processes can influence the distributions of these chemical species in polar regions. For HNO3 we also find that small concentrations of liquid aerosols above 30 km could play a major role in determining the abundance in polar winter at high latitudes. The model is also used to examine the chemical budget of ozone in the midlatitude stratosphere and to set constraints on the effectiveness of bromine relative to chlorine for ozone loss and the role of the HO2 + BrO reaction. Recent laboratory data used in this modeling study suggest that this process greatly enhances the effectiveness of bromine for ozone destruction, making bromine-catalyzed chemistry second only to HO(x)-catalyzed ozone destruction in the contemporary stratosphere at midlatitudes below about 18 km. The calculated vertical distribution of ozone in the lower stratosphere agrees well with observations, as does the total column ozone during most seasons and latitudes, with the important exception of southern hemisphere winter and spring.

Garcia, Rolando R.↗

NASA NTRS1995

A numerical model of nonmigrating diurnal tides between the surface and 65 km

Observations of surface pressure and middle atmosphere temperatures and winds indicate that a substantial nonmigrating component is present in the diurnal tide. The nonmigrating tides, which propagate with a zonal phase speed that is different from the earth's rotation, are attributed to the diurnal heating of geographically fixed sources. In this study we utilize a classical tidal model to examine the propagation characteristics of diurnal tides. The global fields of tropospheric sensible, radiative, and latent heating used to drive the model are supplied from summer and winter diurnal climatologies of the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM2). A novel aspect of this study is the focus on the relative importance of the nonmigrating components. The classical model successfully reproduces many observed features of the low-latitude diurnal surface pressure tides. In the middle atmosphere, the simulated migrating (or sun-synchronous) tide shows qualitative agreement with November-March Limb Infrared Monitor of the Stratosphere (LIMS) observations. Tropospheric solar heating is clearly the dominant driving force for the migrating tide, with secondary contributions from boundary-layer sensible heating and tropospheric latent heat release. The leading modes of the zonal mean tide are also driven chiefly by tropospheric solar heating. The higher-order modes of the zonal mean and eastward propagating tides may be attributed to the joint effects of tropospheric solar heating, sensible heating, and latent heat release. The LIMS and other data reveal features that cannot be explained or examined within the context of the classical model used in the present study. These include upward phase propagation, vertical attenuation, and temporal variations in the migrating diurnal tide.

Lieberman, Ruth S.↗

NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗

NASA NTRS1996

Numerical Model Studies of the Martian Mesoscale Circulations

Studies concerning mesoscale topographical effects on Martian flows examined low-level jets in the near equatorial latitudes and the dynamical intensification of flow by steep terrain. Continuation of work from previous years included evaluating the dissipation of cold air mass outbreaks due to enhanced sensible heat flux, further sensitivity and scaling evaluations for generalization of the characteristics of Martian mesoscale circulation caused by horizontal sensible heat-flux gradients, and evaluations of the significance that non-uniform surface would have on enhancing the polar CO2 ice sublimation during the spring. The sensitivity of maximum and minimum atmospheric temperatures to changes in wind speed, surface albedo, and deep soil temperature was investigated.

Segal, M.↗

NASA NTRS1997

Numerical Modeling for Combustion of Thermoplastic Materials in Microgravity

A time-dependent, three-dimensional model is under development to predict the temperature field, burning rate, and bubble bursting characteristics of burning thermoplastic materials in microgravity. Model results will be compared with experiments performed under microgravity and normal gravity conditions. The model will then be used to study the effects of variations in material properties and combustion conditions on burning rate and combustion behavior.

Butler, Kathryn M.↗

NASA NTRS1998

Numerical Modeling Studies of Wake Vortex Transport and Evolution Within the Planetary Boundary Layer

In support of the wake vortex effect of the Terminal Area Productivity program, we have put forward four tasks to be accomplished in our proposal. The first task is validation of two-dimensional wake vortex-turbulence interaction. The second task is investigation of three-dimensional interaction between wake vortices and atmospheric boundary layer (ABL) turbulence. The third task is ABL studies. The, fourth task is addition of a Klemp-Durran condition at the top boundary for TASS model. The accomplishment of these tasks will increase our understanding of the dynamics of wake vortex and improve forecasting systems responsible for air safety and efficiency. The first two tasks include following three parts: (a) Determine significant length scale for vortex decay and transport, especially the length scales associated with the onset of Crow instability (Crow, 1970); (b) Study the effects of atmospheric turbulence on the decay of the wake vortices; and (c) Determine the relationships between decay rate, transport properties and atmospheric parameters based on large eddy simulation (LES) results and the observational data. These parameters may include turbulence kinetic energy, dissipation rate, wind shear and atmospheric stratification. The ABL studies cover LES modeling of turbulence structure within planetary boundary layer under transition and stable stratification conditions. Evidences have shown that the turbulence in the stable boundary layer can be highly intermittent and the length scales of eddies are very small compared to those in convective case. We proposed to develop a nesting grid mesh scheme and a modified Klemp-Durran conditions (Klemp and Wilhelmson, 1978) at the top boundary for TASS model to simulate planetary boundary layer under stable stratification conditions. During the past year, our group has made great efforts to carry out the above mentioned four tasks simultaneously. The work accomplished in the last year will be described in the next section.

Lin, Yuh-Lang↗