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.

Steady State Analysis of Towed Marine Cables

Efficient numerical schemes were presented for the steady state solutions of towed marine cables. For most of towed systems,the steady state problem can be resolved into two-point boundary-value problem,or initial value problem in some special cases where the initial values are available directly.A new technique was proposed and attempted to solve the two-point boundary-value problem rather than the conventional shooting method due to its algorithm complexity and low efficiency.First,the boundary conditions are transformed into a set of nonlinear governing equations about the initial values,then bisection method is employed to solve these nonlinear equations with the aid of 4th order Runge-Kutta method.In common sense,non-uniform (sheared) current is assumed,which varies in magnitude and direction with depth.The schemes are validated through the DE Zoysa's example,then several numerical examples are also presented to illustrate the numerical schemes.

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.

An Improved Time Domain Approach for Analysis of Floating Bridges Based on Dynamic Finite Element Method and State-Space Model

The floating bridge bears the dead weight and live load with buoyancy,and has wide application prospect in deep-water transportation infrastructure.The structural analysis of floating bridge is challenging due to the complicated fluid-solid coupling effects of wind and wave.In this research,a novel time domain approach combining dynamic finite element method and state-space model(SSM)is established for the refined analysis of floating bridges.The dynamic coupled effects induced by wave excitation load,radiation load and buffeting load are carefully simulated.High-precision fitted SSMs for pontoons are established to enhance the calculation efficiency of hydrodynamic radiation forces in time domain.The dispersion relation is also introduced in the analysis model to appropriately consider the phase differences of wave loads on pontoons.The proposed approach is then employed to simulate the dynamic responses of a scaled floating bridge model which has been tested under real wind and wave loads in laboratory.The numerical results are found to agree well with the test data regarding the structural responses of floating bridge under the considered environmental conditions.The proposed time domain approach is considered to be accurate and effective in simulating the structural behaviors of floating bridge under typical environmental conditions.

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.

Determination of Scheduled Travel Time for a Fixed Transit Route Based on Multistate Model

In order to provide the guideline for bus drivers to adjust speed to minimize scheduled deviation,the method for setting bus scheduled travel time is proposed. Firstly,multistate model is introduced to fit historical travel time data and identify different service states. Based on the calibrated travel time distribution parameters,an optimization model is proposed,followed by a Monte Carlo( MC) simulation based genetic algorithm( GA)procedure to obtain the optimal scheduled time. A case study from a fixed bus route from Shenzhen is used to demonstrate the model applicability. The sensitivity analysis is conducted to study the effects of parameters setting on optimal slack time for each segment. The results show that multistate model fits travel time under peak hours better than Lognormal distribution,and the length of scheduled travel time basically reflects travel time reliability.

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

Steady-State Configuration and Tension Calculations of Marine Cables Under Complex Currents via Separated Particle Swarm Optimization?

Under complex currents, the motion governing equations of marine cables are complex and nonlinear, and the calculations of cable configuration and tension become difficult compared with those under the uniform or simple currents. To obtain the numerical results, the usual Newton-Raphson iteration is often adopted, but its stability depends on the initial guessed solution to the governing equations. To improve the stability of numerical calculation, this paper proposed separated the particle swarm optimization, in which the variables are separated into several groups, and the dimension of search space is reduced to facilitate the particle swarm optimization. Via the separated particle swarm optimization, these governing nonlinear equations can be solved successfully with any initial solution, and the process of numerical calculation is very stable. For the calculations of cable configuration and tension of marine cables under complex currents, the proposed separated swarm particle optimization is more effective than the other particle swarm optimizations.

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.

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.

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.

Experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering boundary effect

The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.

Experimental Investigation on Vertical Hydraulic Transport of Ores in Deepsea Mining

Deepsea mining has been proposed since the 1960s to alleviate the lack of resources on land.Vertical hydraulic transport of collected ores from the seabed to the sea surface is considered the most promising method for industrial applications.In the present study,an indoor model test of the vertical hydraulic transport of particles was conducted.A noncontact optical method has been proposed to measure the local characteristics of the particles inside a vertical pipe,including the local concentration and particle velocity.The hydraulic gradient of ore transport was evaluated with various particle size distributions,particle densities,feeding concentrations and mixture flow velocities.During transport,the local concentration is larger than the feeding concentration,whereas the particle velocity is less than the mixture velocity.The qualitative effects of the local concentration and local fluid velocity on the particle velocity and slip velocity were investigated.The local fluid velocity contributes significantly to particle velocity and slip velocity,whereas the effect of the local concentration is marginal.A higher feeding concentration and mixture flow velocity result in an increased hydraulic gradient.The effect of the particle size gradation is slight,whereas the particle density plays a crucial role in the transport.

Optimal Structural Design of the Midship of a VLCC Based on the Strategy Integrating SVM and GA

In this paper a hybrid process of modeling and optimization, which integrates a support vector machine （SVM） and genetic algorithm （GA）, was introduced to reduce the high time cost in structural optimization of ships. SVM, which is rooted in statistical learning theory and an approximate implementation of the method of structural risk minimization, can provide a good generalization performance in metamodeling the input-output relationship of real problems and consequently cuts down on high time cost in the analysis of real problems, such as FEM analysis. The GA, as a powerful optimization technique, possesses remarkable advantages for the problems that can hardly be optimized with common gradient-based optimization methods, which makes it suitable for optimizing models built by SVM. Based on the SVM-GA strategy, optimization of structural scantlings in the midship of a very large crude carrier （VLCC） ship was carried out according to the direct strength assessment method in common structural rules （CSR）, which eventually demonstrates the high efficiency of SVM-GA in optimizing the ship structural scantlings under heavy computational complexity. The time cost of this optimization with SVM-GA has been sharply reduced, many more loops have been processed within a small amount of time and the design has been improved remarkably.

Hydrodynamic Interactions of Three Barges in Close Proximity in A Floatover Installation

The hydrodynamics of side-by-side barges are much more complex than those of a single barge in waves because of wave shielding, viscous effects and water resonance in the gap. In the present study, hydrodynamic coefficients in the frequency domain were calculated for both the system of multiple bodies and the isolated body using both low-order and higher-order boundary-element methods with different element numbers. In these calculations, the damping-lid method was used to modify the free-surface boundary conditions in the gap and to make the hydrodynamic results more reasonable. Then far-field, mid-field and near-field methods were used to calculate wave-drift forces for both the multi-body system and the isolated body. The results show that the higher-order method has faster convergence speed than the low-order method for the multi-body case. Comparison of different methods of computing drift force showed that mid-field and far-field methods have better convergence than the near-field method. In addition, corresponding model tests were performed in the Deepwater Offshore Basin at Shanghai Jiao Tong University. Comparison between numerical and experimental results showed good agreement.

CFD Investigations of Ship Maneuvering in Waves Using naoe-FOAM-SJTU Solver

Ship maneuvering in waves includes the performance of ship resistance, seakeeping, propulsion, and maneuverability. It is a complex hydrodynamic problem with the interaction of many factors. With the purpose of directly predicting the behavior of ship maneuvering in waves, a CFD solver named naoe-FOAM-SJTU is developed by the Computational Marine Hydrodynamics Lab(CMHL) in Shanghai Jiao Tong University. The solver is based on open source platform OpenFOAM and has introduced dynamic overset grid technology to handle complex ship hull-propeller-rudder motion system. Maneuvering control module based on feedback control mechanism is also developed to accurately simulate corresponding motion behavior of free running ship maneuver. Inlet boundary wavemaker and relaxation zone technique is used to generate desired waves. Based on the developed modules, unsteady Reynolds-averaged Navier-Stokes(RANS) computations are carried out for several validation cases of free running ship maneuver in waves including zigzag, turning circle, and course keeping maneuvers. The simulation results are compared with available benchmark data. Ship motions, trajectories, and other maneuvering parameters are consistent with available experimental data, which indicate that the present solver can be suitable and reliable in predicting the performance of ship maneuvering in waves. Flow visualizations, such as free surface elevation, wake flow, vortical structures, are presented to explain the hydrodynamic performance of ship maneuvering in waves. Large flow separation can be observed around propellers and rudders. It is concluded that RANS approach is not accurate enough for predicting ship maneuvering in waves with large flow separations and detached eddy simulation(DES) or large eddy simulation(LES) computations are required to improve the prediction accuracy.