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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

SWIM: A Semi-Analytical Ocean Color Inversion Algorithm for Optically Shallow Waters

Ocean color remote sensing provides synoptic-scale, near-daily observations of marine inherent optical properties (IOPs). Whilst contemporary ocean color algorithms are known to perform well in deep oceanic waters, they have difficulty operating in optically clear, shallow marine environments where light reflected from the seafloor contributes to the water-leaving radiance. The effect of benthic reflectance in optically shallow waters is known to adversely affect algorithms developed for optically deep waters [1, 2]. Whilst adapted versions of optically deep ocean color algorithms have been applied to optically shallow regions with reasonable success [3], there is presently no approach that directly corrects for bottom reflectance using existing knowledge of bathymetry and benthic albedo.To address the issue of optically shallow waters, we have developed a semi-analytical ocean color inversion algorithm: the Shallow Water Inversion Model (SWIM). SWIM uses existing bathymetry and a derived benthic albedo map to correct for bottom reflectance using the semi-analytical model of Lee et al [4]. The algorithm was incorporated into the NASA Ocean Biology Processing Groups L2GEN program and tested in optically shallow waters of the Great Barrier Reef, Australia. In-lieu of readily available in situ matchup data, we present a comparison between SWIM and two contemporary ocean color algorithms, the Generalized Inherent Optical Property Algorithm (GIOP) and the Quasi-Analytical Algorithm (QAA).

QAA↗
NASA NTRS2012

Assessing Tsunami Vulnerabilities of Geographies with Shallow Water Equations

Tsunami preparedness is crucial for saving human lives in case of disasters that involve massive water movement. In this work, we develop a framework for visual assessment of tsunami preparedness of geographies. Shallow water equations (also called Saint Venant equations) are a set of hyperbolic partial differential equations that are derived by depth-integrating the Navier-Stokes equations and provide a great abstraction of water masses that have lower depths compared to their free surface area. Our specific contribution in this study is to use Microsoft's XNA Game Studio to import underwater and shore line geographies, create different tsunami scenarios, and visualize the propagation of the waves and their impact on the shore line geography. Most importantly, we utilized the computational power of graphical processing units (GPUs) as HLSL based shader files and delegated all of the heavy computations to the GPU. Finally, we also conducted a validation study, in which we have tested our model against a controlled shallow water experiment. We believe that such a framework with an easy to use interface that is based on readily available software libraries, which are widely available and easily distributable, would encourage not only researchers, but also educators to showcase ideas.

Aras, Rifat↗
DOE OSTI2024

A high-order explicit Runge-Kutta approximation technique for the shallow water equations

Here, we introduce a high-order space–time approximation of the Shallow Water Equations with sources that is invariant-domain preserving (IDP), well-balanced with respect to rest states, and employs a novel explicit Runge–Kutta (ERK) introduced in Ern and Guermond (SIAM J. Sci. Comput. 44(5), A3366–A3392, 2022) for systems of non-linear conservation equations. The resulting method is then numerically illustrated through verification and validation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Energy conserving and potential-enstrophy dissipating schemes for the shallow water equations

To incorporate potential enstrophy dissipation into discrete shallow water equations with no or arbitrarily small energy dissipation, a family of finite-difference schemes have been derived with which potential enstrophy is guaranteed to decrease while energy is conserved (when the mass flux is nondivergent and time is continuous). Among this family of schemes, there is a member that minimizes the spurious impact of infinite potential vorticities associated with infinitesimal fluid depth. The scheme is, therefore, useful for problems in which the free surface may intersect with the lower boundary.

Arakawa, Akio↗
DOE OSTI2023

Shallow-Water Habitat in the Lower Columbia River Estuary: A Highly Altered System

Abstract Decreases in shallow-water habitat area (SWHA) in the Lower Columbia River and Estuary (LCRE) have adversely affected salmonid populations. We investigate the causes by hindcasting SWHA from 1928 to 2004, system-wide, based on daily higher high water (HHW) and system hypsometry. Physics-based regression models are used to represent HHW along the system as a function of river inflow, tides, and coastal processes, and hypsometry is used to estimate the associated SWHA. Scenario modeling is employed to attribute SWHA losses to levees, flow regulation, diversion, navigational development, and climate-induced hydrologic change, for subsidence scenarios of up to 2 m, and for 0.5 m fill. For zero subsidence, the system-wide annual-average loss of SWHA is 55 ± 5%, or 51 × 10 5 ha/year; levees have caused the largest decrease ( $${54}_{-14}^{+5}$$ 54 - 14 + 5 %, or ~ 50 × 10 5 ha/year). The loss in SWHA due to operation of the hydropower system is small, but spatially and seasonally variable. During the spring freshet critical to juvenile salmonids, the total SWHA loss was $${63}_{-3}^{+2}$$ 63 - 3 + 2 %, with the hydropower system causing losses of 5–16% (depending on subsidence). Climate change and navigation have caused SWHA losses of $${5}_{-5}^{+16}$$ 5 - 5 + 16 % and $${4}_{-6}^{+14}$$ 4 - 6 + 14 %, respectively, but with high spatial variability; irrigation impacts have been small. Uncertain subsidence causes most of the uncertainty in estimates; the sum of the individual factors exceeds the total loss, because factors interact. Any factor that reduces mean or peak flows (reservoirs, diversion, and climate change) or alters tides and along-channel slope (navigation) becomes more impactful as assumed historical elevations are increased to account for subsidence, while levees matter less.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2003

Strangelove Ocean and Deposition of Unusual Shallow-Water Carbonates After the End-Permian Mass Extinction

The severe mass extinction of marine and terrestrial organisms at the end of the Permian Period (approx. 251 Ma) was accompanied by a rapid negative excursion of approx. 3 to 4 per mil in the carbon-isotope ratio of the global surface oceans and atmosphere that persisted for some 500,000 into the Early Triassic. Simulations with an ocean-atmosphere/carbon-cycle model suggest that the isotope excursion can be explained by collapse of ocean primary productivity (a Strangelove Ocean) and changes in the delivery and cycling of carbon in the ocean and on land. Model results also suggest that perturbations of the global carbon cycle resulting from the extinctions led to short-term fluctuations in atmospheric pCO2 and ocean carbonate deposition, and to a long-term (>1 Ma) decrease in sedimentary burial of organic carbon in the Triassic. Deposition of calcium carbonate is a major sink of river-derived ocean alkalinity and for CO2 from the ocean/atmosphere system. The end of the Permian was marked by extinction of most calcium carbonate secreting organisms. Therefore, the reduction of carbonate accumulation made the oceans vulnerable to a build-up of alkalinity and related fluctuations in atmospheric CO2. Our model results suggest that an increase in ocean carbonate-ion concentration should cause increased carbonate accumulation rates in shallow-water settings. After the end-Permian extinctions, early Triassic shallow-water sediments show an abundance of abiogenic and microbial carbonates that removed CaCO3 from the ocean and may have prevented a full 'ocean-alkalinity crisis' from developing.

Rampino, Michael R.↗
DOE OSTI2022

A mass conservative, well balanced, tangency preserving and energy decaying method for the shallow water equations on a sphere

Here a fully discrete surface finite element method is proposed for solving the viscous shallow water equations in a bounded Lipschitz domain on the sphere based on a general triangular mesh. The method consists of a modified Crank–Nicolson method in time and a Galerkin surface finite element method in space for the fluid thickness H and the fluid velocity u. A finite element space tangential to the sphere at all finite element nodes is proposed to approximate the fluid velocity u. The proposed method has second-order accuracy in time and first-order accuracy in space, and preserves mass conservation, well balancedness, tangency of velocity to the sphere, and energy decay. Numerical experiments are presented to illustrate the accuracy of the proposed method and the preservation of the physical properties, including mass conservation, well balancedness, and energy decay. A numerical simulation of ocean mesoscale activity on a circular basin with a continental shelf is provided.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Generalized energy and potential enstrophy conserving finite difference schemes for the shallow water equations

The conditions under which finite difference schemes for the shallow water equations can conserve both total energy and potential enstrophy are considered. A method of deriving such schemes using operator formalism is developed. Several such schemes are derived for the A-, B- and C-grids. The derived schemes include second-order schemes and pseudo-fourth-order schemes. The simplest B-grid pseudo-fourth-order schemes are presented.

Abramopoulos, Frank↗
DOE OSTI2023

SERGHEI (SERGHEI-SWE) v1.0: a performance-portable high-performance parallel-computing shallow-water solver for hydrology and environmental hydraulics

The Simulation EnviRonment for Geomorphology, Hydrodynamics, and Ecohydrology in Integrated form (SERGHEI) is a multi-dimensional, multi-domain, and multi-physics model framework for environmental and landscape simulation, designed with an outlook towards Earth system modelling. At the core of SERGHEI's innovation is its performance-portable high-performance parallel-computing (HPC) implementation, built from scratch on the Kokkos portability layer, allowing SERGHEI to be deployed, in a performance-portable fashion, in graphics processing unit (GPU)-based heterogeneous systems. In this work, we explore combinations of MPI and Kokkos using OpenMP and CUDA backends. In this contribution, we introduce the SERGHEI model framework and present with detail its first operational module for solving shallow-water equations (SERGHEI-SWE) and its HPC implementation. This module is designed to be applicable to hydrological and environmental problems including flooding and runoff generation, with an outlook towards Earth system modelling. Its applicability is demonstrated by testing several well-known benchmarks and large-scale problems, for which SERGHEI-SWE achieves excellent results for the different types of shallow-water problems. Finally, SERGHEI-SWE scalability and performance portability is demonstrated and evaluated on several TOP500 HPC systems, with very good scaling in the range of over 20 000 CPUs and up to 256 state-of-the art GPUs.

58 GEOSCIENCES↗
DOE OSTI2021

Local time stepping for the shallow water equations in MPAS

In this work we assess the performance of a set of local time-stepping (LTS) schemes for the shallow water equations implemented in the Model for Prediction Across Scales (MPAS). The goal of LTS is to speed up the simulation by allowing different time-steps on different regions of the computational grid. The LTS schemes considered here were originally introduced by Hoang et al. (2019) [26], who laid out the mathematical foundation of the methods. Here, the authors take on the task of presenting a fast, efficient and scalable parallel implementation of these LTS methods on high performance computing machines, with the aim to provide a recipe for other climate modeling groups that may be interested in employing LTS algorithms in their codes. As a matter of fact, even if MPAS is our framework of choice, our approach is general enough and could be of interest to other groups beyond the MPAS community. Due to their nature, LTS methods possess an inherent load imbalance that needs to be carefully addressed in order to obtain efficient scalability. Even more important is the far from trivial task of computing the right-hand side terms only on specific LTS regions during the time-stepping procedure. An inefficient handling of this task causes a drastic decay of the CPU time performance, making the LTS algorithms practically of no use. The emphasis of the present work is therefore on the computational and parallel aspects of the LTS methods, whose proper treatment is crucial to make the methods run faster against existing strategies, such as for instance high-order explicit global time-stepping schemes. This is in fact the ultimate goal of using an LTS procedure and it is the one to which we direct all our optimization efforts.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

An energy and potential enstrophy conserving scheme for the shallow water equations

A three-dimensional finite difference scheme for the solution of the shallow water momentum equations which accounts for the conservation of potential enstrophy in the flow of a homogeneous incompressible shallow atmosphere over steep topography as well as for total energy conservation is presented. The scheme is derived to be consistent with a reasonable scheme for potential vorticity advection in a long-term integration for a general flow with divergent mass flux. Numerical comparisons of the characteristics of the present potential enstrophy-conserving scheme with those of a scheme that conserves potential enstrophy only for purely horizontal nondivergent flow are presented which demonstrate the reduction of computational noise in the wind field with the enstrophy-conserving scheme and its convergence even in relatively coarse grids.

Arakawa, A.↗
NASA NTRS1986

Using Hough harmonics to validate and assess nonlinear shallow-water models

The implementation of a technique for locating programming errors in shallow-water codes, establishing the correctness of the code, and assessing the performance of the numerical model under various flow conditions is described. The right-hand side of the differential equations is modified in such a way that the exact solution of the nonlinear initial-value problem is known, so that the truncation errors of the numerical scheme can be studied in detail. The exact solution is prescribed to be any linear combination of Hough harmonics which propagate in time according to their natural frequencies.

Dee, Dick P.↗
NASA NTRS1989

Solving the shallow water equations on the Cray X-MP/48 and the connection machine 2

The shallow water equations in Cartesian coordinates and 2-D are solved on the Connection Machine 2 (CM-2) using both the spectral and finite difference methods. A description of these implementations is presented together with a brief discussion of the CM-2 as it relates to these specific computations. The finite difference code was written both in C* and *LISP and the spectral code was written in *LISP. The performance of the codes is compared with a FORTRAN version that was optimized for the Cray X-MP/48.

Swarztrauber, Paul N.↗
NASA NTRS1986

Documentation of the Goddard Laboratory for atmospheres fourth-order two-layer shallow water model

The theory and numerical treatment used in the 2-level GLA fourth-order shallow water model are described. This model was designed to emulate the horizontal finite differences used by the GLA Fourth-Order General Circulation Model (Kalnay et al., 1983) in addition to its grid structure, form of high-latitude and global filtering, and time-integration schemes. A user's guide is also provided instructing the user on how to create initial conditions, execute the model, and post-process the data history.

Takacs, L. L.↗
NASA NTRS1998

Shallow Water Optical Water Quality Buoy

This NASA grant was funded as a result of an unsolicited proposal submission to Kennedy Space Center. The proposal proposed the development and testing of a shallow water optical water quality buoy. The buoy is meant to work in shallow aquatic systems (ponds, rivers, lagoons, and semi-enclosed water areas where strong wind wave action is not a major environmental During the project period of three years, a demonstration of the buoy was conducted. The last demonstration during the project period was held in November, 1996 when the buoy was demonstrated as being totally operational with no tethered communications line. During the last year of the project the buoy was made to be solar operated by large gel cell batteries. Fund limitations did not permit the batteries in metal enclosures as hoped for higher wind conditions, however the system used to date has worked continuously for in- situ operation of over 18 months continuous deployment. The system needs to have maintenance and somewhat continuous operational attention since various components have limited lifetime ages. For example, within the last six months the onboard computer has had to be repaired as it did approximately 6 months after deployment. The spectrograph had to be repaired and costs for repairs was covered by KB Science since no ftmds were available for this purpose after the grant expired. Most recently the computer web page server failed and it is currently being repaired by KB Science. In addition, the cell phone operation is currently being ftmded by Dr. Bostater in order to maintain the system's operation. The above points need to be made to allow NASA to understand that like any sophisticated measuring system in a lab or in the field, necessary funding and maintenance is needed to insure the system's operational state and to obtain quality factor. The proposal stated that the project was based upon the integration of a proprietary and confidential sensor and probe design that was developed by KB Science and Engineering and is currently patented by KB Science. The buoy's purpose was to collected hyperspectral optical signatures for analysis and resulting estimation of water quality parameters such as chlorophyll-a, seston and dissolved organic matter (DOC). The ultimate goal of the project was to develop a buoy that would integrate a probe to measure upwelling light from a source and thus relate this backscattered light to water quality parameters.

Bostater, Charles↗
NASA NTRS2017

Comparison of Damage from Hydrocode Simulations of an Asteroid Airburst or Impact on Land, in Deep, or in Shallow Water

If an asteroid is discovered to be on a collision course with Earth and there is insufficient time for a deflection effort to make it miss Earth completely, should it be redirected to a land or ocean impact? While distance from densely populated areas should obviously be maximized, the differing ability of air blast, seismic waves, and tsunami waves to cause damage at distance does not make the choice between land and ocean impacts an immediately obvious one. More broadly this work is a step towards improving damage models from asteroid impacts. This extended abstract follows the hypothetical scenario of the 2017 IAA Planetary Defense Conference where a 100-250m diameter asteroid is on a potential impact course with Earth. A hydrocode was used to simulate impacts into the most sparsely populated areas along the eastern end of the hypothetical impact corridor- specifically in the Gobi Desert, in the shallow waters of the Sea of Japan, and in the deep waters of the Japan Trench in the Pacific Ocean.

Shallow Water↗