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

Diurnal forcing of planetary atmospheres

The utility of the Mars Planetary Boundary Layer Model (MPBL) for calculations in support of the Mars 94 balloon mission was substantially enhanced by the introduction of a balloon equation of motion into the model. Both vertical and horizontal excursions of the balloon are calculated along with its volume, temperature, and pressure. The simulations reproduce the expected 5-min vertical oscillations of a constant density balloon at altitude on Mars. The results of these calculations are presented for the nominal target location of the balloon. A nonlinear balanced model was developed for the Martian atmosphere. It was used to initialize a primitive equation model for the simulations of the Earth's atmosphere at the time of the El Chichon eruption in 1982. It is also used as an assimilation model to update the temperature and wind fields at frequent intervals.

Houben, Howard C.↗
NASA NTRS1992

Results of a zonally truncated three-dimensional model of the Venus middle atmosphere

Although the equatorial rotational speed of the solid surface of Venus is only 4 m s(exp-1), the atmospheric rotational speed reaches a maximum of approximately 100 m s(exp-1) near the equatorial cloud top level (65 to 70 km). This phenomenon, known as superrotation, is the central dynamical problem of the Venus atmosphere. We report here the results of numerical simulations aimed at clarifying the mechanism for maintaining the equatorial cloud top rotation. Maintenance of an equatorial rotational speed maximum above the surface requires waves or eddies that systematically transport angular momentum against its zonal mean gradient. The zonally symmetric Hadley circulation is driven thermally and acts to reduce the rotational speed at the equatorial cloud top level; thus wave or eddy transport must counter this tendency as well as friction. Planetary waves arising from horizontal shear instability of the zonal flow (barotropic instability) could maintain the equatorial rotation by transporting angular momentum horizontally from midlatitudes toward the equator. Alternatively, vertically propagating waves could provide the required momentum source. The relative motion between the rotating atmosphere and the pattern of solar heating, which as a maximum where solar radiation is absorbed near the cloud tops, drives diurnal and semidiurnal thermal tides that propagate vertically away from the cloud top level. The effect of this wave propagation is to transport momentum toward the cloud top level at low latitudes and accelerate the mean zonal flow there. We employ a semispectral primitive equation model with a zonal mean flow and zonal wavenumbers 1 and 2. These waves correspond to the diurnal and semidiurnal tides, but they can also be excited by barotropic or baroclinic instability. Waves of higher wavenumbers and interactions between the waves are neglected. Symmetry about the equator is assumed, so the model applies to one hemisphere and covers the altitude range 30 to 110 km. Horizontal resolution is 1.5 deg latitude, and vertical resolution is 1.5 km. Solar and thermal infrared heating, based on Venus observations and calculations drive the model flow. Dissipation is accomplished mainly by Rayleigh friction, chosen to produce strong dissipation above 85 km in order to absorb upward propagating waves and limit extreme flow velocities there, yet to give very weak Rayleigh friction below 70 km; results in the cloud layer do not appear to be sensitive to the Rayleigh friction. The model also has weak vertical diffusion, and very weak horizontal diffusion, which has a smoothing effect on the flow only at the two grid points nearest the pole.

Newman, M.↗
NASA NTRS1983

Theoretical study of multiple equilibria in simple axisymmetric tropical circulations

The possibility that the asymmetric part of the atmospheric circulation can possess multiple equilibrium states is examined using a two-layer axisymmetric model involving balance equations on an equatorial beta plane. Mountains are excluded from consideration and a Newtonian cooling formulation represents thermal forcing. A temperature maximum at 25 deg N is selected to simulate summer conditions in the Northern Hemisphere. Steady-state solutions obtained are investigated for stability with regard to first and second y-mode perturbations. A single stable mode is found, together with two other quasi-stable states. Attention is given to numerically modeling multiple equilibria in symmetric circulations, and one steady-state is determined for the two-layer model. A model employing primitive equations with the Boussinesq approximation is also examined, and it also furnishes only one steady state. The reasons for the lack of multiple steady-states as derived by the models are discussed.

Goswami, B. N.↗
NASA NTRS1987

Finite element techniques for the Navier-Stokes equations in the primitive variable formulation and the vorticity stream-function formulation

Finite element procedures for the Navier-Stokes equations in the primitive variable formulation and the vorticity stream-function formulation have been implemented. For both formulations, streamline-upwind/Petrov-Galerkin techniques are used for the discretization of the transport equations. The main problem associated with the vorticity stream-function formulation is the lack of boundary conditions for vorticity at solid surfaces. Here an implicit treatment of the vorticity at no-slip boundaries is incorporated in a predictor-multicorrector time integration scheme. For the primitive variable formulation, mixed finite-element approximations are used. A nine-node element and a four-node + bubble element have been implemented. The latter is shown to exhibit a checkerboard pressure mode and a numerical treatment for this spurious pressure mode is proposed. The two methods are compared from the points of view of simulating internal and external flows and the possibilities of extensions to three dimensions.

Glaisner, F.↗
NASA NTRS1973

Numerical simulation of the world ocean circulation

A multi-level model, based on the primitive equations, is developed for simulating the temperature and velocity fields produced in the world ocean by differential heating and surface wind stress. The model ocean has constant depth, free slip at the lower boundary, and neglects momentum advection; so that there is no energy exchange between the barotropic and baroclinic components of the motion, although the former influences the latter through temperature advection. The ocean model was designed to be coupled to the UCLA atmospheric general circulation model, for the study of the dynamics of climate and climate changes. But here, the model is tested by prescribing the observed seasonally varying surface wind stress and the incident solar radiation, the surface air temperature and humidity, cloudiness and the surface wind speed, which, together with the predicted ocean surface temperature, determine the surface flux of radiant energy, sensible heat and latent heat.

Takano, K.↗
NASA NTRS1974

The GISS model of the global atmosphere

A description and numerical results are presented for a global atmospheric circulation model developed at the Goddard Institute for Space Studies (GISS). The model version described is a 9-level primitive-equation model in sigma coordinates. It includes a realistic distribution of continents, oceans and topography. Detailed calculations of energy transfer by solar and terrestrial radiation make use of cloud and water vapor fields calculated by the model. The model hydrologic cycle includes two precipitation mechanisms: large-scale supersaturation and a parameterization of subgrid-scale cumulus convection. Results are presented both from a comparison of the 13th to the 43rd days (January) of one integration with climatological statistics, and from five short-range forecasting experiments. In the extended integration, the near-equilibrium January-mean model atmosphere exhibits an energy cycle in good agreement with observational estimates, together with generally realistic zonal mean fields of winds, temperature, humidity, transports, diabatic heating, evaporation, precipitation, and cloud cover.

Somerville, R. C. J.↗
NASA NTRS1974

Short-range forecasts with the GISS model of the global atmosphere

A nine-layer, primitive equation (PE) model of the global atmosphere developed at the Goddard Institute for Space Studies (GISS) has been used to generate six 48-hr forecasts during December 1972 and January 1973. Operational analyses north of 18 N and experimental global analyses made available by the National Meteorological Center (NMC), NOAA, were used as the initial conditions; the operational analyses were used to verify the forecasts at 12-hr intervals over the Northern Hemisphere north of 22 N. The combined analyses were used to verify the forecasts in the global domain. RMS errors of the sea-level pressure, 1000-mb heights, and vector geostrophic winds, and 500-mb heights and vector geostrophic winds indicate that the GISS forecasts have skill comparable to those made by operational PE models.

Druyan, L. M.↗
NASA NTRS1975

Three-dimensional time dependent computation of turbulent flow

The three-dimensional, primitive equations of motion are solved numerically for the case of isotropic box turbulence and the distortion of homogeneous turbulence by irrotational plane strain at large Reynolds numbers. A Gaussian filter is applied to governing equations to define the large scale field. This gives rise to additional second order computed scale stresses (Leonard stresses). The residual stresses are simulated through an eddy viscosity. Uniform grids are used, with a fourth order differencing scheme in space and a second order Adams-Bashforth predictor for explicit time stepping. The results are compared to the experiments and statistical information extracted from the computer generated data.

Kwak, D.↗
NASA NTRS1974

Short-range forecasts with the GISS model of the global atmosphere

Results of tests carried out during the period from December 1972 through January 1973 to determine the short-term forecasting accuracy of a multilevel numerical primitive-equation (PE) model of the global atmosphere developed at the Goddard Institute for Space Studies (GISS). Six 48-hr forecasts were carried out with the aid of this model, using nine vertical levels and a horizontal grid spacing of 4 deg in latitude and 5 deg in longitude for an effective grid point separation averaging slightly more than 400 km. Verification of forecast sea-level pressures, 1000-mb heights, and 500-mb heights, as well as 1000-mb and 500-mb vector geostrophic winds, shows that the model has forecast skill comparable to that of operational PE models. Based on the 36-hr evolution of 18 extratropical cyclones, the model forecasts exhibit a tendency toward underestimating their propagation speeds and overestimating their central pressures. Both deficiencies are attributed to inadequate horizontal grid resolution. Quantitative verification of forecast surface temperatures over the eastern United States shows a forecast skill equal to that achieved by combined dynamical-statistical procedures.

Druyan, L. M.↗
NASA NTRS1975

A numerical investigation of severe thunderstorm gust fronts

A numerical model was developed to simulate the evolution and structure of severe thunderstorm gust fronts. The model is a non-hydrostatic, fine resolution, cross-sectional primitive equation model. Two-dimensional horizontal and vertical equations of motion, the continuity equation, and the thermodynamic energy equation were utilized. It was shown that two dominant factors influencing gust front configuration are surface friction and the solenoidal field coincident with the front. It is suggested that solenoidal accelerations oppose the deceleration of surface friction. After a downdraft is initiated in the model, these opposing tendencies soon reach a balance and the gust front achieves a quasi-steady configuration. Thus, the experiments indicate that surface friction does not induce a cycle of front formation and collapse. In addition, the effect of evaporative cooling in producing a vigorous downdraft was parameterized by a local cooling function. Greater cooling in the downdraft results in a more intense gust front that exhibits stronger wind maximums and greater shears. The ambient air stability was shown to be an important factor influencing the depth of the cold outflow.

Mitchell, K. E.↗
NASA NTRS1976

Calculation of wind-driven surface currents in the North Atlantic Ocean

Calculations to simulate the wind driven near surface currents of the North Atlantic Ocean are described. The primitive equations were integrated on a finite difference grid with a horizontal resolution of 2.5 deg in longitude and latitude. The model ocean was homogeneous with a uniform depth of 100 m and with five levels in the vertical direction. A form of the rigid-lid approximation was applied. Generally, the computed surface current patterns agreed with observed currents. The development of a subsurface equatorial countercurrent was observed.

Rees, T. H.↗
NASA NTRS1976

Atmospheric analysis and prediction model development, volume 1

A set of hemispheric atmospheric analysis and prediction models was designed and tested. All programs were executed on either a 63 x 63 or 187 x 187 polar stereographic grid of the Northern Hemisphere. Parameters for objective analysis included sea surface temperature, sea level pressure, and twelve levels (from 1,000 to 100 millibars) of temperatures, heights, and winds. Stratospheric extensions (up to 10 millibars) were also provided. Four versions of a complex atmospheric prediction model, based on primitive equations, were programmed and tested. These models were executed on either the 63 x 63 or 187 x 187 grid, using either five or ten computational layers. The coarse-mesh (63 x 63) models were tested using real data for the period 21-23 April 1976. The fine-mesh (187 x 187) models were debugged, but insufficient computer resources precluded production tests. Preliminary test results for the 63 x 63 models are provided. Problem areas and proposed solutions are discussed.

Kesel, P. G.↗
NASA NTRS1976

Meteorological analysis models, volume 2

As part of the SEASAT program, two sets of analysis programs were developed. One set of programs produce 63 x 63 horizontal mesh analyses on a polar stereographic grid. The other set produces 187 x 187 third mesh analyses. The parameters analyzed include sea surface temperature, sea level pressure and twelve levels of upper air temperature, height and wind analyses. Both sets use operational data provided by a weather bureau. The analysis output is used to initialize the primitive equation forecast models also included.

Langland, R. A.↗
NASA NTRS1976

Prediction of convective activity using a system of parasitic-nested numerical models

A limited area, three dimensional, moist, primitive equation (PE) model is developed to test the sensitivity of quantitative precipitation forecasts to the initial relative humidity distribution. Special emphasis is placed on the squall-line region. To accomplish the desired goal, time dependent lateral boundaries and a general convective parameterization scheme suitable for mid-latitude systems were developed. The sequential plume convective parameterization scheme presented is designed to have the versatility necessary in mid-latitudes and to be applicable for short-range forecasts. The results indicate that the scheme is able to function in the frontally forced squallline region, in the gently rising altostratus region ahead of the approaching low center, and in the over-riding region ahead of the warm front. Three experiments are discussed.

Perkey, D. J.↗
NASA NTRS1976

A numerical study of the radiative-dynamic interaction in the stratosphere by a time-dependent zonal mean circulation model

The response of a zonal mean circulation model was compared to different heating functions used to simulate the interannual climate cycle. Two heating functions were used: A simple Newtonian heating function and a sophisticated heating function (Ramanathan, 1976). The emphasis was on the coupling between radiative transfer and the circulation in the stratosphere. The model was formulated in terms of primitive equations with 10 deg latitudinal resolution and 18 levels in pressure. Eddy transport of heat in the stratosphere was treated by the Reed and German (1965) formulation, while Stone's (1973) formulation is used in the troposphere. Momentum eddy fluxes are specified by climate statistics after Harwood and Pyle (1975). The model dynamics resulting from the two runs are compared to demonstrate the benefits of using detailed heating functions.

Turner, R. E.↗
NASA NTRS1977

Numerical simulation of the Gulf Stream and Mid-Ocean eddies

Oceanographic field measurements have revealed intense, transient mesoscale motions in many parts of the world ocean. The circulation of the western North Atlantic, considered in the present study as a rectangular basin, is simulated with a primitive equation model that has five levels and a horizontal grid size of 37 km. The model ocean is driven by a 2.5 gyre pattern of steady zonal wind stress and by a Newtonian-type surface heating. Two cases are considered: the first uses a Laplacian formulation for the subgrid-scale lateral diffusions of heat and momentum; the second uses a highly scale-selective biharmonic formulation for these diffusions. An analysis of the heat transport, in the biharmonic experiment, shows that the horizontal transport of heat by eddies is much larger than the subgrid-scale horizontal heat diffusion. In the Gulf Stream region, the eddy heat transport is comparable to the effect of a lateral diffusion coefficient of 10 to the 7th sq cm/s.

Semtner, A. J., Jr.↗
NASA NTRS1977

Updating prediction models by dynamical relaxation - An examination of the technique

A dynamical relaxation technique for updating prediction models is analyzed with the help of the linear and nonlinear barotropic primitive equations. It is assumed that a complete four-dimensional time history of some prescribed subset of the meteorological variables is known. The rate of adaptation of the flow variables toward the true state is determined for a linearized f-model, and for mid-latitude and equatorial beta-plane models. The results of the analysis are corroborated by numerical experiments with the nonlinear shallow-water equations.

Davies, H. C.↗
NASA NTRS1977

Atmospheric model development in support of SEASAT. Volume 2: Analysis models

As part of the SEASAT program of NASA, two sets of analysis programs were developed for the Jet Propulsion Laboratory. One set of programs produce 63 x 63 horizontal mesh analyses on a polar stereographic grid. The other set produces 187 x 187 third mesh analyses. The parameters analyzed include sea surface temperature, sea level pressure and twelve levels of upper air temperature, height and wind analyses. The analysis output is used to initialize the primitive equation forecast models.

Langland, R. A.↗