Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Primitive equations”

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

Numerical simulations of a transverse indirect circulation and low-level jet in the exit region of an upper-level jet

A numerical study was performed of a severe weather event (tornado) which occurred on May 10, 1973 in the Ohio region. The situation was modeled with a primitive equation mesoscale dynamic formulation. Account was taken of precipitation, the planetary boundary layer parameters as bulk quantities, the vertical pressure gradient, and lateral boundary conditions based on radiosonde data. Two 12-hr simulations, adiabatic and nondivergent, respectively, were analyzed for relationships between upper and lower level jets. In the adiabatic formulation, a transverse circulation with a low level jet formed at the exit region of the upper level jet. The nondivergent situation led to similar, but weaker, phenomena. Both forms suggest that indirect circulation in the exit zone of an upper level jet is strongly influenced by the initial structure of the jet.

Brill, K. F.↗
NASA NTRS1985

Zonally averaged model of dynamics, chemistry and radiation for the atmosphere

A nongeostrophic theory of zonally averaged circulation is formulated using the nonlinear primitive equations on a sphere, taking advantage of the more direct relationship between the mean meridional circulation and diabatic heating rate which is available in isentropic coordinates. Possible differences between results of nongeostrophic theory and the commonly used geostrophic formulation are discussed concerning: (1) the role of eddy forcing of the diabatic circulation, and (2) the nonlinear nearly inviscid limit vs the geostrophic limit. Problems associated with the traditional Rossby number scaling in quasi-geostrophic formulations are pointed out and an alternate, more general scaling based on the smallness of mean meridional to zonal velocities for a rotating planet is suggested. Such a scaling recovers the geostrophic balanced wind relationship for the mean zonal flow but reveals that the mean meridional velocity is in general ageostrophic.

Tung, K. K.↗
NASA NTRS1985

Quantitative non-destructive testing

The work undertaken during this period included two primary efforts. The first is a continuation of theoretical development from the previous year of models and data analyses for NDE using the Optical Thermal Infra-Red Measurement System (OPTITHIRMS) system, which involves heat injection with a laser and observation of the resulting thermal pattern with an infrared imaging system. The second is an investigation into the use of the thermoelastic effect as an effective tool for NDE. As in the past, the effort is aimed towards NDE techniques applicable to composite materials in structural applications. The theoretical development described produced several models of temperature patterns over several geometries and material types. Agreement between model data and temperature observations was obtained. A model study with one of these models investigated some fundamental difficulties with the proposed method (the primitive equation method) for obtaining diffusivity values in plates of thickness and supplied guidelines for avoiding these difficulties. A wide range of computing speeds was found among the various models, with a one-dimensional model based on Laplace's integral solution being both very fast and very accurate.

Welch, C. S.↗
NASA NTRS1985

Rossby wavetrains in the stratosphere forced by localised disturbances in the troposphere

Studies have linked elements of the stratospheric circulation with particular localized features in the tropospheric circulation. This suggests a study of the response of the stratosphere to forcing by localized disturbances in the troposphere. A multilevel, primitive equation model is used of the stratosphere and mesosphere the height of whose lower boundary at 300 mb can be prescribed. Localized height disturbances which grow to steady amplitude are applied at this lower boundary, and the response of the initially axially symmetric stratosphere is studied. The perturbation, centered at 45 deg N, has a Gaussian distribution with a half width of 15 deg., corresponding roughly to the size of persistent troughs and blocking ridges in the troposphere. The forcing is nearly at fully amplitude 10 days after being switched on, and thereafter remains steady. Two types of experiments are conducted: in one the forcing is of small amplitude (100 gpm) and in the other it is of large amplitude (600 gpm). These pairs of experiments are compared to determine how nonlinear processes affect the perturbation fields. This is done locally by defining the perturbation response to the forcing as an anomaly, i.e., as a departure from the response in a controlled experiment in which no asymmetric forcing was applied at the lower boundary of the model. Experiments have been conducted for a number of atmospheric states obtained as zonal means of observations made by a stratospheric sounding unit (SSU). In this summary, results for the zonal mean flow on January 19, 1982 are outlined.

Marks, C. J.↗
NASA NTRS1985

Mean zonal winds and planetary waves induced by internal gravity wave packets

The generation of mean zonal winds and planetary waves induced by internal gravity wave packets are investigated by numerical experiments. The numerical model is a 3-dimensional primitive equation model on the equatorial beta-plane with boundaries at the equator and North Pole. Because of equatorial boundary conditions, equatorial symmetric modes only are included in the model. Internal gravity wave packets are produced by thermal heating at the bottom boundary. Because of restriction of grid size, convective instability is parameterized by time-independent vertical and horizontal eddy viscosity. Other dissipation mechanisms (e.g., Newtonian cooling, ion drag, molecular viscosity etc.) are also included in the model.

Takahashi, M.↗
NASA NTRS1985

A comparison of two- and three-dimensional tracer transport within a stratospheric circulation model

Use of the residual circulation for stratospheric tracer transport has been compared to a fully three-dimensional calculation. The wind fields used in this study were obtained from a global, semispectral, primitive equation model, extending from 10 to 100 km in altitude. Comparisons were done with a passive tracer and an ozone-like substance over a two-month period corresponding to a Northern Hemisphere winter. It was found that the use of the residual circulation can lead to errors in the tracer concentrations of about a factor of 2. The error is made up of two components. One is fluctuating with a period of approximately one month and reflects directly the wave transience that occurs on that time-scale. The second part is increasing steadily over the integration period and results from an overestimate of the vertical transport by the residual circulation. Furthermore, the equatorward and upward mixing that occurs with transport by the three-dimensional circulation at low latitudes is not well reproduced when the residual circulation is used.

Schneider, H.-R.↗
NASA NTRS1986

Synoptic impacts of NOAA 6/Tiros N data on numerical circulation forecasts

Profiles of retrieved satellite observed temperatures provided to the Israel Meteorological Service (IMS) by the National Environmental Satellite Data and Information Service were used to enhance the specification of atmospheric conditions during three winter periods. Introduction of the satellite data was via three alternative schemes, and the enhanced specifications were used to initialize 48-hour numerical forecasts with the IMS five-layer primitive equation model; control forecasts were run in parallel from initial conditions without the satellite data. The study selects one forecast from each period for which discernible impacts on forecast synoptic fields can be attributed to the different initializations. The first two case studies examine the results of only one assimilation scheme and the third relates to all three schemes. Verification of the prognostics show the positive impact of including satellite temperatures in the initialization, but there were also examples of negative impact.

Druyan, L. M.↗
NASA NTRS1986

The application of the quasi-geostrophic Eliassen-Palm flux to the analysis of stratospheric data

The quasi-geostrophic form of the Eliassen-Palm flux is calculated for steady linear planetary wave 1 in a realistic zonal flow with dissipation. This flux shows a region of divergence in the polar stratosphere that is similar to that found in time averages of observations. Comparison with the full primitive equation form of the Eliassen-Palm flux indicates that in the model this divergent region is spurious and results from the overestimation of the momentum flux divergence when it is calculated from geostrophic winds. An approximate expression for the Eliassen-Palm flux in terms of balanced winds is suggested as an alternative to the quasi-geostrophic form.

Robinson, W. A.↗
NASA NTRS1987

Application of satellite data in variational analysis for global cyclonic systems

The research goal was a variational data assimilation method that incorporates as dynamical constraints, the primitive equations for a moist, convectively unstable atmosphere and the radiative transfer equation. Variables to be adjusted include the three-dimensional vector wind, height, temperature, and moisture from rawinsonde data, and cloud-wind vectors, moisture, and radiance from satellite data. This presents a formidable mathematical problem. In order to facilitate thorough analysis of each of the model components, four variational models that divide the problem naturally according to increasing complexity are defined. Each model is summarized.

Achtemeier, G. L.↗
NASA NTRS1987

The terminal area simulation system. Volume 1: Theoretical formulation

A three-dimensional numerical cloud model was developed for the general purpose of studying convective phenomena. The model utilizes a time splitting integration procedure in the numerical solution of the compressible nonhydrostatic primitive equations. Turbulence closure is achieved by a conventional first-order diagnostic approximation. Open lateral boundaries are incorporated which minimize wave reflection and which do not induce domain-wide mass trends. Microphysical processes are governed by prognostic equations for potential temperature water vapor, cloud droplets, ice crystals, rain, snow, and hail. Microphysical interactions are computed by numerous Orville-type parameterizations. A diagnostic surface boundary layer is parameterized assuming Monin-Obukhov similarity theory. The governing equation set is approximated on a staggered three-dimensional grid with quadratic-conservative central space differencing. Time differencing is approximated by the second-order Adams-Bashforth method. The vertical grid spacing may be either linear or stretched. The model domain may translate along with a convective cell, even at variable speeds.

Proctor, F. H.↗
NASA NTRS1987

Orbital atmospheric physics and dynamics

There are two ways of modeling the upper atmosphere. One is the empirical model that makes use of experimental data on means and excursions from the mean and fits the data in a self-consistent manner. The other approach is to deal directly with the physical processes. This is difficult since what is happening is extremely complex. Data measured using an interferometer to give Doppler shifts of airglow lines showed 300 to 800 m/sec winds with a complex structure in the upper region of the thermosphere at high latitudes. Ionospheric electric fields, strongly influenced by interaction with the solar wind, drive the ionized component and large neutral winds result due to momentum transfer between the charged particles and the neutrals. Frictional heating results from movement of ions through the neutrals, which also influences the compositional structure. These are examples of the complex interactions involved. The NCAR General Circulation Model (tropospheric) was adapted for use at thermospheric altitudes: the Thermospheric General Circulation Model (TGCM). The model makes use partly of primitive equations and partly of empirical data for some quantities such as electron density, magnetic field, and ion drift.

Roble, Raymond↗
NASA NTRS1987

A study of the adequacy of quasi-geostrophic dynamics for modeling the effect of frontal cyclones on the larger scale flow

The evolution of individual cyclone waves is studied in order to see how well quasi-geostrophic (QG) dynamics can simulate the behavior of primitive equations (PE) dynamics. This work is an extension of a similar study (Mudrick, 1982); emphasis is placed here on adding a frontal zone and other more diverse features to the basic states used. In addition, sets of PE integrations, with and without friction, are used to study the formation of surface occluded fronts within the evolving cyclones. Results of the study are summarized at the beginning of the report.

Mudrick, Stephen↗
NASA NTRS1986

Numerical studies of frontal dynamics

Efforts concentrated on the development of a two dimensional primitive equation (PE) model of frontogenesis that simultaneously incorporates the frontagenetical mechanisms of confluence and horizontal shear. Applying this model to study the effects of upper level frontogenesis, it appeared to be dominated by tilting effects associated with cross front variation of vertical motion, in which subsidence is maximized within and to the warm side of the frontal zone. Results suggest that aspects characteristic of three-dimensional baroclinic waves may be abstracted to a significant extent in a two dimensional framework. They also show that upper-level frontogenesis and tropopause folding can occur in the absence of three-dimensional curvature effects, commonly believed to be necessary for realistic upper-level frontogenesis. An implication of the dominant tilting effects is that they may have to be adequately resolved by numerical weather prediction models, thus requiring better horizontal and vertical resolution.

Keyser, Daniel↗
NASA NTRS1985

A fully implicit scheme for global numerical weather prediction

A fast-slow factored scheme is presented for use with shallow-water primitive equation numerical weather prediction models. The technique was developed to reduce the rotational mode errors which arise when the fast and slow terms of the governing differential equations are treated simultaneously. The method factors out the fast and slow terms along the coordinate directions by means of a modified Crank-Nicolson scheme. A finite-difference spatial discretization is carried out in the zonal and meridional directions to reduce the factorization error to near-zero, and that time steps of 60-90 min can be used to obtain acceptably accurate results, even in the presence of fine spatial structures in the flow.

Augenbaum, J. M.↗
NASA NTRS1985

Sensitivity of a mesoscale model to initial specification of relative humidity, liquid water and vertical motion

The influence of synoptic scale initial conditions on the accuracy of mesoscale precipitation modeling is investigated. Attention is focused on the relative importance of the water vapor, cloud water, rain water, and vertical motion, with the analysis carried out using the Limited Area Mesoscale Prediction System (LAMPS). The fully moist primitive equation model has 15 levels and a terrain-following sigma coordinate system. A K-theory approach was implemented to model the planetary boundary layer. A total of 15 sensitivity simulations were run to investigate the effects of the synoptic initial conditions of the four atmospheric variables. The absence of synoptic cloud and rain water amounts in the initialization caused a 2 hr delay in the onset of precipitation. The delay was increased if synoptic-scale vertical motion was used instead of mesoscale values. Both the delays and a choice of a smoothed moisture field resulted in underestimations of the total rainfall.

Kalb, M. W.↗
NASA NTRS1986

Evaluation of a multivariate variational assimilation of conventional and satellite data for the diagnosis of cyclone systems

A variational data assimilation method for the study of cyclone-scale weather systems is described. The variational data assimilation method is to incorporate primitive equations for a moist, convectively unstable atmosphere and the radiative transfer equation. The variables to be adjusted include the three-dimensional vector wind, height, temperature, and moisture from rawinsonde data, and cloud-wind vectors, moisture, and radiance from satellite data. The development of variational model 1 which contains two nonlinear horizontal momentum equations, an integrated continuity equation, and a hydrostatic equation is examined. Examples applying the assimilation model to rawinsonde and satellite data are presented.

Achtemeier, Gary L.↗
NASA NTRS1986

Transport of atmospheric tracers by planetary waves during a winter stratospheric warming event - A three-dimensional model simulation

In a simple three-dimensional primitive equation model, a wave number 1 major stratospheric warming is simulated. With the aid of two idealized tracers it is shown that the transport during a major warming event is characterized by a small, well-organized tongue of subtropical air flowing around the displaced winter vortex into the polar cap and by a wide area with strong quasi-horizontal mixing (surf zone). The description of these dynamical processes requires a full three-dimensional space resolution.

Kouker, W.↗
NASA NTRS1986

Nongeostrophic theory of zonally averaged circulation. I - Formulation

A nongeostrophic theory of zonally averaged circulation is formulated using the nonlinear primitive equations (mass conservation, thermodynamics, and zonal momentum) on a sphere. The relationship between the mean meridional circulation and diabatic heating rate is studied. Differences between results of nongeostropic theory and the geostrophic formulation concerning the role of eddy forcing of the diabatic circulation and the nonlinear nearly inviscid limit versus the geostrophic limit are discussed. Consideration is given to the Eliassen-Palm flux divergence, the Eliassen-Palm pseudodivergence, the nonacceleration theorem, and the nonlinear nongeostrophic Taylor relationship.

Tung, Ka Kit↗