Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting...Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.展开更多
A new scour countermeasure using solidified slurry for offshore foundation has been proposed recently.Fluidized solidified slurry is pumped to seabed area around foundation for scour protection or pumped into the deve...A new scour countermeasure using solidified slurry for offshore foundation has been proposed recently.Fluidized solidified slurry is pumped to seabed area around foundation for scour protection or pumped into the developed scour holes for scour repair as the fluidized material solidifies gradually.In the pumping operation and solidification,the engineering behaviors of solidified slurry require to be considered synthetically for the reliable application in scour repair and protection of ocean engineering such as the pumpability related flow value,flow diffusion behavior related rheological property,anti-scour performance related retention rate in solidification and bearing capacity related strength property after solidification.In this study,a series of laboratory tests are conducted to investigate the effects of mix proportion(initial water content and binder content)on the flow value,rheological properties,density,retention rate of solidified slurry and unconfined compressive strength(UCS).The results reveal that the flow value increases with the water content and decreases with the binder amount.All the solidified slurry exhibits Bingham plastic behavior when the shear rate is larger than 5 s^(-1).The Bingham model has been employed to fit the rheology test results,and empirical formulas for obtaining the density,yield stress and viscosity are established,providing scientific support for the numerical assessment of flow and diffusion of solidified slurry.Retention rate of solidified slurry decreases with the water flow velocity and flow value,which means the pumpability of solidified slurry is contrary to anti-scour performance.The unconfined compressive strength after solidification reduces as the water content increases and binder content decreases.A design and application procedure of solidified soil for scour repair and protection is also proposed for engineering reference.展开更多
Ice-going ships play a crucial role in polar transportation and resource extraction.Different from the existing modeling approach which assumes that ships remain stationary,dynamic overset grid technology and DFBI(Dyn...Ice-going ships play a crucial role in polar transportation and resource extraction.Different from the existing modeling approach which assumes that ships remain stationary,dynamic overset grid technology and DFBI(Dynamic Fluid-Body Interaction)method are employed in this paper to enable the free-running motion of the ship in modeling.A numerical model capable of simulating a ship navigating through pack ice area is proposed,which uses Computational Fluid Dynamics(CFD)method to solve the flow field and applies the Discrete Element Method(DEM)to simulate ship-ice and ice-ice interactions.Besides,the proposed high-precision method for generating pack ice area can be used in conjunction with the proposed numerical model.By comparing the numerical results with the available model test data and experimental observations,the effectiveness of the numerical model is validated,demonstrating its strong capability of predicting resistance and simulating ship navigation in pack ice,as well as its significant potential and applicability for further studies.展开更多
To study the rolling motion of a ship in the presence of water on its deck,a linear-plus-quadratic damping term was incorporated into its equation of motion.Ship model tests indicates that the key dynamics of the phys...To study the rolling motion of a ship in the presence of water on its deck,a linear-plus-quadratic damping term was incorporated into its equation of motion.Ship model tests indicates that the key dynamics of the physical system are preserved in the ship rolling equation with the linear-plus-quadratic type damping term.To take into account the presence of randomness in the excitation and the response,a new method was developed and a Melnikov criterion was obtained to provide an upper bound on the domain of the potential chaotic rolling motion(erratic rocking).Additionally,the Melnikov criterion proposed in this study was verified by the utilization of phase plane diagrams and Poincare maps.Furthermore,this research has made the initial endeavor to systematically modify the system parameters in the rolling equation of motion for ship stability analysis.展开更多
Floating breakwaters(FBs)are commonly employed for the protection of coastal installations.In this work,a convextype floating breakwater(FB)is proposed,and its hydrodynamic characteristics are studied through systemat...Floating breakwaters(FBs)are commonly employed for the protection of coastal installations.In this work,a convextype floating breakwater(FB)is proposed,and its hydrodynamic characteristics are studied through systematic laboratory experiments.Two different deck widths and two different mooring systems are set in the experiment.The transmission coefficients,reflection coefficients,motion responses and mooring forces of convex-type FBs are obtained in experiments.The influences of the deck width and mooring system on the hydrodynamic characteristics of the proposed FB are analyzed and compared.The experimental results show that the reflection coefficient and mooring force of the convex-type FB with a cross-mooring system are significantly larger than those of the convextype FB with a parallel-mooring system.A convex-type FB with a larger deck width has a higher reflection coefficient.The convex-type FBs with cross-and parallel-mooring systems have similar surge and heave motions,but the cross-mooring results in small roll motion.In addition,reliable prediction formulas for the transmission coefficient of convex-type FBs with different mooring systems have been developed,which are important for engineering design.展开更多
This study investigates the hydrodynamic process of a cylinder ejected from the water’s surface through high-speed camera experiments.Using digital image processing methods,the images obtained through experiments are...This study investigates the hydrodynamic process of a cylinder ejected from the water’s surface through high-speed camera experiments.Using digital image processing methods,the images obtained through experiments are processed and analyzed.Although the dynamics of rising buoyant cylinders have been thoroughly investigated,the pop-up height of the cylinders has not been extensively explored.Statistical analysis of the kinematic and dynamic data of cylinders is conducted.Research has shown that after the cylinder rises,it pops out of the water’s surface.Within the experimental range,the pop-up height of the cylinder is related to the release depth.Furthermore,the pop-up height and release depth of the cylinder vary linearly under vertical release conditions.Under horizontal release conditions,the relationship between pop-up height and release depth shows irregular changes mainly because of the unstable shedding of the wake vortex.展开更多
This paper investigated the use of magnesium phosphate cement (MPC) for solidifying sludge with different humic acid (HA) content (ranging from 0 to 4.5%) and explored the solidification mechanism. Fluidity, setting t...This paper investigated the use of magnesium phosphate cement (MPC) for solidifying sludge with different humic acid (HA) content (ranging from 0 to 4.5%) and explored the solidification mechanism. Fluidity, setting time, unconfined compressive strength (UCS), the strength formation mechanism, and the spontaneous imbibition process of solidified sludge (SS) were studied. The results indicate that MPC can be used as a low-alkalinity curing agent. As the HA content increases, fluidity and setting time also increase, while hydration temperature and strength decrease. Additionally, the failure mode of SS transitions from brittleness to ductility. The strength of SS is composed of the cementation strength provided by MPC hydration products, matric suction, osmotic suction, and the structural strength of the sludge. MPC reduces the structural strength caused by the shrinkage of pure sludge under the action of matric suction, but the incorporation of MPC significantly improved the strength when the sludge is eroded by water. X-ray diffraction (XRD) and scanning electron microscopy (SEM) show that the sludge and MPC can form a dense solid body, forming various hydration products, and synergistically improve the mechanical properties of the sludge.展开更多
Air blasts induced by rock-ice avalanches are common natural phenomena known for their far-field destructive impact.In this study,remote sensing images,eyewitness videos and numerical modelling were comprehensively ap...Air blasts induced by rock-ice avalanches are common natural phenomena known for their far-field destructive impact.In this study,remote sensing images,eyewitness videos and numerical modelling were comprehensively applied to analyze the initiation and propagation of the 2021 Chamoli avalancheinduced air blast.Our findings indicate that air blasts are observed from the avalanche source area to the Rishiganga valley,but nearly disappear in the Dhauliganga valley.The most intense air blast is concentrated on the left side of Ronti Gad valley,with maximum velocity and pressure estimated at over 70 m/s and 20 kPa,respectively.Such high pressure results in widespread tree breakage in the area.Based on the analysis of the Chamoli event,we further discussed the potential contribution of the avalanche flow regime,avalanche dynamics and geomorphology to the destructive potential of air blasts.Rapidly moved sliding mass can impart the air blast a high initial momentum,and this process will be exaggerated when the avalanche impacts valley walls at bends.However,when the rock-ice avalanche transforms into a debris-enriched flash flood,free water within the flowing mass can displace air,inhibiting the generation of air blasts.Our work offers new insights into the generation and propagation of rock-ice avalanche-induced air blasts,underscoring the importance of including this type of hazard during avalanche risk assessment in high-altitude glacial regions.展开更多
The integrated systems of unmanned surface vehicles(USVs) and remotely operated vehicles(ROVs) have been extensively applied in marine exploration and seabed coverage. However, the simultaneous navigation of USV-ROV s...The integrated systems of unmanned surface vehicles(USVs) and remotely operated vehicles(ROVs) have been extensively applied in marine exploration and seabed coverage. However, the simultaneous navigation of USV-ROV systems is frequently limited by strong disturbances induced by waves or currents. This paper develops a novel rigidflexible coupling multibody dynamic model that incorporates disturbances of variable-length marine cables with geometrically nonlinear motion. A hybrid Lagrangian-Eulerian absolute nodal coordinate formulation(ANCF) element is developed to accurately model subsea cables which undergo significant overall motion, substantial deformation,and mass flow during the deployment of underwater equipment. Furthermore, the governing equations of the coupled USV-umbilical-ROV system are derived, considering wave-induced forces and current disturbances. A numerical solver based on the Newmark-beta method is proposed, along with an adaptive meshing technique near the release point. After validating three experimental cases, the cable disturbances at both the USV and ROV ends—caused by ocean currents, heave motion, and simultaneous navigation—are comprehensively compared and evaluated. Finally,it is demonstrated that a PD controller with disturbance compensation can enhance the simultaneous navigation performance of USV-ROV systems.展开更多
The oscillating hydrofoil represents a promising technology for harvesting energy from tidal currents.While previous research has primarily focused on oscillating hydrofoils utilizing a fully activated control strateg...The oscillating hydrofoil represents a promising technology for harvesting energy from tidal currents.While previous research has primarily focused on oscillating hydrofoils utilizing a fully activated control strategy,the industry predominantly employs a semi-activated control strategy in existing tidal current energy converters.It is essential to identify the differences in predicted energy-harvesting performance between these two controlling strategies through experimental modeling or numerical studies.Furthermore,the suitability of the fully activated control strategy in predicting the energy-harvesting capabilities of oscillating hydrofoils is evaluated.The 2D numerical models of hydrofoil based on fully activated and semi-activated control strategies have been developed and validated.The amplitudes of heaving and pitching movements for the fully activated hydrofoil are determined to match those of the semi-activated hydrofoil.The results show that the main difference between the two control strategies lies in the phase shift occurring between the pitching and heaving motions.This phase shift affects the lift force and its coordination with the heaving velocity,which in turn affects the power output.Notably,the maximum relative efficiency difference obtained between the fully activated and semi-activated control strategies can reach 191%.展开更多
Comprehensive investigations have been conducted to study the structure and overconsolidation of upper Shanghai clays, i.e. Layers 2–6 clays, typically located at depths of 30–40 m. However, limited information is a...Comprehensive investigations have been conducted to study the structure and overconsolidation of upper Shanghai clays, i.e. Layers 2–6 clays, typically located at depths of 30–40 m. However, limited information is available on their anisotropy, and even less is known about the correlation between structure, overconsolidation, and anisotropy. In this study, the undrained anisotropy characteristics of shear strength and small-strain shear stiffness in upper Shanghai Layers 2–6 clays were thoroughly assessed using a series of K0-consolidated undrained triaxial compression (TC) and triaxial extension (TE) tests (K0 is the coefficient of lateral earth pressure at rest). The effective stress paths, shear strength, and small-strain shear stiffness from the undrained TC and TE tests demonstrate the anisotropic behaviors in upper Shanghai clays. Analyses of data from upper Shanghai clays and other clays worldwide indicate that the shear strength anisotropy ratio (Ks) converges at 0.8 as the overconsolidation ratio (OCR) and plasticity index (Ip) increase, while the small-strain shear stiffness anisotropy ratio (Re) converges at 1.0. The influence of OCR on Ks and Re is more pronounced than that of Ip and sensitivity (St). Nevertheless, no clear correlation between Ks and Re is observed in upper Shanghai clays.展开更多
Veno-arterial extracorporeal membrane oxygenation(VA-ECMO)is a life support intervention for patients with refractory cardiogenic shock or severe cardiopulmonary failure.However,the choice of cannulation strategy rema...Veno-arterial extracorporeal membrane oxygenation(VA-ECMO)is a life support intervention for patients with refractory cardiogenic shock or severe cardiopulmonary failure.However,the choice of cannulation strategy remains contentious,partly due to insufficient understanding of hemodynamic characteristics associated with the site of arterial cannulation.In this study,a geometrical multiscale model was built to offer a mathematical tool for addressing the issue.The outflow cannula of ECMO was inserted into the ascending aorta in the case of central cannulation,whereas it was inserted into the right subclavian artery(RSA)or the left iliac artery(LIA)in the case of peripheral cannulation.Numerical simulations conducted on three patient-specific aortas demonstrated that the central cannulation outperformed the two types of peripheral cannulation in evenly delivering ECMO flow to branch arteries.Both the central and RSA cannulations could maintain an approximately normal hemodynamic state in the aortas,although the area of aortic walls exposed to abnormal wall shear stress(WSS)was considerably enlarged in comparison with the normal physiological condition.In contrast,the LIA cannulation not only led to insufficient delivery of ECMO flow to the right upper body(with ECMO flow fractions<0.5),but also induced marked flow disturbance in the aorta,causing about 40%of the abdominal aortic wall and over 65%of the resting aortic wall to suffer from high time-averaged WSS(>5 Pa)and low time-averaged WSS(<0.4 Pa),respectively.The LIA cannulation also resulted in significantly prolonged blood residence time(>40 s)in the ascending aorta,which,along with abnormal WSS,may considerably increase the risk of thrombosis.In summary,our numerical study elucidated the impact of arterial cannulation site in VA-ECMO intervention on aortic hemodynamics and ECMO flow distribution.The findings provide compensatory biomechanical information for traditional clinical studies and may serve as a theoretical reference for guiding the evaluation and selection of cannulation strategies in clinical practice.展开更多
Knowledge-Based Engineering （KBE） is introduced into the ship structural design in this paper. From the implementation of KBE, the design solutions for both Rules Design Method （RDM） and Interpolation Design Meth...Knowledge-Based Engineering （KBE） is introduced into the ship structural design in this paper. From the implementation of KBE, the design solutions for both Rules Design Method （RDM） and Interpolation Design Method （IDM） are generated. The corresponding Finite Element （FE） models are generated. Topological design of the longitudinal structures is studied where the Gaussian Process （GP） is employed to build the surrogate model for FE analysis. Multi-objective optimization methods inspired by Pareto Front are used to reduce the design tank weight and outer surface area simultaneously. Additionally, an enhanced Level Set Method （LSM） which employs implicit algorithm is applied to the topological design of typical bracket plate which is used extensively in ship structures. Two different sets of boundary conditions are considered. The proposed methods show satisfactory efficiency and accuracy.展开更多
The phase change characteristic of the power source of an underwater glider propelled by the ocean＇s thermal energy is the key factor in glider attitude control. A numerical model has been established based on the en...The phase change characteristic of the power source of an underwater glider propelled by the ocean＇s thermal energy is the key factor in glider attitude control. A numerical model has been established based on the enthalpy method to analyze the phase change heat transfer process under convective boundary conditions. Phase change is not an isothermal process, but one that occurs at a range of temperature. The total melting time of the material is very sensitive to the surrounding temperature. When the temperature of the surroundings decreases 8 degrees, the total melting time increases 1.8 times. But variations in surrounding temperature have little effect on the initial temperature of phase change, and the slope of the temperature time history curve remains the same. However, the temperature at which phase change is completed decreases significantly. Our research shows that the phase change process is also affected by container size, boundary conditions, and the power source＇s cross sectional area. Materials stored in 3 cylindrical containers with a diameter of 38ram needed the shortest phase change time. Our conclusions should be helpful in effective design of underwater glider power systems.展开更多
The size and shape of the effective test area are crucial to consider when short-crested waves are created by segmented wavemakers. The range of the effective test area of short-crested waves simulated by two-sided se...The size and shape of the effective test area are crucial to consider when short-crested waves are created by segmented wavemakers. The range of the effective test area of short-crested waves simulated by two-sided segmented wavemakers is analyzed in this paper. The experimental investigation on the wave field distribution of short-crested waves generated by two-sided segmented wavemakers is conducted by using an array of wave gauges. Wave spectra and directional spreading function are analyzed and the results show that when the main direction is at a certain angle with the normal line of wave generators, the wave field of 3D short-crested waves generated by two-sided segmented wavemakers has good spatial uniformity within the model test area. The effective test area can provide good wave environments for seakeeping model tests of various ocean engineering structures in the deep ocean engineering basin.展开更多
Research on the diffusion characteristics of swells contributes positively to wave energy forecasting, swell monitoring, and early warning. In this work, the South Indian Ocean westerly index(SIWI) and Indian Ocean sw...Research on the diffusion characteristics of swells contributes positively to wave energy forecasting, swell monitoring, and early warning. In this work, the South Indian Ocean westerly index(SIWI) and Indian Ocean swell diffusion effect index(IOSDEI) are defined on the basis of the 45-year(September 1957–August 2002) ERA-40 wave reanalysis data from the European Centre for Medium-Range Weather Forecasts(ECMWF) to analyze the impact of the South Indian Ocean westerlies on the propagation of swell acreage. The following results were obtained: 1) The South Indian Ocean swell mainly propagates from southwest to northeast. The swell also spreads to the Arabian Sea upon reaching low-latitude waters. The 2.0-meter contour of the swell can reach northward to Sri Lankan waters. 2) The size of the IOSDEI is determined by the SIWI strength. The IOSDEI requires approximately 2–3.5 days to fully respond to the SIWI. The correlations between SIWI and IOSDEI show obvious seasonal differences, with the highest correlations found in December–January–February(DJF) and the lowest correlations observed in June–July–August(JJA). 3) The SIWI and IOSDEI have a common period of approximately 1 week in JJA and DJF. The SIWI leads by approximately 2–3 days in this common period.展开更多
Short-term forecasts of wave energy play a key role in the daily operation,maintenance planning,and electrical grid operation of power farms.In this study,we propose a short-term wave energy forecast scheme and use th...Short-term forecasts of wave energy play a key role in the daily operation,maintenance planning,and electrical grid operation of power farms.In this study,we propose a short-term wave energy forecast scheme and use the North Indian Ocean(NIO)as a case study.Compared with the traditional forecast scheme,our proposed scheme considers more forecast elements.In addition to the traditional short-term forecast factors related to wave energy(wave power,significant wave height(SWH),wave period),our scheme emphasizes the forecast of a series of key factors that are closely related to the effectiveness of the energy output,capture efficiency,and conversion efficiency.These factors include the available rate,total storage,effective storage,co-occurrence of wave power-wave direction,co-occurrence of the SWH-wave period,and the wave energy at key points.In the regional nesting of nu-merical simulations of wave energy in the NIO,the selection of the southern boundary is found to have a significant impact on the simulation precision,especially during periods of strong swell processes of the South Indian Ocean(SIO)westerly.During tropical cyclone‘VARDAH’in the NIO,as compared with the simulation precision obtained with no expansion of the southern boundary(scheme-1),when the southern boundary is extended to the tropical SIO(scheme-2),the improvement in simulation precision is significant,with an obvious increase in the correlation coefficient and decrease in error.In addition,the improvement is much more significant when the southern boundary extends to the SIO westerly(scheme-3).In the case of strong swell processes generated by the SIO westerly,the improvement obtained by scheme-3 is even more significant.展开更多
To ensure the safe performance of deep-sea mining vehicles(DSMVs),it is necessary to study the mechanical characteristics of the interaction between the seabed soil and the track plate.The rotation and digging motions...To ensure the safe performance of deep-sea mining vehicles(DSMVs),it is necessary to study the mechanical characteristics of the interaction between the seabed soil and the track plate.The rotation and digging motions of the track plate are important links in the contact between the driving mechanism of the DSMV and seabed soil.In this study,a numerical simulation is conducted using the coupled Eulerian–Lagrangian(CEL)large deformation numerical method to investigate the interaction between the track plate of the DSMV and the seabed soil under two working conditions:rotating condition and digging condition.First,a soil numerical model is established based on the elastoplastic mechanical characterization using the basic physical and mechanical properties of the seabed soil obtained by in situ sampling.Subsequently,the soil disturbance mechanism and the dynamic mechanical response of the track plate under rotating and digging conditions are obtained through the analysis of the sensitivity of the motion parameters,the grouser structure,the layered soil features and the soil heterogeneity.The results indicate that the above parameters remarkably influence the interaction between the DSMV and the seabed soil.Therefore,it is important to consider the rotating and digging motion of the DSMV in practical engineering to develop a detailed optimization design of the track plate.展开更多
The Voronoi grain-based breakable block model(VGBBM)based on the combined finite-discrete element method(FDEM)was proposed to explicitly characterize the failure mechanism and predict the deformation behavior of hard-...The Voronoi grain-based breakable block model(VGBBM)based on the combined finite-discrete element method(FDEM)was proposed to explicitly characterize the failure mechanism and predict the deformation behavior of hard-rock mine pillars.The influence of the microscopic parameters on the macroscopic mechanical behavior was investigated using laboratory-scale models.The field-scale pillar models(width-to-height,W/H=1,2 and 3)were calibrated based on the empirically predicted stress-strain curves of Creighton mine pillars.The results indicated that as the W/H ratios increased,the VGBBM effectively predicted the transition from strain-softening to pseudo-ductile behavior in pillars,and explicitly captured the separated rock slabs and the V-shaped damage zones on both sides of pillars and conjugate shear bands in core zones of pillars.The volumetric strain field revealed significant compressional deformation in core zones of pillars.While the peak strains of W/H=1 and 2 pillars were relatively consistent,there were significant differences in the strain energy storage and release mechanism.W/H was the primary factor influencing the deformation and strain energy in the pillar core.The friction coefficient of the structural plane was also an important factor affecting the pillar strength and the weakest discontinuity angle.The fracture surface was controlled by the discontinuity angle and the friction coefficient.This study demonstrated the capability of the VGBBM in predicting the strengths and deformation behavior of hard-rock pillars in deep mine design.展开更多
The collapse pressure is a key parameter when RTPs are applied in harsh deep-water environments.To investigate the collapse of RTPs,numerical simulations and hydrostatic pressure tests are conducted.For the numerical ...The collapse pressure is a key parameter when RTPs are applied in harsh deep-water environments.To investigate the collapse of RTPs,numerical simulations and hydrostatic pressure tests are conducted.For the numerical simulations,the eigenvalue analysis and Riks analysis are combined,in which the Hashin failure criterion and fracture energy stiffness degradation model are used to simulate the progressive failure of composites,and the“infinite”boundary conditions are applied to eliminate the boundary effects.As for the hydrostatic pressure tests,RTP specimens were placed in a hydrostatic chamber after filled with water.It has been observed that the cross-section of the middle part collapses when it reaches the maximum pressure.The collapse pressure obtained from the numerical simulations agrees well with that in the experiment.Meanwhile,the applicability of NASA SP-8007 formula on the collapse pressure prediction was also discussed.It has a relatively greater difference because of the ignorance of the progressive failure of composites.For the parametric study,it is found that RTPs have much higher first-ply-failure pressure when the winding angles are between 50°and 70°.Besides,the effect of debonding and initial ovality,and the contribution of the liner and coating are also discussed.展开更多
基金supported by the National Natural Science Foundation of China(51809169,51879159)Chang Jiang Scholars Program(T2014099)+2 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(2016-23/09)National Key Research and Development Program of China(2019YFB1704203,2019YFC0312400).
文摘Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.
基金financially supported by the Science and Technology Commission Foundation of Shanghai(Grant Nos.22DZ1208903,20DZ2251900)the National Natural Science Foundation of China(Grant No.51679134)。
文摘A new scour countermeasure using solidified slurry for offshore foundation has been proposed recently.Fluidized solidified slurry is pumped to seabed area around foundation for scour protection or pumped into the developed scour holes for scour repair as the fluidized material solidifies gradually.In the pumping operation and solidification,the engineering behaviors of solidified slurry require to be considered synthetically for the reliable application in scour repair and protection of ocean engineering such as the pumpability related flow value,flow diffusion behavior related rheological property,anti-scour performance related retention rate in solidification and bearing capacity related strength property after solidification.In this study,a series of laboratory tests are conducted to investigate the effects of mix proportion(initial water content and binder content)on the flow value,rheological properties,density,retention rate of solidified slurry and unconfined compressive strength(UCS).The results reveal that the flow value increases with the water content and decreases with the binder amount.All the solidified slurry exhibits Bingham plastic behavior when the shear rate is larger than 5 s^(-1).The Bingham model has been employed to fit the rheology test results,and empirical formulas for obtaining the density,yield stress and viscosity are established,providing scientific support for the numerical assessment of flow and diffusion of solidified slurry.Retention rate of solidified slurry decreases with the water flow velocity and flow value,which means the pumpability of solidified slurry is contrary to anti-scour performance.The unconfined compressive strength after solidification reduces as the water content increases and binder content decreases.A design and application procedure of solidified soil for scour repair and protection is also proposed for engineering reference.
文摘Ice-going ships play a crucial role in polar transportation and resource extraction.Different from the existing modeling approach which assumes that ships remain stationary,dynamic overset grid technology and DFBI(Dynamic Fluid-Body Interaction)method are employed in this paper to enable the free-running motion of the ship in modeling.A numerical model capable of simulating a ship navigating through pack ice area is proposed,which uses Computational Fluid Dynamics(CFD)method to solve the flow field and applies the Discrete Element Method(DEM)to simulate ship-ice and ice-ice interactions.Besides,the proposed high-precision method for generating pack ice area can be used in conjunction with the proposed numerical model.By comparing the numerical results with the available model test data and experimental observations,the effectiveness of the numerical model is validated,demonstrating its strong capability of predicting resistance and simulating ship navigation in pack ice,as well as its significant potential and applicability for further studies.
文摘To study the rolling motion of a ship in the presence of water on its deck,a linear-plus-quadratic damping term was incorporated into its equation of motion.Ship model tests indicates that the key dynamics of the physical system are preserved in the ship rolling equation with the linear-plus-quadratic type damping term.To take into account the presence of randomness in the excitation and the response,a new method was developed and a Melnikov criterion was obtained to provide an upper bound on the domain of the potential chaotic rolling motion(erratic rocking).Additionally,the Melnikov criterion proposed in this study was verified by the utilization of phase plane diagrams and Poincare maps.Furthermore,this research has made the initial endeavor to systematically modify the system parameters in the rolling equation of motion for ship stability analysis.
基金financially supported by the National Natural Science Foundation of China(Grant No.52088102)New Cornerstone Science Foundation through the XPLORER PRIZE.
文摘Floating breakwaters(FBs)are commonly employed for the protection of coastal installations.In this work,a convextype floating breakwater(FB)is proposed,and its hydrodynamic characteristics are studied through systematic laboratory experiments.Two different deck widths and two different mooring systems are set in the experiment.The transmission coefficients,reflection coefficients,motion responses and mooring forces of convex-type FBs are obtained in experiments.The influences of the deck width and mooring system on the hydrodynamic characteristics of the proposed FB are analyzed and compared.The experimental results show that the reflection coefficient and mooring force of the convex-type FB with a cross-mooring system are significantly larger than those of the convextype FB with a parallel-mooring system.A convex-type FB with a larger deck width has a higher reflection coefficient.The convex-type FBs with cross-and parallel-mooring systems have similar surge and heave motions,but the cross-mooring results in small roll motion.In addition,reliable prediction formulas for the transmission coefficient of convex-type FBs with different mooring systems have been developed,which are important for engineering design.
基金supported in part by the National Natural Science Foundation of China under Grant No.41706108the Science and Technology Commission of Shanghai Municipality Project 20dz1206600+2 种基金the Natural Science Foundation of Shanghai under Grant No.20ZR1424800the Oceanic Interdisciplinary Program of Shanghai Jiao Tong University under Grant No.SL2022ZD106the Natural Science Foundation of Chongqing under Grant No.cstc2021jcyj-msxmX0650.
文摘This study investigates the hydrodynamic process of a cylinder ejected from the water’s surface through high-speed camera experiments.Using digital image processing methods,the images obtained through experiments are processed and analyzed.Although the dynamics of rising buoyant cylinders have been thoroughly investigated,the pop-up height of the cylinders has not been extensively explored.Statistical analysis of the kinematic and dynamic data of cylinders is conducted.Research has shown that after the cylinder rises,it pops out of the water’s surface.Within the experimental range,the pop-up height of the cylinder is related to the release depth.Furthermore,the pop-up height and release depth of the cylinder vary linearly under vertical release conditions.Under horizontal release conditions,the relationship between pop-up height and release depth shows irregular changes mainly because of the unstable shedding of the wake vortex.
基金This research work was financially supported by the National Natural Science Foundation of China(Grant No.51972209).
文摘This paper investigated the use of magnesium phosphate cement (MPC) for solidifying sludge with different humic acid (HA) content (ranging from 0 to 4.5%) and explored the solidification mechanism. Fluidity, setting time, unconfined compressive strength (UCS), the strength formation mechanism, and the spontaneous imbibition process of solidified sludge (SS) were studied. The results indicate that MPC can be used as a low-alkalinity curing agent. As the HA content increases, fluidity and setting time also increase, while hydration temperature and strength decrease. Additionally, the failure mode of SS transitions from brittleness to ductility. The strength of SS is composed of the cementation strength provided by MPC hydration products, matric suction, osmotic suction, and the structural strength of the sludge. MPC reduces the structural strength caused by the shrinkage of pure sludge under the action of matric suction, but the incorporation of MPC significantly improved the strength when the sludge is eroded by water. X-ray diffraction (XRD) and scanning electron microscopy (SEM) show that the sludge and MPC can form a dense solid body, forming various hydration products, and synergistically improve the mechanical properties of the sludge.
基金supported by the National Natural Science Foundation of China(Grant Nos.U2244227,42277126 and 41977215).
文摘Air blasts induced by rock-ice avalanches are common natural phenomena known for their far-field destructive impact.In this study,remote sensing images,eyewitness videos and numerical modelling were comprehensively applied to analyze the initiation and propagation of the 2021 Chamoli avalancheinduced air blast.Our findings indicate that air blasts are observed from the avalanche source area to the Rishiganga valley,but nearly disappear in the Dhauliganga valley.The most intense air blast is concentrated on the left side of Ronti Gad valley,with maximum velocity and pressure estimated at over 70 m/s and 20 kPa,respectively.Such high pressure results in widespread tree breakage in the area.Based on the analysis of the Chamoli event,we further discussed the potential contribution of the avalanche flow regime,avalanche dynamics and geomorphology to the destructive potential of air blasts.Rapidly moved sliding mass can impart the air blast a high initial momentum,and this process will be exaggerated when the avalanche impacts valley walls at bends.However,when the rock-ice avalanche transforms into a debris-enriched flash flood,free water within the flowing mass can displace air,inhibiting the generation of air blasts.Our work offers new insights into the generation and propagation of rock-ice avalanche-induced air blasts,underscoring the importance of including this type of hazard during avalanche risk assessment in high-altitude glacial regions.
基金financially supported in part by the General Program of the National Natural Science Foundation of China (Grant No.12272221)the State Key Laboratory of Ocean Engineering (Shanghai Jiao Tong University)(Grant No. GKZD010087)。
文摘The integrated systems of unmanned surface vehicles(USVs) and remotely operated vehicles(ROVs) have been extensively applied in marine exploration and seabed coverage. However, the simultaneous navigation of USV-ROV systems is frequently limited by strong disturbances induced by waves or currents. This paper develops a novel rigidflexible coupling multibody dynamic model that incorporates disturbances of variable-length marine cables with geometrically nonlinear motion. A hybrid Lagrangian-Eulerian absolute nodal coordinate formulation(ANCF) element is developed to accurately model subsea cables which undergo significant overall motion, substantial deformation,and mass flow during the deployment of underwater equipment. Furthermore, the governing equations of the coupled USV-umbilical-ROV system are derived, considering wave-induced forces and current disturbances. A numerical solver based on the Newmark-beta method is proposed, along with an adaptive meshing technique near the release point. After validating three experimental cases, the cable disturbances at both the USV and ROV ends—caused by ocean currents, heave motion, and simultaneous navigation—are comprehensively compared and evaluated. Finally,it is demonstrated that a PD controller with disturbance compensation can enhance the simultaneous navigation performance of USV-ROV systems.
基金supported by the Shandong Natural Science Foundation Youth Project(No.ZR2023QE075)the Open Project Program of Shandong Marine Aerospace Equipment Technological Innovation Center(No.MAETIC202210)the Shandong Natural Science Foundation Project(No.ZR2022ME145)。
文摘The oscillating hydrofoil represents a promising technology for harvesting energy from tidal currents.While previous research has primarily focused on oscillating hydrofoils utilizing a fully activated control strategy,the industry predominantly employs a semi-activated control strategy in existing tidal current energy converters.It is essential to identify the differences in predicted energy-harvesting performance between these two controlling strategies through experimental modeling or numerical studies.Furthermore,the suitability of the fully activated control strategy in predicting the energy-harvesting capabilities of oscillating hydrofoils is evaluated.The 2D numerical models of hydrofoil based on fully activated and semi-activated control strategies have been developed and validated.The amplitudes of heaving and pitching movements for the fully activated hydrofoil are determined to match those of the semi-activated hydrofoil.The results show that the main difference between the two control strategies lies in the phase shift occurring between the pitching and heaving motions.This phase shift affects the lift force and its coordination with the heaving velocity,which in turn affects the power output.Notably,the maximum relative efficiency difference obtained between the fully activated and semi-activated control strategies can reach 191%.
基金the financial support received from the National Natural Science Foundation of China(Grant Nos.42072317 and 41727802)the Science and Technology Commission of Shanghai Municipality(Funding No.21DZ1204300).
文摘Comprehensive investigations have been conducted to study the structure and overconsolidation of upper Shanghai clays, i.e. Layers 2–6 clays, typically located at depths of 30–40 m. However, limited information is available on their anisotropy, and even less is known about the correlation between structure, overconsolidation, and anisotropy. In this study, the undrained anisotropy characteristics of shear strength and small-strain shear stiffness in upper Shanghai Layers 2–6 clays were thoroughly assessed using a series of K0-consolidated undrained triaxial compression (TC) and triaxial extension (TE) tests (K0 is the coefficient of lateral earth pressure at rest). The effective stress paths, shear strength, and small-strain shear stiffness from the undrained TC and TE tests demonstrate the anisotropic behaviors in upper Shanghai clays. Analyses of data from upper Shanghai clays and other clays worldwide indicate that the shear strength anisotropy ratio (Ks) converges at 0.8 as the overconsolidation ratio (OCR) and plasticity index (Ip) increase, while the small-strain shear stiffness anisotropy ratio (Re) converges at 1.0. The influence of OCR on Ks and Re is more pronounced than that of Ip and sensitivity (St). Nevertheless, no clear correlation between Ks and Re is observed in upper Shanghai clays.
基金supported by the National Natural Science Foundation of China(Grant Nos.12372309,12061131015).
文摘Veno-arterial extracorporeal membrane oxygenation(VA-ECMO)is a life support intervention for patients with refractory cardiogenic shock or severe cardiopulmonary failure.However,the choice of cannulation strategy remains contentious,partly due to insufficient understanding of hemodynamic characteristics associated with the site of arterial cannulation.In this study,a geometrical multiscale model was built to offer a mathematical tool for addressing the issue.The outflow cannula of ECMO was inserted into the ascending aorta in the case of central cannulation,whereas it was inserted into the right subclavian artery(RSA)or the left iliac artery(LIA)in the case of peripheral cannulation.Numerical simulations conducted on three patient-specific aortas demonstrated that the central cannulation outperformed the two types of peripheral cannulation in evenly delivering ECMO flow to branch arteries.Both the central and RSA cannulations could maintain an approximately normal hemodynamic state in the aortas,although the area of aortic walls exposed to abnormal wall shear stress(WSS)was considerably enlarged in comparison with the normal physiological condition.In contrast,the LIA cannulation not only led to insufficient delivery of ECMO flow to the right upper body(with ECMO flow fractions<0.5),but also induced marked flow disturbance in the aorta,causing about 40%of the abdominal aortic wall and over 65%of the resting aortic wall to suffer from high time-averaged WSS(>5 Pa)and low time-averaged WSS(<0.4 Pa),respectively.The LIA cannulation also resulted in significantly prolonged blood residence time(>40 s)in the ascending aorta,which,along with abnormal WSS,may considerably increase the risk of thrombosis.In summary,our numerical study elucidated the impact of arterial cannulation site in VA-ECMO intervention on aortic hemodynamics and ECMO flow distribution.The findings provide compensatory biomechanical information for traditional clinical studies and may serve as a theoretical reference for guiding the evaluation and selection of cannulation strategies in clinical practice.
基金financially supported by the Project of Ministry of Education and Finance of China(Grant Nos.200512 and 201335)the Project of the State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University(Grant No.GKZD010053-10)
文摘Knowledge-Based Engineering （KBE） is introduced into the ship structural design in this paper. From the implementation of KBE, the design solutions for both Rules Design Method （RDM） and Interpolation Design Method （IDM） are generated. The corresponding Finite Element （FE） models are generated. Topological design of the longitudinal structures is studied where the Gaussian Process （GP） is employed to build the surrogate model for FE analysis. Multi-objective optimization methods inspired by Pareto Front are used to reduce the design tank weight and outer surface area simultaneously. Additionally, an enhanced Level Set Method （LSM） which employs implicit algorithm is applied to the topological design of typical bracket plate which is used extensively in ship structures. Two different sets of boundary conditions are considered. The proposed methods show satisfactory efficiency and accuracy.
基金Supported by the Sustainable Energy Propulsion System Program of 211’s Engineering Foundation, Shanghai Jiaotong University.
文摘The phase change characteristic of the power source of an underwater glider propelled by the ocean＇s thermal energy is the key factor in glider attitude control. A numerical model has been established based on the enthalpy method to analyze the phase change heat transfer process under convective boundary conditions. Phase change is not an isothermal process, but one that occurs at a range of temperature. The total melting time of the material is very sensitive to the surrounding temperature. When the temperature of the surroundings decreases 8 degrees, the total melting time increases 1.8 times. But variations in surrounding temperature have little effect on the initial temperature of phase change, and the slope of the temperature time history curve remains the same. However, the temperature at which phase change is completed decreases significantly. Our research shows that the phase change process is also affected by container size, boundary conditions, and the power source＇s cross sectional area. Materials stored in 3 cylindrical containers with a diameter of 38ram needed the shortest phase change time. Our conclusions should be helpful in effective design of underwater glider power systems.
基金financially supported by the National Natural Science Foundation of China(Grant No.51239007)
文摘The size and shape of the effective test area are crucial to consider when short-crested waves are created by segmented wavemakers. The range of the effective test area of short-crested waves simulated by two-sided segmented wavemakers is analyzed in this paper. The experimental investigation on the wave field distribution of short-crested waves generated by two-sided segmented wavemakers is conducted by using an array of wave gauges. Wave spectra and directional spreading function are analyzed and the results show that when the main direction is at a certain angle with the normal line of wave generators, the wave field of 3D short-crested waves generated by two-sided segmented wavemakers has good spatial uniformity within the model test area. The effective test area can provide good wave environments for seakeeping model tests of various ocean engineering structures in the deep ocean engineering basin.
基金supported by the National Key R&D Program (No.2017YFC1405103)the Joint Funds of the National Natural Science Foundation of China (No.U170 6220)+1 种基金the National Natural Science Foundation of China (Nos.41901006, 41471005, and 41271016)the Natural Science Foundation of Shandong Province (No.ZR 2019BD005)。
文摘Research on the diffusion characteristics of swells contributes positively to wave energy forecasting, swell monitoring, and early warning. In this work, the South Indian Ocean westerly index(SIWI) and Indian Ocean swell diffusion effect index(IOSDEI) are defined on the basis of the 45-year(September 1957–August 2002) ERA-40 wave reanalysis data from the European Centre for Medium-Range Weather Forecasts(ECMWF) to analyze the impact of the South Indian Ocean westerlies on the propagation of swell acreage. The following results were obtained: 1) The South Indian Ocean swell mainly propagates from southwest to northeast. The swell also spreads to the Arabian Sea upon reaching low-latitude waters. The 2.0-meter contour of the swell can reach northward to Sri Lankan waters. 2) The size of the IOSDEI is determined by the SIWI strength. The IOSDEI requires approximately 2–3.5 days to fully respond to the SIWI. The correlations between SIWI and IOSDEI show obvious seasonal differences, with the highest correlations found in December–January–February(DJF) and the lowest correlations observed in June–July–August(JJA). 3) The SIWI and IOSDEI have a common period of approximately 1 week in JJA and DJF. The SIWI leads by approximately 2–3 days in this common period.
基金This work was supported by the open fund project of Shandong Provincial Key Laboratory of Ocean Engineer-ing,Ocean University of China(No.kloe201901)the Major International(Regional)Joint Research Project of the National Science Foundation of China(No.41520104008).
文摘Short-term forecasts of wave energy play a key role in the daily operation,maintenance planning,and electrical grid operation of power farms.In this study,we propose a short-term wave energy forecast scheme and use the North Indian Ocean(NIO)as a case study.Compared with the traditional forecast scheme,our proposed scheme considers more forecast elements.In addition to the traditional short-term forecast factors related to wave energy(wave power,significant wave height(SWH),wave period),our scheme emphasizes the forecast of a series of key factors that are closely related to the effectiveness of the energy output,capture efficiency,and conversion efficiency.These factors include the available rate,total storage,effective storage,co-occurrence of wave power-wave direction,co-occurrence of the SWH-wave period,and the wave energy at key points.In the regional nesting of nu-merical simulations of wave energy in the NIO,the selection of the southern boundary is found to have a significant impact on the simulation precision,especially during periods of strong swell processes of the South Indian Ocean(SIO)westerly.During tropical cyclone‘VARDAH’in the NIO,as compared with the simulation precision obtained with no expansion of the southern boundary(scheme-1),when the southern boundary is extended to the tropical SIO(scheme-2),the improvement in simulation precision is significant,with an obvious increase in the correlation coefficient and decrease in error.In addition,the improvement is much more significant when the southern boundary extends to the SIO westerly(scheme-3).In the case of strong swell processes generated by the SIO westerly,the improvement obtained by scheme-3 is even more significant.
基金supported by the Natural Science Foundation of Hainan Province(Grant No.520LH015)the Fundamental Research Funds for the Central Universities and the Major Projects of Strategic Emerging Industries in Shanghai(Grant No.BH3230001).
文摘To ensure the safe performance of deep-sea mining vehicles(DSMVs),it is necessary to study the mechanical characteristics of the interaction between the seabed soil and the track plate.The rotation and digging motions of the track plate are important links in the contact between the driving mechanism of the DSMV and seabed soil.In this study,a numerical simulation is conducted using the coupled Eulerian–Lagrangian(CEL)large deformation numerical method to investigate the interaction between the track plate of the DSMV and the seabed soil under two working conditions:rotating condition and digging condition.First,a soil numerical model is established based on the elastoplastic mechanical characterization using the basic physical and mechanical properties of the seabed soil obtained by in situ sampling.Subsequently,the soil disturbance mechanism and the dynamic mechanical response of the track plate under rotating and digging conditions are obtained through the analysis of the sensitivity of the motion parameters,the grouser structure,the layered soil features and the soil heterogeneity.The results indicate that the above parameters remarkably influence the interaction between the DSMV and the seabed soil.Therefore,it is important to consider the rotating and digging motion of the DSMV in practical engineering to develop a detailed optimization design of the track plate.
基金the National Natural Science Foundation of China(No.42377172)the National Key Research and Development Plan Project of China(No.2023YFC2907204).
文摘The Voronoi grain-based breakable block model(VGBBM)based on the combined finite-discrete element method(FDEM)was proposed to explicitly characterize the failure mechanism and predict the deformation behavior of hard-rock mine pillars.The influence of the microscopic parameters on the macroscopic mechanical behavior was investigated using laboratory-scale models.The field-scale pillar models(width-to-height,W/H=1,2 and 3)were calibrated based on the empirically predicted stress-strain curves of Creighton mine pillars.The results indicated that as the W/H ratios increased,the VGBBM effectively predicted the transition from strain-softening to pseudo-ductile behavior in pillars,and explicitly captured the separated rock slabs and the V-shaped damage zones on both sides of pillars and conjugate shear bands in core zones of pillars.The volumetric strain field revealed significant compressional deformation in core zones of pillars.While the peak strains of W/H=1 and 2 pillars were relatively consistent,there were significant differences in the strain energy storage and release mechanism.W/H was the primary factor influencing the deformation and strain energy in the pillar core.The friction coefficient of the structural plane was also an important factor affecting the pillar strength and the weakest discontinuity angle.The fracture surface was controlled by the discontinuity angle and the friction coefficient.This study demonstrated the capability of the VGBBM in predicting the strengths and deformation behavior of hard-rock pillars in deep mine design.
基金financially supported by National Natural Science Foundation of China(Grant Nos.52088102,51879249)Fundamental Research Funds for the Central Universities(Grant No.202261055)。
文摘The collapse pressure is a key parameter when RTPs are applied in harsh deep-water environments.To investigate the collapse of RTPs,numerical simulations and hydrostatic pressure tests are conducted.For the numerical simulations,the eigenvalue analysis and Riks analysis are combined,in which the Hashin failure criterion and fracture energy stiffness degradation model are used to simulate the progressive failure of composites,and the“infinite”boundary conditions are applied to eliminate the boundary effects.As for the hydrostatic pressure tests,RTP specimens were placed in a hydrostatic chamber after filled with water.It has been observed that the cross-section of the middle part collapses when it reaches the maximum pressure.The collapse pressure obtained from the numerical simulations agrees well with that in the experiment.Meanwhile,the applicability of NASA SP-8007 formula on the collapse pressure prediction was also discussed.It has a relatively greater difference because of the ignorance of the progressive failure of composites.For the parametric study,it is found that RTPs have much higher first-ply-failure pressure when the winding angles are between 50°and 70°.Besides,the effect of debonding and initial ovality,and the contribution of the liner and coating are also discussed.
