SEARCH · Engineering Papers

Results for “Shallow water 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1983

The effects of the Asselin time filter on numerical solutions to the linearized shallow-water wave equations

In the present investigation, a one-dimensional linearized analysis is used to determine the effect of Asselin's (1972) time filter on both the computational stability and phase error of numerical solutions for the shallow water wave equations, in cases with diffusion but without rotation. An attempt has been made to establish the approximate optimal values of the filtering parameter nu for each of the 'lagged', Dufort-Frankel, and Crank-Nicholson diffusion schemes, suppressing the computational wave mode without materially altering the physical wave mode. It is determined that in the presence of diffusion, the optimum filter length depends on whether waves are undergoing significant propagation. When moderate propagation is present, with or without diffusion, the Asselin filter has little effect on the spatial phase lag of the physical mode for the leapfrog advection scheme of the three diffusion schemes considered.

Schlesinger, R. E.↗

NASA NTRS1981

Determining solutions of hyperbolic systems from incomplete data

An investigation is conducted regarding first-order hyperbolic systems of partial differential equations, taking into account problems for which complete initial data are not available. Problems of the considered kind arise in geophysical applications where satellites are used to collect data. In global weather prediction, it is possible to derive atmospheric temperature and pressure reasonably well over the whole globe from satellite measurements; obtaining the wind field globally is more difficult. It is pointed out that a simple model of atmospheric flow investigated in numerical weather prediction is governed by the shallow water equations. The effect of the Coriolis term on the linearized shallow-water equations is studied.

Bube, K. P.↗

NASA NTRS1994

On the transition towards slow manifold in shallow-water and 3D Euler equations in a rotating frame

The long-time, asymptotic state of rotating homogeneous shallow-water equations is investigated. Our analysis is based on long-time averaged rotating shallow-water equations describing interactions of large-scale, horizontal, two-dimensional motions with surface inertial-gravity waves field for a shallow, uniformly rotating fluid layer. These equations are obtained in two steps: first by introducing a Poincare/Kelvin linear propagator directly into classical shallow-water equations, then by averaging. The averaged equations describe interaction of wave fields with large-scale motions on time scales long compared to the time scale 1/f(sub o) introduced by rotation (f(sub o)/2-angular velocity of background rotation). The present analysis is similar to the one presented by Waleffe (1991) for 3D Euler equations in a rotating frame. However, since three-wave interactions in rotating shallow-water equations are forbidden, the final equations describing the asymptotic state are simplified considerably. Special emphasis is given to a new conservation law found in the asymptotic state and decoupling of the dynamics of the divergence free part of the velocity field. The possible rising of a decoupled dynamics in the asymptotic state is also investigated for homogeneous turbulence subjected to a background rotation. In our analysis we use long-time expansion, where the velocity field is decomposed into the 'slow manifold' part (the manifold which is unaffected by the linear 'rapid' effects of rotation or the inertial waves) and a formal 3D disturbance. We derive the physical space version of the long-time averaged equations and consider an invariant, basis-free derivation. This formulation can be used to generalize Waleffe's (1991) helical decomposition to viscous inhomogeneous flows (e.g. problems in cylindrical geometry with no-slip boundary conditions on the cylinder surface and homogeneous in the vertical direction).

Mahalov, A.↗

NASA NTRS1988

Effects of using a posteriori methods for the conservation of integral invariants

The nature and effect of using a posteriori adjustments to nonconservative finite-difference schemes to enforce integral invariants of the corresponding analytic system are examined. The method of a posteriori integral constraint restoration is analyzed for the case of linear advection, and the harmonic response associated with the a posteriori adjustments is examined in detail. The conservative properties of the shallow water system are reviewed, and the constraint restoration algorithm applied to the shallow water equations are described. A comparison is made between forecasts obtained using implicit and a posteriori methods for the conservation of mass, energy, and potential enstrophy in the complete nonlinear shallow-water system.

Takacs, Lawrence L.↗

NASA NTRS1975

Initialization procedures for primitive equation models

A linear analysis and comparison of the damping properties of six dynamic initialization schemes is presented, indicating that the Okamura-Rivas scheme has the most efficient damping properties over the whole frequency range, and suggesting that it should be faster than the other methods and given more stable results. The results obtained with a nonlinear shallow water equations model agree well with the linear analysis. The Okamura-Rivas scheme attains complete balance in the equivalent of 5 to 6 hours of leapfrog forecasting, and requires in this model an order of magnitude less computation than the balance equation solution.

Grant, W. K. F.↗

NASA NTRS1985

Midlatitude tropical interactions during winter

Pre-FGGE and FGGE/MONEX data are used to identify short term midlatitude tropical and longitudinal interactions during the winter monsoon. These interactions occur as cold surges, which develop over the East Asian continent and penetrate deep into the tropics with fast gravity wave speed. The observed interactions that occur after a surge include cyclogenesis and enhanced convection in the equatorial region, feedback from equatorial convection to midlatitude circulation systems, tropical east-west (Walker) circulations, and cross-equatorial influence. These interactions are also studied theoretically by analytical solutions of linearized shallow water equations. Response to transient forcing (monsoon surges) are mainly in Rossby and Kelvin modes. When the forcing time scale is short, significant gravity modes are also excited. The responses closely resemble observed winter monsoon flow. Responses to stationary forcing show that deep (barotropic) motions propagate energy away into high latitudes and that shallow (baroclinic) motions are trapped around the equator. It is shown that the barotropic teleconnection-type response to tropical sources found in previous numerical studies was due to the specified vertical wind shear and surface friction.

Chang, C. P.↗

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 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 NTRS2017

Atmospheric Waves and Dynamics Beneath Jupiters Clouds from Radio Wavelength Observations

We observed Jupiter at wavelengths near 2 cm with the Karl G. Jansky Very Large Array (VLA) in February 2015. These frequencies are mostly sensitive to variations in ammonia abundance and probe between approx. 0.5- 2.0 bars of pressure in Jupiters atmosphere; within and below the visible cloud deck which has its base near 0.7 bars. The resultant observed data were projected into a cylindrical map of the planet with spatial resolution of approx. 1500 km at the equator. We have examined the data for atmospheric waves and observed a prominent bright belt of radio hotspot features near 10 N, likely connected to the same equatorial wave associated with the 5-m hotspots. We conducted a passive tracer power spectral wave analysis for the entire map and latitude regions corresponding to eastward and westward jets and compare our results to previous studies. The power spectra analysis revealed that the atmosphere sampled in our observation (excluding the NEB region) is in a 2-D turbulent regime and its dynamics are predominately governed by the shallow water equations. The Great Red Spot (GRS) is also very prominent and has a noticeable meridional asymmetry and we compare it, and nearby storms, with optical images. We find that the meridional radio profile has a global north-south hemisphere distinction and find correlations of it to optical intensity banding and to shear zones of the zonal wind profile over select regions of latitude. Amateur optical images taken before and after our observation complemented the radio wave- length map to investigate dynamics of the equatorial region in Jupiters atmosphere. We find that two radio hotspots at 2 cm are well correlated with optical plumes in the NEB, additionally revealing they are not the same 5 m hotspot features correlated with optical dark patches between adjacent plumes. This analysis exploits the VLAs upgraded sensitivity and explores the opportunities now possible when studying gas giants, especially atmospheric dynamics of layers beneath upper level clouds.

radio observations↗

NASA NTRS2001

Estimates of M2 Tidal Energy Dissipation from TOPEX/Poseidon Altimeter Data

Most of the tidal energy dissipation in the ocean occurs in shallow seas, as has long been recognized. However, recent work has suggested that a significant fraction of the dissipation, perhaps 1 TW or more, occurs in the deep ocean. This paper builds further evidence for that conclusion. More than 6 years of data from the TOPEX/Poseidon satellite altimeter are used to map the tidal dissipation rate throughout the world ocean. The dissipation rate is estimated as a balance between the rate of working by tidal forces and the energy flux divergence, computed using currents derived by least squares fitting of the altimeter data and the shallow water equations. Such calculations require dynamical assumptions, in particular about the nature of dissipation. To assess sensitivity of dissipation estimates to input assumptions, a large suite of tidal inversions based on a wide range of drag parameterizations and employing both real and synthetic altimeter data are compared. These experiments and Monte Carlo error fields from a generalized inverse model are used to establish error uncertainties for the dissipation estimates. Owing to the tight constraints on tidal elevation fields provided by the altimeter, area integrals of the energy balance are remarkably insensitive to required dynamical assumptions. Tidal energy dissipation is estimated for all major shallow seas (excluding individual polar seas) and compared with previous model and data-based estimates. Dissipation in the open ocean is significantly tnhanced around major bathymetric features, in a manner consistent with simple theories the generation of baroclinic tides.

Egbert, Gary D.↗

NASA NTRS2001

Orbital, Rotational and Climatic Interactions: Energy Dissipation and Angular Momentum Exchange in the Earth-Moon System

A numerical ocean tide model has been developed and tested using highly accurate TOPEX/Poseidon (T/P) tidal solutions. The hydrodynamic model is based on time stepping a finite difference approximation to the non-linear shallow water equations. Two novel features of our implementation are a rigorous treatment of self attraction and loading (SAL), and a physically based parameterization for internal tide (IT) radiation drag. The model was run for a range of grid resolutions, and with variations in model parameters and bathymetry. For a rational treatment of SAL and IT drag, the model run at high resolution (1/12 degree) fits the T/P solutions to within 5 cm RMS in the open ocean. Both the rigorous SAL treatment and the IT drag parameterization are required to obtain solutions of this quality. The sensitivity of the solution to perturbations in bathymetry suggest that the fit to T/P is probably now limited by errors in this critical input. Since the model is not constrained by any data, we can test the effect of dropping sea-level to match estimated bathymetry from the last glacial maximum (LGM). Our results suggest that the 100 m drop in sea-level in the LGM would have significantly increased tidal amplitudes in the North Atlantic, and increased overall tidal dissipation by about 40%. However, details in tidal solutions for the past 20 ka are sensitive to the assumed stratification. IT drag accounts for a significant fraction of dissipation, especially in the LGM when large areas of present day shallow sea were exposed, and this parameter is poorly constrained at present.

Egbert, Gary D.↗

NASA NTRS1979

On the influence of orography on large-scale atmospheric flow

The steady response to orography as described by shallow-water equations on the sphere is examined in an attempt to provide insight into the dynamical effects of large-scale orographic features on atmospheric motion. The model equations and the zonal flows and orography used in the study are described. The results for simple mountains and for the earth orography are given. The two-dimensional nature of the horizontal propagation on the sphere is emphasized. The results give interesting indications of the regions of influence of mountains and suggest that quantitative theories of the stationary waves must involve a full representation of the spherical domain.

Grose, W. L.↗

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 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 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.↗