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.展开更多
A tunneled planing hull has unique hybrid hydrodynamic and aerodynamic characteristics due to the presence of a tunnel.In this paper,experimental and numerical investigations on hydrody namic analysis of a tunneled pl...A tunneled planing hull has unique hybrid hydrodynamic and aerodynamic characteristics due to the presence of a tunnel.In this paper,experimental and numerical investigations on hydrody namic analysis of a tunneled planing hull are carried out.The resistance tests of models with three dif ferent masses(127.4 kg,159.5 kg,202.9 kg)are conducted for the Froude number in the range of 0.761≤Fn≤1.925.The results of resistance measured by towing tank imply that the tunneled planing hull with a larger displacement has a superior resistance performance.The numerical simulation of Reynolds Average Navier Stokes(RANS)equations based on the finite volume method is performed to analyze the hull characteristics in calm water(M=159.5 kg)with two degrees of freedom(sinkage and trim).The numerical results are compared with the experimental data,which shows good agreement.Pressure distribution,wave profiles and lift forces obtained by SST k-ωand Realizable k-εturbulence models are compared and discussed.Finally,the local fluid flow of streamline around the hull can be divided into four regions due to the presence of a tunnel,which is different from the behaviors of the conventional planing monohull with prismatic form.展开更多
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 LN...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.展开更多
In the present work, the CFD-based method coupled with the dynamic overset grid technique is applied to investigate the hydrodynamic performance of the fully appended ONR tumblehome ship model under self-propulsion co...In the present work, the CFD-based method coupled with the dynamic overset grid technique is applied to investigate the hydrodynamic performance of the fully appended ONR tumblehome ship model under self-propulsion condition in head waves. All the computations are carried out by our in-house CFD solver naoe-FOAM-SJTU and the overset grid module is used to update the ship motions with rotating propellers while a self-developed 3D wave tank module is applied to generate desired wave environment. The ship model is advancing at its model point obtained with previous CFD results in calm water and the simulation is according to the bench- mark case from the Tokyo 2015 CFD Workshop in ship hydrodynamics. The predicted results, i.e. ship motions and instantaneous advancing speeds are presented and compared with the availa- ble experimental data. Propulsion coefficients, Kr and KQ, as well as detailed information of the flow field are also given to explain the hydrodynamic performance during the self-propulsion in waves. Good agreements are achieved which indicate that the present approach is applicable for the direct simulation of self-propulsion in waves.展开更多
In this paper, numerical simulations of FPSO ship motion coupled with LNG tank sloshing with low-filling ratios are conducted. The fully coupled problem is addressed with our own unsteady RANS solver: naoe-FOAM-SJTU d...In this paper, numerical simulations of FPSO ship motion coupled with LNG tank sloshing with low-filling ratios are conducted. The fully coupled problem is addressed with our own unsteady RANS solver: naoe-FOAM-SJTU developed based on the open source tool librar- ies of OpenFOAM. The internal tank sloshing and external wave flow are solved simultaneously. The FPSO model includes 2 LNG tanks. For the ship 3-DOFs are released in the regular beam waves. The filling ratios of the 2 tanks are 20% - 20%, lower than the external free surface. This kind of low-filling condition reduces ship roll motion significantly, and produces complex free surface shapes in tanks. 4 different incident wave frequencies are considered in the simulation in comparison with the existing experimental data. The comparison shows that the numerical re- sults are in good agreement with the experimental data, proving the reliability of the proposed method. The filling conditions with large wave amplitudes are studied further, and due to the coupling effect, violent sloshing occurs in tanks and impulsive pressure forms on bulkhead.展开更多
A numerical triangulation and transformation into the time domain of a Kirchhoff approximation(KA)method is proposed for the modeling of bistatic scattering from an underwater non-penetrable target.The time domain sol...A numerical triangulation and transformation into the time domain of a Kirchhoff approximation(KA)method is proposed for the modeling of bistatic scattering from an underwater non-penetrable target.The time domain solution in this approximation can be split up into two parts:the solution of reflected field,contributing around the specular direction,and the solution of shadow radiation,contributing around the forward direction.An average solution in the time domain satisfying the reciprocity principle is presented.The solution is expressed in terms of non-singular functions.The proposed method is validated against a normal mode method for bistatic scattering from a rigid sphere.Moreover,the reflected and shadow highlights on the surface of the sphere are shown to verify the integration surface of the reflected field and shadow radiation.It is also tested against a finite element method and an experiment involving a scaled Benchmark Target Strength Simulation Submarine model.The time-angle bistatic spectra for the model are evaluated by the direct and transformed average solutions of KA,and the former accelerates its speed of calculation.The results are good,and show that this method can be used to predict the bistatic scattered field of a non-penetrable target.展开更多
In this paper, 2 detached-eddy simulation (DES) approaches, namely SST-DES and SST-DDES are implemented, integrated in to the naoe-FOAM-SJTU solver which is developed based on the open source platform OpenFOAM. Flow p...In this paper, 2 detached-eddy simulation (DES) approaches, namely SST-DES and SST-DDES are implemented, integrated in to the naoe-FOAM-SJTU solver which is developed based on the open source platform OpenFOAM. Flow past 2 cylinders in tandem arrangement is selected as the benchmark case for the validation of the SST-DES and SST-DDES approaches. The experiment was previously conducted in 2 different wind tunnels at the NASA Langley Re- search Center. Time-averaged flow fields and some quantities of computational results are com- pared with experiments. In addition, the 3D instantaneous flow structures are also given and discussed. It is shown that the current implementation of SST-DES and SST-DDES is able to re- solve some characteristics for massively separated complex turbulent flows.展开更多
Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a l...Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a liquid sloshing tank potentially affect the structural safety of cargo tanks or vessels. The challenges of numerical study on FSI problems involve not only multidisciplinary features, but also accurate description of non-linear free surface. A fully Lagrangian particle-based method , the moving particle semi-implicit and finite element coupled method ( MPS-FEM), is developed to numerically study the FSI problems. Taking into account the advantage of the Lagrangian method for large deformations of both fluid and solid bounda- ties, the MPS method is used to simulate the fluid field while the finite element method(FEM) to calculate the structure field. Besides, the partitioning strategy is employed to couple the MPS and FEM modules. To validate accuracy of the proposed algorithm, a benchmark case is numer- ically investigated. Both the patterns of free surface and the deflections of the elastic structures are in good agreement with the experimental data. Then, the present FSI solver is applied to the comparative study of the mitigating effects of rigid baffles and elastic baffles on the sloshing motions and impact loads.展开更多
文摘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.
文摘A tunneled planing hull has unique hybrid hydrodynamic and aerodynamic characteristics due to the presence of a tunnel.In this paper,experimental and numerical investigations on hydrody namic analysis of a tunneled planing hull are carried out.The resistance tests of models with three dif ferent masses(127.4 kg,159.5 kg,202.9 kg)are conducted for the Froude number in the range of 0.761≤Fn≤1.925.The results of resistance measured by towing tank imply that the tunneled planing hull with a larger displacement has a superior resistance performance.The numerical simulation of Reynolds Average Navier Stokes(RANS)equations based on the finite volume method is performed to analyze the hull characteristics in calm water(M=159.5 kg)with two degrees of freedom(sinkage and trim).The numerical results are compared with the experimental data,which shows good agreement.Pressure distribution,wave profiles and lift forces obtained by SST k-ωand Realizable k-εturbulence models are compared and discussed.Finally,the local fluid flow of streamline around the hull can be divided into four regions due to the presence of a tunnel,which is different from the behaviors of the conventional planing monohull with prismatic form.
文摘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.
文摘In the present work, the CFD-based method coupled with the dynamic overset grid technique is applied to investigate the hydrodynamic performance of the fully appended ONR tumblehome ship model under self-propulsion condition in head waves. All the computations are carried out by our in-house CFD solver naoe-FOAM-SJTU and the overset grid module is used to update the ship motions with rotating propellers while a self-developed 3D wave tank module is applied to generate desired wave environment. The ship model is advancing at its model point obtained with previous CFD results in calm water and the simulation is according to the bench- mark case from the Tokyo 2015 CFD Workshop in ship hydrodynamics. The predicted results, i.e. ship motions and instantaneous advancing speeds are presented and compared with the availa- ble experimental data. Propulsion coefficients, Kr and KQ, as well as detailed information of the flow field are also given to explain the hydrodynamic performance during the self-propulsion in waves. Good agreements are achieved which indicate that the present approach is applicable for the direct simulation of self-propulsion in waves.
文摘In this paper, numerical simulations of FPSO ship motion coupled with LNG tank sloshing with low-filling ratios are conducted. The fully coupled problem is addressed with our own unsteady RANS solver: naoe-FOAM-SJTU developed based on the open source tool librar- ies of OpenFOAM. The internal tank sloshing and external wave flow are solved simultaneously. The FPSO model includes 2 LNG tanks. For the ship 3-DOFs are released in the regular beam waves. The filling ratios of the 2 tanks are 20% - 20%, lower than the external free surface. This kind of low-filling condition reduces ship roll motion significantly, and produces complex free surface shapes in tanks. 4 different incident wave frequencies are considered in the simulation in comparison with the existing experimental data. The comparison shows that the numerical re- sults are in good agreement with the experimental data, proving the reliability of the proposed method. The filling conditions with large wave amplitudes are studied further, and due to the coupling effect, violent sloshing occurs in tanks and impulsive pressure forms on bulkhead.
基金the Oceanic Interdisciplinary Program of Shanghai Jiao Tong University(project number SL2021PT108)。
文摘A numerical triangulation and transformation into the time domain of a Kirchhoff approximation(KA)method is proposed for the modeling of bistatic scattering from an underwater non-penetrable target.The time domain solution in this approximation can be split up into two parts:the solution of reflected field,contributing around the specular direction,and the solution of shadow radiation,contributing around the forward direction.An average solution in the time domain satisfying the reciprocity principle is presented.The solution is expressed in terms of non-singular functions.The proposed method is validated against a normal mode method for bistatic scattering from a rigid sphere.Moreover,the reflected and shadow highlights on the surface of the sphere are shown to verify the integration surface of the reflected field and shadow radiation.It is also tested against a finite element method and an experiment involving a scaled Benchmark Target Strength Simulation Submarine model.The time-angle bistatic spectra for the model are evaluated by the direct and transformed average solutions of KA,and the former accelerates its speed of calculation.The results are good,and show that this method can be used to predict the bistatic scattered field of a non-penetrable target.
文摘In this paper, 2 detached-eddy simulation (DES) approaches, namely SST-DES and SST-DDES are implemented, integrated in to the naoe-FOAM-SJTU solver which is developed based on the open source platform OpenFOAM. Flow past 2 cylinders in tandem arrangement is selected as the benchmark case for the validation of the SST-DES and SST-DDES approaches. The experiment was previously conducted in 2 different wind tunnels at the NASA Langley Re- search Center. Time-averaged flow fields and some quantities of computational results are com- pared with experiments. In addition, the 3D instantaneous flow structures are also given and discussed. It is shown that the current implementation of SST-DES and SST-DDES is able to re- solve some characteristics for massively separated complex turbulent flows.
文摘Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a liquid sloshing tank potentially affect the structural safety of cargo tanks or vessels. The challenges of numerical study on FSI problems involve not only multidisciplinary features, but also accurate description of non-linear free surface. A fully Lagrangian particle-based method , the moving particle semi-implicit and finite element coupled method ( MPS-FEM), is developed to numerically study the FSI problems. Taking into account the advantage of the Lagrangian method for large deformations of both fluid and solid bounda- ties, the MPS method is used to simulate the fluid field while the finite element method(FEM) to calculate the structure field. Besides, the partitioning strategy is employed to couple the MPS and FEM modules. To validate accuracy of the proposed algorithm, a benchmark case is numer- ically investigated. Both the patterns of free surface and the deflections of the elastic structures are in good agreement with the experimental data. Then, the present FSI solver is applied to the comparative study of the mitigating effects of rigid baffles and elastic baffles on the sloshing motions and impact loads.
