Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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

17 records
All sourcesNASA NTRSDOE OSTI
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 NTRS1985

A Lagrangian method for the shallow water equations based on a Voronoi mesh - Flows on a rotating sphere

A Lagrangian scheme using the Voronoi mesh is applied to study shallow water flow on a sphere. Discrete approximations to the shallow water equations are obtained for the surfaces of a nonrotating and a rotating sphere, and discrete differential operators are defined for the gradient and the divergence on the sphere. Dissipation is put into the model, when needed, by merging fluid points when they get too close to each other. The full numerical scheme is described and results of numerical computations on various test cases are given, including zonal flow and the Riemann problem.

Augenbaum, J. M.↗
NASA NTRS1989

Semi-Lagrangian integration of a gridpoint shallow water model on the sphere

A stable, semi-Lagrangian, semi-implicit, two-time-level, gridpoint integration scheme for the shallow water equations on the sphere is presented. A rotated spherical coordinate system is used to integrate the equations of motion at each gridpoint poleward of a certain latitude, thus overcoming problems associated with the polar singularity. The results of medium term integrations of large scale test patterns using a long time step are presented.

Mcdonald, A.↗
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 NTRS1985

A fully implicit scheme for global numerical weather prediction

A fast-slow factored scheme is presented for use with shallow-water primitive equation numerical weather prediction models. The technique was developed to reduce the rotational mode errors which arise when the fast and slow terms of the governing differential equations are treated simultaneously. The method factors out the fast and slow terms along the coordinate directions by means of a modified Crank-Nicolson scheme. A finite-difference spatial discretization is carried out in the zonal and meridional directions to reduce the factorization error to near-zero, and that time steps of 60-90 min can be used to obtain acceptably accurate results, even in the presence of fine spatial structures in the flow.

Augenbaum, J. M.↗
NASA NTRS1988

Semi-Lagrangian integration of a grid-point shallow water model on the sphere

This paper describes a semi-Lagrangian technique for integrating the equations of motion on the global domain. The technique uses an auxiliary spherical coordinate system at each near-polar gridpoint of the latitude-longitude grid; the auxiliary system is obtained by a rotation such that the new equator passes through the gridpoint in question and the new coordinate directions coincide with those of the original system at that point. The technique was applied to the shallow water equations, incorporating a semiimplicit treatment of the adjustment terms on a C-grid, with two-time levels. A five day integration was successfully carried out for a situation involving strong cross-polar flow. No filtering or diffusion was required to maintain stability over a five day period.

Mcdonald, A.↗
NASA NTRS1973

Amplification and decay of long nonlinear waves.

The interaction of weakly nonlinear waves with slowly varying boundaries is considered. Special emphasis is given to rotating fluids, but the analysis applies with minor modifications to waves in stratified fluids and shallow-water waves. An asymptotic solution of a variant of the Korteweg-de Vries equation with variable coefficients is developed that produces a 'Green's law' for the amplification of waves of finite amplitude. For shallow-water waves in water of variable depth, the result predicts wave growth proportional to the -1/3 power of the depth.

Leibovich, S.↗
NASA NTRS1991

Stability of propagating modons for small-amplitude perturbations

A stability analysis based on Arnol'd-Liapunov arguments in a gauge-variable formalism is applied to a generalized form of propagating modons. The method was previously applied to the special case of Stern's modon in a quiescent background flow. The analysis presented here shows that propagating-modon solutions in the shallow-water equations are stable to small-amplitude perturbations, regardless of the sign of their translational speed, as long as this is much smaller in magnitude than the solid-rotation speed of the earth. Sufficient conditions for stability are that the given flow be quasigeostrophic, i.e., that the Rossby number be small, and that the associated Froude number be one order of magnitude smaller than the Rossby number.

Sakuma, H.↗
NASA NTRS2018

Symmetric Equations on the Surface of a Sphere as Used by Model GISS:IB

Standard vector calculus formulas of Cartesian three space are projected onto the surface of a sphere. This produces symmetric equations with three nonindependent horizontal velocity components. Each orthogonal axis has a velocity component that rotates around its axis (eastward velocity rotates around the north–south axis) and a specific angular momentum component that is the product of the velocity component multiplied by the cosine of axis’ latitude. Angular momentum components align with the fixed axes and simplify several formulas, whereas the rotating velocity components are not orthogonal and vary with location. Three symmetric coordinates allow vector resolution and calculus operations continuously over the whole spherical surface, which is not possible with only two coordinates. The symmetric equations are applied to one-layer shallow water models on cubed-sphere and icosahedral grids, the latter being computationally simple and applicable to an ocean domain. Model results are presented for three different initial conditions and five different resolutions.

Icosahedral grid↗
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 NTRS1990

Integration of the shallow water equations on the sphere using a vector semi-Lagrangian scheme with a multigrid solver

A vector semi-Lagrangian semi-implicit two-time-level finite-difference integration scheme for the shallow water equations on the sphere is presented. A C-grid is used for the spatial differencing. The trajectory-centered discretization of the momentum equation in vector form eliminates pole problems and, at comparable cost, gives greater accuracy than a previous semi-Lagrangian finite-difference scheme which used a rotated spherical coordinate system. In terms of the insensitivity of the results to increasing timestep, the new scheme is as successful as recent spectral semi-Lagrangian schemes. In addition, the use of a multigrid method for solving the elliptic equation for the geopotential allows efficient integration with an operation count which, at high resolution, is of lower order than in the case of the spectral models. The properties of the new scheme should allow finite-difference models to compete with spectral models more effectively than has previously been possible.

Bates, J. R.↗
NASA NTRS1979

Explicit large time-step schemes for the shallow water equations

Modifications to explicit finite difference schemes for solving the shallow water equations for meteorological applications by increasing the time step for the fast gravity waves are analyzed. Terms associated with the gravity waves in the shallow water equations are treated on a coarser grid than those associated with the slow Rossby waves, which contain much more of the available energy and must be treated with higher accuracy, enabling a several-fold increase in time step without degrading the accuracy of the solution. The method is presented in Cartesian and spherical coordinates for a rotating earth, using generalized leapfrog, frozen coefficient, and Fourier filtering finite difference schemes. Computational results verify the numerical stability of the approach.

Turkel, E.↗
NASA NTRS2003

Deviation of Long-Period Tides from Equilibrium: Kinematics and Geostrophy

New empirical estimates of the long-period fortnightly (Mf) tide obtained from TOPEX/Poseidon (T/P) altimeter data confirm significant basin-scale deviations from equilibrium. Elevations in the low-latitude Pacific have reduced amplitude and lag those in the Atlantic by 30 deg or more. These interbasin amplitude and phase variations are robust features that are reproduced by numerical solutions of the shallow-water equations, even for a constant-depth ocean with schematic interconnected rectangular basins. A simplified analytical model for cooscillating connected basins also reproduces the principal features observed in the empirical solutions. This simple model is largely kinematic. Zonally averaged elevations within a simple closed basin would be nearly in equilibrium with the gravitational potential, except for a constant offset required to conserve mass. With connected basins these offsets are mostly eliminated by interbasin mass flux. Because of rotation, this flux occurs mostly in a narrow boundary layer across the mouth and at the western edge of each basin, and geostrophic balance in this zone supports small residual offsets (and phase shifts) between basins. The simple model predicts that this effect should decrease roughly linearly with frequency, a result that is confirmed by numerical modeling and empirical T/P estimates of the monthly (Mm) tidal constituent. This model also explains some aspects of the anomalous nonisostatic response of the ocean to atmospheric pressure forcing at periods of around 5 days.

Egbert, Gary D.↗
NASA NTRS2016

Traveling Weather Disturbances in Mars Southern Extratropics: Sway of the Great Impact Basins

As on Earth, between late autumn and early spring on Mars middle and high latitudes within its atmosphere support strong mean thermal contrasts between the equator and poles (i.e. "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that this strong baroclinicity supports vigorous, large-scale eastward traveling weather systems (i.e. transient synoptic-period waves). Within a rapidly rotating, differentially heated, shallow atmosphere such as on Earth and Mars, such large-scale, extratropical weather disturbances are critical components of the global circulation. These wave-like disturbances act as agents in the transport of heat and momentum, and moreover generalized tracer quantities (e.g., atmospheric dust, water vapor and water-ice clouds) between low and high latitudes of the planet. The character of large-scale, traveling extratropical synoptic-period disturbances in Mars' southern hemisphere during late winter through early spring is investigated using a high-resolution Mars global climate model (Mars GCM). This global circulation model imposes interactively lifted (and radiatively active) dust based on a threshold value of the instantaneous surface stress. Compared to observations, the model exhibits a reasonable "dust cycle" (i.e. globally averaged, a more dusty atmosphere during southern spring and summer occurs). In contrast to their northern-hemisphere counterparts, southern synoptic-period weather disturbances and accompanying frontal waves have smaller meridional and zonal scales, and are far less intense synoptically. Influences of the zonally asymmetric (i.e. east-west varying) topography on southern large-scale weather disturbances are examined. Simulations that adapt Mars' full topography compared to simulations that utilize synthetic topographies emulating essential large-scale features of the southern middle latitudes indicate that Mars' transient barotropic/baroclinic eddies are significantly influenced by the great impact basins of this hemisphere (e.g., Argyre and Hellas). In addition, the occurrence of a southern storm zone in late winter and early spring is keyed particularly to the western hemisphere via orographic influences arising from the Tharsis highlands, and the Argyre and Hellas impact basins. Geographically localized transient-wave activity diagnostics are constructed that illuminate fundamental differences amongst such simulations and these are described.

Weather disturbances↗
NASA NTRS2017

Extratropical Weather Systems on Mars: Radiatively-Active Water Ice Effects

Extratropical, large-scale weather disturbances, namely transient, synoptic-period,baroclinic barotropic eddies - or - low- (high-) pressure cyclones (anticyclones), are components fundamental to global circulation patterns for rapidly rotating, differentially heated, shallow atmospheres such as Earth and Mars. Such "wave-like" disturbances that arise via (geophysical) fluid shear instability develop, mature and decay, and travel west-to-east in the middle and high latitudes within terrestrial-like planetary atmospheres. These disturbances serve as critical agents in the transport of heat and momentum between low and high latitudes of the planet. Moreover, they transport trace species within the atmosphere (e.g., water vapor/ice, other aerosols (dust), chemical species, etc). Between early autumn through early spring, middle and high latitudes on Mars exhibit strong equator-to-pole mean temperature contrasts (i.e., "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that such strong baroclinicity supports vigorous, large-scale eastward traveling weather systems [Banfield et al., 2004; Barnes et al., 1993]. A good example of traveling weather systems, frontal wave activity and sequestered dust activity from MGS/MOC image analyses is provided in Figure 1 (cf. Wang et al. [2005]). Utilizing an upgraded and evolving version of the NASA Ames Research Center (ARC) Mars global climate model, investigated here are key dynamical and physical aspects of simulated northern hemisphere (NH) large-scale extratropica lweather systems,with and without radiatively-active water ice clouds. Mars Climate Model:

cyclones↗
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 NTRS2014

Modeling Mars Cyclogenesis and Frontal Waves: Seasonal Variations and Implications on Dust Activity

Between late autumn through early spring,middle and high latitudes onMars exhibit strong equator-to-polemean temperature contrasts (i.e., "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that such strong baroclinicity supports vigorous, large-scale eastward traveling weather systems (i.e., transient synoptic period waves) [1, 2]. For a rapidly rotating, differentially heated, shallow atmosphere such as on Earth and Mars, these large-scale, extratropical weather disturbances are critical components of the global circulation. The wave-like disturbances serve as agents in the transport of heat and momentum between low and high latitudes of the planet. Through cyclonic/anticyclonic winds, intense shear deformations, contractions-dilatations in temperature and density, and sharp perturbations amongst atmospheric tracers (i.e., dust, volatiles (e.g., water vapor) and condensates (e.g., water-ice cloud particles)), Mars' extratropical weather systems have significant sub-synoptic scale ramifications by supporting atmospheric frontal waves (Fig. 1).

Mars Global Surveyor (MGS)↗