The wind shear theory of the formation of temperature zone sporadic E layers.
Wind shear theory of formation of temperate zone blanketing sporadic E layers, noting motion equation for ions in ionospheric E region
SEARCH · Search NASA
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.
Wind shear theory of formation of temperate zone blanketing sporadic E layers, noting motion equation for ions in ionospheric E region
Solution for shock layer energy equation with mass and convective heat transfer at hypersonic speed
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.
It is found that layered ice and dust deposits accumulate on the smooth and banded terrains of the north and south polar cap deposits of Mars, while erosion exposes these layers on the equator-facing slopes of layered terrain. The stripping of layers occurs in areas of strong and multidirectional wind, forming striped terrain. Also noted is the production of steep scarps with crescentic planform by local wind scour, in frequent association with dune fields positioned immediately downwind. Unconformities, which are common within layered deposits in areas of complex topography at the north pole, indicate that deposition episodes of 5-10 layers alternate with the slight erosion of pole-facing trough walls. It is suggested that poles may migrate poleward, replacing the relief increase by depositional infilling when they near the poles.
Transonic flow fields about a low-aspect-ratio advanced technology wing have been computed using a viscous/inviscid zonal approach. The flow field near the wing where viscous effects are important was solved using the 'Reynolds-Averaged Navier-Stokes Equations' in 'thin-layer' form. The Euler equations were used to determine the flow field in regions away from the wing where viscous effects are insignificant. A zonal grid using an H-H topology was generated around the wing by first solving a set of Poisson's equations for the global grid. This grid was then subdivided into separate zones of viscous or inviscid flow as suggested by the flow physics. A series of flow cases were computed and compared with corresponding sets of experimental data. All cases showed good agreement with experiment in terms of the pressure field. Also, a good correlation between computed separated surface flow and experimental oil flow was obtained.
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.
A comparison study has been conducted using four recently developed parabolized Navier-Stokes (PNS) codes which have the capability of predicting finite-rate, chemically reacting flows over three-dimensional bodies. These are the (1) UPS code, (2) the STUFF code, (3) the TONIC code, and (4) the VRA-PNS code. All of the codes use the same seven-species, single-temperature air chemistry model, but otherwise they are unique, with different capabilities and characteristics. The differences include upwinding vs central differencing, strongly-coupled vs weakly-coupled chemistry, shock capturing vs shock fitting, finite volume vs finite difference, and full PNS vs thin-layer PNS equations. Three test cases were utilized to compare the codes. The comparisons presented indicate a good agreement among the codes tested.
The effect of forcing by monthly wind-stress variability is assessed with respect to the behavior of an isopycnic coordinate general circulation model. The isopycnic model incorporates an embedded mixed layer, and the equations are given for the interior layers, the surface mixed layer, and the atmospheric forcing. The anomalous variability of several variables is examined in a simulation using data from 1970-1985 with specific attention given to the anomalous response related to the 1982-1983 ENSO. The model fields describing sea-surface temperature (SST), sea level, and currents are found to be comparable to those reported from observational campaigns. Anomalous values of SST compare best with observed data when averaged over large spatial regions, and the amplitudes of the mean horizontal currents and anomalous zonal currents are comparable except with respect to amplitude data. The present study permits the simulation of basin-scale tropical Pacific mixed-layer depth and the mean vertical velocity field.
The seasonal formation of a large (500-800 km diameter) anticyclonic eddy in the upper 300-400 m of the eastern Arabian Sea during the northest monsoon period (December-April) is indicated fom hydrographic and satellite altimetry sea level observations, as well as from numerical model experiments. The center of the eddy circulation is approximately 10 deg N, 70 deg E, just to the west of the north-south Laccadive Island chain. In this paper the eddy is called the Laccadive High (LH). In some ways it is like a mirrorlike counterpart to the Great Whirl that develops during the southwest monsoon of the Somali coast (western Arabian Sea). The LH occurs at the same latitude but on the opposite side of the basin during the reversed monsoon. It is different from the Great Whirl, however, in its formation process, its intensity, and its decay. The hydrographic data obtained from surveys all during a single season give sufficiently close station spacing to allow reasonable contouring of the geopotential surfaces and of the properties within and around the LH region with minimum time aliasing. The Geostat altimeter record extends over 4 years, during which the seasonal variability of the LH indicates a dynamic relief of approximately 15-20 cm, which is in good agreement with the hydrographics observations. The altimetry time series also suggests a westward translation of the LH by January with a subsequent dissipation in midbasin. The model used is a wind-forced three-layer primitive equation model which depicts a LH agreement with the timing, position, and amplitude of both the hydrographic and altimetric measurements. The numerical simulation includes a passive tracer located in the Western Bay of Bengal; the western advection of the tracer around the south coasts of Sri Lanka and India in December and January is consistent with the appearance of low-salinity water observed to extend into the Arabian Sea during this period. The modeling studies suggest that both local and remote forcing are important in formation of the LH.
Studies are performed on both the Stratospheric Aerosol and Gas Experiment (SAGE) II (1985-1989) and Solar Backscatter Ultraviolet Spectrometer (SBUV)-Version 6 (1979-1986) global ozone vertical structure satellite data sets to determine the long-term trends in ozone as a function of altitude (pressure) and latitude. SAGE II data are only available during the period of increasing solar activity and show increases in ozone with time in the upper stratosphere which are attributed largely to rising solar activity. Looking at this data set independently, the solar effects and trends are highly coupled and cannot be clearly separated. However, a study of combined SBUV and SAGE II data over the 11-year solar cycle shows a clear response of ozone to 11-year solar variations and allows a decoupling of solar effects, quasibiennial oscillations (QBO), and trends. The detailed pattern of long-term ozone trends become clear using this approach. In the upper stratosphere, ozone depletion increases sharply with latitude. Global trends are fairly symmetric about the equator but are somewhat stronger in the Southern Hemisphere. Near the equator, some layers of ozone are decreasing with time while others appear to be increasing. Near 30 mb, there is evidence of intrusion to mid latitudes of high latitude negative trends. Near 15 mb, trends appear to be very weak. Near the tropopause there appears to be strong ozone depletion on a global scale. Two regions of unexpectedly strong ozone response to 11-year solar variations were detected: the first near 2 mb and the second near 30 mb at low latitudes and near 15 mb at mid latitudes.
A steady-state scheme for data assimilation in the context of a single, short period (relative to a day), sun-synchronous, polar-orbiting satellite is examined. If the satellite takes observations continuously, the gains, which are the weights for blending observations and predictions together, are steady in time. For a linear system forced by random noise, the optimal steady-state gains (Wiener gains) are equivalent to those of a Kalman filter. Computing the Kalman gains increases the computational cost of the model by a large factor, but computing the Wiener gains does not. The latter are computed by iteration using prior estimates of the gains to assimilate simulated observations of one run of the model, termed 'truth' into another run termed 'prediction'. At each stage, the prediction errors form the basis for the next estimate of the gains. Steady state is achieved after three or four iterations. Further simplification is achieved by making the gains depend on longitudinal distance from the observation point, not on absolute longitude. For a single-layer primitive equation model, the scheme works well even if only the mass field is observed but not the velocity field. Although the scheme was developed for Mars Observer, it should be applicable to data retrieved from Earth atmosphere satellites, for example, UARS.
As part of JPL's ocean data assimilation effort to study ocean circulation and seasonal-interannual climate variability, sea level anomaly observed by TOPEX altimeter, together with sea surface temperature and wind stress data, are assimilated into a simple coupled ocean atmosphere model of the tropical Pacific. Model-data consistency is examined. Impact of the assimilation (as initialization) on El Nino Southern Oscillation (ENSO) forecasts is evaluated. The coupled model consists of a shallow water component with two baroclinic modes, an Ekman shear layer, a simplified mixed-layer temperature equation, and a statistical atmosphere based on dominant correlations between historical surface temperature and wind stress anomaly data. The adjoins method is used to fit the coupled model to the data over various six-month periods from late 1996 to early 1998 by optimally adjusting the initial state, model parameters, and basis functions of the statistical atmosphere. On average, the coupled model can be fitted to the data to approximately within the data and representation errors (5 cm, 0.5 C, and 10 sq m/sq m for sea level, surface temperature, and pseudo wind stress anomalies, respectively). The estimated fields resemble observed spatio-temporal structure reasonably well. Hindcasts/forecasts of the 1997/1998 El Nino initialized from forced estimated ocean states and parameters are much more realistic than those simply initialized from ocean states (see figure below). In particular, the ability of the model to produce significant warming beyond the initial state is dramatically improved. Parameter estimation, which compensates for some model errors, is found to be important to obtaining better fits of the model to data and to improving forecasts.
The unsteady, compressible, thin-layer Navier-Stokes equations are used to obtain three-dimensional, asymmetric, vortex-flow solutions around cones and cone-cylinder configurations. The equations are solved using an implicit, upwind, flux-difference splitting, finite-volume scheme. The computational applications cover asymmetric flows around a 5 semi-apex angle cone of unit length at various Reynolds number. Next, a cylindrical afterbody of various length is added to the conical forebody to study the effect of the length of cylindrical afterbody on the flow asymmetry. All the asymmetric flow solutions are obtained by using a short-duration side-slip disturbance.
A comprehensive set of scaling parameters is developed for the ice accretion process by analyzing the energy equations of the dynamic freezing zone and the already frozen ice layer, the continuity equation associated with supercooled liquid droplets entering into and impacting within the dynamic freezing zone, and energy equation of the ice layer. No initial arbitrary judgments are made regarding the relative magnitudes of each of the terms. The method of intrinsic reference variables in employed in order to develop the appropriate scaling parameters and their relative significance in rime icing conditions in an orderly process, rather than utilizing empiricism. The significance of these parameters is examined and the parameters are combined with scaling criteria related to droplet trajectory similitude.
Semigroup solution of time dependent continuity equation for ion density in ionospheric F 2 region
A technique is implemented for computing hypersonic aeroheating, shear stress, and other flow properties on the windward side of a three-dimensional (3D) blunt body. The technique uses a 2D/axisymmetric flow solver modified by scale factors for a, corresponding equivalent axisymmetric body. Examples are given in which a 2D solver is used to calculate the flow at selected meridional planes on elliptic paraboloids in reentry flight. The report describes the equations and the codes used to convert the body surface parameters into input used to scale the 2D viscous shock layer equations in the axisymmetric viscous shock layer code. Very good agreement is obtained with solutions to finite rate chemistry 3D thin viscous shock layer equations for a finite rate catalytic body.
A unified central-difference finite-volume Euler and Navier-Stokes solver with four-stage Runge-Kutta time stepping is presented. The computer code developed for this purpose is capable of solving the standard set and nonstandard sets (zero-total-pressure loss) of Euler equations and the thin-layer and full Navier-Stokes equations. Applications are presented for conical supersonic flows with weak shocks using the standard and nonstandard sets of Euler equations, and the thin-layer and full Navier-Stokes equations for sharp and round-edged delta wings. Applications are also presented for three-dimensional transonic and subsonic flows using the standard set of Euler equations for sharp-edged delta wings. The computational results of the different sets of equations are compared with each other and with the experimental results and conclusions on the validity of these sets to these applications, are presented.
This paper solves an 'incremental' form of the sensitivity equations derived by differentiating the discretized thin-layer Navier Stokes equations with respect to certain design variables of interest. The equations are solved with a parallel, preconditioned Generalized Minimal RESidual (GMRES) solver on a distributed-memory architecture. The 'serial' sensitivity analysis code is parallelized by using the Single Program Multiple Data (SPMD) programming model, domain decomposition techniques, and message-passing tools. Sensitivity derivatives are computed for low and high Reynolds number flows over a NACA 1406 airfoil on a 32-processor Intel Hypercube, and found to be identical to those computed on a single-processor Cray Y-MP. It is estimated that the parallel sensitivity analysis code has to be run on 40-50 processors of the Intel Hypercube in order to match the single-processor processing time of a Cray Y-MP.