A numerical method based on solutions of Euler/Navier-Stokes （N-S） equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advectio...A numerical method based on solutions of Euler/Navier-Stokes （N-S） equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advection upwind splitting method （AUSM） scheme, upwind AUSM/van Leer scheme, AUSM＋ scheme and AUSMDV scheme are implemented for spatial discretization, and van Albada limiter is also applied. For temporal discretization, both explicit Runge-Kutta method and implicit lower-upper symmetric Gauss-Seidel （LU-SGS） method are attempted. Simultaneously, overset grid technique is adopted. In detail, hole-map method is utilized to identify intergrid boundary points （IGBPs）. Furthermore, aimed at identification issue of donor elements, inverse-map method is implemented. Eventually, blade surface pressure distributions derived from numerical simulation are validated compared with experimental data, showing that all the schemes mentioned above have the capability to predict the rotor flow field accurately. At the same time, vorticity contours are illustrated for analysis, and other characteristics are also analyzed.展开更多
With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor co...With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.展开更多
A method of unstructured dynamic overset grids is developed for the numerical simulation of helicopter unsteady rotorairframe aerodynamic interaction.For the effective treatment of the relative motion between the roto...A method of unstructured dynamic overset grids is developed for the numerical simulation of helicopter unsteady rotorairframe aerodynamic interaction.For the effective treatment of the relative motion between the rotor and the airframe,the domain of flowfield is divided into two overset subzones,namely,a rotational subzone containing the blades and a stationary subzone containing the airframe.The overset part of two subzones is used to convect the flow variables of the two zones.The Taylor series expansion is used to obtain a second-order spatial accuracy,and dual-time stepping is adopted to improve the solution accuracy.Mesh deformation from the blade motion in forward flight is treated by using a spring analogy.Validation is made by numerically simulating the flows around a wind tunnel configuration and comparing the predicted time-averaged and instantaneous inflow and airframe surface pressure distributions with the experimental data.It shows that the present method is efficient and robust for the prediction of complicated unsteady rotor-airframe aerodynamic interaction phenomena.展开更多
In the present study, the performance of the NTNU Blind Test 1 wind turbine is analyzed in the computational fluid dynamics (CFD) simulations by using the CFD code FANS with structured overset grids. First, the numeri...In the present study, the performance of the NTNU Blind Test 1 wind turbine is analyzed in the computational fluid dynamics (CFD) simulations by using the CFD code FANS with structured overset grids. First, the numerical methods including the governing equations, the turbulence closure model, and the flow solver are introduced. In addition, the NTNU BT1 wind tunnel experiment is described. Then, structured overset grid blocks are generated in the computational domain with fully resolved wind turbine geometry, including the blades, hub, nacelle, and tower. Afterward, unsteady Reynolds averaged Navier-Stokes (RANS) simulations with the two-layer k - ε turbulence model are performed with an inlet velocity of 10 m/s and a tip-speed ratio (TSR) of 6. The overset-grid capability of FANS is leveraged to handle the rotation of the rotor. Finally, simulations are performed for a range of TSRs and a comparison is made among the present CFD results, other numerical results obtained from representative methods, and the experimental data. It is observed that the CFD-predicted thrust coefficients match the experimental measurement at low TSRs while under-predicting the values at high TSRs, and potential reasons for this deviation are discussed.展开更多
A cell-centred overset unstructured grids approach is developed.In this approach,the intergrid boundary is initially established based on the wall distance from the cell centre,and is then optimized.To accelerate the ...A cell-centred overset unstructured grids approach is developed.In this approach,the intergrid boundary is initially established based on the wall distance from the cell centre,and is then optimized.To accelerate the intergrid-boundary definition much more,a neighbor-toneighbor donor search algorithm based on advancing-front method is modified with the help of minimum cuboid boxes.To simplify the communications between different grid cell types and to obtain second-order spatial accuracy,a new interpolation method is constructed based on linear reconstruction,which employs only one layer of fringe cells along the intergrid boundary.For unsteady flows with relative motion,the intergrid boundary can be redefined fast and automatically.Several numerical results show that the present dynamic overset unstructured grids approach is accurate and reliable.展开更多
We present high order accurate numerical methods for the wave equation that combines efficient Hermite methods with geometrically flexible discontinuous Galerkin methods by using overset grids.Near boundaries we use t...We present high order accurate numerical methods for the wave equation that combines efficient Hermite methods with geometrically flexible discontinuous Galerkin methods by using overset grids.Near boundaries we use thin boundary fitted curvilinear grids and in the volume we use Cartesian grids so that the computational complexity of the solvers approaches a structured Cartesian Hermite method.Unlike many other overset methods we do not need to add artificial dissipation but we find that the built-in dissipation of the Hermite and discontinuous Galerkin methods is sufficient to maintain the stability.By numerical experiments we demonstrate the stability,accuracy,efficiency,and the applicability of the methods to forward and inverse problems.展开更多
In this article, the overset grid assembly method is improved to efficiently solve several critical problems that occur when applying overset grids to the complicated geometries and moving body. First, instead of usin...In this article, the overset grid assembly method is improved to efficiently solve several critical problems that occur when applying overset grids to the complicated geometries and moving body. First, instead of using the two-step searching (i.e. cut and paste), a single-step searching method based on the grid cell size is proposed to modify holes and optimize the grid overlapping automatically. Second, discrepancies between the wall surface representations, where the grids overlap, are handled by introducing the wall's normal directions to the hole-map and projecting the interpolated points on the recipient mesh into the donor mesh. Finally, the dynamic overset method is modified to address the complex moving body problem. At every time step, the initial hole surface of the previous time step is dynamically adjusted to accomplish hole cutting and avoid the time consuming hole-map procedure. Numerical experiments show that the enhanced overset grid assembly method obtains satisfactory results.展开更多
In the present study,an efficient overset grid method by means of parallel implicit hole-cutting is proposed for the sake of simulating unsteady flows in aerospace engineering involving multiple bodies in relative mov...In the present study,an efficient overset grid method by means of parallel implicit hole-cutting is proposed for the sake of simulating unsteady flows in aerospace engineering involving multiple bodies in relative movement.In view of the degraded computational efficiency and robustness for conventional overset grid assembly,several innovative techniques are developed within the overset grid assembly process,viz.,a bookkeeping alternative digital tree method to speed up the donor-cell searching,a fast parallel advancing front algorithm to accelerate the wall-distance calculation and a message-passing strategy with efficient information communication and lower storage expenditure within distributed computational architecture.The contribution of the developed techniques is evidenced by comparison with the existing alternative ways in terms of computing efficiency.Subsequently,the overset grid method is embedded into an inhouse programed URANS solver to examine its capability in predicting the flow field of complex applications such as helicopter,store separation and component deploying.Results show that the developed overset grid methodology is,in practice,able to resolve the aerodynamic characteristics of complex aerospace engineering with a high-fidelity flow topology and accuracy.展开更多
The hydrodynamic performance of high-speed planing hulls has gained considerable interest,with recent advancements in computational fluid dynamics and hull design techniques enhancing the understanding of planing hull...The hydrodynamic performance of high-speed planing hulls has gained considerable interest,with recent advancements in computational fluid dynamics and hull design techniques enhancing the understanding of planing hull hydrodynamics.In this study,we conducted a numerical investigation using the Reynolds-averaged Navier-Stokes approach with overset grids to capture large motions at high speeds.This study aims to improve the hydrodynamic performances of planing hulls,specifically focusing on total resistance,trim,and sinkage.The initial Fridsma hull with a deadrise angle of 20°has been used for validation,demonstrating good agreement with measurements at different Froude numbers.Subsequently,new configurations based on the Fridsma hull have been designed by varying the deadrise angle,number of chines,and transverse steps.Our findings reveal a correlation between the deadrise angle,the number of chines,and the Froude number.As the deadrise angle increases,total resistance also increases.Additionally,a single chine yields superior results at higher Froude numbers,while multiple chines offer advantages at lower values.The introduction of transverse steps consistently increases total resistance,highlighting their role in improving planing hull performance.This research not only offers valuable insights into planing hull design but also leverages state-of-the-art numerical methods to advance the understanding of hydrodynamic behaviors at high ship speeds.展开更多
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.展开更多
Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of th...Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.展开更多
A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or...A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or non-viscous, field equations with nonlinear boundary conditions applied to the free surface and fluid-solid interaction interfaces. An Arbitrary-LagrangianEulerian (ALE) mesh system is used to construct the numerical model. A multi-block numerical scheme of study is adopted allowing for the relative motion between moving overset grids, which are independent of one another. This provides a convenient method to overcome the difficulties in matching fluid meshes with large solid motions. Nonlinear numerical equations describing nonlinear fluid-solid interaction dynamics are derived through a numerical discretization scheme of study. A coupling iteration process is used to solve these numerical equations. Numerical examples are presented to demonstrate applications of the model developed.展开更多
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 direct...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.展开更多
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.展开更多
Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical d...Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors.WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction,and its advantage is to reduce dissipation in smooth regions of flow while preserving high-resolution around discontinuities.Here WENO-K scheme is adopted to reconstruct left and right state values within the Roe Riemann solver updating the inviscid fluxes on a structured dynamic overset grid.To minimize the accuracy loss for high-order reconstruction on artificial boundaries of overset grid,a method of multilayer fringes is proposed to carry out interpolation between background grid and blade grid.Massively parallel computing considering automatic load balance on averagely partitioned overset grid is developed to reduce the wall-clock time of an unsteady simulation.Numerical results for Caradonna-Tung(C-T)rotor in hover at the conditions of subsonic and transonic tip Mach numbers show that the thrust coefficient error for the result of WENO-K scheme is no more than 3%.Compared with WENO-JS scheme,WENO-K scheme achieves about 40%improvement on accuracy of predicting rotor thrust with only 4.1%extra computational cost.More importantly,WENO-K scheme can capture more sophisticated unsteady flow structures and resolve tip vortices to a larger wake age with an increment of about 270°compared to WENO-JS scheme.展开更多
Parametric rolling is one of five types of the ship stability failure modes as proposed by IMO. The periodic change of the metacentric height is often considered as the internal cause of this phenomenon. Parametric ro...Parametric rolling is one of five types of the ship stability failure modes as proposed by IMO. The periodic change of the metacentric height is often considered as the internal cause of this phenomenon. Parametric rolling is a complex nonlinear hydrodynamic problem, often accompanied by large amplitude vertical motions of ships. In recent years,the Reynolds-averaged Navier–Stokes(RANS) equation simulations for viscous flows have made great progress in the field of ship seakeeping. In this paper, the parametric rolling for the C11 containership in regular waves is studied both experimentally and numerically. In the experiments, parametric rolling amplitudes at different drafts, forward speeds and wave steepnesses are analyzed. The differences in the steady amplitudes of parametric rolling are observed for two drafts. The effect of the incident wave steepness(or wave amplitude) is also studied, and this supports previous results obtained on limits of the stability for parametric rolling. In numerical simulations, the ship motions of parametric rolling are analyzed by use of the potential-flow and viscous-flow methods. In the viscousflow method, the Reynolds-averaged Navier–Stokes equations are solved using the overset grid method. The numerical accuracies of the two methods at different wave steepnesses are also discussed.展开更多
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.展开更多
This paper presents 3D （three-dimensional） CFD （computational fluid dynamic） simulation to analyse the FIR （flow-induced response） especially the yaw motion of a buoyancy can. The numerical cases are conducted w...This paper presents 3D （three-dimensional） CFD （computational fluid dynamic） simulation to analyse the FIR （flow-induced response） especially the yaw motion of a buoyancy can. The numerical cases are conducted with a buoyancy can under different reduced velocities utilizing our in-house code naoe-FOAM-SJTU, a solver based on the open source toolkit OpenFOAM. SST-DDES （shear-stress transport-delayed detached-eddy simulation） model is applied to handle the flowseparation and overset grid method is utilized to solve a large amplitude 6-DOF （6 degrees of freedom） motions. Free decay test and VIM （vortex-induced motion） test are built numerically. In VIM cases, the responses of trajectory, amplitude, frequency are calculated in a series of reduced velocities. With the increase of reduced velocity, yaw frequency is increased, which is similar to surge and sway frequency. And yaw frequency is equal to the sway frequency, which is consistent with experimental results. Furthermore, comparing two cases, one fixed in rotation and the other one free in rotation, it can be concluded that release in the degree of rotation can decrease the sway amplitude but make no difference in the surge amplitude.展开更多
High-speed planing crafts have successfully evolved through developments in the last several decades.Classical approaches such as inviscid potential flow–based methods and the empirically based Savitsky method provid...High-speed planing crafts have successfully evolved through developments in the last several decades.Classical approaches such as inviscid potential flow–based methods and the empirically based Savitsky method provide general understanding for practical design.However,sometimes such analyses suffer inaccuracies since the air–water interface effects,especially in the transition phase,are not fully accounted for.Hence,understanding the behaviour at the transition speed is of fundamental importance for the designer.The fluid forces in planing hulls are dominated by phenomena such as flow separation at various discontinuities viz.,knuckles,chines and transom,with resultant spray generation.In such cases,the application of potential theory at high speeds introduces limitations.This paper investigates the simulation of modelling of the pre-planing behaviour with a view to capturing the air–water interface effects,with validations through experiments to compare the drag,dynamic trim and wetted surface area.The paper also brings out the merits of gridding strategies to obtain reliable results especially with regard to spray generation due to the air–water interface effects.The verification and validation studies serve to authenticate the use of the multi-gridding strategies on the basis of comparisons with simulations using model tests.It emerges from the study that overset/chimera grids give better results compared with single unstructured hexahedral grids.Two overset methods are investigated to obtain reliable estimation of the dynamic trim and drag,and their ability to capture the spray resulting from the air–water interaction.The results demonstrate very close simulation of the actual flow kinematics at steady-speed conditions in terms of spray at the air–water interface,drag at the pre-planing and full planing range and dynamic trim angles.展开更多
High-resolution numerical simulations for wake vortical flows have long been a challenge in rotor aerodynamics.A novel spectrum-optimized sixth-order Weighted Essentially NonOscillatory(WENO)scheme is proposed to disc...High-resolution numerical simulations for wake vortical flows have long been a challenge in rotor aerodynamics.A novel spectrum-optimized sixth-order Weighted Essentially NonOscillatory(WENO)scheme is proposed to discretize inviscid fluxes on moving overset grids,and the Improved Delayed Detached Eddy Simulation(IDDES)is employed to resolve turbulent vortices.The integration of these methods facilitates a comprehensive numerical investigation into the unsteady vortical flows over coaxial rotors in hover.The results highlight the substantial improvement in numerical resolution,in terms of both spatial structure and temporal evolution of unsteady multiscale wake vortices.Coaxial rotors in hover manifest three primary scales of wake vortex structures:(A)the helical evolution of primary blade tip vortices and the periodic occurrence of strong Blade-Vortex Interactions(BVI);(B)the continuous shedding of small-scale horseshoeshaped vortices from the trailing edges of rotor blades,forming the vortex sheets;(C)the emergence of small-scale secondary vortex braids induced by interactions between rotor tip vortices and the vortex sheets.These vortex structures and their interactions cause high-frequency oscillations in rotor disk loads and induce unsteady perturbations in the local flow field.Interactions among these primary vortices,coupled with the generation of secondary vortices,result in the dissipation,distortion,and breakup of the rotor tip vortices,ultimately forming a vortex soup.Notably,a substantial quantity of seemingly weak small-scale secondary vortex braids significantly contribute to energy dissipation during the evolution of wake vortices for coaxial rotors in hover.ó2024 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).展开更多
基金Innovation Foundation of BUAA for PhD Graduates (300522)
文摘A numerical method based on solutions of Euler/Navier-Stokes （N-S） equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advection upwind splitting method （AUSM） scheme, upwind AUSM/van Leer scheme, AUSM＋ scheme and AUSMDV scheme are implemented for spatial discretization, and van Albada limiter is also applied. For temporal discretization, both explicit Runge-Kutta method and implicit lower-upper symmetric Gauss-Seidel （LU-SGS） method are attempted. Simultaneously, overset grid technique is adopted. In detail, hole-map method is utilized to identify intergrid boundary points （IGBPs）. Furthermore, aimed at identification issue of donor elements, inverse-map method is implemented. Eventually, blade surface pressure distributions derived from numerical simulation are validated compared with experimental data, showing that all the schemes mentioned above have the capability to predict the rotor flow field accurately. At the same time, vorticity contours are illustrated for analysis, and other characteristics are also analyzed.
基金Dr.X.L.Li at the China Meteorological Administration.This study was supported by grants from the National Natural Science Foundation of China(Nos.41575103 and 91637210).
文摘With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.
基金supported by the Basic Research Program of Northwestern Polytechnical University (Grant No. JC201219)the Postdoctoral Science Foundation of China (Grant No. 20100481368)
文摘A method of unstructured dynamic overset grids is developed for the numerical simulation of helicopter unsteady rotorairframe aerodynamic interaction.For the effective treatment of the relative motion between the rotor and the airframe,the domain of flowfield is divided into two overset subzones,namely,a rotational subzone containing the blades and a stationary subzone containing the airframe.The overset part of two subzones is used to convect the flow variables of the two zones.The Taylor series expansion is used to obtain a second-order spatial accuracy,and dual-time stepping is adopted to improve the solution accuracy.Mesh deformation from the blade motion in forward flight is treated by using a spring analogy.Validation is made by numerically simulating the flows around a wind tunnel configuration and comparing the predicted time-averaged and instantaneous inflow and airframe surface pressure distributions with the experimental data.It shows that the present method is efficient and robust for the prediction of complicated unsteady rotor-airframe aerodynamic interaction phenomena.
基金Project supported by the National Natural Science Foundation of China(Grant No.52131102).
文摘In the present study, the performance of the NTNU Blind Test 1 wind turbine is analyzed in the computational fluid dynamics (CFD) simulations by using the CFD code FANS with structured overset grids. First, the numerical methods including the governing equations, the turbulence closure model, and the flow solver are introduced. In addition, the NTNU BT1 wind tunnel experiment is described. Then, structured overset grid blocks are generated in the computational domain with fully resolved wind turbine geometry, including the blades, hub, nacelle, and tower. Afterward, unsteady Reynolds averaged Navier-Stokes (RANS) simulations with the two-layer k - ε turbulence model are performed with an inlet velocity of 10 m/s and a tip-speed ratio (TSR) of 6. The overset-grid capability of FANS is leveraged to handle the rotation of the rotor. Finally, simulations are performed for a range of TSRs and a comparison is made among the present CFD results, other numerical results obtained from representative methods, and the experimental data. It is observed that the CFD-predicted thrust coefficients match the experimental measurement at low TSRs while under-predicting the values at high TSRs, and potential reasons for this deviation are discussed.
基金supported by the National Basic Research Program of China (2009CB724104)
文摘A cell-centred overset unstructured grids approach is developed.In this approach,the intergrid boundary is initially established based on the wall distance from the cell centre,and is then optimized.To accelerate the intergrid-boundary definition much more,a neighbor-toneighbor donor search algorithm based on advancing-front method is modified with the help of minimum cuboid boxes.To simplify the communications between different grid cell types and to obtain second-order spatial accuracy,a new interpolation method is constructed based on linear reconstruction,which employs only one layer of fringe cells along the intergrid boundary.For unsteady flows with relative motion,the intergrid boundary can be redefined fast and automatically.Several numerical results show that the present dynamic overset unstructured grids approach is accurate and reliable.
基金This work was supported in part by the National Science Foundation under Grant NSF-1913076.Any opinions,findings,and conclusions or recommendations expressed in this material are those of the author(s)and do not necessarily reflect the views of the National Science Foundation.
文摘We present high order accurate numerical methods for the wave equation that combines efficient Hermite methods with geometrically flexible discontinuous Galerkin methods by using overset grids.Near boundaries we use thin boundary fitted curvilinear grids and in the volume we use Cartesian grids so that the computational complexity of the solvers approaches a structured Cartesian Hermite method.Unlike many other overset methods we do not need to add artificial dissipation but we find that the built-in dissipation of the Hermite and discontinuous Galerkin methods is sufficient to maintain the stability.By numerical experiments we demonstrate the stability,accuracy,efficiency,and the applicability of the methods to forward and inverse problems.
基金National Natural Science Foundation of China (90716010)
文摘In this article, the overset grid assembly method is improved to efficiently solve several critical problems that occur when applying overset grids to the complicated geometries and moving body. First, instead of using the two-step searching (i.e. cut and paste), a single-step searching method based on the grid cell size is proposed to modify holes and optimize the grid overlapping automatically. Second, discrepancies between the wall surface representations, where the grids overlap, are handled by introducing the wall's normal directions to the hole-map and projecting the interpolated points on the recipient mesh into the donor mesh. Finally, the dynamic overset method is modified to address the complex moving body problem. At every time step, the initial hole surface of the previous time step is dynamically adjusted to accomplish hole cutting and avoid the time consuming hole-map procedure. Numerical experiments show that the enhanced overset grid assembly method obtains satisfactory results.
基金supported by the National Natural Science Foundation of China(Nos.11672133,12002161)Open Foundations of EDL Laboratory,China(No.EDL19092111)+2 种基金supports from National Science Foundation of Shaanxi Province,China(No.2021JQ-078)Fundamental Research Fund of Zhuhai,China(No.ZH22017003210011PWC)Aeronautical Science Foundation of China(No.F2021110)are acknowledged as well。
文摘In the present study,an efficient overset grid method by means of parallel implicit hole-cutting is proposed for the sake of simulating unsteady flows in aerospace engineering involving multiple bodies in relative movement.In view of the degraded computational efficiency and robustness for conventional overset grid assembly,several innovative techniques are developed within the overset grid assembly process,viz.,a bookkeeping alternative digital tree method to speed up the donor-cell searching,a fast parallel advancing front algorithm to accelerate the wall-distance calculation and a message-passing strategy with efficient information communication and lower storage expenditure within distributed computational architecture.The contribution of the developed techniques is evidenced by comparison with the existing alternative ways in terms of computing efficiency.Subsequently,the overset grid method is embedded into an inhouse programed URANS solver to examine its capability in predicting the flow field of complex applications such as helicopter,store separation and component deploying.Results show that the developed overset grid methodology is,in practice,able to resolve the aerodynamic characteristics of complex aerospace engineering with a high-fidelity flow topology and accuracy.
基金Supported by the UK Department for Transport,as part of the UK Shipping Office for Reducing Emissions(UK SHORE)Programme and the UK Engineering and Physical Sciences Research Council(EPSRC)[grant number EP/Y024605/1].
文摘The hydrodynamic performance of high-speed planing hulls has gained considerable interest,with recent advancements in computational fluid dynamics and hull design techniques enhancing the understanding of planing hull hydrodynamics.In this study,we conducted a numerical investigation using the Reynolds-averaged Navier-Stokes approach with overset grids to capture large motions at high speeds.This study aims to improve the hydrodynamic performances of planing hulls,specifically focusing on total resistance,trim,and sinkage.The initial Fridsma hull with a deadrise angle of 20°has been used for validation,demonstrating good agreement with measurements at different Froude numbers.Subsequently,new configurations based on the Fridsma hull have been designed by varying the deadrise angle,number of chines,and transverse steps.Our findings reveal a correlation between the deadrise angle,the number of chines,and the Froude number.As the deadrise angle increases,total resistance also increases.Additionally,a single chine yields superior results at higher Froude numbers,while multiple chines offer advantages at lower values.The introduction of transverse steps consistently increases total resistance,highlighting their role in improving planing hull performance.This research not only offers valuable insights into planing hull design but also leverages state-of-the-art numerical methods to advance the understanding of hydrodynamic behaviors at high ship speeds.
文摘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.
基金China Postdoctoral Science Foundation (20100481368)National Key Laboratory Foundation of China
文摘Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.
文摘A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or non-viscous, field equations with nonlinear boundary conditions applied to the free surface and fluid-solid interaction interfaces. An Arbitrary-LagrangianEulerian (ALE) mesh system is used to construct the numerical model. A multi-block numerical scheme of study is adopted allowing for the relative motion between moving overset grids, which are independent of one another. This provides a convenient method to overcome the difficulties in matching fluid meshes with large solid motions. Nonlinear numerical equations describing nonlinear fluid-solid interaction dynamics are derived through a numerical discretization scheme of study. A coupling iteration process is used to solve these numerical equations. Numerical examples are presented to demonstrate applications of the model developed.
基金the National Natural Science Foundation of China (51809169,51879159,51490675,11432009, 51579145)Chang Jiang Scholars Program (T2014099)+2 种基金Shanghai Excellent Academic Leaders Program (17XD1402300)Program for Professor of Special Appointment (Eastern Scholar)at Shanghai Institutions of Higher Learning (2013022)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China (2016-23/09).
文摘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.
文摘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.
基金co-supported by the National Natural Science Foundation of China(No.12072285)Shaanxi Science foundation for Distinguished Young Scholars,China(No.2020JC-13)。
文摘Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors.WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction,and its advantage is to reduce dissipation in smooth regions of flow while preserving high-resolution around discontinuities.Here WENO-K scheme is adopted to reconstruct left and right state values within the Roe Riemann solver updating the inviscid fluxes on a structured dynamic overset grid.To minimize the accuracy loss for high-order reconstruction on artificial boundaries of overset grid,a method of multilayer fringes is proposed to carry out interpolation between background grid and blade grid.Massively parallel computing considering automatic load balance on averagely partitioned overset grid is developed to reduce the wall-clock time of an unsteady simulation.Numerical results for Caradonna-Tung(C-T)rotor in hover at the conditions of subsonic and transonic tip Mach numbers show that the thrust coefficient error for the result of WENO-K scheme is no more than 3%.Compared with WENO-JS scheme,WENO-K scheme achieves about 40%improvement on accuracy of predicting rotor thrust with only 4.1%extra computational cost.More importantly,WENO-K scheme can capture more sophisticated unsteady flow structures and resolve tip vortices to a larger wake age with an increment of about 270°compared to WENO-JS scheme.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51379045,51679043,and 51679053)the National Defense Basic Scientific Research Foundation(Grant No.B2420132001)
文摘Parametric rolling is one of five types of the ship stability failure modes as proposed by IMO. The periodic change of the metacentric height is often considered as the internal cause of this phenomenon. Parametric rolling is a complex nonlinear hydrodynamic problem, often accompanied by large amplitude vertical motions of ships. In recent years,the Reynolds-averaged Navier–Stokes(RANS) equation simulations for viscous flows have made great progress in the field of ship seakeeping. In this paper, the parametric rolling for the C11 containership in regular waves is studied both experimentally and numerically. In the experiments, parametric rolling amplitudes at different drafts, forward speeds and wave steepnesses are analyzed. The differences in the steady amplitudes of parametric rolling are observed for two drafts. The effect of the incident wave steepness(or wave amplitude) is also studied, and this supports previous results obtained on limits of the stability for parametric rolling. In numerical simulations, the ship motions of parametric rolling are analyzed by use of the potential-flow and viscous-flow methods. In the viscousflow method, the Reynolds-averaged Navier–Stokes equations are solved using the overset grid method. The numerical accuracies of the two methods at different wave steepnesses 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.
基金Acknowledgements This work is supported by the National Natural Science Foundation of China （51490675, 51379125, 11432009, 51579145）, Chang Jiang Scholars Program （T2014099）, Shanghai Excellent Academic Leaders Program （17XD1402300）, Shanghai Key Laboratory of Marine Engineering （K2015-11）, Program for Professor of Special Appointment （Eastern Scholar） at Shanghai Institutions of Higher Learning （2013022）, Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China （2016-23/09） and Lloyd＇s Register Foundation for doctoral student, to which the authors are most grateful.
文摘This paper presents 3D （three-dimensional） CFD （computational fluid dynamic） simulation to analyse the FIR （flow-induced response） especially the yaw motion of a buoyancy can. The numerical cases are conducted with a buoyancy can under different reduced velocities utilizing our in-house code naoe-FOAM-SJTU, a solver based on the open source toolkit OpenFOAM. SST-DDES （shear-stress transport-delayed detached-eddy simulation） model is applied to handle the flowseparation and overset grid method is utilized to solve a large amplitude 6-DOF （6 degrees of freedom） motions. Free decay test and VIM （vortex-induced motion） test are built numerically. In VIM cases, the responses of trajectory, amplitude, frequency are calculated in a series of reduced velocities. With the increase of reduced velocity, yaw frequency is increased, which is similar to surge and sway frequency. And yaw frequency is equal to the sway frequency, which is consistent with experimental results. Furthermore, comparing two cases, one fixed in rotation and the other one free in rotation, it can be concluded that release in the degree of rotation can decrease the sway amplitude but make no difference in the surge amplitude.
文摘High-speed planing crafts have successfully evolved through developments in the last several decades.Classical approaches such as inviscid potential flow–based methods and the empirically based Savitsky method provide general understanding for practical design.However,sometimes such analyses suffer inaccuracies since the air–water interface effects,especially in the transition phase,are not fully accounted for.Hence,understanding the behaviour at the transition speed is of fundamental importance for the designer.The fluid forces in planing hulls are dominated by phenomena such as flow separation at various discontinuities viz.,knuckles,chines and transom,with resultant spray generation.In such cases,the application of potential theory at high speeds introduces limitations.This paper investigates the simulation of modelling of the pre-planing behaviour with a view to capturing the air–water interface effects,with validations through experiments to compare the drag,dynamic trim and wetted surface area.The paper also brings out the merits of gridding strategies to obtain reliable results especially with regard to spray generation due to the air–water interface effects.The verification and validation studies serve to authenticate the use of the multi-gridding strategies on the basis of comparisons with simulations using model tests.It emerges from the study that overset/chimera grids give better results compared with single unstructured hexahedral grids.Two overset methods are investigated to obtain reliable estimation of the dynamic trim and drag,and their ability to capture the spray resulting from the air–water interaction.The results demonstrate very close simulation of the actual flow kinematics at steady-speed conditions in terms of spray at the air–water interface,drag at the pre-planing and full planing range and dynamic trim angles.
文摘High-resolution numerical simulations for wake vortical flows have long been a challenge in rotor aerodynamics.A novel spectrum-optimized sixth-order Weighted Essentially NonOscillatory(WENO)scheme is proposed to discretize inviscid fluxes on moving overset grids,and the Improved Delayed Detached Eddy Simulation(IDDES)is employed to resolve turbulent vortices.The integration of these methods facilitates a comprehensive numerical investigation into the unsteady vortical flows over coaxial rotors in hover.The results highlight the substantial improvement in numerical resolution,in terms of both spatial structure and temporal evolution of unsteady multiscale wake vortices.Coaxial rotors in hover manifest three primary scales of wake vortex structures:(A)the helical evolution of primary blade tip vortices and the periodic occurrence of strong Blade-Vortex Interactions(BVI);(B)the continuous shedding of small-scale horseshoeshaped vortices from the trailing edges of rotor blades,forming the vortex sheets;(C)the emergence of small-scale secondary vortex braids induced by interactions between rotor tip vortices and the vortex sheets.These vortex structures and their interactions cause high-frequency oscillations in rotor disk loads and induce unsteady perturbations in the local flow field.Interactions among these primary vortices,coupled with the generation of secondary vortices,result in the dissipation,distortion,and breakup of the rotor tip vortices,ultimately forming a vortex soup.Notably,a substantial quantity of seemingly weak small-scale secondary vortex braids significantly contribute to energy dissipation during the evolution of wake vortices for coaxial rotors in hover.ó2024 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).
