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 55 records · Page 3

Monsoon Depression Amplification by Horizontal Shear and Humidity Gradients: A Shallow Water Perspective

Transient, synoptic-scale vortices produce a large fraction of total rainfall in most monsoon regions and are often associated with extreme precipitation. However, the mechanism of their amplification remains a topic of active research. For monsoon depressions, which are the most prominent synoptic-scale vortex in the Asian–Australian monsoon, recent work has suggested that meridional gradients in zonal wind in the vortex environment may produce growth through barotropic instability, while meridional gradients in environmental humidity have also been proposed to cause amplification through coupling with precipitating convection. Here, a two-dimensional shallow water model on a sphere with parameterized precipitation is used to examine the relative role played by these two environmental gradients. By systematically varying the meridional moisture gradient and meridional wind shear for both weak, quasi-linear waves and finite-amplitude isolated vortices, we show that rotational winds in the initial vortex are amplified most strongly by meridional shear of the environmental zonal wind, while vortex precipitation rates are most sensitive to environmental moisture gradients. The growth rate in the presence of both gradients is less than the sum of growth rates in the presence of isolated gradients, as the phase relation between moisture and vorticity anomalies becomes distorted with increasing shear. In conclusion, these results suggest that background meridional gradients in both zonal wind and environmental humidity can contribute to the amplification of vortices to monsoon depression strength, but with some degree of decoupling of the dry rotational flow and the moist convection.

54 ENVIRONMENTAL SCIENCES↗

Design and Analysis of a Floating-Wind Shallow-Water Mooring System Featuring Polymer Springs

This paper presents a mooring system design featuring polymer springs for the VolturnUS-S 15-MW reference floating wind turbine in site conditions for the New York Bight at a 50-m water depth. Polymer springs have a nonlinear stress-strain curve that allows a stiffer response at low loads and a more flexible response at higher loads, potentially reducing peak mooring line tensions. The mooring dynamics model MoorDyn has been extended to model springs with nonlinear tension-strain curves from a user-inputted look-up table. This MoorDyn modeling advancement is verified against OrcaFlex simulation results. Using MoorDyn's updated capabilities, a spring-equipped catenary mooring system is designed for the 15-MW floating system, along with a baseline catenary mooring system that does not use polymer springs. The floating wind turbine simulator OpenFAST is used to simulate the mooring systems in design-driving load cases to show the effect of polymer springs on key dynamic behaviours. The results show that the spring-equipped design reduces peak tensions by up to 60%, whereas the turbine offsets stay within a maximum of 7.2 m, which is still a reasonable offset limit for cable considerations. The reduction in peak tensions results in a significant decrease in damage equivalent loads - on the order of 50% for upwind lines in fully loaded conditions. These results show that polymer springs can effectively reduce peak tensions and fatigue loads in mooring systems at shallow water depths.

floating offshore wind↗

Design and Analysis of a Floating-Wind Shallow-Water Mooring System Featuring Polymer Springs: Preprint

In this paper, a mooring system featuring polymer springs is designed for the VolturnUS-S 15 MW reference floating wind turbine in site conditions for the New York Bight at a 50-m water depth. Polymer springs have a nonlinear stress-strain curve that allows a stiffer response at low loads and a more flexible response at higher loads, potentially reducing peak mooring line tensions. The mooring dynamics model MoorDyn has been extended to model springs with nonlinear tension-strain curves from a user-inputted look-up table. This MoorDyn modelling advancement is verified against OrcaFlex simulation results. Using MoorDyn's updated capabilities, a spring-equipped catenary mooring system is designed for the 15 MW floating system, along with a baseline catenary mooring system that does not use polymer springs. The floating wind turbine simulator OpenFAST is used to simulate the mooring systems in design-driving load cases to show the effect of polymer springs on key dynamic behaviours. The results show that the spring-equipped design reduces peak tensions by up to 60%, while the turbine offsets stay within a maximum of 7.2 m, which is still a reasonable offset limit for cable considerations. The reduction in peak tensions results in a significant decrease in damage equivalent loads, on the order of 50% for upwind lines in fully loaded conditions. These results show that polymer springs can effectively reduce peak tensions and fatigue loads in mooring systems at shallow water depths.

floating offshore wind↗

Uncertainties in Simulating Flooding During Hurricane Harvey Using 2D Shallow Water Equations

Abstract Flooding is one of the most impactful weather‐related natural hazards. Numerical models that solve the two dimensional (2D) shallow water equations (SWE) represent the first‐principles approach to simulate all types of spatial flooding, such as pluvial, fluvial, and coastal flooding, and their compound dynamics. High spatial resolution (e.g., () m) is needed in 2D SWE simulations to capture flood dynamics accurately, resulting in formidable computational challenges. Thus, relatively coarser spatial resolutions are used for large‐scale simulations of flooding, which introduce uncertainties in the results. It is unclear how the uncertainty associated with the model resolution compares to the uncertainties in precipitation data sets and assumptions regarding boundary conditions when channelized flows interact with other water bodies. In this study, we compare these three sources of uncertainties in 2D SWE simulations for the 2017 Houston flooding event. Our results show that precipitation uncertainty and mesh resolution have more significant impacts on the simulated streamflow and inundation dynamics than the choice of the downstream boundary condition at the watershed outlet. We point out the viability to confine the uncertainty of coarsening mesh resolution by using the variable resolution mesh (VRM) which refines critical topographic features with far fewer grid cells. Specifically, in simulations with VRM, the simulated inundation depths over the refined region are comparable to that use the finest uniform mesh. This study contributes to understanding the challenges and pathways for applying 2D SWE models to improve the realism of flood simulations over large scales.

54 ENVIRONMENTAL SCIENCES↗

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

Balanced Central Schemes for the Shallow Water Equations on Unstructured Grids

We present a two-dimensional, well-balanced, central-upwind scheme for approximating solutions of the shallow water equations in the presence of a stationary bottom topography on triangular meshes. Our starting point is the recent central scheme of Kurganov and Petrova (KP) for approximating solutions of conservation laws on triangular meshes. In order to extend this scheme from systems of conservation laws to systems of balance laws one has to find an appropriate discretization of the source terms. We first show that for general triangulations there is no discretization of the source terms that corresponds to a well-balanced form of the KP scheme. We then derive a new variant of a central scheme that can be balanced on triangular meshes. We note in passing that it is straightforward to extend the KP scheme to general unstructured conformal meshes. This extension allows us to recover our previous well-balanced scheme on Cartesian grids. We conclude with several simulations, verifying the second-order accuracy of our scheme as well as its well-balanced properties.

Bryson, Steve↗

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

To improve the simulation of nonlinear aspects of the flow over steep topography, a potential enstrophy and energy conserving scheme for the shallow water equations is derived. It is pointed out that a family of schemes can conserve total energy for general flow and potential enstrophy for flow with no mass flux divergence. The newly derived scheme is a unique member of this family, that conserves both potential enstrophy and energy for general flow. Comparison by means of numerical experiment with a scheme that conserves (potential) enstrophy for purely horizontal nondivergent flow demonstrated the considerable superiority of the newly derived potential enstrophy and energy conserving scheme, not only in suppressing a spurious energy cascade but also in determining the overall flow regime. The potential enstrophy and energy conserving scheme for a spherical grid is also presented.

Arakawa, A.↗

Variational data assimilation with a semi-Lagrangian semi-implicit global shallow-water equation model and its adjoint

An adjoint model is developed for variational data assimilation using the 2D semi-Lagrangian semi-implicit (SLSI) shallow-water equation global model of Bates et al. with special attention being paid to the linearization of the interpolation routines. It is demonstrated that with larger time steps the limit of the validity of the tangent linear model will be curtailed due to the interpolations, especially in regions where sharp gradients in the interpolated variables coupled with strong advective wind occur, a synoptic situation common in the high latitudes. This effect is particularly evident near the pole in the Northern Hemisphere during the winter season. Variational data assimilation experiments of 'identical twin' type with observations available only at the end of the assimilation period perform well with this adjoint model. It is confirmed that the computational efficiency of the semi-Lagrangian scheme is preserved during the minimization process, related to the variational data assimilation procedure.

Li, Y.↗

A Shallow-Water Model for Convective Self-Aggregation

Randomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow-water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. Here we show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping time scale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

54 ENVIRONMENTAL SCIENCES↗

Comparison of horizontal difference schemes for the shallow water equations on a sphere

The accuracy of horizontal difference schemes used in the hydrodynamics parts of General Circulation Models are compared by means of numerical experiments for the shallow water equations on a sphere. As expected, the phase lag of moving waves decreases as the order of accuracy of a scheme increases or as the grid resolution increases. Overall, Takano and Wurtele's partial fourth order energy and potential enstrophy conserving scheme on the C grid is most accurate. It is clearly superior to the other schemes for the Rossby-Haurwitz wave number 6 initial conditions for coarse grid resolution.

Russell, Gary L.↗

A multigrid solver for semi-implicit global shallow-water models

A multigrid solver is developed for the discretized two-dimensional elliptic equation on the sphere that arises from a semiimplicit time discretization of the global shallow-water equations. Different formulations of the semiimplicit scheme result in variable-coefficient Helmholtz-type equations for which no fast direct solvers are available. The efficiency of the multigrid solver is optimal, in the sense that the total operation count is proportional to the number of unknowns. Numerical experiments using initial data derived from actual 300-mb height and wind velocity fields indicate that the present model has very good accuracy and stability properties.

Barros, Saulo R. M.↗

Shallow water experiments using an improved Kreiss-Oliger 2-4 scheme

The Abarbanel-Gottleib (AG) modification of the Kreiss-Oliger (KO) scheme is applied to a global 1-level shallow water model and the results compared to results obtained using the standard KO scheme. The standard KO scheme consists of the second order leapfrog scheme for time derivatives and fourth order space differencing for spatial derivatives. The AG scheme modifies the KO scheme in a manner similar to that done by Wilson to the multidimensional leapfrog scheme, obtaining a 35-40% improvement of the maximum time step. Application of the AG scheme expands the domain of influence of the KO scheme by including off-axis grid points.

Takacs, L. L.↗

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↗

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

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

Assessing the fidelity of shallow-water carbonates as records of the Ni isotope composition of surface seawater

Nickel is a bioessential metal that is used in enzymes important to the C, N, and O cycles, and changes in its marine abundance and bioavailability may have affected the evolutionary trajectory of early life. Changes over time in the Ni isotope composition (δ 60 Ni) of surface seawater, which reflects biological demand for Ni, could allow for the reconstruction of the dynamics of Ni demand over Earth’s history, but this approach would require geologic records of surface seawater. Here, we investigate the fidelity of shallow-water carbonates as a record of the Ni isotope composition of surface seawater by determining how Ni is first partitioned into natural carbonates and then how post-depositional processes influence the Ni signal. Our samples come from the Great Bahama Bank, which is a well-studied, modern carbonate platform often used to study ancient platforms. We found that Ni is fractionated from seawater upon incorporation into carbonates capturing shallow (<18 cm), recent deposition (0.1 ‰–0.4 ‰ lighter than seawater). Variation among these [Ni] and δ 60 Ni values may be controlled by variation in mineral proportions. Meteoric diagenesis shifts δ 60 Ni to lower values, which we attribute to isotopically light meteoric fluids. In contrast, carbonates that experienced sediment-buffered marine diagenesis with respect to Ca isotopes and Sr/Ca ratios do not appear to differ in δ 60 Ni values from sediments generally representative of their initial deposition. The sensitivity of δ 60 Ni to diagenetic reset in these samples appears comparable to the sensitivities of Ca isotopes and Sr/Ca ratios, to first order. Thus, in general, carbonates that experienced sediment-buffered marine diagenesis with respect to these elements may hold the most promise as a record of the δ 60 Ni of coeval surface seawater. Additionally, we use our results to infer that the fraction of Ni removed from seawater into carbonates is less than 10 % of the total Ni output from the global oceans and incorporation of this Ni sink into global biogeochemical models will only have a minor impact on the modeled modern Ni budget.

Great Bahama Bank carbonates↗

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