SEARCH · Engineering Papers
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.
As part of the SEASAT Satellite program of NASA, a survey of representative primitive equation atmospheric prediction models that exist in the world today was written for the Jet Propulsion Laboratory. Seventeen models developed by eleven different operational and research centers throughout the world are included in the survey. The surveys are tutorial in nature describing the features of the various models in a systematic manner.
The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.
The spin-up flow in a cylinder of homogeneous fluid has been examined both experimentally and numerically. A series of laser-Doppler measurements was made of the zonal flow over a range of Ekman numbers and Rossby numbers at various locations in the interior of the flow. These measurements exceed previous ones in accuracy. The weak inertial modes excited by the impulsive start are detectable. The numerical simulations used the primitive equations in axisymmetric form and employed finite-difference techniques on both constant and variable grids. The number of grid points necessary to resolve the Ekman layers was determined. A thorough comparison of the simulations and the experimental measurements is made which includes the details of the amplitude and frequency of the inertial modes. Agreement to within the experimental tolerance is achieved. Analytical results for conditions identical to those in the experiments are not available but some similar linear and nonlinear theories are also compared with the experiments.
A numerical forecast model is described which uses three-dimensional primitive equations to predict the evolution of the wind, temperature and moisture fields along with the subsynoptic-scale convective activity. There are 15 vertical levels set on a longitude-latitude horizontal grid. In this study, the horizontal grid distance is 70 km (0.625 deg latitude and 0.816 deg longitude at 40 deg N) and the time step is 90 sec. Orography is included and terrain coordinates are adopted. The forecast model is used to test two tornado-producing storms with very different synoptic structures. A 12 hr forecast is carried out for both cases. The rms errors of the forecast indicate much smaller errors in the middle troposphere than in the lower troposphere. The results demonstrate that a rather fine-mesh limited-area numerical model is capable of providing a valuable short-range forecast of the circulation and the mesoscale precipitation features.
The complete model is a semispectral model in which the longitudinal dependence is represented by expansion in zonal harmonics while the latitude and height dependencies are represented by a finite difference grid. The model is based on the primitive equations in the log pressure coordinate system. The lower boundary of the model domain is set at the 100 mb level (i.e., near the tropopause) and the effects of tropospheric forcing are included in the lower boundary condition. The upper boundary is at approximately 96 km, and the latitudinal extent is either global or hemispheric. The basic differential equations and boundary conditions are outlined. The finite difference equations are described. The initial conditions are discussed and a sample calculation is presented. The FORTRAN code is given in the appendix.
A solution is presented for matching boundary conditions across the interface of an isentropic and sigma coordinate hybrid model. A hybrid model based on the flux form of the primitive equations is developed which allows direct vertical exchange between the model domains, satisfies conservation principles with respect to transport processes, and maintains a smooth transition across the interface without need for artificial adjustment or parameterization schemes. The initial hybrid model simulations of a jet streak propagating in a zonal channel are used to test the feasibility of the hybrid model approach. High efficiency of the hybrid model is demonstrated.
A time-dependent, primitive equation, beta plane model that is two-dimensional in the horizontal has been developed to model wind-driven equatorial ocean circulation. A simple vertical structure consisting of two layers above the thermocline with the same constant density permits a steady-state undercurrent in the model. An analytical study of the linear dynamics of the model suggests that the addition of inertial effects is needed to simulate the undercurrent properly. Also, both linear and nonlinear dynamics of the model are investigated numerically. Such nonlinear response to wind stress as a strong eastward equatorial undercurrent and an intense eastward 'countercurrent' at three deg N are noted in the numerical results.
A time-dependent primitive equation model for an equatorial channel is used to assess the interaction of equatorial Kelvin and mixed Rossby-gravity waves with the mean flow. The proposed model involves a semiimplicit time-differencing scheme and a finite-difference grid in the meridional plane. It is shown that forced equatorial waves interact with mean flow to produce equatorial jets characterized by downward-moving westerly (Kelvin wave forcing) and easterly (mixed Rossby-gravity wave forcing) shear zones, respectively. For parameters characteristic of the observed waves in the equatorial stratosphere, the wave-mean flow interaction process always reduces the amplitude of any initial cross-equatorial mean wind shear. The mean flow profile tends to become symmetric about the equator as the interaction process continues.
The sea state in the North Atlantic Ocean has been monitored in February, 1977 to obtain data for verification of Geos 3 radio altimeter wave form data. The altimeter sea state data were found to be in good agreement with in situ sea state measurements and with aircraft remote sensing measurements of wave height. The average difference between the supporting data and Geos 3 was 0.34 m, and the standard deviation 0.61 m. Contour maps of sea state constructed for large regions of the earth's oceans on the basis of Geos 3 altimeter measurements exhibited systematic differences from operationally available computerized maps. It is thought that the underestimation of significant wave height maxima in the operational maps and the underestimation of the speed of movement of oceanic cyclones may be due to basic errors in the output of atmospheric primitive-equation models. Because of the narrow swath width of the Geos-3 sensor, data from several satellite orbits must be combined to construct contour maps, leading to low values of significant wave heights.
An investigation of the Martian polar cap winds and their response to a variety of factors is carried out by numerical experiments based on a zonally symmetric primitive equation model. The seasonal thermal forcing, mass exchange between polar caps and atmosphere, large-scale topography, and polar cap size are discussed, noting that topography has a small effect, but the circulation intensity increases with cap size. The model results show that surface winds near the edge of a retreating polar cap are enhanced, and that the surficial wind indicators near the south pole are formed during spring and those near the north pole during winter. It is suggested that the high-latitude dune fields in the northern hemisphere are formed when the terrain is covered by frost, and that the saltating particles are 'snowflakes' which formed by the mechanism proposed by Pollack (1976). The model results for the winter simulation compare favorably with general circulation model (GCM) calculations.
A general circulation model is presented which uses quadratically conservative, fourth order horizontal space differences on an unstaggered grid and second order vertical space differences with a forward-backward or a smooth leap frog time scheme to solve the primitive equations of motion. The dynamic equations for motion, finite difference equations, a discussion of the structure and flow chart of the program code, a program listing, and three relevent papers are given.
A mathematical model package for thermal pollution analyses and prediction is presented. These models, intended as user's manuals, are three dimensional and time dependent using the primitive equation approach. Although they have sufficient generality for application at sites with diverse topographical features; they also present specific instructions regarding data preparation for program execution and sample problems. The mathematical formulation of these models is presented including assumptions, approximations, governing equations, boundary and initial conditions, numerical method of solution, and same results.
The paper presents a simulation of the zonally averaged circulation in the middle atmosphere using a numerical model based on the primitive equations in log pressure coordinates. The circulation is driven radiatively by heating due to solar ultraviolet absorption by ozone and infrared cooling due to carbon dioxide and ozone; Rayleigh friction with a a short time constant above 70 km is included to simulate the strong mechanical dissipation which is hypothesized to exist in the vicinity of the mesopause due to turbulence associated with gravity waves and tides near the mesopause.
A unified theoretical approach to both the four-dimensional assimilation of asynoptic data and the initialization problem is attempted. This approach relies on the derivation of certain relationships between geopotential tendencies and tendencies of the horizontal velocity field in primitive-equation models of atmospheric flow. The approach is worked out and analyzed in detail for some simple barotropic models. Certain independent results of numerical experiments for the time-continuous assimilation of real asynoptic meteorological data into a complex, baroclinic weather prediction model are discussed in the context of the present approach. Tentative inferences are drawn for practical assimilation procedures.
The monthly mean sea surface temperature anomalies over tropical Altantic and rainfall anomalies over two selected stations for 25 years (1948-1972) were examined. It is found that the most severe drought events are associated with the simultaneous occurrence of warm sea surface temperature anomalies over north and cold sea surface temperature anomalies over south tropical Atlantic. Simultaneous occurrences of warm sea surface temperature anomaly at 15 deg N, 45 deg W and cold sea surface temperature anomaly at 15 deg S, 5 deg W were always associated with negative anomalies of rainfall, and vice versa. A simple primitive equation model is used to calculate the frictionally controlled and thermally driven circulation due to a prescribed heating function in a resting atmosphere.
The tidal heating components for the dusty Martian atmosphere are computed based on dust optical parameters estimated from Viking Lander imaging data, and used to compute the variation of the tidal surface pressure components at the Viking Lander sites as a function of season and the total vertical extinction optical depth of the atmosphere. An atmospheric tidal model is used which is based on the inviscid, hydrostatic primitive equations linearized about a motionless basic state the temperature of which varies only with height, and the profiles of the tidal forcing components are computed using a delta-Eddington approximation to the radiative transfer equations. Comparison of the model results with the observed variations of surface pressure and overhead dust opacity at the Viking Lander 1 site reveal that the dust opacities and optical parameters derived from imaging data are roughly representative of the global dust haze necessary to reproduce the observed surface pressure amplitudes, with the exception of the model-inferred asymmetry parameter, which is smaller during the onset of a great storm. The observed preferential enhancement of the semidiurnal tide with respect to the diurnal tide during dust storm onset is shown to be due primarily to the elevation of the tidal heating source in a very dusty atmosphere.
Namias' hypothesis, that anomalous snowcover on the eastern side of the North American continent can generate an anomalous east coast low pressure system and an anomalous inland high pressure system, is consistent with the time-averaged anomalous response from a nonlinear, primitive equation channel model with an idealized, flat land-sea arrangement. An attempt to understand and describe this anomalous response in the nonlinear model as a linear response to anomalous diabatic heating was largely unsuccessful, primarily because the anomalous eddy fluxes were also important. This unsuccessful attempt to describe the nonlinear model's time averages by linear theory then motivated several comparisons between linear and nonlinear severely truncated quasi-geostrophic models. It was also found in these models that the eddy fluxes were extremely important for forcing or dissipating the stationary eddies.
The influence on the stratospheric mean circulation of planetary wavenumber 2 disturbances excited by steady forcing at the 100 mb level is investigated using a global semi-spectral primitive equation model. There exists a critical forcing amplitude below which the waves have little effect on the mean flow, while above which the waves produce subseasonal time scale vacillations in the winter hemisphere, including both major and minor warmings. Wave transience induced by the evolving mean flow distribution is responsible for generating the vacillations, while thermal damping serves to reduce the wave-driven changes in the mean flow.