This paper presents a comparative study of a meshless level-set method in the simulation of sloshing flows. The numerical moving particle semi-implicit （MPS） method and a grid based schemes of the MPS and level-set ...This paper presents a comparative study of a meshless level-set method in the simulation of sloshing flows. The numerical moving particle semi-implicit （MPS） method and a grid based schemes of the MPS and level-set methods are outlined and two violent sloshing cases are considered. The computed results are compared with the corresponding experimental data for validation. The impact pressure and the deformations of free surface induced by sloshing are comparatively analyzed, and are in good agreement with experimental ones. Results show that both the MPS and level-set methods are good tools for simulation of violent sloshing flows. However, the second pressure peaks as well as breaking and splashing of free surface by the MPS method are captured better than by the level-set method.展开更多
Numerical weather prediction(NWP) is a core technology in weather forecast and disaster mitigation. China’s NWP research and operational applications have been attached great importance by the meteorological communit...Numerical weather prediction(NWP) is a core technology in weather forecast and disaster mitigation. China’s NWP research and operational applications have been attached great importance by the meteorological community.Fundamental achievements have been made in the theories, methods, and NWP model development in China, which are of certain international impacts. In this paper, the scientific and technological progress of NWP in China since1949 is summarized. The current status and recent progress of the domestically developed NWP system-GRAPES(Global/Regional Assimilation and Pr Ediction System) are presented. Through independent research and development in the past 10 years, the operational GRAPES system has been established, which includes both regional and global deterministic and ensemble prediction models, with resolutions of 3-10 km for regional and 25-50 km for global forecasts. Major improvements include establishment of a new non-hydrostatic dynamic core, setup of four-dimensional variational data assimilation, and development of associated satellite application. As members of the GRAPES system, prediction models for atmospheric chemistry and air pollution, tropical cyclones, and ocean waves have also been developed and put into operational use. The GRAPES system has been an important milestone in NWP science and technology in China.展开更多
A 3D dynamic core of the non-hydrostatic model GRAPES（Global/Regional Assimilation and Prediction System） is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. T...A 3D dynamic core of the non-hydrostatic model GRAPES（Global/Regional Assimilation and Prediction System） is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. Three-dimensional Coriolis forcing is introduced to the new core, and full representation of the Coriolis forcing makes it straightforward to share code between the Yin and Yang subdomains. Similar to that in the original GRAPES model, a semi-implicit semi-Lagrangian scheme is adopted for temporal integration and advection with additional arrangement for cross-boundary transport. Under a non-centered second-order temporal and spatial discretization, the dry nonhydrostatic frame is summarized as the solution of an elliptical problem. The resulting Helmholtz equation is solved with the Generalized Conjugate Residual solver in cooperation with the classic Schwarz method. Even though the coefficients of the equation are quite different from those in the original model, the computational procedure of the new core is just the same. The bi-cubic Lagrangian interpolation serves to provide Dirichlet-type boundary conditions with data transfer between the subdomains. The dry core is evaluated with several benchmark test cases, and all the tests display reasonable numerical stability and computing performance. Persistency of the balanced flow and development of both the mountain-induced Rossby wave and Rossby–Haurwitz wave confirms the appropriate installation of the 3D Coriolis terms in the semi-implicit semi-Lagrangian dynamic core on the Yin-Yang grid.展开更多
A meshless numerical simulation method, the moving-particle semi-implicit method （MPS） is presented in this paper to study the sloshing phenomenon in ocean and naval engineering. As a meshless method, MPS uses parti...A meshless numerical simulation method, the moving-particle semi-implicit method （MPS） is presented in this paper to study the sloshing phenomenon in ocean and naval engineering. As a meshless method, MPS uses particles to replace the mesh in traditional methods, the governing equations are discretized by virtue of the relationship of particles, and the Poisson equation of pressure is solved by incomplete Cholesky conjugate gradient method （ICCG）, the free surface is tracked by the change of numerical density. A numerical experiment of viscous liquid sloshing tank was presented and compared with the result got by the difference method with the VOF, and an additional modification step was added to make the simulation more stable. The results show that the MPS method is suitable for the simulation of viscous liquid sloshing, with the advantage in arranging the particles easily, especially on some complex curved surface.展开更多
A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 a...A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 are chosen for the simulation. It shows that the result of MPS is very close to results of experiments or mesh-numerical simulations. In the simulation of MPS, vortices are found periodically in bilges of ship sections. In section S. S. 5.0 and section S. S. 7.0, which are close to the middle ship, two little vortices are found at different bilges of the section, in section S. S. 0.5, which is close to the bow, only one big vortex is found at the bottom of the section, these vortices patterns are consistent with the theory of Ikeda. The distribution of shear stress and pressure on the rolling hull of ship section is calculated. When vortices are in bilges of the section, the sign clmnge of pressure can be found, but in section S. S. 0.5, there is no sign change of pressure because only one vortex in the bottom of the section. With shear stress distribution, it can be found the shear stress in bilges is bigger than that at other part of the ship section. As the free surface is considered, the shear stress of both sides near the free surface is close to zero and even sign changed.展开更多
The recent progresses in the research and development of (NWP) in China are reviewed in this paper. The most impressive achievements are the development of direct assimilation of satellite irradiances with a 3DVAR (th...The recent progresses in the research and development of (NWP) in China are reviewed in this paper. The most impressive achievements are the development of direct assimilation of satellite irradiances with a 3DVAR (three-dimentional variational) data assimilation system and a non-hydrostatic modei with a semi-Lagrangian semi-implicit scheme. Progresses have also been made in modei physics and modei application to precipitation and environmental forecasts. Some scientific issues of great importance for further development are discussed.展开更多
The Global/Regional Assimilation and PrEdiction System(GRAPES) is a newly developed global non-hydrostatic numerical prediction model,which will become the next generation medium-range opera-tional model at China Mete...The Global/Regional Assimilation and PrEdiction System(GRAPES) is a newly developed global non-hydrostatic numerical prediction model,which will become the next generation medium-range opera-tional model at China Meteorological Administration(CMA).The dynamic framework of GRAPES is featuring with fully compressible equations,nonhydrostatic or hydrostatic optionally,two-level time semi-Lagrangian and semi-implicit time integration,Charney-Phillips vertical staggering,and complex three-dimensional pre-conditioned Helmholtz solver,etc.Concerning the singularity of horizontal momentum equations at the poles,the polar discretization schemes are described,which include adoption of Arakawa C horizontal grid with ν at poles,incorporation of polar filtering to maintain the computational stability,the correction to Helmholtz equation near the poles,as well as the treatment of semi-Lagrangian interpolation to improve the departure point accuracy,etc.The balanced flow tests validate the rationality of the treatment of semi-Lagrangian departure point calculation and the polar discretization during long time integration.Held and Suarez tests show that the conservation proper-ties of GRAPES model are quite good.展开更多
A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal d...A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal direction and z-level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohal Sea. The results are in good agreement both with the analytical solutions and measurement results.展开更多
Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model（moving particle semi-implicit（MPS）） is c...Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model（moving particle semi-implicit（MPS）） is coupled with either a Newtonian or a Hamiltonian Lagrangian structure model（MPS or HMPS） in a mathematically-physically consistent manner. The fluid model is founded on the solution of Navier-Stokes and continuity equations. The structure models are configured either in the framework of Newtonian mechanics on the basis of conservation of linear and angular momenta, or Hamiltonian mechanics on the basis of variational principle for incompressible elastodynamics. A set of enhanced schemes are incorporated for projection-based fluid model（Enhanced MPS）, thus, the developed coupled solvers for fluid structure interaction（FSI） are referred to as Enhanced MPS-MPS and Enhanced MPS-HMPS. Besides, two smoothed particle hydrodynamics（SPH）-based FSI solvers, being developed by the authors, are considered and their potential applicability and comparable performance are briefly discussed in comparison with MPS-based FSI solvers. The SPH-based FSI solvers are established through coupling of projection-based incompressible SPH（ISPH） fluid model and SPH-based Newtonian/Hamiltonian structure models, leading to Enhanced ISPH-SPH and Enhanced ISPH-HSPH. A comparative study is carried out on the performances of the FSI solvers through a set of benchmark tests, including hydrostatic water column on an elastic plate,high speed impact of an elastic aluminum beam, hydroelastic slamming of a marine panel and dam break with elastic gate.展开更多
This paper proposes an approach for rendering breaking waves out of large-scale of particle-based simulation. Moving particle semi-implicit (MPS) is used to solve the governing equation, and 2D simulation is expanded ...This paper proposes an approach for rendering breaking waves out of large-scale of particle-based simulation. Moving particle semi-implicit (MPS) is used to solve the governing equation, and 2D simulation is expanded to 3D representation by giving motion variation using fractional Brownian motion (fBm). The waterbody surface is reconstructed from the outlines of 2D simulation. The splashing effect is computed according to the properties of the particles. Realistic features of the wave are ren-dered on GPU, including the reflective and refractive effect and the effect of splash. Experiments showed that the proposed method can simulate large scale breaking waves efficiently.展开更多
An overlapping moving particle semi-implicit （MPS） method is applied for 3-D free surface flows based on our in-house particle solver MLParticle-SJTU. In this method, the coarse particles are distributed in the whol...An overlapping moving particle semi-implicit （MPS） method is applied for 3-D free surface flows based on our in-house particle solver MLParticle-SJTU. In this method, the coarse particles are distributed in the whole domain and the fine particles are distributed in the local region of interest at the same time. With the fine particles being generated and removed dynamically, an algorithm of generating particles based on the 3-D overlapping volume is developed. Then, a 3-D dam break flow with an obstacle is simulated to validate the overlapping MPS. The qualitative comparison among experimental data and the results obtained by the VOF and the MPS shows that the shape of the free surface obtained by the overlapping MPS is more accurate than that obtained by the UNI-coarse and close to that obtained by the UNI-fine in the overlapping domain. In addition, the water height and the impact pressure at Pi are also in an overall agreement with experimental data. Finally, the CPU time required by the overlapping MPS is about half of that required by the UNl-fine.展开更多
The main river, the Dongting Lake and river networks in the Jingjiang reach of the Yangtze River constitute a complex water system, for which a full 2-D hydrodynamic model is established instead of the traditional 1-D...The main river, the Dongting Lake and river networks in the Jingjiang reach of the Yangtze River constitute a complex water system, for which a full 2-D hydrodynamic model is established instead of the traditional 1-D or compound models for simulation of such complex systems, based on the latest developments of computer technologies and numerical methods. To better handle irregular boundaries and keep the computation cost well in a reasonable limit, unstructured grids of moderate scale are used. In addition, a dynamic boundary tracking method is proposed to simulate variable flow domains at different floods, especially, when the moderate scale gird can not describe flows in narrow river-network channels at low water levels. The t9 semi-implicit method and the Eulerian-Lagrangian Method （ELM） are adopted, which make the model unconditionally stable with respect to the gravity wave speed and Courant number restrictions. Properties and efficiency of the model are discussed, and it is concluded that the new model is robust and efficient enough for the simulation of a big, complex water system. Validation tests show that the simulation results agree well with field data. It takes about 0.96 h to complete the computation of a 76 d flood, which indicates that the model is efficient enough for engineering applications.展开更多
A semi-implicit and Eulerian - Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-moment turbulence closure model and tr...A semi-implicit and Eulerian - Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-moment turbulence closure model and transport equations of salinity and temperature. The simulation for flooding and drying of mudflats has been improved. The model is applied to Xiamen waters. Based on extensive survey data, water level elevation, temperature and salinity field along the eastern open boundary and at the Jiulong River inlets and runoffs are analyzed, specified and calibrated. The computed results show good agreement with the measured data, reproduce flooding, emergence of large and complex mudflat region.展开更多
Semi-implicit direct kinetics(SIDK)is an innovative method for the temporal discretization of neutronic equations proposed by J.Banfield.The key approximation of the SIDK method is to substitute a timeaveraged quantit...Semi-implicit direct kinetics(SIDK)is an innovative method for the temporal discretization of neutronic equations proposed by J.Banfield.The key approximation of the SIDK method is to substitute a timeaveraged quantity for the fission source term in the delayed neutron differential equations.Hence,these equations are decoupled from prompt neutron equations and an explicit analytical representation of precursor groups is obtained,which leads to a significant reduction in computational cost.As the fission source is not known in a time step,the original study suggested using a constant quantity pertaining to the previous time step for this purpose,and a reduction in the size of the time step was proposed to lessen the imposed errors.However,this remedy notably diminishes the main advantage of the SIDK method.We discerned that if the original method is properly introduced into the algorithm of the point-implicit solver along with some modifications,the mentioned drawbacks will be mitigated adequately.To test this idea,a novel multigroup,multi-dimensional diffusion code using the finitevolume method and a point-implicit solver is developed which works in both transient and steady states.In addition to the SIDK,two other kinetic methods,i.e.,direct kinetics and higher-order backward discretization,are programmed into the diffusion code for comparison with the proposed model.The final code is tested at different conditions of two well-known transient benchmark problems.Results indicate that while the accuracy of the improved SIDK is closely comparable with the best available kinetic methods,it reduces the total time required for computation by up to 24%.展开更多
Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering application...Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accur展开更多
A 3-D numerical formulation is proposed on the horizontal Cartesian, vertical sigma-coordinate grid for modeling non-hydrostatic pressure flee-surface flows. The pressure decomposition technique and 0 semi-implicit me...A 3-D numerical formulation is proposed on the horizontal Cartesian, vertical sigma-coordinate grid for modeling non-hydrostatic pressure flee-surface flows. The pressure decomposition technique and 0 semi-implicit method are used, with the solution procedure being split into two steps. First, with the implicit parts of non-hydrostatic pressures excluded, the provisional velocity field and free surface are obtained by solving a 2-D Poisson equation. Second, the theory of the differential operator is employed to derive the 3-D Poisson equation for non-hydrostatic pressures, which is solved to obtain the non-hydrostatic pressures and to update the provisional velocity field. When the non-orthogonal sigma-coordinate transformation is introduced, additional terms come into being, resulting in a 15-diagonal, diagonally dominant but unsymmetric linear system in the 3-D Poisson equation for non-hydrostatic pressures. The Biconjugate Gradient Stabilized （BiCGstab） method is used to solve the resulting 3-D unsymmetric linear system instead of the conjugate gradient method, which can only be used for symmetric, positive-definite linear systems. Three test cases are used for validations. The successful simulations of the small-amplitude wave, a supercritical flow over a ramp and a turbulent flow in the open channel indicate that the new model can simulate well non-hydrostatic flows, supercritical flows and turbulent flows.展开更多
A numerical simulation method for parachute Fluid-Structure Interaction （FSI） problem using Semi-Implicit Method for Pres- sure-Linked Equations （SIMPLE） algorithm is proposed. This method could be used in both co...A numerical simulation method for parachute Fluid-Structure Interaction （FSI） problem using Semi-Implicit Method for Pres- sure-Linked Equations （SIMPLE） algorithm is proposed. This method could be used in both coupling computation of para- chute FSI and flow field analysis. Both fiat circular parachute and conical parachute are modeled and simulated by this new method. Flow field characteristics at various angles of attack are further simulated for the conical parachute model. Compari- son with the space-time FSI technique shows that this method also provides similar and reasonable results.展开更多
The definition of a reference state close to the realistic atmosphere in an atmospheric model is essential for deriving prognostic deviations and improving numerical accuracy.In this study,a new dynamical framework al...The definition of a reference state close to the realistic atmosphere in an atmospheric model is essential for deriving prognostic deviations and improving numerical accuracy.In this study,a new dynamical framework allowing easy switching between a one-dimensional(1D)and a three-dimensional(3D)time-independent reference state is developed for the semi-implicit semi-Lagrangian solver in a global non-hydrostatic atmospheric model on Yin–Yang grids.The 3D reference state is introduced with consideration of additional horizontal gradient terms of referencestate terms,which is different from the 1D reference state.It is characterized by reduced magnitude of deviations,more accurate pressure gradient force,as well as alleviated numerical noise.Four idealized benchmark tests and multiple full-physics real-case forecasts are carried out to assess the impact of the 3D and 1D reference states.The 3D reference state shows significant advantages in the simulation of atmospheric transport and wave propagation in the idealized experiments.In the real-case forecasts,batched forecasts from June to August 2021 show a comprehensive improvement in medium-range prediction by using the 3D reference state.The new scheme achieves an enhanced prediction skill for large-scale circulation and extends the effective forecast period by 0.8 days in the Northern Hemisphere.展开更多
In this paper,two formulation theorems of time-difference fidelity schemes for general quadratic and cubic physical conservation laws are respectively constructed and proved,with earlier major conserving time-discreti...In this paper,two formulation theorems of time-difference fidelity schemes for general quadratic and cubic physical conservation laws are respectively constructed and proved,with earlier major conserving time-discretized schemes given as special cases.These two theorems can provide new mathematical basis for solving basic formulation problems of more types of conservative time- discrete fidelity schemes,and even for formulating conservative temporal-spatial discrete fidelity schemes by combining existing instantly conserving space-discretized schemes.Besides.the two theorems can also solve two large categories of problems about linear and nonlinear computational instability. The traditional global spectral-vertical finite-difference semi-implicit model for baroclinic primitive equations is currently used in many countries in the world for operational weather forecast and numerical simulations of general circulation.The present work,however,based on Theorem 2 formulated in this paper,develops and realizes a high-order total energy conserving semi-implicit time-difference fidelity scheme for global spectral-vertical finite-difference model of baroclinic primitive equations.Prior to this,such a basic formulation problem remains unsolved for long,whether in terms of theory or practice.The total energy conserving semi-implicit scheme formulated here is applicable to real data long-term numerical integration. The experiment of thirteen FGGE data 30-day numerical integration indicates that the new type of total energy conserving semi-implicit fidelity scheme can surely modify the systematic deviation of energy and mass conserving of the traditional scheme.It should be particularly noted that,under the experiment conditions of the present work,the systematic errors induced by the violation of physical laws of conservation in the time-discretized process regarding the traditional scheme designs(called type Z errors for short)can contribute up to one-third of the total systematic root-mean-square(RMS)error at the end of se展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.51379125,51411130131 and 11272120)the National Key Basic Research Development of China(973 Program,Grant No.2013CB036103)+1 种基金the High Te-chnology of Marine Research Project of the Ministry of Indu-stry and the Information Technology of Chinathe Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(Grant No.2013022)
文摘This paper presents a comparative study of a meshless level-set method in the simulation of sloshing flows. The numerical moving particle semi-implicit （MPS） method and a grid based schemes of the MPS and level-set methods are outlined and two violent sloshing cases are considered. The computed results are compared with the corresponding experimental data for validation. The impact pressure and the deformations of free surface induced by sloshing are comparatively analyzed, and are in good agreement with experimental ones. Results show that both the MPS and level-set methods are good tools for simulation of violent sloshing flows. However, the second pressure peaks as well as breaking and splashing of free surface by the MPS method are captured better than by the level-set method.
基金Supported by the National Key Research and Development Program of China(2017YFC1501900)Middle-and Long-term Development Strategic Research Project of the Chinese Academy of Engineering(2019-ZCQ-06)。
文摘Numerical weather prediction(NWP) is a core technology in weather forecast and disaster mitigation. China’s NWP research and operational applications have been attached great importance by the meteorological community.Fundamental achievements have been made in the theories, methods, and NWP model development in China, which are of certain international impacts. In this paper, the scientific and technological progress of NWP in China since1949 is summarized. The current status and recent progress of the domestically developed NWP system-GRAPES(Global/Regional Assimilation and Pr Ediction System) are presented. Through independent research and development in the past 10 years, the operational GRAPES system has been established, which includes both regional and global deterministic and ensemble prediction models, with resolutions of 3-10 km for regional and 25-50 km for global forecasts. Major improvements include establishment of a new non-hydrostatic dynamic core, setup of four-dimensional variational data assimilation, and development of associated satellite application. As members of the GRAPES system, prediction models for atmospheric chemistry and air pollution, tropical cyclones, and ocean waves have also been developed and put into operational use. The GRAPES system has been an important milestone in NWP science and technology in China.
基金supported by the National Natural Science Foundation of China (Grant No. 41175095)the National Key Technology R&D Program(Grant No. 2012BAC22B01)a research project of the Chinese Academy of Meteorological Sciences (Grant No. 2014Z001)
文摘A 3D dynamic core of the non-hydrostatic model GRAPES（Global/Regional Assimilation and Prediction System） is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. Three-dimensional Coriolis forcing is introduced to the new core, and full representation of the Coriolis forcing makes it straightforward to share code between the Yin and Yang subdomains. Similar to that in the original GRAPES model, a semi-implicit semi-Lagrangian scheme is adopted for temporal integration and advection with additional arrangement for cross-boundary transport. Under a non-centered second-order temporal and spatial discretization, the dry nonhydrostatic frame is summarized as the solution of an elliptical problem. The resulting Helmholtz equation is solved with the Generalized Conjugate Residual solver in cooperation with the classic Schwarz method. Even though the coefficients of the equation are quite different from those in the original model, the computational procedure of the new core is just the same. The bi-cubic Lagrangian interpolation serves to provide Dirichlet-type boundary conditions with data transfer between the subdomains. The dry core is evaluated with several benchmark test cases, and all the tests display reasonable numerical stability and computing performance. Persistency of the balanced flow and development of both the mountain-induced Rossby wave and Rossby–Haurwitz wave confirms the appropriate installation of the 3D Coriolis terms in the semi-implicit semi-Lagrangian dynamic core on the Yin-Yang grid.
基金the National Natural Science Foundation under Grant No.50579035
文摘A meshless numerical simulation method, the moving-particle semi-implicit method （MPS） is presented in this paper to study the sloshing phenomenon in ocean and naval engineering. As a meshless method, MPS uses particles to replace the mesh in traditional methods, the governing equations are discretized by virtue of the relationship of particles, and the Poisson equation of pressure is solved by incomplete Cholesky conjugate gradient method （ICCG）, the free surface is tracked by the change of numerical density. A numerical experiment of viscous liquid sloshing tank was presented and compared with the result got by the difference method with the VOF, and an additional modification step was added to make the simulation more stable. The results show that the MPS method is suitable for the simulation of viscous liquid sloshing, with the advantage in arranging the particles easily, especially on some complex curved surface.
基金the National Natural Science Foundation of China (Grant No.50579035)
文摘A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 are chosen for the simulation. It shows that the result of MPS is very close to results of experiments or mesh-numerical simulations. In the simulation of MPS, vortices are found periodically in bilges of ship sections. In section S. S. 5.0 and section S. S. 7.0, which are close to the middle ship, two little vortices are found at different bilges of the section, in section S. S. 0.5, which is close to the bow, only one big vortex is found at the bottom of the section, these vortices patterns are consistent with the theory of Ikeda. The distribution of shear stress and pressure on the rolling hull of ship section is calculated. When vortices are in bilges of the section, the sign clmnge of pressure can be found, but in section S. S. 0.5, there is no sign change of pressure because only one vortex in the bottom of the section. With shear stress distribution, it can be found the shear stress in bilges is bigger than that at other part of the ship section. As the free surface is considered, the shear stress of both sides near the free surface is close to zero and even sign changed.
文摘The recent progresses in the research and development of (NWP) in China are reviewed in this paper. The most impressive achievements are the development of direct assimilation of satellite irradiances with a 3DVAR (three-dimentional variational) data assimilation system and a non-hydrostatic modei with a semi-Lagrangian semi-implicit scheme. Progresses have also been made in modei physics and modei application to precipitation and environmental forecasts. Some scientific issues of great importance for further development are discussed.
基金Supported by the Ministry of Science and Technology of China (Grant Nos 2006BAC02B01 and 2006BAC03B03)the National High Technology Research and Development Program of China (863 Program) (Grant No 2006AA01A123)
文摘The Global/Regional Assimilation and PrEdiction System(GRAPES) is a newly developed global non-hydrostatic numerical prediction model,which will become the next generation medium-range opera-tional model at China Meteorological Administration(CMA).The dynamic framework of GRAPES is featuring with fully compressible equations,nonhydrostatic or hydrostatic optionally,two-level time semi-Lagrangian and semi-implicit time integration,Charney-Phillips vertical staggering,and complex three-dimensional pre-conditioned Helmholtz solver,etc.Concerning the singularity of horizontal momentum equations at the poles,the polar discretization schemes are described,which include adoption of Arakawa C horizontal grid with ν at poles,incorporation of polar filtering to maintain the computational stability,the correction to Helmholtz equation near the poles,as well as the treatment of semi-Lagrangian interpolation to improve the departure point accuracy,etc.The balanced flow tests validate the rationality of the treatment of semi-Lagrangian departure point calculation and the polar discretization during long time integration.Held and Suarez tests show that the conservation proper-ties of GRAPES model are quite good.
基金The Major State Basic Research Program of China under contract No. 2012CB417002the National Natural Science Foundation of China under contract Nos 50909065 and 51109039
文摘A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal direction and z-level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohal Sea. The results are in good agreement both with the analytical solutions and measurement results.
文摘Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model（moving particle semi-implicit（MPS）） is coupled with either a Newtonian or a Hamiltonian Lagrangian structure model（MPS or HMPS） in a mathematically-physically consistent manner. The fluid model is founded on the solution of Navier-Stokes and continuity equations. The structure models are configured either in the framework of Newtonian mechanics on the basis of conservation of linear and angular momenta, or Hamiltonian mechanics on the basis of variational principle for incompressible elastodynamics. A set of enhanced schemes are incorporated for projection-based fluid model（Enhanced MPS）, thus, the developed coupled solvers for fluid structure interaction（FSI） are referred to as Enhanced MPS-MPS and Enhanced MPS-HMPS. Besides, two smoothed particle hydrodynamics（SPH）-based FSI solvers, being developed by the authors, are considered and their potential applicability and comparable performance are briefly discussed in comparison with MPS-based FSI solvers. The SPH-based FSI solvers are established through coupling of projection-based incompressible SPH（ISPH） fluid model and SPH-based Newtonian/Hamiltonian structure models, leading to Enhanced ISPH-SPH and Enhanced ISPH-HSPH. A comparative study is carried out on the performances of the FSI solvers through a set of benchmark tests, including hydrostatic water column on an elastic plate,high speed impact of an elastic aluminum beam, hydroelastic slamming of a marine panel and dam break with elastic gate.
基金Project partly supported by the National Institute of Information andCommunication Technology (NICT), Japan
文摘This paper proposes an approach for rendering breaking waves out of large-scale of particle-based simulation. Moving particle semi-implicit (MPS) is used to solve the governing equation, and 2D simulation is expanded to 3D representation by giving motion variation using fractional Brownian motion (fBm). The waterbody surface is reconstructed from the outlines of 2D simulation. The splashing effect is computed according to the properties of the particles. Realistic features of the wave are ren-dered on GPU, including the reflective and refractive effect and the effect of splash. Experiments showed that the proposed method can simulate large scale breaking waves efficiently.
基金supported by the National Natural Science Foundation of China(Grant Nos.51379125,51490675,11432009 and 51579145)
文摘An overlapping moving particle semi-implicit （MPS） method is applied for 3-D free surface flows based on our in-house particle solver MLParticle-SJTU. In this method, the coarse particles are distributed in the whole domain and the fine particles are distributed in the local region of interest at the same time. With the fine particles being generated and removed dynamically, an algorithm of generating particles based on the 3-D overlapping volume is developed. Then, a 3-D dam break flow with an obstacle is simulated to validate the overlapping MPS. The qualitative comparison among experimental data and the results obtained by the VOF and the MPS shows that the shape of the free surface obtained by the overlapping MPS is more accurate than that obtained by the UNI-coarse and close to that obtained by the UNI-fine in the overlapping domain. In addition, the water height and the impact pressure at Pi are also in an overall agreement with experimental data. Finally, the CPU time required by the overlapping MPS is about half of that required by the UNl-fine.
基金supported by the Eleventh"Five-Year Plan" Science and Technology Program of China(Grant No. 2008BAB29B08)the National Key Basic Research Program of China(973 Program,Grant No.2007CB714100)supported by the Yangtze River Scientific Research Institute project(Grant No.CKSQ2010075)
文摘The main river, the Dongting Lake and river networks in the Jingjiang reach of the Yangtze River constitute a complex water system, for which a full 2-D hydrodynamic model is established instead of the traditional 1-D or compound models for simulation of such complex systems, based on the latest developments of computer technologies and numerical methods. To better handle irregular boundaries and keep the computation cost well in a reasonable limit, unstructured grids of moderate scale are used. In addition, a dynamic boundary tracking method is proposed to simulate variable flow domains at different floods, especially, when the moderate scale gird can not describe flows in narrow river-network channels at low water levels. The t9 semi-implicit method and the Eulerian-Lagrangian Method （ELM） are adopted, which make the model unconditionally stable with respect to the gravity wave speed and Courant number restrictions. Properties and efficiency of the model are discussed, and it is concluded that the new model is robust and efficient enough for the simulation of a big, complex water system. Validation tests show that the simulation results agree well with field data. It takes about 0.96 h to complete the computation of a 76 d flood, which indicates that the model is efficient enough for engineering applications.
文摘A semi-implicit and Eulerian - Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-moment turbulence closure model and transport equations of salinity and temperature. The simulation for flooding and drying of mudflats has been improved. The model is applied to Xiamen waters. Based on extensive survey data, water level elevation, temperature and salinity field along the eastern open boundary and at the Jiulong River inlets and runoffs are analyzed, specified and calibrated. The computed results show good agreement with the measured data, reproduce flooding, emergence of large and complex mudflat region.
文摘Semi-implicit direct kinetics(SIDK)is an innovative method for the temporal discretization of neutronic equations proposed by J.Banfield.The key approximation of the SIDK method is to substitute a timeaveraged quantity for the fission source term in the delayed neutron differential equations.Hence,these equations are decoupled from prompt neutron equations and an explicit analytical representation of precursor groups is obtained,which leads to a significant reduction in computational cost.As the fission source is not known in a time step,the original study suggested using a constant quantity pertaining to the previous time step for this purpose,and a reduction in the size of the time step was proposed to lessen the imposed errors.However,this remedy notably diminishes the main advantage of the SIDK method.We discerned that if the original method is properly introduced into the algorithm of the point-implicit solver along with some modifications,the mentioned drawbacks will be mitigated adequately.To test this idea,a novel multigroup,multi-dimensional diffusion code using the finitevolume method and a point-implicit solver is developed which works in both transient and steady states.In addition to the SIDK,two other kinetic methods,i.e.,direct kinetics and higher-order backward discretization,are programmed into the diffusion code for comparison with the proposed model.The final code is tested at different conditions of two well-known transient benchmark problems.Results indicate that while the accuracy of the improved SIDK is closely comparable with the best available kinetic methods,it reduces the total time required for computation by up to 24%.
基金supported by the National Natural Science Foundation of China(Grant No.50979059)
文摘Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accur
基金Project supported by the National Key Basic Research Program of China (973 Program, Grant No. 2007CB714100)the Key Project of National Natural Science Foundation of China (Grant No. 51039003)supported by the Yangtze River Scientific Research Institute project (Grant No. CKSQ2010075)
文摘A 3-D numerical formulation is proposed on the horizontal Cartesian, vertical sigma-coordinate grid for modeling non-hydrostatic pressure flee-surface flows. The pressure decomposition technique and 0 semi-implicit method are used, with the solution procedure being split into two steps. First, with the implicit parts of non-hydrostatic pressures excluded, the provisional velocity field and free surface are obtained by solving a 2-D Poisson equation. Second, the theory of the differential operator is employed to derive the 3-D Poisson equation for non-hydrostatic pressures, which is solved to obtain the non-hydrostatic pressures and to update the provisional velocity field. When the non-orthogonal sigma-coordinate transformation is introduced, additional terms come into being, resulting in a 15-diagonal, diagonally dominant but unsymmetric linear system in the 3-D Poisson equation for non-hydrostatic pressures. The Biconjugate Gradient Stabilized （BiCGstab） method is used to solve the resulting 3-D unsymmetric linear system instead of the conjugate gradient method, which can only be used for symmetric, positive-definite linear systems. Three test cases are used for validations. The successful simulations of the small-amplitude wave, a supercritical flow over a ramp and a turbulent flow in the open channel indicate that the new model can simulate well non-hydrostatic flows, supercritical flows and turbulent flows.
基金Supported by the Fundamental Research Funds(No.HEUCF140115)the Research Funds for State Key Laboratory of Ocean Engineering in Shanghai Jiao Tong University(No.1310)
基金supported by the National Natural Science Foundation of China (Grant No. 10577003)Monash University of Australia
文摘A numerical simulation method for parachute Fluid-Structure Interaction （FSI） problem using Semi-Implicit Method for Pres- sure-Linked Equations （SIMPLE） algorithm is proposed. This method could be used in both coupling computation of para- chute FSI and flow field analysis. Both fiat circular parachute and conical parachute are modeled and simulated by this new method. Flow field characteristics at various angles of attack are further simulated for the conical parachute model. Compari- son with the space-time FSI technique shows that this method also provides similar and reasonable results.
基金Supported by the National Natural Science Foundation of China(42375153,42075151,and 42205157).
文摘The definition of a reference state close to the realistic atmosphere in an atmospheric model is essential for deriving prognostic deviations and improving numerical accuracy.In this study,a new dynamical framework allowing easy switching between a one-dimensional(1D)and a three-dimensional(3D)time-independent reference state is developed for the semi-implicit semi-Lagrangian solver in a global non-hydrostatic atmospheric model on Yin–Yang grids.The 3D reference state is introduced with consideration of additional horizontal gradient terms of referencestate terms,which is different from the 1D reference state.It is characterized by reduced magnitude of deviations,more accurate pressure gradient force,as well as alleviated numerical noise.Four idealized benchmark tests and multiple full-physics real-case forecasts are carried out to assess the impact of the 3D and 1D reference states.The 3D reference state shows significant advantages in the simulation of atmospheric transport and wave propagation in the idealized experiments.In the real-case forecasts,batched forecasts from June to August 2021 show a comprehensive improvement in medium-range prediction by using the 3D reference state.The new scheme achieves an enhanced prediction skill for large-scale circulation and extends the effective forecast period by 0.8 days in the Northern Hemisphere.
基金The work is supported by the National Natural Science Foundation of China(49675267).
文摘In this paper,two formulation theorems of time-difference fidelity schemes for general quadratic and cubic physical conservation laws are respectively constructed and proved,with earlier major conserving time-discretized schemes given as special cases.These two theorems can provide new mathematical basis for solving basic formulation problems of more types of conservative time- discrete fidelity schemes,and even for formulating conservative temporal-spatial discrete fidelity schemes by combining existing instantly conserving space-discretized schemes.Besides.the two theorems can also solve two large categories of problems about linear and nonlinear computational instability. The traditional global spectral-vertical finite-difference semi-implicit model for baroclinic primitive equations is currently used in many countries in the world for operational weather forecast and numerical simulations of general circulation.The present work,however,based on Theorem 2 formulated in this paper,develops and realizes a high-order total energy conserving semi-implicit time-difference fidelity scheme for global spectral-vertical finite-difference model of baroclinic primitive equations.Prior to this,such a basic formulation problem remains unsolved for long,whether in terms of theory or practice.The total energy conserving semi-implicit scheme formulated here is applicable to real data long-term numerical integration. The experiment of thirteen FGGE data 30-day numerical integration indicates that the new type of total energy conserving semi-implicit fidelity scheme can surely modify the systematic deviation of energy and mass conserving of the traditional scheme.It should be particularly noted that,under the experiment conditions of the present work,the systematic errors induced by the violation of physical laws of conservation in the time-discretized process regarding the traditional scheme designs(called type Z errors for short)can contribute up to one-third of the total systematic root-mean-square(RMS)error at the end of se
