Simulation of Free Surface Flow over Spillway

A numerical approach is proposed to simulate and study the effect of geometry on the free surface flow over a tunnel spillway. A three-step solution procedure is proposed to speed up the solution. The first step is to obtain an approximate free surface profile and mean velocity distribution, assuming 1D steady flow. Next, the 3D turbulent flow field is computed while the water surface profile is kept fixed. Finally, the water surface is set free to move and generate waves. The governing equations for weakly compressible flow (compressible hydrodynamic flow) are solved with an explicit finite volume method. A boundary fitted grid system is used to accurately resolve the flow near the free surface with steep waves. A mixed Lagrangian-Eulerian approach is proposed to calculate the new free surface position. The numerical results of a time-averaged free surface profile as well as pressure and velocity distribution have been compared with some experimental data.     

Spectral Collocation Model for Solitary Wave Attenuation and Mass Transport over Viscous Mud

In this paper, wave attenuation and mass transport of a water-mud system due to a solitary wave on the free surface is modeled by using the Chebyshev-Chebyshev collocation spectral method for spatial discretization and a fourth-order multistage scheme for time integration. The governing equations are formulated in Lagrangian coordinates and perturbation equations for shallow water waves are derived. An iteration-by-subdomain technique is introduced to tackle the interface in the two-layer system. The numerical model is tested against available analytical solutions and good agreement has been found. Numerical simulations of the water-mud system with different layer thicknesses suggest that the accuracy of the existing boundary layer theory for fluid-mud interaction is limited when the mud layer is thin because the assumption of irrotational core may not be valid. Although the paper is focused on solitary waves and Newtonian fluid-mud, the methodology can be extended to oscillatory, nonlinear water waves over a non-Newtonian mud bottom.     

Saltatory propagation of Ca2+ waves by Ca2+ sparks

Punctate releases of Ca2+, called Ca2+ sparks, originate at the regular array of t-tubules in cardiac myocytes and skeletal muscle. During Ca2+ overload sparks serve as sites for the initiation and propagation of Ca2+ waves in myocytes. Computer simulations of spark-mediated waves are performed with model release sites that reproduce the adaptive Ca2+ release observed for the ryanodine receptor. The speed of these waves is proportional to the diffusion constant of Ca2+, D, rather than D, as is true for reaction-diffusion equations in a continuous excitable medium. A simplified "fire-diffuse-fire" model that mimics the properties of Ca2+-induced Ca2+ release (CICR) from isolated sites is used to explain this saltatory mode of wave propagation. Saltatory and continuous wave propagation can be differentiated by the temperature and Ca2+ buffer dependence of wave speed.     

Flood Simulation Using a Well-Balanced Shallow Flow Model

This work extends and improves a one-dimensional shallow flow model to two-dimensional (2D) for real-world flood simulations. The model solves a prebalanced formulation of the fully 2D shallow water equations, including friction source terms using a finite volume Godunov-type numerical scheme. A reconstruction method ensuring nonnegative depth is used along with a Harten, Lax, and van Leer approximate Riemann solver with the contact wave restored for calculation of interface fluxes. A local bed modification method is proposed to maintain the well-balanced property of the algorithm for simulations involving wetting and drying. Second-order accurate scheme is achieved by using the slope limited linear reconstruction together with a Runge-Kutta time integration method. The model is applicable to calculate different types of flood wave ranging from slow-varying inundations to extreme and violent floods, propagating over complex domains including natural terrains and dense urban areas. After validating against an analytical case of flow sloshing in a domain with a parabolic bed profile, the model is applied to simulate an inundation event in a 36?km2 floodplain in Thamesmead near London. The numerical predictions are compared with analytical solutions and alternative numerical results.     

Agonist-induced calcium waves in oscillatory cells: a biological example of Burgers' equation

Oscillations in cytosolic Ca2+ concentrations in living cells are often a manifestation of propagating waves of Ca2+. Numerical simulations with a realistic model of inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ wave trains lead to wave speeds that increase linearly at long times when (a) IP3 levels are in the range for Ca2+ oscillations, (b) a gradient of phase is established by either an initial ramp or pulse of IP3, and (c) IP3 concentrations asymptotically become uniform. We explore this phenomenon with analytical and numerical methods using a simple two-variable reduction of the De Young-Keizer model of the IP3 receptor that includes the influence of Ca2+ buffers. For concentrations of IP3 in the oscillatory regime, numerical solution of the resulting reaction diffusion equations produces nonlinear wave trains that shows the same asymptotic growth of wave speed. Due to buffering, diffusion of Ca2+ is quite slow and, as previously noted, these waves occur without appreciable bulk movement of Ca2+. Thus, following Neu and Murray, we explore the behavior of these waves using an asymptotic expansion based on the small size of the buffered diffusion constant for Ca2+. We find that the gradient in phase of the wave obeys Burgers' equation asymptotically in time. This result is used to explain the linear increase of the wave speed observed in the simulations.     

Interaction of Nonlinear Progressive Waves with Two Serially Arranged Submerged Obstacles

The purpose of the present study is to develop a numerical model for the investigation of water waves propagating over a pair of impermeable submerged obstacles. The mathematic model is formulated by coupling solutions of the Navier–Stokes equations and transport equations for the surface elevation using the volume of fluid method. Based on a staggered computational mesh, an explicit numerical algorithm is employed with a predictor–corrector procedure of pressure and velocity field. The proposed model provides good agreement with other experimental results and validates its good performance. Regarding the spatial harmonic evolutions of various cases, it is noted that the present fluctuating mode of harmonic amplitudes exists upstream and at the gap between obstacles. The results show that the nonlinearity of propagating waves becomes stronger than the initial wave in such areas, and reveals much steeper wave profiles compared to the initial ones. The fluctuating harmonic amplitudes vary with the gap width and form two hydrodynamic cycles. The vortices play an important role in the wave reflection as they form a water column wall to reflect the incoming waves. The reflection ratio depends on the extent of vortex development near the upstream obstacle. The maximum wave reflection occurs in cases with dimensionless gap width S/L equal to 3/8 and 7/8 in this study.     

Lagrangian Modeling of Weakly Nonlinear Nonhydrostatic Shallow Water Waves in Open Channels

A Lagrangian, nonhydrostatic, Boussinesq model for weakly nonlinear and weakly dispersive flow is presented. The model is an extension of the hydrostatic model—dynamic river model. The model uses a second-order, staggered grid, predictor-corrector scheme with a fractional step method for the computation of the nonhydrostatic pressure. Numerical results for solitary waves and undular bores are compared with Korteweg-de Vries analytical solutions and published numerical, laboratory, and theoretical results. The model reproduced well known features of solitary waves, such as wave speed, wave height, balance between nonlinear steepening and wave dispersion, nonlinear interactions, and phase shifting when waves interact. It is shown that the Lagrangian moving grid is dynamically adaptive in that it ensures a compression of the grid size under the wave to provide higher resolution in this region. Also the model successfully reproduced a train of undular waves (short waves) from a long wave such that the predicted amplitude of the leading wave in the train agreed well with published numerical and experimental results. For prismatic channels, the method has no numerical diffusion and it is demonstrated that a simple second-order scheme suffices to provide an efficient and economical solution for predicting nonhydrostatic shallow water flows.     

Using Finite-Element Method to Simulate Wave Transformations in Surf Zone

In this study, a finite element method proposed by Hsu et al. in 2003 is extended to develop a numerical model for the simulation of wave transformation in the surf zone. The governing equation is the elliptic mild-slope equation including the energy dissipation of wave breaking. At the open boundaries with varying depth, the reflected waves caused by shoaling are adopted to the radiation boundary conditions. The rationality of the present numerical model is examined through the cases of offshore parallel breakwater problems. The results of calculation are in good agreement with experimental results.     

Kinematic and Diffusion Waves: Analytical and Numerical Solutions to Overland and Channel Flow

Flood wave propagation is the unifying concept in representing open channel and overland flow. Therefore, understanding flood wave routing theory and solving the governing equations accurately is an important issue in hydrology and hydraulics. In an attempt to contribute to the understanding of this subject, in this study: (1) an analytical solution is derived for diffusion waves with constant wave celerity and hydraulic diffusivity applied to overland flow problems; and (2) an algorithm is developed using the MacCormack explicit finite difference method to solve the kinematic and diffusion wave governing equations for both overland and open channel flow. The MacCormack method is particularly well suited to approximate nonlinear differential equations. The analytical solutions provide the practicing engineer with computational speed in obtaining results for overland flow problems, and a means to check the validity of the numerical models. On the other hand, for larger scale catchment-stream problems, the verified numerical methods provide efficient and accurate algorithms to obtain solutions. Both the analytical approaches and the MacCormack algorithm are used to solve the same synthetic examples. Comparison of results shows that the numerical and analytical solutions are in close agreement. Furthermore, the MacCormack algorithm is applied to a real catchment: a segment of the Duke University West Campus storm water drainage system. In order to check the accuracy of the results obtained by the MacCormack method, the results are compared to predictions of the Environmental Protection Agency storm water management model (SWMM) as calibrated with measured rainfall and surface runoff flow data. The results obtained from SWMM are in good agreement with the results obtained from applying the MacCormack algorithm.     

RANS∕Laplace Calculations of Nonlinear Waves Induced by Surface-Piercing Bodies

An interactive zonal numerical method has been developed for the prediction of free surface flows around surface-piercing bodies, including both viscous and nonlinear wave effects. In this study, a Laplace solver for potential flow body-wave problems is used in conjunction with a Reynolds-averaged Navier-Stokes (RANS) method for accurate resolution of viscous, nonlinear free surface flows around a vertical strut and a series 60 ship hull. The Laplace equation for potential flow is solved in the far field to provide the nonlinear waves generated by the body. The RANS method is used in the near field to resolve the turbulent boundary layers, wakes, and nonlinear waves around the body. Both the kinematic and dynamic boundary conditions are satisfied on the exact free surface to ensure accurate resolution of the divergent and transverse waves. The viscous-inviscid interaction between the potential flow and viscous flow regions is captured through a direct matching of the velocity and pressure fields in an overlapping RANS and potential flow computational region. The numerical results demonstrate the capability of an interactive RANS∕Laplace coupling method for accurate and efficient resolution of the body boundary layer, the viscous wake, and the nonlinear waves induced by surface-piercing bodies.     

Comparison between Computed and Experimentally Generated Impulse Waves

Large water waves caused by massive slide impacts are a potential hazard along waterways, coastal areas and Alpine regions. Experimental research has been conducted at the Swiss Laboratory of Hydraulics to assess the risk from landslide-generated impulse waves. Analogously, the Centro Elettrotecnico Sperimentale Italiano performed numerical simulations of initial landslide and consequent impulse wave propagation using two mathematical models based on the conservative shallow-water equations. This paper presents the experimental test results and numerical predictions of impulse waves in a flume for a range of stillwater depths, landslide volumes, and impact velocities at laboratory scale. The comparison between the measured and predicted wave free surface profiles generally produced corresponding wave heights, although the initial wave peak is too steep and arrives too early. Excluding spurious random effects, the relative differences between measured and numerically computed maximum wave heights ranged within ±20%, which can be considered satisfactory from the engineering point of view.     

Stochastic Analysis of a Single-Degree-of-Freedom Nonlinear Experimental Moored System Using an Independent-Flow-Field Model

Stochastic characteristics of the surge response of a nonlinear single-degree-of-freedom moored structure subjected to random wave excitations are examined in this paper. Sources of nonlinearity of the system include a complex geometric configuration and wave-induced quadratic drag. A Morison-type model with an independent-flow-field formulation and a three-term-polynomial approximation of the nonlinear restoring force is employed for its proven excellent prediction capability for the experimental results investigated. Wave excitations considered in this study include nearly periodic waves, which take into account the presence of tank noise, noisy periodic waves that have predominant periodic components with designed additive random perturbations, and narrow-band random waves. A unified wave excitation model is used to describe all the wave conditions. A modulating factor governing the degree of randomness in the wave excitations is introduced. The corresponding Fokker–Planck formulation is applied and numerically solved for the response probability density functions (PDFs). Experimental results and simulations are compared in detail via the PDFs in phase space. The PDFs portray coexisting multiple response attractors and indicate their relative strengths, and experimental response behaviors, including transitions and interactions, are accordingly interpreted from the ensemble perspective. Using time-averaged probability density functions as an invariant measure, probability distributions of large excursions in experimental and simulated responses to various random wave excitations are demonstrated and compared. Asymptotic long-term behaviors of the experimental responses are then inferred.     

Applicability of Sediment Transport Capacity Formulas to Dam-Break Flows over Movable Beds

In this paper, we investigate the extent to which well-known sediment transport capacity formulas can be used in one-dimensional (1D) numerical modeling of dam-break waves over movable beds. The 1D model considered here is a one-layer model based on the shallow-water equations, a bed update (Exner) equation, a space-lag equation for the nonequilibrium sediment transport and an empirical formula calculating the sediment transport capacity of the flow. The model incorporates a variety of sediment transport capacity formulas proposed by Meyer-Peter and Müller, Bagnold, Engelund and Hansen, Ackers and White, Smart and Jaeggi, van Rijn, Rickenmann, Cheng, Abrahams and Camenen, and Larson. We examine the performance of each formula by simulating four idealized laboratory cases on dam-break waves over sandy beds. Comparisons between numerical results and measurements show that for each case better predictions are obtained using a particular formula, but overall, formulas proposed by Meyer-Peter and Müller (with the factor 8 being replaced by 12), Smart and J?ggi, Cheng, Abrahams and Camenen, and Larson rank as the best predictors for the entire range of conditions studied here. Moreover, results show that in the cases where a bed step exists, implementing a mass failure mechanism in the numerical modeling plays an important role in reproducing the bed and water profiles.     

Hydrodynamic Investigation of Coastal Bridge Collapse during Hurricane Katrina

A number of U.S. coastal bridges have been destroyed by hurricanes, including three highway bridges in Mississippi and Louisiana during Hurricane Katrina (2005). This paper addresses three fundamental questions on the coastal bridge failures: (1) what were the hydrodynamic conditions near the failed bridge during the hurricane; (2) what was the cause of the bridge collapse; and (3) what was the magnitude of the hydrodynamic loading on the bridge under the extreme hurricane conditions. Guided by field observations of winds, waves, and water levels, two numerical models for storm surges and water waves are coupled to hindcast the hydrodynamic conditions. Fairly good agreement between the modeled and measured high watermarks and offshore wave heights is found, allowing an estimate of the surge and wave conditions near the bridges in nested domains with higher resolutions. The output of the coupled wave-surge models is utilized to determine the static buoyant force and wave forces on the bridge superstructure based on empirical equations derived from small-scale hydraulic tests for elevated decks used in the coastal and offshore industry. It is inferred that the bridge failure was caused by the wind waves accompanied by the storm surge, which raised the water level to an elevation where surface waves generated by strong winds over a relatively short fetch were able to strike the bridge superstructure. The storm waves produced both an uplift force and a horizontal force on the bridge decks. The magnitude of wave uplift force from individual waves exceeded the weight of the simple span bridge decks and the horizontal force overcame the resistance provided by the connections of the bridge decks to the pilings. The methodology for determining the hydrodynamic forcing on bridge decks can be used to produce a preliminary assessment of the vulnerability of existing coastal bridges in hurricane-prone areas.     

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer,and is illustrated by establishing second order (M=l) recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.     

Hydrodynamic Performance of Wave-Driven Artificial Upwelling Device

Artificial upwelling and mixing is a new technology for enhancing open ocean mariculture using nutrient-rich deep ocean water. A wave-driven artificial upwelling device was developed based on mathematical modeling and hydraulic experiments. Results of the numerical modeling are in good agreement with hydraulic experiments. The mathematical model was applied to simulate the operation of a prototype device in typical Hawaiian waves. The prototype device consists of a buoy with a water chamber 4.0 m in diameter, a flow-controlling valve, and a long tail pipe 300 m in length and 1.2 m in diameter. This device can produce an upwelling flow of 0.62 m3∕s in regular Hawaiian waves with a period of 12 s and a wave height of 1.9 m. Two analytical solutions to the simplified governing equations were also derived that provide a basis for the derivation of predictive equations, useful for preliminary engineering design and analysis.     

Generation of Three-Dimensional Fully Nonlinear Water Waves by a Submerged Moving Object

This paper describes a numerical investigation on the generation of three-dimensional (3D) fully nonlinear water waves by a submerged object moving at speeds varied from subcritical to supercritical conditions in an unbounded fluid domain. Considering a semispheroid as the moving object, simulations of the time evolutions of 3D free-surface elevation and flow field are performed. The present 3D model results are found to agree reasonably well with other published vertical two-dimensional (2D) and quasi-3D numerical solutions using Boussinesq-type models. Different from the 2D cases with near critical moving speeds, the 3D long-term wave pattern suggests that in addition to the circularly expanded upstream advancing solitonlike waves, a sequence of divergent and transverse waves are also developed behind the moving object. The velocity distributions and associated fluid-particle trajectories at the free-surface and middle layers are presented to show the 3D feature of the motion. The results under various vertical positions (referred as gap) of a moving object are also compared. It is found the gap has shown a substantial influence on the generated waves, especially in the wake region, when an object moves at a near critical or subcritical speed. However, the results under the case with a high supercritical moving speed indicate the gap has a negligible effect on the generated upstream and downstream waves.     

Analysis of Debris Wave Development with One-Dimensional Shallow-Water Equations

The objective of this contribution is to analyze the formation of debris waves in natural channels. Numerical simulations are carried out with a 1D code, based on shallow-water equations and on the weighted averaged flux method. The numerical code represents the incised channel geometry with a power-law relation between local width and flow depth and accounts for all source terms in the momentum equation. The debris mixture is treated as a homogeneous fluid over a fixed bottom, whose rheological behavior alternatively follows Herschel-Bulkley, Bingham, or generalized viscoplastic models. The code is first validated by applying it to dam-break tests on mudflows down a laboratory chute and verifying its efficiency in the simulation of rapid transients. Then, following the analytical method developed by Trowbridge, the stability of a uniform flow for a generalized viscoplastic fluid is examined, showing that debris flows become unstable for Froude numbers well below 1. Applications of the code to real debris flow events in the Cortina d’Ampezzo area (Dolomites) are presented and compared with available measured hydrographs. A statistical analysis of debris waves shows that a good representation of wave statistics can be obtained with a proper calibration of rheological parameters. Finally, it is shown that a minimum duration of debris event and channel length are required for waves showing up, and an explanation, confirmed both by field data and numerical simulations, is provided.     

Viscous Effect on the Roll Motion of a Rectangular Structure

The viscous effect on the roll motion of a rectangular structure was investigated in a two-dimensional wave tank. The structure was used to simulate a simplified barge in the beam sea condition. The structure with a draft one-half of its height was hinged at the center of gravity and free to roll (1 degree of freedom) by waves. The dynamic characteristics of the structure were first identified, including its roll natural period. The dynamic response of the barge-like structure under wave actions was then tested with regular waves with a range of wave periods that are shorter, equal to, and longer than its roll natural period. Particle image velocimetry was used to obtain the velocity field in the vicinity of the structure. The coupled interactions between the incident waves and the structure were demonstrated by examining the vortical flow fields to elucidate the effect of viscous damping (also called the eddy making damping) to the roll motion of the structure over wave periods. For incoming waves with a wave period same as the roll natural period, the structure roll motion was, as expected, greatly reduced by the viscous-damping effect. At wave periods shorter than the roll natural period, the structure roll motion was slightly reduced by the viscous effect. However, at wave periods longer than the roll natural period, the viscous effect due to flow separation at structure corners indeed amplified the roll motion. This indicates that not only can the viscous effect damp out the roll motion, it can also amplify the roll motion.     

Computational Fluid Dynamics Modeling for the Design of Large Primary Settling Tanks

Settling tanks are used to remove solids at wastewater treatment plants. Many numerical models have been proposed to simulate the settling process and to improve tank efficiency. In this research, a three-dimensional (3D) numerical model is developed to simulate large primary settling tanks. In the proposed model, the non-Newtonian properties of the sludge flow in the settling tank are described by a Bingham plastic rheological model. To eliminate the singularity inherited in the rheological model, a modified constitutive relation is used in both the yielded and unyielded regions. Hindered settling of particles in the settling tank is also modeled. Tracer study, where a massless scalar is injected and transported, is done to investigate the tank’s residence time. This numerical model is used to improve the design of the primary settling tanks, which will be built in Chicago. The Metropolitan Water Reclamation District of Greater Chicago (MWRDGC) is in the process of building new preliminary treatment facilities at their Calumet Water Reclamation Plant (CWRP), including twelve 155-ft-diameter primary settling tanks (PSTs) designed to treat flows up to (480?million?gal./day (MGD). The computational fluid dynamics (CFD) model simulated solids removal efficiencies based on a particle size distribution similar to the one observed in the CWRP influent. The results were used to establish the design basis for tank side-water depth, inlet feedwell dimensions, etc., resulting in improved performance and substantial reduction in construction costs.