In this study,the potential Kadomtsev-Petviashvili(pKP)equation,which describes the oblique interaction of surface waves in shallow waters,is solved by the new extended direct algebraic method.The results of the study...In this study,the potential Kadomtsev-Petviashvili(pKP)equation,which describes the oblique interaction of surface waves in shallow waters,is solved by the new extended direct algebraic method.The results of the study show that by applying the new direct algebraic method to the pKP equation,the behavior of the obliquely interacting surface waves in two dimensions can be analyzed.This article fairly clarifies the behaviors of surface waves in shallow waters.In the literature,several mathematical models have been developed in attempt to study these behaviors,with nonlinear mathematics being one of the most important steps;however,the investigations are still at a level that can be called‘baby steps’.Therefore,every study to be carried out in this context is of great importance.Thus,this study will serve as a reference to guide scientists working in this field.展开更多
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.I...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 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 e...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(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 t...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.展开更多
The geological conditions of shallow offshore delta oil reservoirs are complex. Under the condition of less well data and larger well spacing, the traditional reservoir configuration method is difficult to solve the d...The geological conditions of shallow offshore delta oil reservoirs are complex. Under the condition of less well data and larger well spacing, the traditional reservoir configuration method is difficult to solve the detailed study of such reservoirs in offshore oil fields. Based on the comprehensive analysis of the seismic phase, data of well log. The paper identifies criteria of the quaternary configuration boundary in shallow water delta of different types with distributary sand dam is established. At the same time, this paper used sensitive factor to construct the edge detection operator based on the amplitude attribute, characterizing the boundary of sand body thickness mutation or physical property mutation quantitatively, realizing the quantitative characterization of four-stage configuration boundary in the region with no wells or few wells, guiding the efficient development of offshore shallow water delta oilfield, and realizing the increase of storage and production of Bohai oilfield.展开更多
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 fl...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.展开更多
The geological conditions of offshore shallow water delta oil reservoirs were complex, with limited well data and large well spacing. Taking A Oilfield in the Bohai Sea Area, China as an example, the target sand body ...The geological conditions of offshore shallow water delta oil reservoirs were complex, with limited well data and large well spacing. Taking A Oilfield in the Bohai Sea Area, China as an example, the target sand body was formed in a shallow water delta sedimentary environment, with well-developed underwater distributary channels and frequent branching and diversion. The reservoir was strong non-uniformity and uneven plane water cut pressure. To this end, based on the existing work of predecessors, combined with seismic, logging, and production dynamics data, and based on the genesis mechanism of shallow water delta reservoirs, the boundary of composite river channels was identified through seismic facies, and logging facies were used to subdivide them into single river levels within the composite river channels. Then, seismic waveform characteristics were applied to track and characterize the plane distribution of single river channels, guiding the efficient development of offshore shallow water delta oil fields and achieving increased storage and production in Bohai Oilfield, China.展开更多
Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes...Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.展开更多
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 ener...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.展开更多
Much attention should be paid to a large FPSO moored permanently in an oil field with water depth of only about 20 m, since shallow water effects on the hydrodynamics may bring about collision and damage. A 160kDWT FP...Much attention should be paid to a large FPSO moored permanently in an oil field with water depth of only about 20 m, since shallow water effects on the hydrodynamics may bring about collision and damage. A 160kDWT FPSO with a permanent soft yoke mooring system is investigated with various shallow water depths and focuses are the low frequency surge motion and mooring load. Computation for the FPSO system is made based on linear 3-D potential fluid theory and time-domain numerical simulation method. Corresponding model test is carried out in the ocean engineering basin of Shanghai Jiao Tong University. It is shown that, in the surge natural period, low frequency surge motion and mooring force increase remarkably with the decrease of water depth. Especially, the smaller the ratio of water depth and draught is, the quicker the increase is. The shallow water effects should be taken into account carefully for determining the design load of a single point mooring system.展开更多
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 porosi...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.展开更多
A semi-implicit and Eulerian - Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-moment turbulence closure model and tr...A semi-implicit and Eulerian - Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-moment turbulence closure model and transport equations of salinity and temperature. The simulation for flooding and drying of mudflats has been improved. The model is applied to Xiamen waters. Based on extensive survey data, water level elevation, temperature and salinity field along the eastern open boundary and at the Jiulong River inlets and runoffs are analyzed, specified and calibrated. The computed results show good agreement with the measured data, reproduce flooding, emergence of large and complex mudflat region.展开更多
The hybrid finite analytic(HFA) method is a kind of numerical scheme in rectangular element. In order to simulate the shallow circulation in irregular bathymetry by HFA scheme, the model in sigma coordinate system was...The hybrid finite analytic(HFA) method is a kind of numerical scheme in rectangular element. In order to simulate the shallow circulation in irregular bathymetry by HFA scheme, the model in sigma coordinate system was obtained. The model has been tested against three cases: 1) Wind induced circulation; 2) Density driven circulation and 3) Seiche oscillation. The results obtained in the present study compare well with those obtained from the corresponding analytical solutions under idealized for the above three cases. The hybrid finite analytic method and the circulation model in sigma coordinate system can be used calculate the flow and water quality in estuaries and coastal waters.展开更多
Vertical errors often present in multibeam swath bathymetric data. They are mainly sourced by sound refraction, internal wave disturbance, imperfect tide correction, transducer mounting, long period heave, static draf...Vertical errors often present in multibeam swath bathymetric data. They are mainly sourced by sound refraction, internal wave disturbance, imperfect tide correction, transducer mounting, long period heave, static draft change, dynamic squat and dynamic motion residuals, etc. Although they can be partly removed or reduced by specific algorithms, the synthesized depth biases are unavoidable and sometimes have an important influence on high precise utilization of the final bathymetric data. In order to. confidently identify the decimeter-level changes in seabed morphology by MBES, we must remove or weaken depth biases and improve the precision of multibeam bathymetry further. The fixed-interval profiles that are perpendicular to the vessel track are generated to adjust depth biases between swaths. We present a kind of postprocessing method to minimize the depth biases by the histogram of cumulative depth biases. The datum line in each profile can be obtained by the maximum value of histogram. The corrections of depth biases can be calculated according to the datum line. And then the quality of final bathymetry can be improved by the corrections. The method is verified by a field test.展开更多
Recently, a new (2+1)-dimensional shallow water wave system, the (2+1)-dlmenslonal displacement shallow water wave system (2DDSWWS), was constructed by applying the variational principle of the analytic mechan...Recently, a new (2+1)-dimensional shallow water wave system, the (2+1)-dlmenslonal displacement shallow water wave system (2DDSWWS), was constructed by applying the variational principle of the analytic mechanics in the Lagrange coordinates. The disadvantage is that fluid viscidity is not considered in the 2DDSWWS, which is the same as the famous Kadomtsev-Petviashvili equation and Korteweg-de Vries equation. Applying dimensional analysis, we modify the 2DDSWWS and add the term related to the fluid viscidity to the 2DDSWWS. The approximate similarity solutions of the modified 2DDSWWS (M2DDSWWS) is studied and four similarity solutions are obtained. For the perfect fluids, the coefficient of kinematic viscosity is zero, then the M2DDSWWS will degenerate to the 2DDSWWS.展开更多
In this paper, a nonlinear model is presented to describe wave transformation in shallow water with the zero- vorticity equation of wave- number vector and energy conservation equation. The nonlinear effect due to an ...In this paper, a nonlinear model is presented to describe wave transformation in shallow water with the zero- vorticity equation of wave- number vector and energy conservation equation. The nonlinear effect due to an empirical dispersion relation (by Hedges) is compared with that of Dalrymple's dispersion relation. The model is tested against the laboratory measurements for the case of a submerged elliptical shoal on a slope beach, where both refraction and diffraction are significant. The computation results, compared with those obtained through linear dispersion relation, show that the nonlinear effect of wave transformation in shallow water is important. And the empirical dispersion relation is suitable for researching the nonlinearity of wave in shallow water.展开更多
The habitat suitability index(HSI) model was used to identify potential sites for sustainable restoration of ark shell, Scapharca subcrenata(Lischke), in the shallow water of Xiaoheishan Island, using a geographic inf...The habitat suitability index(HSI) model was used to identify potential sites for sustainable restoration of ark shell, Scapharca subcrenata(Lischke), in the shallow water of Xiaoheishan Island, using a geographic information system framework. The seven input variables of the HSI model were sediment composition, water temperature, salinity, dissolved oxygen, water depth, p H, and ammonia. A non-linear suitability function for each variable factor was used to transform the value into a normalized quality index ranging from 0(nonsuitability) to 1(best suitability). In present study, the analysis of habitat suitability was conducted for four seasons respectively. The majority of the study area has a high HSI value(>0.6) year round, which implies a strong suitability for restoration, with the optimal habitat located on the eastern side of the island.Correspondence analysis indicated that water temperature was the main factor causing seasonal variation,while sediment composition and water depth were the two major reasons for the differences in sites. The results of this work could provide support for restoration decision making through identification of potential sites for sustainable establishment of S. subcrenata.展开更多
To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying (WD) method is proposed. In the new method, a "predicted water depth" is evaluated e...To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying (WD) method is proposed. In the new method, a "predicted water depth" is evaluated explicitly based on the simplified shallow water equations and used to determine the status (wet or dry) together with the direction of flow. Compared with previous WD method, besides the water elevation, more factors, such as the flow velocity and the surface shear stress, are taken into account in the new method to determine the moving boundary. In addition, a formula is deduced to determine the threshold, as critical water depth, which needs to be preset before simulations. The new WD method is tested with five cases including three 1D ones and two 2D ones. The results show that the new WD method can simulate the wetting and drying process, in beth typical and practical cases, with smooth manner and achieves effective estimation of the retention volume at shallow water body.展开更多
A particular porosity method named "slot method" is implemented in a depth-integrated shallow water flow model (DIVAST) to simulate wetting and drying processes. Discussed is the relationship between the shape fac...A particular porosity method named "slot method" is implemented in a depth-integrated shallow water flow model (DIVAST) to simulate wetting and drying processes. Discussed is the relationship between the shape factors of the "slot" and the preset depth used in "wetting-drying" algorithm. Two typical tests are conducted to examine the performance of the method with the effect of the shape factors of the "slot" being checked in detail in the first test. Numerical results demonstrate that: 1 ) no additional effort to improve the finite difference scheme is needed to implement "slot method" in DIVAST, and 2) "slot method" will simulate wetting and diying processes correctly if the shape factors of the "slot" being selected properly.展开更多
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...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.展开更多
文摘In this study,the potential Kadomtsev-Petviashvili(pKP)equation,which describes the oblique interaction of surface waves in shallow waters,is solved by the new extended direct algebraic method.The results of the study show that by applying the new direct algebraic method to the pKP equation,the behavior of the obliquely interacting surface waves in two dimensions can be analyzed.This article fairly clarifies the behaviors of surface waves in shallow waters.In the literature,several mathematical models have been developed in attempt to study these behaviors,with nonlinear mathematics being one of the most important steps;however,the investigations are still at a level that can be called‘baby steps’.Therefore,every study to be carried out in this context is of great importance.Thus,this study will serve as a reference to guide scientists working in this field.
文摘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.
基金the National Natural Science Foundation of China(11901555,11871448,12001009).
文摘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.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11975156 and 12175148)。
文摘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.
文摘The geological conditions of shallow offshore delta oil reservoirs are complex. Under the condition of less well data and larger well spacing, the traditional reservoir configuration method is difficult to solve the detailed study of such reservoirs in offshore oil fields. Based on the comprehensive analysis of the seismic phase, data of well log. The paper identifies criteria of the quaternary configuration boundary in shallow water delta of different types with distributary sand dam is established. At the same time, this paper used sensitive factor to construct the edge detection operator based on the amplitude attribute, characterizing the boundary of sand body thickness mutation or physical property mutation quantitatively, realizing the quantitative characterization of four-stage configuration boundary in the region with no wells or few wells, guiding the efficient development of offshore shallow water delta oilfield, and realizing the increase of storage and production of Bohai oilfield.
文摘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.
文摘The geological conditions of offshore shallow water delta oil reservoirs were complex, with limited well data and large well spacing. Taking A Oilfield in the Bohai Sea Area, China as an example, the target sand body was formed in a shallow water delta sedimentary environment, with well-developed underwater distributary channels and frequent branching and diversion. The reservoir was strong non-uniformity and uneven plane water cut pressure. To this end, based on the existing work of predecessors, combined with seismic, logging, and production dynamics data, and based on the genesis mechanism of shallow water delta reservoirs, the boundary of composite river channels was identified through seismic facies, and logging facies were used to subdivide them into single river levels within the composite river channels. Then, seismic waveform characteristics were applied to track and characterize the plane distribution of single river channels, guiding the efficient development of offshore shallow water delta oil fields and achieving increased storage and production in Bohai Oilfield, China.
文摘Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.
基金"333"Project Scientific Research Foundation of Jiangsu ProvinceScience Fundation of Hohai University(3853)
文摘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.
基金This work was financially supported bythe National High Technology and Development Programof China (Grant No.2004AA616180)
文摘Much attention should be paid to a large FPSO moored permanently in an oil field with water depth of only about 20 m, since shallow water effects on the hydrodynamics may bring about collision and damage. A 160kDWT FPSO with a permanent soft yoke mooring system is investigated with various shallow water depths and focuses are the low frequency surge motion and mooring load. Computation for the FPSO system is made based on linear 3-D potential fluid theory and time-domain numerical simulation method. Corresponding model test is carried out in the ocean engineering basin of Shanghai Jiao Tong University. It is shown that, in the surge natural period, low frequency surge motion and mooring force increase remarkably with the decrease of water depth. Especially, the smaller the ratio of water depth and draught is, the quicker the increase is. The shallow water effects should be taken into account carefully for determining the design load of a single point mooring system.
基金supported by the National Natural Science Foundation of China (Grants No. 50909065 and 51109039)the National Basic Research Program of China (973 Program, Grant No. 2012CB417002)
文摘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.
文摘A semi-implicit and Eulerian - Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-moment turbulence closure model and transport equations of salinity and temperature. The simulation for flooding and drying of mudflats has been improved. The model is applied to Xiamen waters. Based on extensive survey data, water level elevation, temperature and salinity field along the eastern open boundary and at the Jiulong River inlets and runoffs are analyzed, specified and calibrated. The computed results show good agreement with the measured data, reproduce flooding, emergence of large and complex mudflat region.
文摘The hybrid finite analytic(HFA) method is a kind of numerical scheme in rectangular element. In order to simulate the shallow circulation in irregular bathymetry by HFA scheme, the model in sigma coordinate system was obtained. The model has been tested against three cases: 1) Wind induced circulation; 2) Density driven circulation and 3) Seiche oscillation. The results obtained in the present study compare well with those obtained from the corresponding analytical solutions under idealized for the above three cases. The hybrid finite analytic method and the circulation model in sigma coordinate system can be used calculate the flow and water quality in estuaries and coastal waters.
基金supported by the project of Shandong Provincial National Science Foundation for Distinguished Young Scholars(Grant No.JQ201113)SDUST's National Science Foundation for Distinguished Young Scholars(Grant No.2010KYJQ102)
文摘Vertical errors often present in multibeam swath bathymetric data. They are mainly sourced by sound refraction, internal wave disturbance, imperfect tide correction, transducer mounting, long period heave, static draft change, dynamic squat and dynamic motion residuals, etc. Although they can be partly removed or reduced by specific algorithms, the synthesized depth biases are unavoidable and sometimes have an important influence on high precise utilization of the final bathymetric data. In order to. confidently identify the decimeter-level changes in seabed morphology by MBES, we must remove or weaken depth biases and improve the precision of multibeam bathymetry further. The fixed-interval profiles that are perpendicular to the vessel track are generated to adjust depth biases between swaths. We present a kind of postprocessing method to minimize the depth biases by the histogram of cumulative depth biases. The datum line in each profile can be obtained by the maximum value of histogram. The corrections of depth biases can be calculated according to the datum line. And then the quality of final bathymetry can be improved by the corrections. The method is verified by a field test.
基金Project supported by the Natural Science Foundation of Guangdong Province of China (Grant No.10452840301004616)the National Natural Science Foundation of China (Grant No.61001018)the Scientific Research Foundation for the Doctors of University of Electronic Science and Technology of China Zhongshan Institute (Grant No.408YKQ09)
文摘Recently, a new (2+1)-dimensional shallow water wave system, the (2+1)-dlmenslonal displacement shallow water wave system (2DDSWWS), was constructed by applying the variational principle of the analytic mechanics in the Lagrange coordinates. The disadvantage is that fluid viscidity is not considered in the 2DDSWWS, which is the same as the famous Kadomtsev-Petviashvili equation and Korteweg-de Vries equation. Applying dimensional analysis, we modify the 2DDSWWS and add the term related to the fluid viscidity to the 2DDSWWS. The approximate similarity solutions of the modified 2DDSWWS (M2DDSWWS) is studied and four similarity solutions are obtained. For the perfect fluids, the coefficient of kinematic viscosity is zero, then the M2DDSWWS will degenerate to the 2DDSWWS.
文摘In this paper, a nonlinear model is presented to describe wave transformation in shallow water with the zero- vorticity equation of wave- number vector and energy conservation equation. The nonlinear effect due to an empirical dispersion relation (by Hedges) is compared with that of Dalrymple's dispersion relation. The model is tested against the laboratory measurements for the case of a submerged elliptical shoal on a slope beach, where both refraction and diffraction are significant. The computation results, compared with those obtained through linear dispersion relation, show that the nonlinear effect of wave transformation in shallow water is important. And the empirical dispersion relation is suitable for researching the nonlinearity of wave in shallow water.
基金The National Natural Science Foundation of China under contract No.41206102the Program of the Key Laboratory of Marine Ecology and Environmental Science and Engineering,State Oceanic Administration under contract No.MESE-2013-01
文摘The habitat suitability index(HSI) model was used to identify potential sites for sustainable restoration of ark shell, Scapharca subcrenata(Lischke), in the shallow water of Xiaoheishan Island, using a geographic information system framework. The seven input variables of the HSI model were sediment composition, water temperature, salinity, dissolved oxygen, water depth, p H, and ammonia. A non-linear suitability function for each variable factor was used to transform the value into a normalized quality index ranging from 0(nonsuitability) to 1(best suitability). In present study, the analysis of habitat suitability was conducted for four seasons respectively. The majority of the study area has a high HSI value(>0.6) year round, which implies a strong suitability for restoration, with the optimal habitat located on the eastern side of the island.Correspondence analysis indicated that water temperature was the main factor causing seasonal variation,while sediment composition and water depth were the two major reasons for the differences in sites. The results of this work could provide support for restoration decision making through identification of potential sites for sustainable establishment of S. subcrenata.
基金supported by the National Natural Science Foundation of China (Grant No.10872144)the Support Plan of Science and Technology of Tianjin (Grant No.07ZCGYSH01700)
文摘To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying (WD) method is proposed. In the new method, a "predicted water depth" is evaluated explicitly based on the simplified shallow water equations and used to determine the status (wet or dry) together with the direction of flow. Compared with previous WD method, besides the water elevation, more factors, such as the flow velocity and the surface shear stress, are taken into account in the new method to determine the moving boundary. In addition, a formula is deduced to determine the threshold, as critical water depth, which needs to be preset before simulations. The new WD method is tested with five cases including three 1D ones and two 2D ones. The results show that the new WD method can simulate the wetting and drying process, in beth typical and practical cases, with smooth manner and achieves effective estimation of the retention volume at shallow water body.
基金the National Natural Science Foundation of China (Grant No.10702050)the Natural Science Foundation of Tianjin (Grant No.07JCYBJC07500)the Support Plan of Science and Technology of Tianjin (Grant No.07ZCGYSH01700)
文摘A particular porosity method named "slot method" is implemented in a depth-integrated shallow water flow model (DIVAST) to simulate wetting and drying processes. Discussed is the relationship between the shape factors of the "slot" and the preset depth used in "wetting-drying" algorithm. Two typical tests are conducted to examine the performance of the method with the effect of the shape factors of the "slot" being checked in detail in the first test. Numerical results demonstrate that: 1 ) no additional effort to improve the finite difference scheme is needed to implement "slot method" in DIVAST, and 2) "slot method" will simulate wetting and diying processes correctly if the shape factors of the "slot" being selected properly.
文摘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.