A CFD Model to Evaluate Near-Surface Oil Spill from a Broken Loading Pipe in Shallow Coastal Waters

Oil spills continue to generate various issues and concerns regarding their effect and behavior in the marine environment,owing to the related potential for detrimental environmental,economic and social implications.It is essential to have a solid understanding of the ways in which oil interacts with the water and the coastal ecosystems that are located nearby.This study proposes a simplified model for predicting the plume-like transport behavior of heavy Bunker C fuel oil discharging downward from an acutely-angled broken pipeline located on the water surface.The results show that the spill overall profile is articulated in three major flow areas.The first,is the source field,i.e.,a region near the origin of the initial jet,followed by the intermediate or transport field,namely,the region where the jet oil flow transitions into an underwater oil plume flow and starts to move horizontally,and finally,the far-field,where the oil re-surface and spreads onto the shore at a significant distance from the spill site.The behavior of the oil in the intermediate field is investigated using a simplified injection-type oil spill model capable of mimicking the undersea trapping and lateral migration of an oil plume originating from a negatively buoyant jet spill.A rectangular domain with proper boundary conditions is used to implement the model.The Projection approach is used to discretize a modified version of the Navier-Stokes equations in two dimensions.A benchmark fluid flow issue is used to verify the model and the results indicate a reasonable relationship between specific gravity and depth as well as agreement with the aerial data and a vertical temperature profile plot.

Residual symmetry, CRE integrability and interaction solutions of two higher-dimensional shallow water wave equations

Two(3+1)-dimensional shallow water wave equations are studied by using residual symmetry and the consistent Riccati expansion(CRE) method. Through localization of residual symmetries, symmetry reduction solutions of the two equations are obtained. The CRE method is applied to the two equations to obtain new B?cklund transformations from which a type of interesting interaction solution between solitons and periodic waves is generated.

High Order Finite Difference WENO Methods for Shallow Water Equations on Curvilinear Meshes

A high order finite difference numerical scheme is developed for the shallow water equations on curvilinear meshes based on an alternative flux formulation of the weighted essentially non-oscillatory(WENO)scheme.The exact C-property is investigated,and comparison with the standard finite difference WENO scheme is made.Theoretical derivation and numerical results show that the proposed finite difference WENO scheme can maintain the exact C-property on both stationarily and dynamically generalized coordinate systems.The Harten-Lax-van Leer type flux is developed on general curvilinear meshes in two dimensions and verified on a number of benchmark problems,indicating smaller errors compared with the Lax-Friedrichs solver.In addition,we propose a positivity-preserving limiter on stationary meshes such that the scheme can preserve the non-negativity of the water height without loss of mass conservation.

Two-dimensional shallow water equations with porosity and their numerical scheme on unstructured grids

In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosity source term appears in the momentum equation. The numerical model of the shallow water equations with porosity is presented with the finite volume method on unstructured grids and the modified Roe-type approximate Riemann solver. The source terms of the bed slope and porosity are both decomposed in the characteristic direction so that the numerical scheme can exactly satisfy the conservative property. The present model was tested with a dam break with discontinuous porosity and a flash flood in the Toce River Valley. The results show that the model can simulate the influence of obstructions, and the numerical scheme can maintain the flux balance at the interface with high efficiency and resolution.

An Explicit High Resolution Scheme for Nonlinear Shallow Water Equations

The present study develops a numerical model of the two-dimensional fully nonlinear shallow water equations （NSWE） for the wave run-up on a beach. The finite volume method （FVM） is used to solve the equations, and a second-order explicit scheme is developed to improve the computation efficiency. The numerical fluxes are obtained by the two dimensional Roe＇ s flux function to overcome the errors caused by the use of one dimensional fluxes in dimension splitting methods. The high-resolution Godunov-type TVD upwind scheme is employed and a second-order accuracy is achieved based on monotonic upstream schemes for conservation laws （MUSCL） variable extrapolation; a nonlinear limiter is applied to prevent unwanted spurious oscillation. A simple but efficient technique is adopted to deal with the moving shoreline boundary. The verification of the solution technique is carried out by comparing the model output with documented results and it shows that the solution technique is robust.

Exact Traveling Wave Solutions for the System of Shallow Water Wave Equations and Modified Liouville Equation Using Extended Jacobian Elliptic Function Expansion Method

In this work, an extended Jacobian elliptic function expansion method is proposed for constructing the exact solutions of nonlinear evolution equations. The validity and reliability of the method are tested by its applications to the system of shallow water wave equations and modified Liouville equation which play an important role in mathematical physics.

Moving-Water Equilibria Preserving HLL-Type Schemes for the Shallow Water Equations

We construct new HLL-type moving-water equilibria preserving upwind schemes for the one-dimensional Saint-Venant system of shallow water equations with nonflat bottom topography.The designed first-and secondorder schemes are tested on a number of numerical examples,in which we verify the well-balanced property as well as the ability of the proposed schemes to accurately capture small perturbations of moving-water steady states.

Solution of shallow-water equations using least-squares finite-element method

A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.

Least-squares finite-element method for shallow-water equations with source terms

Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.

High-resolution central difference scheme for the shallow water equations

A two-dimensional nonoscillatory central difference scheme was extended to the shallow water equations. A high-resolution numerical method for solving the shallow water equations was presented. In order to prevent oscillation, the nonlinear limiter is employed to approximate the discrete slopes. The main advantage of the presented method is simplicity comparable with the upwind schemes. This method does not require Riemann solvers or some form of flux difference splitting methods. Furthermore, the discrete derivatives of flux can be approximated by the component-wise approach and thus the computation of Jacobian can be avoided. The method retains high resolution and high accuracy similar to the upwind results. It is applied to simulating several tests, including circular dam-break problem, shock focusing problem and partial dam-break problem. The results are in good agreement with the numerical results obtained by other methods. The simulated results also demonstrate that the presented method is stable and efficient.

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION (Ⅱ)——THE DISCRETE-TIME CASE ALONG CHARACTERISTICS

The mixed finite element(MFE) methods for a shallow water equation system consisting of water dynamics equations,silt transport equation,and the equation of bottom topography change were derived.A fully discrete MFE scheme for the discrete_time along characteristics is presented and error estimates are established.The existence and convergence of MFE solution of the discrete current velocity,elevation of the bottom topography,thickness of fluid column,and mass rate of sediment is demonstrated.

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTA-TION (Ⅰ)-THE CONTINUOUS-TIME CASE

An initial-boundary value problem for shallow equation system consisting of water dynamics equations,silt transport equation, the equation of bottom topography change,and of some boundary and initial conditions is studied, the existence of its generalized solution and semidiscrete mixed finite element(MFE) solution was discussed, and the error estimates of the semidiscrete MFE solution was derived.The error estimates are optimal.

Analytical Solution of Two Extended Model Equations for Shallow Water Waves By Adomian’S Decomposition Method

In this paper, we consider two extended model equations for shallow water waves. We use Adomian’s decomposition method (ADM) to solve them. It is proved that this method is a very good tool for shallow water wave equations and the obtained solutions are shown graphically.

THE CAUCHY PROBLEM FOR THE TWO LAYER VISCOUS SHALLOW WATER EQUATIONS

In this paper,the Cauchy problem for the two layer viscous shallow water equations is investigated with third-order surface-tension terms and a low regularity assumption on the initial data.The global existence and uniqueness of the strong solution in a hybrid Besov space are proved by using the Littlewood-Paley decomposition and Friedrichs'regularization method.

Explicit Kinetic Flux Vector Splitting Scheme for the 2-D Shallow Water Wave Equations

Originally, the kinetic flux vector splitting (KFVS) scheme was developed as a numerical method to solve gas dynamic problems. The main idea in the approach is to construct the flux based on the microscopical description of the gas. In this paper, based on the analogy between the shallow water wave equations and the gas dynamic equations, we develop an explicit KFVS method for simulating the shallow water wave equations. A 1D steady flow and a 2D unsteady flow are presented to show the robust and accuracy of the KFVS scheme.

Analytical Solution of Two Extended Model Equations for Shallow Water Waves by He’s Variational Iteration Method

In this paper, we consider two extended model equations for shallow water waves. We use He’s variational iteration method (VIM) to solve them. It is proved that this method is a very good tool for shallow water wave equations and the obtained solutions are shown graphically.

A second-order accurate fluid-in-cell (FLIC) method for the 2D shallow water equations with topography

The fluid-in-cell （FLIC） approach of Gentry et al. （1966） is extended to second-order accuracy in space and applied to solve the 2D shallow water equations with topography. The FLIC method can be interpreted in a finite volume sense, it therefore conserves both water mass and momentum. Like the original FLIC method the second-order FLIC method presented here is able to handle wetting-drying fronts without any special treatment. Moreover, the resulting method is shock capturing and well-balanced, satisfying both the C- and extended C-properties exactly. Published by Elsevier Ltd on behalf of The Chinese Society of Theoretical and Applied Mechanics.

A Similarity Technique for Solving Two-Layer Shallow-Water Equations

This paper is devoted to the analysis of the two-layer shallow-water equations representing gravity currents. A similarity technique which is the characteristic function method is applied for this study. The application of the characteristic function method makes it possible to obtain the similarity forms depending on a group of infinitesimal transformations. Thus, the number of independent variables is reduced by one and the governing partial differential equations with the auxiliary conditions reduce to a system of ordinary differential equations with the appropriate auxiliary conditions. Numeric solutions are presented and discussed.

Wave Numerical Model for Shallow Water

The history of forecasting wind waves by wave energy conservation equation Is briefly described. Several currently used wave numerical models for shallow water based on different wave theories are discussed. Wave energy conservation models for the simulation of shallow water waves are introduced, with emphasis placed on the SWAN model, which takes use of the most advanced wave research achievements and has been applied to several theoretical and field conditions. The characteristics and applicability of the model, the finite difference numerical scheme of the action balance equation and its source terms computing methods are described in detail. The model has been verified with the propagation refraction numerical experiments for waves propagating in following and opposing currents; finally, the model is applied to the Haian Gulf area to simulate the wave height and wave period field there, and the results are compared with observed data.