Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “total variation”

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

Response of Tropical Clouds to the Interannual Variation of Sea Surface Temperature

Connections between the large-scale interannual variations of clouds, deep convection, atmospheric winds. vertical thermodynamic structure, and SSTs over global tropical oceans are examined over the period July 1983 - December 1990. The SST warming associated with El Nino had a significant impact on the global tropical cloud field, although the warming itself was confined to the equatorial central and eastern Pacific. Extensive variations of the total cloud field occurred in the northeastern Indian, western and central Pacific, and western Atlantic Oceans. The changes of high and middle clouds dominated the total cloud variation in these regions. Total cloud variation was relatively weak in the eastern Pacific and the Atlantic because of the cancellation between the changes of high and low clouds. The variation of low clouds dominated the total cloud change in those areas. The destabilization of the lapse rate between 900 and 750 mb was more important for enhancing convective instability than was the change of local SSTs in the equatorial central Pacific during the 1997 El Nino. This destabilization is associated with anomalous rising motion in that region. As a result. convection and high and middle clouds increased in the equatorial central Pacific, In the subtropical Pacific, both the change of lapse rate between 900 and 750 mb associated A,ith anomalous subsidence and the decrease of boundary-layer buoyancy due to a decrease of temperature and moisture played an important role in enhancing convective stability. Consequently, convection, as well its high and middle clouds, decreased in these areas. The change ot'low clouds in the equatorial and southeastern Atlantic was correlated to both local SSTs and the SST changes in the equatorial eastern Pacific. In this area. the increase of low clouds was consistent with the sharper inversion during the 1987 El Nino, The strengthening of the inversion was not caused by a local SST change. although the local SST change appeared to he correlated to the change of low clouds. The coherence between clouds and SST tendency shows that SST tendency leads cloud variation in the equatorial Pacific. Thus, the change of clouds does not dominate the sign of SST tendency even though the cloud change was maximum during the 1987 El Nino. In some ideas of the Indian, subtropical Pacific, and North Atlantic Oceans, cloud change leads SST tendency. Cloud change might affect SST tendency in these regions.

Fu, Rong↗

Downscaling Satellite Precipitation with Emphasis on Extremes: A Variational 1-Norm Regularization in the Derivative Domain

The increasing availability of precipitation observations from space, e.g., from the Tropical Rainfall Measuring Mission (TRMM) and the forthcoming Global Precipitation Measuring (GPM) Mission, has fueled renewed interest in developing frameworks for downscaling and multi-sensor data fusion that can handle large data sets in computationally efficient ways while optimally reproducing desired properties of the underlying rainfall fields. Of special interest is the reproduction of extreme precipitation intensities and gradients, as these are directly relevant to hazard prediction. In this paper, we present a new formalism for downscaling satellite precipitation observations, which explicitly allows for the preservation of some key geometrical and statistical properties of spatial precipitation. These include sharp intensity gradients (due to high-intensity regions embedded within lower-intensity areas), coherent spatial structures (due to regions of slowly varying rainfall),and thicker-than-Gaussian tails of precipitation gradients and intensities. Specifically, we pose the downscaling problem as a discrete inverse problem and solve it via a regularized variational approach (variational downscaling) where the regularization term is selected to impose the desired smoothness in the solution while allowing for some steep gradients(called 1-norm or total variation regularization). We demonstrate the duality between this geometrically inspired solution and its Bayesian statistical interpretation, which is equivalent to assuming a Laplace prior distribution for the precipitation intensities in the derivative (wavelet) space. When the observation operator is not known, we discuss the effect of its misspecification and explore a previously proposed dictionary-based sparse inverse downscaling methodology to indirectly learn the observation operator from a database of coincidental high- and low-resolution observations. The proposed method and ideas are illustrated in case studies featuring the downscaling of a hurricane precipitation field.

Hurricanes↗

Nonlinear evolution of tidally distorted accretion disks: Two-dimensional simulations

According to a previously published linear analysis, the tidal distortion of accretion disks in binary star systems produces a local hydrodynamic instability to m = 1 internal waves, which may have arbitrarily small wavelengths in the absence of viscosity. The instability is three-dimensional and approximately incrompressible. To explore the nonlinear outcome of this instability, we develop a shearing-sheet approximation on scales comparable to the disk thickness. The large-scale azimuthal variation of the disk is represented by varying the local metric with time (local orbital phase). The hydrodynamic equations can then be posed two-dimensionally on local meridional planes. We solve these equations with a second-order gasdynamical code based on the Total-Variation-Diminishing scheme. Our simulations confirm the predicted linear growth rate. The modes saturate chaotically at velocities scaling as the product of the linear growth rate and the wavelength. If the wavelength is small compared with the disk thickness, the modes remain nearly incompressible even when nonlinear. The two-dimensional power spectrum of velocities after saturation is roughly isotropic and extends over a broad range of scales in an approximately power-law fashion. We measure the heating rate associated with the nonlinear dissipation of the modes. The dissipation implies a secular torque on the disk and a return of angular momentum to the secondary star via the tidal potential. The estimated torque is somewhat larger than the tidal torque produced by maximal disk viscosity (alpha approximately 1). At least in these two-dimensional simulations, however, there is no significant angular momentum flux within the disk.

Ryu, Dongsu↗

Constraints on mantle viscosity from relative sea level variations in Hudson Bay

Frechet kernels for the RSL data in Hunson Bay are computed to determine the detailed depth-dependent sensitivity of the data to variations in a viscosity profile which is consistent with the past inferences. The RSL data provide a robust constraint on the average viscosity in the top half of the lower mantle only (the average must be near 10 exp 21 Pa s). The data admit models whose average viscosity in the deep mantle (below 1800 km depth) and in the upper mantle can differ significantly from the value (near 10 exp 21 Pa s) estimated for the top half of the lower mantle. The total variation in the viscosity from the surface to the CMB can exceed an order of magnitude or more and still satisfy the constraint provided by the Hudson Bay RSL data set. The necessity of invoking an isoviscous mantle model in previous studies is a consequence of the limited class of viscosity model solutions employed in those studies.

Mitrovica, J. X.↗

The oceanic contribution to the Earth's seasonal angular momentum budget

Seasonal variations in the speed of the Earth's rotation manifest themselves as fluctuations in the length of the day (LOD) with an amplitude of about 1000 microseconds. We know from previous work that at least 95% of these variations can be accounted for in terms of angular momentum exchanged between the atmosphere and the solid Earth. Here we examine the respective contributions of the Antarctic Circumpolar Current (ACC) and the global oceans to the Earth's seasonal angular momentum budget, using in situ data from the Drake Passage and results from both the oceanic regional model (Fine Resolution Antarctic Model -- FRAM) of Webb et al. (1991) and the global ocanic model of Maier-Reimer et al. (1993) as analyzed by Brosche et al. (1990). The estimated annual contribution of the ACC (2-4 microsec) is much smaller than the total variation in the oceanic models or the existing LOD-AAM residual (both approximately 15-20 microsec). The estimated semi-annual ACC contribution (3-8 microsec) is offset by counter-current further north in both oceanic models, which exhibit larger semi-annual variations in planetary angular momentum. Further refinements in the Earth's seasonal angular momentum budget, therefore, will require the full (planetary plus relative) contribution of the global oceans in addition to that of the ACC.

Dickey, J. O.↗

Atmospheric pressure (surface)

The total variation of pressure from day to day is relatively small. Rapid but slightly greater variations occur as the result of the passage of frontal systems, while the passage of a hurricane can cause somewhat larger, but still not significant changes for pressure environment design of space vehicles. Surface pressure extremes for various locations and their extreme ranges are given. These data use the results of a study of pressure extremes.

Daniels, G. E.↗

Crop identification and area estimation over large geographic areas using LANDSAT MSS data

The author has identified the following significant results. LANDSAT MSS data was adequate to accurately identify wheat in Kansas; corn and soybean estimates in Indiana were less accurate. Computer-aided analysis techniques were effectively used to extract crop identification information from LANDSAT data. Systematic sampling of entire counties made possible by computer classification methods resulted in very precise area estimates at county, district, and state levels. Training statistics were successfully extended from one county to other counties having similar crops and soils if the training areas sampled the total variation of the area to be classified.

Bauer, M. E.↗

HR 7275 - A new variable star

Three years of photometry in V and B of the UBV system are presented to confirm the suspicion of Herbst (1973) that HR 7275 is a variable star. The photometry is used to derive the photometric period, which proves to be about 3% shorter than the spectroscopically determined optical period of 28.59 d. Total variation observed during the three years was 0.22 m in the V, and the light curve was always asymmetrical.

Fried, R. E.↗

TVD finite difference schemes and artificial viscosity

The total variation diminishing (TVD) finite difference scheme can be interpreted as a Lax-Wendroff scheme plus an upwind weighted artificial dissipation term. If a particular flux limiter is chosen and the requirement for upwind weighting is removed, an artificial dissipation term which is based on the theory of TVD schemes is obtained which does not contain any problem dependent parameters and which can be added to existing MacCormack method codes. Numerical experiments to examine the performance of this new method are discussed.

Davis, S. F.↗

Second- and third-order upwind difference schemes for hyperbolic conservation laws

Second- and third-order two time-level five-point explicit upwind-difference schemes are described for the numerical solution of hyperbolic systems of conservation laws and applied to the Euler equations of inviscid gas dynamics. Nonliner smoothing techniques are used to make the schemes total variation diminishing. In the method both hyperbolicity and conservation properties of the hyperbolic conservation laws are combined in a very natural way by introducing a normalized Jacobian matrix of the hyperbolic system. Entropy satisfying shock transition operators which are consistent with the upwind differencing are locally introduced when transonic shock transition is detected. Schemes thus constructed are suitable for shockcapturing calculations. The stability and the global order of accuracy of the proposed schemes are examined. Numerical experiments for the inviscid Burgers equation and the compressible Euler equations in one and two space dimensions involving various situations of aerodynamic interest are included and compared.

Yang, J. Y.↗

Generalized formulation of TVD Lax-Wendroff schemes

The work of Davis which imports the concept of total variation diminution (TVD) into non-upwinded, Lax-Wendroff type schemes is reformulated in a way which is easier to analyze. The analysis reveals a class of TVD schemes not observed by Davis. Only the case of one dimensional linear advection is treated.

Roe, P. L.↗

A geometric approach to high resolution TVD schemes

A geometric approach, similar to Van Leer's MUSCL schemes, is used to construct a second-order accurate generalization of Godunov's method for solving scalar conservation laws. By making suitable approximations, a scheme is obtained which is easy to implement and total variation diminishing. The entropy condition is also investigated from the standpoint of the spreading of rarefaction waves. Quantitative information is obtained for Godunov's method on the rate of spreading which explain the kinks in rarefaction waves often observed at the sonic point.

Goodman, J. B.↗

Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Harten, A.↗

Shock capturing

Recent developments which have improved the understanding of how finite difference methods resolve discontinuous solutions to hyperbolic partial differential equations are discussed. As a result of this understanding improved shock capturing methods are currently being developed and tested. Some of these methods are described and numerical results are presented showing their performance on problems containing shocks in one and two dimensions. A conservative difference scheme is defined. Conservation implies that, except in very special circumstances, shocks must be spread over at least two grid intervals. These two interval shocks are actually attained in one dimension if the shock is steady and an upwind scheme is used. By analyzing this case, the reason for this excellent shock resolution can be determined. This result is used to provide a mechanism for improving the resolution of two dimensional steady shocks. Unfortunately, this same analysis shows that these results cannot be extended to shocks which move relative to the computing grid. Total variation diminishing (TVD) finite difference schemes and flux limiters are introduced to deal with money shocks and contact discontinuities.

Davis, S. F.↗

Generalized formulation of a class of explicit and implicit TVD schemes

A one parameter family of second order explicit and implicit total variation diminishing (TVD) schemes is reformulated so that a simpler and wider group of limiters is included. The resulting scheme can be viewed as a symmetrical algorithm with a variety of numerical dissipation terms that are designed for weak solutions of hyperbolic problems. This is a generalization of Roe and Davis's recent works to a wider class of symmetric schemes other than Lax-Wendroff. The main properties of the present class of schemes are that they can be implicit, and when steady state calculations are sought, the numerical solution is independent of the time step.

Yee, H. C.↗

A new class of high accuracy TVD schemes for hyperbolic conservation laws

A new family of high accuracy Total Variation Diminishing (TVD) schemes has been developed. Members of the family include the conventional second-order TVD upwind scheme, various other second-order accurate TVD schemes with lower truncation error, and even a third-order accurate TVD approximation. All the schemes are defined with a five-point grid bandwidth. In this paper, the new algorithms are described for scalar equations, systems, and arbitrary coordinates. Selected numerical results are provided to illustrate the new algorithms and their properties.

Chakravarthy, S. R.↗

Application of TVD schemes for the Euler equations of gas dynamics

First-order, second-order, and implicit total variation diminishing (TVD) schemes are reviewed using the modified flux approach. Some transient and steady-state calculations are then carried out to illustrate the applicability of these schemes to the Euler equations. It is shown that the second-order explicit TVD schemes generate good shock resolution for both transient and steady-state one-dimensional and two-dimensional problems. Numerical experiments for a quasi-one-dimensional nozzle problem show that the second-order implicit TVD scheme produces a fairly rapid convergence rate and remains stable even when running with a Courant number of 10 to the 6th.

Yee, H. C.↗

Computing with high-resolution upwind schemes for hyperbolic equations

Computational aspects of modern high-resolution upwind finite-difference schemes for hyperbolic systems of conservation laws are examined. An operational unification is demonstrated for constructing a wide class of flux-difference-split and flux-split schemes based on the design principles underlying total variation diminishing (TVD) schemes. Consideration is also given to TVD scheme design by preprocessing, the extension of preprocessing and postprocessing approaches to general control volumes, the removal of expansion shocks and 'glitches', relaxation methods for implicit TVD schemes, and a new family of high-accuracy TVD schemes.

Chakravarthy, S. R.↗