Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shallow water”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Forced baroclinic ocean motions. III - The linear equatorial basin case

The linear response to simple wind stress of an equatorial ocean described by baroclinic shallow water equations is studied. Three types of basin are considered in treating the linear spin-up of the equatorial ocean: a symmetric basin with zonal walls distant from the equator (compared to the equatorial radius of deformation); a symmetric basin with zonal walls near the equator; and an asymmetric basin with one wall near the equator and one distant. The approach to the steady (Sverdrup) solutions is analyzed, with special attention given to the fast planetary response. Numerical treatments of spin-up for various winds in each type of basin are also presented.

Cane, M. A.↗
NASA NTRS1980

Mountain winds (revisited)

The prediction of extremely high wind speeds, at ground level on the downstream side of a mountain range, is possible by solving the initial value problem for a two-layered nonlinear shallow water model of the atmosphere. Three different numerical methods are described to find the solutions which may involve shocks: (1) the vonNeumann-Richtmyer artificial viscosity method, (2) a filtering scheme, and (3) a hybrid method.

Isaacson, E.↗
NASA NTRS1980

SEASAT-SAR data analysis in the US: An update

A graph is shown which compares SEASAT synthetic aperture radar (SAR) wave length measurements and Krasman's shallow water dispersion relationship for a swell system of deep water length 210 m and period 11.7 sec. The data were taken during the Duck experiment. The status of the production of optically correlated SAR data is given along with a tabulation of both digitally and optically processed SAR image data available to the user community as of Aug. 1979. Results from several SAR workshops are mentioned.

Dunne, J. A.↗
NASA NTRS1980

On time discretizations for spectral methods

New methods are introduced for the time integration of the Fourier and Chebyshev methods of solution for dynamic differential equations. These methods are unconditionally stable, even though no matrix inversions are required. Time steps are chosen by accuracy requirements alone. For the Fourier method both leapfrog and Runge-Kutta methods are considered. For the Chebyshev method only Runge-Kutta schemes are tested. Numerical calculations are presented to verify the analytic results. Applications to the shallow water equations are presented.

Gottlieb, D.↗
NASA NTRS1981

Refraction of coastal ocean waves

Refraction of gravity waves in the coastal area off Cape Hatteras, NC as documented by synthetic aperture radar (SAR) imagery from Seasat orbit 974 (collected on September 3, 1978) is discussed. An analysis of optical Fourier transforms (OFTs) from more than 70 geographical positions yields estimates of wavelength and wave direction for each position. In addition, independent estimates of the same two quantities are calculated using two simple theoretical wave-refraction models. The OFT results are then compared with the theoretical results. A statistical analysis shows a significant degree of linear correlation between the data sets. This is considered to indicate that the Seasat SAR produces imagery whose clarity is sufficient to show the refraction of gravity waves in shallow water.

Shuchman, R. A.↗
NASA NTRS1981

On the computational noise of finite-difference schemes used in ocean models

Different distributions of variables over the horizontal array of grid points in an ocean circulation model are investigated, using the shallow water equations as a guide in the choice of finite-difference schemes for use in ocean modeling. It is shown that the scheme with diffusive dissipation, in which the horizontal velocity is carried at the center and the height field is carried at each corner of a rectangular grid, successively suppresses numerical noise in a coarse (greater than 100 km) grid ocean model. For resolutions smaller than 50 km, it is shown that the scheme in which zonal velocity is carried at points to the east and west of the point of a rectangular grid where the height is carried, with meridional velocity carried to the north and south of the height point, can be free of noise for the gravest mode.

Batteen, M. L.↗
NASA NTRS1981

A simple atmospheric model on the sphere with 100% parallelism

A simple shallow water equations model used to test numerical weather prediction schemes is presented. The spherical geometry is transformed into a formally doubly periodic configuration which can be implemented on a massively parallel computer at 100% parallelism. A nonstaggered grid with all variables defined at the grid points (circles) is used. The parallel structure was simulated by performing step 1 in a FORTRAN array, on an Amdahl serial machine, changing the signs of the winds on the dark side.

Kalnay, E.↗
NASA NTRS1981

The response of a linear baroclinic equatorial ocean to periodic forcing

An investigation is conducted regarding the periodic response of the linear inviscid shallow water equations in a meridionally unbounded basin to zonal forcings at a single low frequency omega. A general solution in the long wave approximation and on an equatorial beta-plane is obtained by summing the Kelvin mode and the finite sum of Rossby modes whose turning points lie equatorward of the turning latitude at frequency omega. The results of the investigation suggest that even if the low frequency forcing has a simple structure, considerable spatial inhomogeneity in the deep ocean response would have to be expected. On the basis of linear inviscid theory, some conclusions are drawn about the causes of the differences between equatorial thermocline response in the Atlantic and Pacific.

Cane, M. A.↗
NASA NTRS1982

Initial conditions and Korteweg-de Vries solitons

The effects of rectangular initial data on the evolution of solitons governed by the Korteweg-de Vries equation is studied. Both isolated and separated disturbances are considered, providing some general insight into how the nature of the initial condition influences the appearance of solitons in the asymptotic state. The analytic approach is based on the inverse scattering transform which relates the initial condition to Schroedinger's equation. The results are used to model the initial shallow-water disturbances, and the results can be summarized by stating that the number of evolved solitons depends on the strength of each rectangular disturbance, the relative amplitudes of the rectangular disturbances, and the relative proximity of the disturbances.

Weidman, P. D.↗
NASA NTRS1982

Description and evaluation of the Acoustic Profiling of Ocean Currents (APOC) system used on R. V. Oceanus cruise 96 on 11-22 May 1981

The underway current profiling system which consists of a microprocessor controlled data logger that collects and formats data from a four beam Ametek-Straza 300 kHz acoustic Doppler current profiler, heading from the ship's gyrocompass, and navigation information from a Loran-C receiver and a satellite navigation unit is discussed. Data are recorded on magnetic tape and real time is calculated. Time averaging is required to remove effects of ship motion. An intercomparison is made with a moored vector measuring current meter (VMCM). The mean difference in hourly averaged APOC and VMCM currents over the four hour intercomparison is a few mm s minus including: two Gulf Stream crossings, a warm core ring survey, and shallow water in a frontal zone to the east of Nantucket Shoals.

Joyce, T. M.↗
NASA NTRS1981

Further SEASAT SAR coastal ocean wave analysis

Analysis techniques used to exploit SEASAT synthetic aperture radar (SAR) data of gravity waves are discussed and the SEASAT SAR's ability to monitor large scale variations in gravity wave fields in both deep and shallow water is evaluated. The SAR analysis techniques investigated included motion compensation adjustments and the semicausal model for spectral analysis of SAR wave data. It was determined that spectra generated from fast Fourier transform analysis (FFT) of SAR wave data were not significantly altered when either range telerotation adjustments or azimuth focus shifts were used during processing of the SAR signal histories, indicating that SEASAT imagery of gravity waves is not significantly improved or degraded by motion compensation adjustments. Evaluation of the semicausal (SC) model using SEASAT SAR data from Rev. 974 indicates that the SC spectral estimates were not significantly better than the FFT results.

Kasischke, E. S.↗
NASA NTRS1981

Optimal interpolation and the Kalman filter

The estimation theory of stochastic-dynamic systems is described and used in a numerical study of optimal interpolation. The general form of data assimilation methods is reviewed. The Kalman-Bucy, KB filter, and optimal interpolation (OI) filters are examined for effectiveness in performance as gain matrices using a one-dimensional form of the shallow-water equations. Control runs in the numerical analyses were performed for a ten-day forecast in concert with the OI method. The effects of optimality, initialization, and assimilation were studied. It was found that correct initialization is necessary in order to localize errors, especially near boundary points. Also, the use of small forecast error growth rates over data-sparse areas was determined to offset inaccurate modeling of correlation functions near boundaries.

Cohn, S.↗
NASA NTRS1981

An explicit mixed numerical method for mesoscale model

A mixed numerical method has been developed for mesoscale models. The technique consists of a forward difference scheme for time tendency terms, an upstream scheme for advective terms, and a central scheme for the other terms in a physical system. It is shown that the mixed method is conditionally stable and highly accurate for approximating the system of either shallow-water equations in one dimension or primitive equations in three dimensions. Since the technique is explicit and two time level, it conserves computer and programming resources.

Hsu, H.-M.↗
NASA NTRS1982

Numerical study of terrain-induced mesoscale motions in a mixed layer

Numerical integrations using a potential enstrophy-conserving scheme are presented for the flow within a mixed layer over hilly terrain using the hydrostatic shallow-water equations with a quadratic drag law. The mesoscale area treated is 150 km on a side; cyclic lateral boundary conditions are used. It is found that for the idealized conditions treated (no surface heating, no entrainment and no pressure adjustments aloft), the topography quickly induces a steady state flow pattern by means of surface friction. Unsteadiness does not occur unless a surface-friction Reynolds number is greater than approximately 100. Effects of varying the Rossby number, Froude number and terrain-height parameter are examined.

Han, Y.-J.↗
NASA NTRS1983

An analytic solution for the response of the neutral atmosphere to the high-latitude convection pattern

The analytic solutions to a two-dimensional, shallow-water model of the high-latitude thermosphere are developed. The equations are a linearized representation of the response of the neutral atmosphere to the momentum source due to the plasma convection pattern. The obtained results show how the neutral winds depend on the Pedersen and Hall ion drag coefficients and the horizontal wave number. The response depends on the dimensionless scale as it appears in the product of the Rossby radius of deformation and the wave number. The flow pattern is shifted in local time from the plasma flow pattern by an amount that depends on the average electron density. At the largest scales there is a 180 deg phase shift between the rotational flow components in the E and F regions. The model applies to quiet magnetic conditions and below 200 km of altitude.

Mikkelsen, I. S.↗
NASA NTRS1984

Computational aspects of the nonlinear normal mode initialization of the GLAS 4th order GCM

Using the normal modes of the GLAS 4th Order Model, a Machenhauer nonlinear normal mode initialization (NLNMI) was carried out for the external vertical mode using the GLAS 4th Order shallow water equations model for an equivalent depth corresponding to that associated with the external vertical mode. A simple procedure was devised which was directed at identifying computational modes by following the rate of increase of BAL sub M, the partial (with respect to the zonal wavenumber m) sum of squares of the time change of the normal mode coefficients (for fixed vertical mode index) varying over the latitude index L of symmetric or antisymmetric gravity waves. A working algorithm is presented which speeds up the convergence of the iterative Machenhauer NLNMI. A 24 h integration using the NLNMI state was carried out using both Matsuno and leap-frog time-integration schemes; these runs were then compared to a 24 h integration starting from a non-initialized state. The maximal impact of the nonlinear normal mode initialization was found to occur 6-10 hours after the initial time.

Navon, I. M.↗
NASA NTRS1984

On the effect of using the Shapiro filter to smooth winds on a sphere

Spatial differencing schemes which are not enstrophy conserving nor implicitly damping require global filtering of short waves to eliminate the build-up of energy in the shortest wavelengths due to aliasing. Takacs and Balgovind (1983) have shown that filtering on a sphere with a latitude dependent damping function will cause spurious vorticity and divergence source terms to occur if care is not taken to ensure the irrotationality of the gradients of the stream function and velocity potential. Using a shallow water model with fourth-order energy-conserving spatial differencing, it is found that using a 16th-order Shapiro (1979) filter on the winds and heights to control nonlinear instability also creates spurious source terms when the winds are filtered in the meridional direction.

Takacs, L. L.↗
NASA NTRS1984

Application of augmented-Lagrangian methods in meteorology: Comparison of different conjugate-gradient codes for large-scale minimization

A Lagrange multiplier method using techniques developed by Bertsekas (1982) was applied to solving the problem of enforcing simultaneous conservation of the nonlinear integral invariants of the shallow water equations on a limited area domain. This application of nonlinear constrained optimization is of the large dimensional type and the conjugate gradient method was found to be the only computationally viable method for the unconstrained minimization. Several conjugate-gradient codes were tested and compared for increasing accuracy requirements. Robustness and computational efficiency were the principal criteria.

Navon, I. M.↗