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.

A THIRD GENERATION SHALLOW WATER WAVE NUMERICAL MODEL-YE-WAM

This paper pnesents a third gneration shallow Whter disode spedtal wave nbotal medeIYE-WAM based on the spedtal action balance equation. The mode accounts for all edevan effectsof currents on waves, incuding tmpotally and spatialy varying depth and current inded refraction,sttalning and fequency shift and also explidtly takeS into aanunt all source terms, speclally adePth-limited breaking dheipation. In addition, an energy forcing scheme is propond and applied to themode’s open boundaries to areUn for the propagution of sedIs into the study spstem The upwinddiffeIenng scheme and a standard hybrid diffdrencing scheme for the propagaion terrn and a simpleEuler method for the source teme are employed.

Sedimentary facies and depositional model of shallow water delta dominated by fluvial for Chang 8 oil-bearing group of Yanchang Formation in southwestern Ordos Basin,China

A systematic analysis of southwestern Ordos Basin's sedimentary characteristics,internal architectural element association styles and depositional model was illustrated through core statistics,well logging data and outcrop observations in Chang 8 oil-bearing group.This analysis indicates that shallow water delta sediments dominated by a fluvial system is the primary sedimentary system of the Chang 8 oil-bearing group of the Yanchang Formation in southwestern Ordos Basin.Four microfacies with fine grain sizes are identified: distributary channels,sheet sandstone,mouth bar and interdistributary fines.According to the sandbody's spatial distribution and internal architecture,two types of sandbody architectural element associations are identified: amalgamated distributary channels and thin-layer lobate sandstone.In this sedimentary system,net-like distributary channels at the delta with a narrow ribbon shape compose the skeleton of the sandbody that extends further into the delta front and shades into contiguous lobate distribution sheet sandstone in the distal delta front.The mouth bar is largely absent in this system.By analyzing the palaeogeomorphology,the palaeostructure background,sedimentary characteristics,sedimentary facies types and spatial distribution of sedimentary facies during the Chang 8 period,a distinctive depositional model of the Chang 8 shallow water fluvial-dominated delta was established,which primarily consists of straight multi-phase amalgamated distributary channels in the delta plain,net-like distributary channels frequently diverting and converging in the proximal delta front,sheet sandstones with dispersing contiguous lobate shapes in the distal delta front,and prodelta or shallow lake mudstones.

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.

Model Tests Research on A Float-Over Barge in Shallow Water Under the Undocking Conditions

In this study,the Jinzhou 9-3 CEPD float-over installation project was investigated.During the undocking condition,the water depth of the motion path of the working barge gradually changed from 10.31 m to 9.41 m.The undocking clearance of the HYSY 228 is smaller than 1 m;therefore,the barge shows highly nonlinear hydrodynamic characteristics,and it is difficult to be accurately simulated by numerical analysis.Thus,it is necessary to obtain the hydrodynamic characteristics and laws of the float-over barge at different water depths by using tank model test,to provide some reference and guidance for float-over operations in shallow water.

Performances of Different Turbulence Models for Simulating Shallow Water Sloshing in Rectangular Tank

Liquid sloshing is a common phenomenon in the transportation of liquid-cargo tanks.Liquid waves lead to fluctuating forces on the tank walls.If these fluctuations are not predicted or controlled,for example,by using baffles,they can lead to large forces and momentums.The volume of fluid(VOF)two-phase numerical model in Open FOAM open-source software has been widely used to model the liquid sloshing.However,a big challenge for modeling the sloshing phenomenon is selecting a suitable turbulence model.Therefore,in the present study,different turbulence models were studied to determine their sloshing phenomenon prediction accuracies.The predictions of these models were validated using experimental data.The turbulence models were ranked by their mean error in predicting the free surface behaviors.The renormalization group(RNG)k-ε and the standard k–ω models were found to be the best and worst turbulence models for modeling the sloshing phenomena,respectively;moreover,the SST k-ω model and v2-f k-ε results were very close to the RNG k-εmodel result.

Effects of gravel on the water absorption characteristics and hydraulic parameters of stony soil

The eastern foothills of the Helan Mountains in China are a typical mountainous region of soil and gravel,where gravel could affect the water movement process in the soil.This study focused on the effects of different gravel contents on the water absorption characteristics and hydraulic parameters of stony soil.The stony soil samples were collected from the eastern foothills of the Helan Mountains in April 2023 and used as the experimental materials to conduct a one-dimensional horizontal soil column absorption experiment.Six experimental groups with gravel contents of 0%,10%,20%,30%,40%,and 50%were established to determine the saturated hydraulic conductivity(K_(s)),saturated water content(θ_(s)),initial water content(θ_(i)),and retention water content(θ_(r)),and explore the changes in the wetting front depth and cumulative absorption volume during the absorption experiment.The Philip model was used to fit the soil absorption process and determine the soil water absorption rate.Then the length of the characteristic wetting front depth,shape coefficient,empirical parameter,inverse intake suction and soil water suction were derived from the van Genuchten model.Finally,the hydraulic parameters mentioned above were used to fit the soil water characteristic curves,unsaturated hydraulic conductivity(K_(θ))and specific water capacity(C(h)).The results showed that the wetting front depth and cumulative absorption volume of each treatment gradually decreased with increasing gravel content.Compared with control check treatment with gravel content of 0%,soil water absorption rates in the treatments with gravel contents of 10%,20%,30%,40%,and 50%decreased by 11.47%,17.97%,25.24%,29.83%,and 42.45%,respectively.As the gravel content increased,inverse intake suction gradually increased,and shape coefficient,K_(s),θ_(s),andθ_(r)gradually decreased.For the same soil water content,soil water suction and K_(θ)gradually decreased with increasing gravel content.At the same soil water suction,C(h)decreased with increasing gravel content,and the water use efficiency worsened.Overall,the water holding capacity,hydraulic conductivity,and water use efficiency of stony soil in the eastern foothills of the Helan Mountains decreased with increasing gravel content.This study could provide data support for improving soil water use efficiency in the eastern foothills of the Helan Mountains and other similar rocky mountainous areas.

Numerical three-dimensional modeling of earthen dam piping failure

A physically-based numerical three-dimensional earthen dam piping failure model is developed for homogeneous and zoned soil dams.This model is an erosion model,coupled with force/moment equilibrium analyses.Orifice flow and two-dimensional(2D)shallow water equations(SWE)are solved to simulate dam break flows at different breaching stages.Erosion rates of different soils with different construction compaction efforts are calculated using corresponding erosion formulae.The dam's real shape,soil properties,and surrounding area are programmed.Large outer 2D-SWE grids are used to control upstream and downstream hydraulic conditions and control the boundary conditions of orifice flow,and inner 2D-SWE flow is used to scour soil and perform force/moment equilibrium analyses.This model is validated using the European Commission IMPACT(Investigation of Extreme Flood Processes and Uncertainty)Test#5 in Norway,Teton Dam failure in Idaho,USA,and Quail Creek Dike failure in Utah,USA.All calculated peak outflows are within 10%errors of observed values.Simulation results show that,for a V-shaped dam like Teton Dam,a piping breach location at the abutment tends to result in a smaller peak breach outflow than the piping breach location at the dam's center;and if Teton Dam had broken from its center for internal erosion,a peak outflow of 117851 m'/s,which is 81%larger than the peak outflow of 65120 m3/s released from its right abutment,would have been released from Teton Dam.A lower piping inlet elevation tends to cause a faster/earlier piping breach than a higher piping inlet elevation.

Evaluation of flow regime of turbidity currents entering Dez Reservoir using extended shallow water model

In this study, the performance of the extended shallow water model （ESWM） in evaluation of the flow regime of turbidity currents entering the Dez Reservoir was investigated. The continuity equations for fluid and particles and the Navier-Stokes equations govern the entire flow of turbidity currents. The shallow water equations governing the flow of the depositing phase of turbidity currents are derived from these equations. A case study was conducted on the flow regime of turbidity currents entering the Dez Reservoir in Iran from January 2002 to July 2003. Facing a serious sedimentation problem, the dead storage of the Dez Reservoir will be full in the coming 10 years, and the inflowing water in the hydropower conduit system is now becoming turbid. Based on the values of the dimensionless friction number （ Nf ≤1 ） and dimensionless entrainment number （ NE≤ 1 ） of turbidity currents, and the coefficient of determination between the observed and predicted deposit depths （R2 = 0.86） for the flow regime of negligible friction and negligible entrainment （NFNE）, the flow regime of turbidity currents coming into the Dez Reservoir is considered to be NFNE. The results suggest that the ESWM is an appropriate approach for evaluation of the flow regime of turbidity currents in dam reservoirs where the characteristics of turbidity currents, such as the deposit depth, must be evaluated.

Applications of sediment sudden deposition model based on the third-generation numerical model for shallow water wave

The existing numerical models for nearshore waves are briefly introduced, and the third-generation numerical model for shallow water wave, which makes use of the most advanced productions of wave research and has been adapted well to be used in the environment of seacoast, lake and estuary area, is particularly discussed. The applied model realizes the significant wave height distribution at different wind directions. To integrate the model into the coastal area sediment, sudden deposition mechanism, the distribution of average silt content and the change of sediment sudden deposition thickness over time in the nearshore area are simulated. The academic productions can give some theoretical guidance to the applications of sediment sudden deposition mechanism for stormy waves in the coastal area. And the advancing directions of sediment sudden deposition model are prospected.

Modeling of Tsunami Generation and Propagation by a Spreading Curvilinear Seismic Faulting in Linearized Shallow-Water Wave Theory

The processes of tsunami evolution during its generation in search for possible amplification mechanisms resulting from unilateral spreading of the sea floor uplift is investigated. We study the nature of the tsunami build up and propagation during and after realistic curvilinear source models represented by a slowly uplift faulting and a spreading slip-fault model. The models are used to study the tsunami amplitude amplification as a function of the spreading velocity and rise time. Tsunami waveforms within the frame of the linearized shallow water theory for constant water depth are analyzed analytically by transform methods (Laplace in time and Fourier in space) for the movable source models. We analyzed the normalized peak amplitude as a function of the propagated uplift length, width and the average depth of the ocean along the propagation path.

Numerical Modeling of Shallow Water Table Behavior with Lisse Effect

Air entrapment is an important consideration in environments with shallow water tables and sandy soil, like the condition of highly conductive sandy soils and flat topography in Florida, USA. It causes water table rises in soils, which are significantly faster and higher than those in soils without air entrapment. Two numerical models, Integrated Hydrologic Model (IHM) and HYDRUS-1D (a single-phase, one-dimensional Richards′ equation model) were tested at an area of west central Florida to help further understanding the shallow water table behavior during a long term air entrapment. This investigation employed field data with two modeling approaches to quantify the variation of air pressurization values. It was found that the air pressurization effect was responsible at time up to 40 cm of water table rise being recorded by the observation well for these two models. The values of air pressurization calculated from IHM and HYDRUS-1D match the previously published values. Results also indicated that the two numerical models did not consider air entrapment effect (as the predictive parameters remain uncertain) and thus results of depth to water table from these models did not compare to the observations for these selected periods. Incorporating air entrapment in prediction models is critical to reproduce shallow water table observations.

Combined refraction and diffraction models for sea waves over complicated topography and with currents in shallow water

-Combined refraction and diffraction models in the form of linear parabolic approximation are derived through smallparameter method. More strictly theoretical basis and more accuracy in the models than Lozano's (1980) are obtained. Some theoretical defects in Liu's model (1985) with consideration of current are not only found but also eliminated. More strict and accurate models are, therefore, presented in this paper.The calculation results and analysis in applying the models to actual wave field with consideration of bottom friction will be given in the following paper.

A Regional Spectral Nested Shallow Water Equation Model

A method to expand meteorological elements in terms of finite double Fourier series in a limited-region and a spectral nested shallow water equation model based upon the method with conformal map projection in rectangular coordinates, have been proposed, and computational stability and efficiency of time integration have been discussed.

Application of the Characteristic CIP Method to a Shallow Water Model on the Sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the effciency of the characteristic constrained interpolation profile （CIP） method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

Study on water resisting property of subsurface aquiclude in shallow coal seam mining

According to the characteristic of clay aquiclude of overburden,the proper simulation materials and proportions of mixture for simulating the plastic clay aquiclude layers were developed,and the plastic similarity conditions were setup.Thus,the simili- tude and simulation method in whole stress-strain was progressed.Furthermore,the simulation condition and material proportions in water reaction property and crack closing property were also put forward.Based on systematical tests,the development and distri- bution of mined cracks in roof and subsurface aquiclude was found and the stability of aquiclude was analyzed at all.At last,the key section and key index of aquclude stability was advanced.It is found that the movement of clay aquiclude follows the movement of the underlying bedrock layers.The basic caving mechanism of the overburden roof strata was also presented.

Wave propagation speeds and source term influences in single and integral porosity shallow water equations

In urban flood modeling,so-called porosity shallow water equations(PSWEs),which conceptually account for unresolved structures, e.g.,buildings, are a promising approach to addressing high CPU times associated with state-of-the-art explicit numerical methods. The PSWE can be formulated with a single porosity term, referred to as the single porosity shallow water model(SP model), which accounts for both the reduced storage in the cell and the reduced conveyance, or with two porosity terms: one accounting for the reduced storage in the cell and another accounting for the reduced conveyance. The latter form is referred to as an integral or anisotropic porosity shallow water model(AP model). The aim of this study was to analyze the differences in wave propagation speeds of the SP model and the AP model and the implications of numerical model results. First, augmented Roe-type solutions were used to assess the influence of the source terms appearing in both models. It is shown that different source terms have different influences on the stability of the models. Second, four computational test cases were presented and the numerical models were compared. It is observed in the eigenvalue-based analysis as well as in the computational test cases that the models converge if the conveyance porosity in the AP model is close to the storage porosity. If the porosity values differ significantly, the AP model yields different wave propagation speeds and numerical fluxes from those of the BP model. In this study, the ratio between the conveyance and storage porosities was determined to be the most significant parameter.

Nonlinear Wave Force on Pier Group in Shallow Water

Based on the 1st order cnoidal wave theory, the wave diffraction around the pier group inshallow water is studied in this paper. The formulas for calculating the nonlinear wave forces are also presented here. In order to verify the theoretical results, model tests are conducted in the wave flume in The State Key Laboratory of Coastal and Offshore Engineering located in Dalian University of Technology. The range of the wave parameters in the experiments is characteristic wave period T g/d^(1/2) = 8.08- 22.86, characteristic wave height H/ d= 0.1 ~ 0.45. The results obtained from the experiments agree with the theoretical results quite well. It is shown that, in shallow water the nonlinear wave forces acting on a pier group are greater than those calculated by linear wave theory, the value of increment in wave force increases with the increases of the nonlinearity of the wave. In the wave range studied in this paper, the nonlinear wave force can reach over 4 times the force calculatecd by linear wave theory. Thus, it is suggested that, when Tg / d^(1/2)> 8, the wave force on the piers in the pier group in shallow water should be calculated by using the cnoidal wave theory.

Sedimentary environment and petrological features of organic-rich fine sediments in shallow water overlapping deposits: A case study of Cambrian Yuertus Formation in northwestern Tarim Basin, NW China

Taking the Cambrian Yuertus Formation outcrop profiles in the Aksu-Keping-Wushi areas of northwestern Tarim Basin as examples, the depositional environments of organic rich fine sediment were analyzed by examining the outcrop profiles macroscopically and microscopically. The study reveals that:(1) The lower part of the Yuertus Formation consists of organic-rich fine sediment or thin rhythmic interbeds of organic-rich fine sediment and siliceous sediment, the formation transforms to terrigenous diamictic grain shoal and inverse grading carbonate rocks upward.(2) The thin limestone interbedded with dark shale rhythmically has inverse grading.(3) The thin-bedded siliceous rock has metasomatic residual granular texture, stromatolithic structure and cementation fabric in vugs.(4) There are iron crust layers at the top of the shallowing diamictic grain shoal, beneath which exposed karst signs, such as karrens, dissolved fissures, sack-like vugs, near surface karst(plastic) breccia, breccia inside the karst system and terrigenous clastic fillings, can be seen.(5) Both the outcrops and seismic profiles show that organic-rich fine sediments above the unconformities or exposed surfaces are characterized by overlapping. The organic-rich fine sediment of the Cambrian Yuertus Formation was deposited in the anoxic-suboxidized restricted gulf lagoon environment, and its formation was controlled by high paleoproductivity and poor oxygen exchange jointly, then a shallow-water overlapping sedimentary model has been established. The results will help enrich and improve the sedimentary theory of organic-rich fine sediments.

Numerical Simulation of Nonlinear Surface Gravity Waves Transformation under Shallow-Water Conditions

This work considers the problems of numerical simulation of non-linear surface gravity waves transformation under shallow bay conditions. The discrete model is built from non-linear shallow-water equations. Are resulted boundary and initial conditions. The method of splitting into physical processes receives system from three equations. Then we define the approximation order and investigate stability conditions of the discrete model. The sweep method was used to calculate the system of equations. This work presents surface gravity wave profiles for different propagation phases.