Engineering Papers⌕ Search

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.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗
NASA NTRS1999

Data Assimilation Into a Coupled Ocean Atmosphere Model: Application to the 1997-1998 El Nino

As part of JPL's ocean data assimilation effort to study ocean circulation and seasonal-interannual climate variability, sea level anomaly observed by TOPEX altimeter, together with sea surface temperature and wind stress data, are assimilated into a simple coupled ocean atmosphere model of the tropical Pacific. Model-data consistency is examined. Impact of the assimilation (as initialization) on El Nino Southern Oscillation (ENSO) forecasts is evaluated. The coupled model consists of a shallow water component with two baroclinic modes, an Ekman shear layer, a simplified mixed-layer temperature equation, and a statistical atmosphere based on dominant correlations between historical surface temperature and wind stress anomaly data. The adjoins method is used to fit the coupled model to the data over various six-month periods from late 1996 to early 1998 by optimally adjusting the initial state, model parameters, and basis functions of the statistical atmosphere. On average, the coupled model can be fitted to the data to approximately within the data and representation errors (5 cm, 0.5 C, and 10 sq m/sq m for sea level, surface temperature, and pseudo wind stress anomalies, respectively). The estimated fields resemble observed spatio-temporal structure reasonably well. Hindcasts/forecasts of the 1997/1998 El Nino initialized from forced estimated ocean states and parameters are much more realistic than those simply initialized from ocean states (see figure below). In particular, the ability of the model to produce significant warming beyond the initial state is dramatically improved. Parameter estimation, which compensates for some model errors, is found to be important to obtaining better fits of the model to data and to improving forecasts.

Lee, Tong↗
DOE OSTI2022

Improving the accuracy of discretisations of the vector transport equation on the lowest-order quadrilateral Raviart-Thomas finite elements

Within finite element models of fluids, vector-valued fields such as velocity or momentum variables are commonly discretised using the Raviart-Thomas elements. However, when using the lowest-order quadrilateral Raviart-Thomas elements, standard finite element discretisations of the vector transport equation typically have a low order of spatial accuracy. This paper describes two schemes that improve the accuracy of transporting such vector-valued fields on two-dimensional curved manifolds. The first scheme that is presented reconstructs the transported field in a higher-order function space, where the transport equation is then solved. The second scheme applies a mixed finite element formulation to the vector transport equation, simultaneously solving for the transported field and its vorticity. In this work, an approach to stabilising this mixed vector-vorticity formulation is presented that uses a Streamline Upwind Petrov-Galerkin (SUPG) method. These schemes are then demonstrated, along with their accuracy properties, through some numerical tests. Two new test cases are used to assess the transport of vector-valued fields on curved manifolds, solving the vector transport equation in isolation. The improvement of the schemes is also shown through two standard test cases for rotating shallow-water models.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Evaluation of HCMM data for assessing soil moisture and water table depth

Data were analyzed for variations in eastern South Dakota. Soil moisture in the 0-4 cm layer could be estimated with 1-mm soil temperatures throughout the growing season of a rainfed barley crop (% cover ranging from 30% to 90%) with an r squared = 0.81. Empirical equations were developed to reduce the effect of canopy cover when radiometrically estimating the 1-mm soil temperature, r squared = 0.88. The corrective equations were applied to an aircraft simulation of HCMM data for a diversity of crop types and land cover conditions to estimate the 0-4 cm soil moisture. The average difference between observed and measured soil moisture was 1.6% of field capacity. HCMM data were used to estimate the soil moisture for four dates with an r squared = 0.55 after correction for crop conditions. Location of shallow alluvial aquifers could be accomplished with HCMM predawn data. After correction of HCMM day data for vegetation differences, equations were developed for predicting water table depths within the aquifer (r=0.8).

Moore, D. G.↗
DOE OSTI2024

Potential vorticity conservation for plasma turbulence in an inhomogeneous magnetic field: Theory and implications

The concept and theory of potential vorticity in drift wave turbulence are extended to the case of an inhomogeneous magnetic field. A one-field magnetic potential vorticity conserving equation is derived via the use of conservative gyrokinetics. The similarity between the corresponding systems for drift wave turbulence and shallow water theory is discussed in detail. Zonal flow physics in an inhomogeneous magnetic field is discussed. In particular, a Charney–Drazin type nonacceleration theorem is derived from the novel system, which conserves magnetic potential vorticity. Extensions of the turbulent equipartition theory to the transport of magnetic potential vorticity are proposed.

Physics↗
NASA NTRS1984

A Lagrangian method for a two layer shallow water model

One of the biggest drawbacks of Eulerian methods is the presence of the nonlinear convective terms which leads to inaccurate representations of advection and fronts. An alternative and conceptually simpler approach is the Lagrangian approach in which the fluid particles themselves are tracked and equations are derived based on local spatial interaction. In this formulation the nonlinear convective terms do not appear. The main drawback with the Lagrangian formulation, however, is that the local spatial interactions are time dependent. That is, at each time level, one must know the current neighbors of a given particle to accurately compute the forces acting on the particle. The Voroni Mesh which is continuously deforming and always linking nearest neighbors is needed. The effectiveness of a simple Lagrangian method for simulating the onset of a Kelvin-Helmholtz instability in a 2 layer fluid is demonstrated.

Augenbaum, J. M.↗
NASA NTRS1997

Evolution of Nonlinear Internal Waves in China Seas

Synthetic Aperture Radar (SAR) images from ERS-I have been used to study the characteristics of internal waves of Taiwan in the East China Sea, and east of Hainan Island in the South China Sea. Rank-ordered packets of internal solitons propagating shoreward from the edge of the continental shelf were observed in the SAR images. Based on the assumption of a semidiurnal tidal origin, the wave speed can be estimated and is consistent with the internal wave theory. By using the SAR images and hydrographic data, internal waves of elevation have been identified in shallow water due to a thicker mixed layer as compared with the bottom layer on the continental shelf. The generation mechanism includes the influences of the tide and the Kuroshio intrusion across the continental shelf for the formations of elevation internal waves. The effects of water depth on the evolution of solitons and wave packets are modeled by nonlinear Kortweg-deVries (KdV) type equation and linked to satellite image observations. The numerical calculations of internal wave evolution on the continental shelf have been performed and compared with the SAR observations. For a case of depression waves in deep water, the solitons first disintegrate into dispersive wave trains and then evolve to a packet of elevation waves in the shallow water area after they pass through a turning point of approximately equal layer depths has been observed in the SAR image and simulated by numerical model.

Liu, Antony K.↗
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 NTRS1983

Evaluation of spatial, radiometric and spectral thematic mapper performance for coastal studies

An area along the southeastern shore of the Chesapeake Bay was subsetted from TM imagery. The subsetted image was then enhanced and classified using an ERDAS 400 system. Results obtained were compared with a chart showing the distribution of both Zolsters marina and Rupplia martime in the Vaucluse Shores and which supports a large community of SAV. Radiative transfer models describing the irradiance reflectance of a water column containing SAV are being refined. Radiative transfer theory was used to model upwelling radiance for an orbiting sensor viewing an estuarine environment. Upwelling radiance was calculated for a clear maritime atmosphere, an optically shallow estuary of either clear or turbid water, and one of three bottom types: vegetation, sand, or mud using TM bands 1, 2, and 3 and MSS bands 4 and 5. A spectral quality index was defined similar to the equation for apparent contrast and used to evaluate the relative effectiveness of TM and MSS bands in detecting submerged vegetation.

Klemas, V.↗
NASA NTRS2016

Evidence for a Heterogeneous Distribution of Water in the Martian Interior

The abundance and distribution of H2O within the terrestrial planets, as well as its timing of delivery, is a topic of vital importance for understanding the chemical and physical evolution of planets and their potential for hosting habitable environments. Analysis of planetary materials from Mars, the Moon, and the eucrite parent body (i.e., asteroid 4Vesta) have confirmed the presence of H2O within their interiors. Moreover, H and N isotopic data from these planetary materials suggests H2O was delivered to the inner solar system very early from a common source, similar in composition to the carbonaceous chondrites. Despite the ubiquity of H2O in the inner Solar System, the only destination with any prospects for past or present habitable environments at this time, outside of the Earth, is Mars. Although the presence of H2O within the martian interior has been confirmed, very little is known regarding its abundance and distribution within the martian interior and how the martian water inventory has changed over time. By combining new analyses of martian apatites within a large number of martian meteorite types with previously published volatile data and recently determined mineral-melt partition coefficients for apatite, we report new insights into the abundance and distribution of volatiles in the martian crust and mantle. Using the subset of samples that did not exhibit crustal contamination, we determined that the enriched shergottite mantle source has 36-73 ppm H2O and the depleted shergottite mantle source has 14-23 ppm H2O. This result is consistent with other observed geochemical differences between enriched and depleted shergottites and supports the idea that there are at least two geochemically distinct reservoirs in the martian mantle. We also estimated the H2O content of the martian crust using the revised mantle H2O abundances and known crust-mantle distributions of incompatible lithophile elements. We determined that the bulk martian crust has approximately 1400 ppm H2O, which is likely distributed toward the martian surface. This crustal water abundance would equate to a global equivalent layer (GEL) of water at a depth of-229 m, which can account for at least some of the surface features on Mars attributed to flowing water and may be sufficient to support the past presence of a shallow sea on Mars' surface.

McCubbin, Francis↗
NASA NTRS1999

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗
NASA NTRS1999

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗
NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗
DOE OSTI2023

Operator inference with roll outs for learning reduced models from scarce and low-quality data

Data-driven modeling has become a key building block in computational science and engineering. However, data that are available in science and engineering are typically scarce, often polluted with noise and affected by measurement errors and other perturbations, which makes learning the dynamics of systems challenging. Here, in this work, we propose to combine data-driven modeling via operator inference with the dynamic training via roll outs of neural ordinary differential equations. Operator inference with roll outs inherits interpretability, scalability, and structure preservation of traditional operator inference while leveraging the dynamic training via roll outs over multiple time steps to increase stability and robustness for learning from low-quality and noisy data. Numerical experiments with data describing shallow water waves and surface quasi-geostrophic dynamics demonstrate that operator inference with roll outs provides predictive models from training trajectories even if data are sampled sparsely in time and polluted with noise of up to 10%.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Experiments on strong interactions between solitary waves

Experiments on the interaction between solitary shallow-water waves propagating in the same direction have been performed in a rectangular channel. Two methods were devised to compensate for the dissipation of the waves in order to compare results with Hirota's (1971) solution for the collision of solitons described by the Korteweg-de Vries equation. Both qualitative and quantitative agreement with theory is obtained using the proposed corrections for wave damping.

Weidman, P. D.↗
DOE OSTI2025

The diurnal susceptibility of subtropical clouds to aerosols

The diurnal susceptibility of clouds and their radiative properties to aerosols is examined during their Lagrangian transition from subtropical stratocumulus to shallow cumulus regimes. Using large-eddy simulations, we analyze the six-day evolution of an air mass along a 3800 km observed trajectory from the coast of Peru toward the equator. Pristine and polluted scenarios are simulated with forcing imposed from weather reanalysis. The polluted scenario exhibits stronger diurnal variations in cloud water, cloud fraction, and albedo, with enhanced nighttime entrainment and suppressed precipitation. The overall response of cloud properties and outgoing shortwave radiation to droplet number concentration follows a distinct diurnal pattern: strong positive cloud adjustments dominate at night and in the morning, while weak negative adjustments prevail in the afternoon. This cycle is driven by the competition between precipitation suppression, which enhances cloud water and coverage, and entrainment drying, which depletes them. In polluted conditions, enhanced entrainment leads to a deeper and more decoupled boundary layer that cannot be sustained by surface fluxes in the afternoon, resulting in negative cloud adjustments. The enhanced entrainment rate under polluted conditions is caused by the reduced sedimentation of cloud and precipitation water from the entrainment zone. While the Twomey effect dominates the diurnal average albedo response, the diurnal variation in the competing cloud adjustments lead to a near-neutral net susceptibility in the afternoon, highlighting the critical role of diurnally varying processes in aerosol-cloud interactions.

Geosciences↗
NASA NTRS1995

Assimilation of TOPEX Sea Level Measurements with a Reduced-Gravity, Shallow Water Model of the Tropical Pacific Ocean

Sea surface height variability measured by TOPEX is analyzed in the tropical Pacific Ocean by way of assimilation into a wind-driven, reduced-gravity, shallow water model using an approximate Kalman filter and smoother. The analysis results in an optimal fit of the dynamic model to the observations, providing it dynamically consistent interpolation of sea level and estimation of the circulation. Nearly 80% of the expected signal variance is accounted for by the model within 20 deg of the equator, and estimation uncertainty is substantially reduced by the voluminous observation. Notable features resolved by the analysis include seasonal changes associated with the North Equatorial Countercurrent and equatorial Kelvin and Rossby waves. Significant discrepancies are also found between the estimate and TOPEX measurements, especially near the eastern boundary. Improvements in the estimate made by the assimilation are validated by comparisons with independent tide gauge and current meter observations. The employed filter and smoother are based on approximately computed estimation error covariance matrices, utilizing a spatial transformation and an symptotic approximation. The analysis demonstrates the practical utility of a quasi-optimal filter and smoother.

Fukumori, Ichiro↗
NASA NTRS2021

PBL Height from AIRS, GPS RO, and MERRA-2 Products in NASA GES DISC and Their 10 Year Seasonal Mean Intercomparison

Within the planetary boundary layer (PBL), surface forcing response, drag, turbulence, and vertical mixing are important processes and play a more critical role here than in the overlying “free atmosphere”. The PBL Height (PBLH) is an important parameter in climate models, weather forecasts, and air quality prediction. The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) provides data processing, archiving, and distribution services for numerous Earth science products. PBLH is a parameter in three products served by GES DISC, which are from the Atmospheric Infrared Sounder (AIRS), the Global Positioning System (GPS) radio occultation (RO) experiment, and the NASA reanalysis product Modern-Era Retrospective analysis for Research and Applications – 2 (MERRA-2). These products have different spatial and temporal resolutions and coverages, and their PBLH definitions are also different. To better serve the PBL research community, we have summarized the specifications of these products. A ten-year seasonal mean intercomparison is also conducted to provide further guidance to users. The intercomparison results show that MERRA-2 has a much shallower PBL than AIRS and GPS RO. An experimental study indicates the different PBLH definition in MERRA-2 caused smaller values of PBLH. The improvement of the water vapor retrieval in AIRS version 7 over version 6 results in the version 7 PBLH agreeing better with GPS RO and MERRA-2 than version 6, especially near the equator and low latitudes.

Feng Ding↗