Application Progress of Computational Fluid Dynamic Techniques for Complex Viscous Flows in Ship and Ocean Engineering

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.

Experimental Study on Engineering Behavior of Solidified Soil for Scour Repair and Protection

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.

Structural Topology Design of Container Ship Based on Knowledge-Based Engineering and Level Set Method

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.

Phase change analysis of an underwater glider propelled by the ocean's thermal energy

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.

Laboratory Research on Effective Test Area of Short-Crested Waves Generated by Two-Sided Segmented Wavemakers

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.

Diffusion Characteristics of Swells in the North Indian Ocean

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.

Case Study of a Short-Term Wave Energy Forecasting Scheme:North Indian Ocean

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.

A Comparative Study on the Post-Buckling Behavior of Reinforced Thermoplastic Pipes(RTPs)Under External Pressure Considering Progressive Failure

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.

Prediction of Height and Period Joint Distributions for Stochastic Ocean Waves

This article proposes a new methodology to predict the wave height and period joint distributions by utilizing a transformed linear simulation method. The proposed transformed linear simulation method is based on a Hermite transformation model where the transformation is chosen to be a monotonic cubic polynomial, calibrated such that the first four moments of the transformed model match the moments of the true process. The proposed new approach is applied for calculating the wave height and period joint distributions of a sea state with the surface elevation data measured at an offshore site, and its accuracy and efficiency are favorably validated by using comparisons with the results from an empirical joint distribution model, from a linear simulation model and from a second-order nonlinear simulation model.

Lateral Dynamic Analysis and Motion Control of a Jet Trencher in Ocean Currents

In order to control the lateral motion of a jet trencher which is important for stable trenching operation,the oscillation characteristics of the jet trencher are researched. The jet trencher is simplified into a single degree of freedom model with restoring and damping force. The nonlinear mathematical model of the trencher laterally oscillating in ocean currents is established,and its approximate analytical solution is obtained.Results show that the analytical solution has small differences with numerical solution based on the fourth-order Runge-Kutta method and can effectively describe the underwater oscillation. A double-loop PID controller is designed to control the lateral motion displacement of the trencher to return to the center of the pipeline route which is effective and robust for the propulsion system.

Knowledge-based engineering approach for the support of ship structure design

The paper presents a knowledge-based engineering （KBE） approach for ship node components design. In the ship design process, many design tasks need design experiences to support. Howev- er, a ship design process is a complicated process with many simultaneously repetitive and time-con- suming activities. In this research, the method combines KBE with Tribon system＇s built-in devel- opment language tools of Vitesse, captures and applies design knowledge for achieving standard com- ponents intelligent design modeling. A case study and industry implementation illustrate the feasibili- ty of the proposed methodology. The KBE technique can provide not only proper references, sug- gests and supports but also knowledge integrated in the ship structure design. Especially, these rules related to the design can avoid lots of design mistakes. During the ship design stage, getting more precise and better designs will not only reduce the time of rework and wasting resources but also shorten the construction time_ imnrov~ clilnl;hz ~nA nrnf;t

Runout prediction of potential landslides based on the multi-source data collaboration analysis on historical cases

Long runout landslides involve a massive amount of energy and can be extremely hazardous owing to their long movement distance,high mobility and strong destructive power.Numerical methods have been widely used to predict the landslide runout but a fundamental problem remained is how to determine the reliable numerical parameters.This study proposes a framework to predict the runout of potential landslides through multi-source data collaboration and numerical analysis of historical landslide events.Specifically,for the historical landslide cases,the landslide-induced seismic signal,geophysical surveys,and possible in-situ drone/phone videos(multi-source data collaboration)can validate the numerical results in terms of landslide dynamics and deposit features and help calibrate the numerical(rheological)parameters.Subsequently,the calibrated numerical parameters can be used to numerically predict the runout of potential landslides in the region with a similar geological setting to the recorded events.Application of the runout prediction approach to the 2020 Jiashanying landslide in Guizhou,China gives reasonable results in comparison to the field observations.The numerical parameters are determined from the multi-source data collaboration analysis of a historical case in the region(2019 Shuicheng landslide).The proposed framework for landslide runout prediction can be of great utility for landslide risk assessment and disaster reduction in mountainous regions worldwide.

Experimental investigation of the inhibition of deep-sea mining sediment plumes by polyaluminum chloride

Deep-sea sediment disturbance may occur when collecting polymetallic nodules,resulting in the creation of plumes that could have a negative impact on the ecological environment.This study aims to investigate the potential solution of using polyaluminum chloride(PAC)in the water jet.The effects of PAC are examined through a self-designed simulation system for deep-sea polymetallic nodule collection and sediment samples from a potential deep-sea mining area.The experimental results showed that the optimal PAC dose was found to be 0.75 g/L.Compared with the test conditions without the addition of PAC,the presence of PAC leads to a reduction in volume,lower characteristic turbidity,smaller diffusion velocity,and shorter settling time of the plume.This indicates that PAC inhibits the entire development process of the plume.The addition of PAC leads to the flocculation of mm-sized particles,resulting in the formation of cm-sized flocs.The flocculation of particles decreases the rate of erosion on the seabed by around 30%.This reduction in erosion helps to decrease the formation of plumes.Additionally,when the size of suspended particles increases,it reduces the scale at which they diffuse.Furthermore,the settling velocity of flocs(around 10^(-2) m/s)is much higher that of compared to sediment particles(around 10^(-5) m/s),which effectively reduces the amount of time the plume remains in suspension.

Characteristics of a Quasi-Zero-Stiffness Isolator Using a Double-Curved Beam as a Negative Stiffness Mechanism

In this paper,a new quasi-zero-stiffness(QZS)nonlinear isolation system using a double-curved beam(DCB)as a negative stiffness structure is proposed,and its performance is investigated.The negative stiffness provided by the DCB to the isolator in the equilibrium position reduces the isolator’s overall dynamic stiffness.Static and dynamic characteristics of the system are investigated.The amplitude-frequency characteristics and force transmissibility equation of the system were derived via the harmonic balance method.The effects of damping ratio and excitation force amplitude on amplitude-frequency and force transmissibility curves are examined,and the isolation performance is compared with that of an equivalent linear isolator supporting the same mass with the same static deflection as nonlinear isolators.Furthermore,MATLAB numerical simulation software is used to perform dynamic time analysis of the nonlinear isolation system.The results indicate that the amplitude-frequency curves of the nonlinear isolation system exhibit bending,accompanied by discontinuous jumps in frequency.The appropriate increase in the damping ratio or reduction in the excitation amplitude benefits the vibration isolation performance of the nonlinear vibration isolation system.Compared with the equivalent linear isolation system,the QZS isolation system exhibits a better low-frequency vibration isolation performance,which provides a theoretical basis for the design of low-frequency nonlinear isolators.

Multiparameter Numerical Investigation of Two Types of Moving Interactions Between the Deep-Sea Mining Vehicle Track Plate and Seabed Soil:Digging and Rotating 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.

Numerical Simulations of Wave Impact Forces on the Open-Type Sea Access Road Using A Two-Phase SPH Model

A numerical study based on a two-dimensional two-phase SPH(Smoothed Particle Hydrodynamics)model to analyze the action of water waves on open-type sea access roads is presented.The study is a continuation of the analyses presented by Chen et al.(2022),in which the sea access roads are semi-immersed.In this new configuration,the sea access roads are placed above the still water level,therefore the presence of the air phase becomes a relevant issue in the determination of the wave forces acting on the structures.Indeed,the comparison of wave forces on the open-type sea access roads obtained from the single and two-phase SPH models with the experimental results shows that the latter are in much better agreement.So in the numerical simulations,a two-phaseδ-SPH model is adopted to investigate the dynamical problems.Based on the numerical results,the maximum horizontal and uplifting wave forces acting on the sea access roads are analyzed by considering different wave conditions and geometries of the structures.In particular,the presence of the girder is analyzed and the differences in the wave forces due to the air cushion effects which are created below the structure are highlighted.

Compression Properties and Energy Absorption of A Novel Double Curved Beam Negative Stiffness Honeycomb Structures

This paper presents the design of a novel honeycomb structure with a double curved beam.The purpose of this design is to achieve vibration isolation for the main engine of an offshore platform and reduce impact loads.An analytical formula for the force-displacement relationship of the honeycomb single-cell structure is presented based on the modal superposition method.This formula provides a theoretical basis for predicting the compression performance of honeycomb structures.The effects of structural geometric parameters,series and parallel connection methods on the mechanical and energy absorption properties are investigated through mathematical modeling and experimental methods.Furthermore,the study focuses on the vibration isolation and impact resistance performance of honeycomb panels.The results show that the designed honeycomb structure has good mechanical and energy absorption performance,and its energy absorption effect is related to the geometric parameters and series and parallel connection methods of the structure.The isolation efficiency of the honeycomb with 4 rows and 3 columns reaches 38%.The initial isolation frequency of the isolator is 11.7 Hz.

Investigation of stress-induced progressive failure of mine pillars using a Voronoi grain-based breakable block model

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.

Influence of Bow Design on Ice Resistance for Arctic LNG Carriers

In order to understand the influence of bow shape on ice resistance and provide guidelines for hull line design in the early design stage,an investigation of the impact of bow shape on ice resistance for the Arctic LNG carriers is carried out based on semi-empirical methods.Firstly,some typical semi-empirical formulas developed for ice resistance estimation of cargo carriers in different ice conditions are summarized.Then,formulas appropriate for ice resistance estimation of Arctic LNG carriers under different ice conditions are verified according to the result comparison between semi-empirical formulas and experimental tests.The comparison result indicates that the Lindqvist formula is appropriate for ice resistance estimation in level ice conditions,Zuev and Dobrodeev formula for ice resistance estimation in broken ice conditions,and Dobrodeev formula for ice resistance estimation in brash ice conditions.After that,the parameters considered in the selected formulas are summarized,and the influence of critical parameters on ice resistance is analyzed.Some parameters describing the ship's bow shape characteristic like ship breadth,waterline angle and stem angle greatly influence the ice resistance.Ice resistance increases with both the growth of ship breadth under all ice conditions and the growth of stem angle in level ice and broken ice conditions while ice resistance decreases with the development of waterline angle under all ice conditions.Finally,the optimization of the bow shape is discussed,and an optimized bow shape with both a large waterline angle and low stem angle is proposed.The optimized bow shape can decrease ice resistance by 9.9%in the level ice condition and reduce ice resistance by 11.3%in the brash ice condition.

A positive trend in the stability of global offshore wind energy

The recognition on the trend of wind energy stability is still extremely rare,although it is closely related to acquisition efficiency,grid connection,equipment lifetime,and costs of wind energy utilization.Using the 40-year(1979–2018)ERA-Interim data from the European Center for Medium-Range Weather Forecasts,this study presented the spatial-temporal distribution and climatic trend of the stability of global offshore wind energy as well as the abrupt phenomenon of wind energy stability in key regions over the past 40 years with the climatic analysis method and Mann-Kendall(M-K)test.The results show the following 5 points.(1)According to the coefficient of variation(C_(v))of the wind power density,there are six permanent stable zones of global offshore wind energy:the southeast and northeast trade wind zones in the Indian,Pacific and Atlantic oceans,the Southern Hemisphere westerly,and a semi-permanent stable zone(North Indian Ocean).(2)There are six lowvalue zones for both seasonal variability index(S_(v))and monthly variability index(M_(v))globally,with a similar spatial distribution as that of the six permanent stable zones.M_(v) and S_(v) in the Arabian Sea are the highest in the world.(3)After C_(v),M_(v) and S_(v) are comprehensively considered,the six permanent stable zones have an obvious advantage in the stability of wind energy over other sea areas,with C_(v) below 0.8,M_(v) within 1.0,and S_(v) within 0.7 all the year round.(4)The global stability of offshore wind energy shows a positive climatic trend for the past four decades.C_(v),M_(v) and S_(v) have not changed significantly or decreased in most of the global ocean during 1979 to2018.That is,wind energy is flat or more stable,while the monthly and seasonal variabilities tend to shrink/smooth,which is beneficial for wind energy utilization.(5)C_(v) in the low-latitude Pacific and M_(v) and S_(v) in both the North Indian Ocean and the low-latitude Pacific have an obvious abrupt phenomenon at the end of the20th century.