Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellit...Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellites,inverting the Earth’s three-dimensional conductivity distribution on a global scale becomes attainable.A key requirement in the global conductivity inversion is to have a forward solver with high-accuracy and efficiency.In this study,a finite volume method for global electromagnetic induction forward modeling is developed based on unstructured grids.Arbitrary polyhedral grids are supported in our algorithms to obtain high geometric adaptability.We employ a cell-centered collocated variable arrangement which allows convenient discretization for complex geometries and straightforward implementation of multigrid technique.To validate the method,we test our code with two synthetic models and compare our finite volume results with an analytical solution and a finite element numerical solution.Good agreements are observed between our solution and other results,indicating acceptable accuracy of the proposed method.展开更多
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.展开更多
In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosi...In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosity source term appears in the momentum equation. The numerical model of the shallow water equations with porosity is presented with the finite volume method on unstructured grids and the modified Roe-type approximate Riemann solver. The source terms of the bed slope and porosity are both decomposed in the characteristic direction so that the numerical scheme can exactly satisfy the conservative property. The present model was tested with a dam break with discontinuous porosity and a flash flood in the Toce River Valley. The results show that the model can simulate the influence of obstructions, and the numerical scheme can maintain the flux balance at the interface with high efficiency and resolution.展开更多
A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method o...A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.展开更多
This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance ...This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.展开更多
A second-order mixing difference scheme with a limiting factor is deduced with the reconstruction gradient method and applied to discretizing the Navier-Stokes equation in an unstructured grid.The transform of nonorth...A second-order mixing difference scheme with a limiting factor is deduced with the reconstruction gradient method and applied to discretizing the Navier-Stokes equation in an unstructured grid.The transform of nonorthogonal diffusion items generated by the scheme in discrete equations is provided.The Delaunay triangulation method is improved to generate the unstructured grid.The computing program based on the SIMPLE algorithm in an unstructured grid is compiled and used to solve the discrete equations of two types of incompressible viscous flow.The numerical simulation results of the laminar flow driven by lid in cavity and flow behind a cylinder are compared with the theoretical solution and experimental data respectively.In the former case,a good agreement is achieved in the main velocity and drag coefficient curve.In the latter case,the numerical structure and development of vortex under several Reynolds numbers match well with that of the experiment.It is indicated that the factor difference scheme is of higher accuracy,and feasible to be applied to Navier-Stokes equation.展开更多
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.展开更多
A depth-integrated model for simulating wave-induced longshore current was developed with unstructured grids. Effects of surface roller and horizontal mixing under combined waves and currents were incorporated in the ...A depth-integrated model for simulating wave-induced longshore current was developed with unstructured grids. Effects of surface roller and horizontal mixing under combined waves and currents were incorporated in the numerical model. Recommended values of model coefficients were also proposed based on sensitivity analysis. Field observations and three series of laboratory measurements including two cases conducted on the plane beach and one implemented on the ideal inlet were employed to examine the predictive capability of this model. For the field case and laboratory cases conducted on the plane beach, numerical results were compared favorably with the measured data. For the case with an ideal inlet, simulated circulation pattern is supposed to be reasonable although some deviations between numerical results and measured data still can be detected.展开更多
By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical so...By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.展开更多
A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing t...A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.展开更多
In most TVD schemes, the r-factors were proposed according to the cell-centered（CC） finite volume method（FVM） framework for the numerical approximation to the convective term. However, it is questionable whether t...In most TVD schemes, the r-factors were proposed according to the cell-centered（CC） finite volume method（FVM） framework for the numerical approximation to the convective term. However, it is questionable whether those r-factors would be appropriate and effective for the vertex-centered（VC） FVM. In the paper, we collected five kinds of r-factor formulae and found out that only three of those, respectively by Bruner（1996）, Darwish and Moukalled（2003） and Cassuli and Zanolli（2005） can be formally extended to a context of the VC FVM. Numerical tests indicate that the TVD schemes and r-factors, after being extended and introduced to a context of the VC FVM, maintained their similar characteristics as in a context of the CC FVM. However, when the gradient-based r-factors and the SUPERBEE scheme were applied simultaneously, non-physical oscillations near the sharp step would appear. In the transient case, the oscillations were weaker in a context of the VC FVM than those in a context of the CC FVM, while the effect was reversed in the steady case. To eliminate disadvantages in the gradient-based r-factor formula, a new modification method by limiting values on the virtual node, namely Фu in the paper, was validated by the tests to effectively dissipate spurious oscillations.展开更多
Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important me...Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important means to reduce storm surge-related losses.Storm surge numerical models are important for storm surge forecasting.To further improve the performance of the storm surge forecast models,we developed a numerical storm surge forecast model based on an unstructured spherical centroidal Voronoi tessellation(SCVT)grid.The model is based on shallow water equations in vector-invariant form,and is discretized by Arakawa C grid.The SCVT grid can not only better describe the coastline information but also avoid rigid transitions,and it has a better global consistency by generating high-resolution grids in the key areas through transition refinement.In addition,the simulation speed of the model is accelerated by using the openACC-based GPU acceleration technology to meet the timeliness requirements of operational ensemble forecast.It only takes 37 s to simulate a day in the coastal waters of China.The newly developed storm surge model was applied to simulate typhoon-induced storm surges in the coastal waters of China.The hindcast experiments on the selected representative typhoon-induced storm surge processes indicate that the model can reasonably simulate the distribution characteristics of storm surges.The simulated maximum storm surges and their occurrence times are consistent with the observed data at the representative tide gauge stations,and the mean absolute errors are 3.5 cm and 0.6 h respectively,showing high accuracy and application prospects.展开更多
The three-dimensional Navier-Stokes equations were solved with the fractional step method where the hydrostatic pressure component was determined first, while the non-hydrostatic component of the pressure was computed...The three-dimensional Navier-Stokes equations were solved with the fractional step method where the hydrostatic pressure component was determined first, while the non-hydrostatic component of the pressure was computed from the pressure Poisson equation in which the coefficient matrix is positive definite and symmetric. The eddy viscosity was calculated from the efficient k-ε turbulence model. The resulting model is computationally efficient and unrestricted to the CFL condition. Computations with and without hydrostatic approximation were compared for the same cases to test the validity of the conventional hydrostatic pressure assumption. The model was verified against analytical solutions and experimental data, with excellent agreement.展开更多
The modified QUICK scheme on unstructured grid was used to improve the advection flux approximation, and the depth-averaged κ-ε turbulence model with the scheme based on FVM by SIMPLE series algorithm was establishe...The modified QUICK scheme on unstructured grid was used to improve the advection flux approximation, and the depth-averaged κ-ε turbulence model with the scheme based on FVM by SIMPLE series algorithm was established and applied to spur-dike flow computation. In this model, the over-relaxed approach was adopted to estimate the diffusion flux in view of its advantages in reducing errors and sustaining numerical stability usually encountered in non-orthogonal meshes. Two spur-dike cases with different defection angles （90°and 135°） were analyzed to validate the model. Computed results show that the predicted velocities and recirculation lengths are in good agreement with the observed data. Moreover, the computations on structured and unstructured grids were compared in terms of the approximately equivalent grid numbers. It can be concluded that the precision with unstructured grids is higher than that with structured grids in spite that the CPU time required is slightly more with unstructured grids Thus, it is significant to apply the method to numerical simulation of practical hydraulic engineering.展开更多
In this paper, an unstructured, collocated finite volume method for solvingthe Navier-Stokes equations was developed by virtue of auxiliary points. The derivatives weredetermined by the Gauss theorem. The proposed met...In this paper, an unstructured, collocated finite volume method for solvingthe Navier-Stokes equations was developed by virtue of auxiliary points. The derivatives weredetermined by the Gauss theorem. The proposed method could provide control volumes with arbitrarygeometry and preserve the second-order accuracy even if highly distorted grids are used. Althougharbitrary number of cell faces can be used, the hybrid quadrilateral/triangular grids are moredesirable for the simplicity of implementation and applications to engineering problems. Thepressure-velocity coupling was treated using a SIMPLE-like algorithm. The Generalized MinimumResidual (GMRES) method with the Incomplete LU (ILU) preconditioner was used to solve linearequations. Four test cases were studied for validating the proposed method. In using this method,grid quality is not important. Thus, engineers can pay mostly attention to physical mechanism ofproblems. Turbulence models can be simply integrated and the method can be straightforwardlyextended to treat three-dimensional problems.展开更多
A new method for sloshing simulation in a sway tank is present, in which the two phase interface is treated as a physical discontinuity, which can be captured by a well-designed high order scheme. Based on Normalized ...A new method for sloshing simulation in a sway tank is present, in which the two phase interface is treated as a physical discontinuity, which can be captured by a well-designed high order scheme. Based on Normalized Variable Diagram （NVD）, a high order discretization scheme with unstructured grids is realized, together with a numerical method for free surface flow with a fixed grid. This method is implemented in an in-house code General Transport Equation Analyzer （ GTEA ） which is an unstructured grids finite volume solver. The present method is first validated by available analytical solutions. A simulation for a 2-D rectangular tank at different excitation frequencies of the sway is carried out. A comparison with experimental data in literature and results obtained by commercial software CFX shows that the sloshing load on the monitor points agrees well with the experimental data, with the same grids, and the present method gives better results on the secondary peak. It is shown that the present method can simulate the free surface overturning and breakup phenomena.展开更多
In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstructio...In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.展开更多
The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the ...The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the bow shock.To improve the accuracy of aero-heating predictions,three multi-dimensional approaches on unstructured grids are coupled in our Reynolds-averaged Navier-Stokes(RANS)solver,including multi-dimensional upwind flux reconstruction(MUP),multi-dimensional limiter(MLP-u2)and multi-dimensional gradient reconstruction(MLR).The coupled multi-dimensional RANS solver is validated by several typical verification and validation(V&V)cases,including hypersonic flows over a cylinder,a blunt biconic,and a double-ellipsoid,with commonly used prism/tetrahedral hybrid grids.Finally,the coupled multi-dimensional solver is applied to simulating the heat flux distribution over a 3D engineering configuration,i.e.a Hermes-like space shuttle model.The obtained numerical results are compared with experimental data.The predicted results demonstrate that the coupled multi-dimensional approach has a good prediction capability on aerodynamic heating over a wide range of complex engineering configurations.展开更多
In particle transport simulations, radiation effects are olden described by the discrete ordinates （Sn） form of Boltzmann equation. In each ordinate direction, the solution is computed by sweeping the radiation flux...In particle transport simulations, radiation effects are olden described by the discrete ordinates （Sn） form of Boltzmann equation. In each ordinate direction, the solution is computed by sweeping the radiation flux across the grid. Parallel Sn sweep on an unstructured grid can be explicitly modeled as topological traversal through an equivalent directed acyclic graph （DAG）, which is a data-driven algorithm. Its traditional design using MPI model results in irregular communication of massive short messages which cannot be efficiently handled by MPI runtime. Meanwhile, in high-end HPC cluster systems, multicore has become the standard processor configuration of a single node. The traditional data-driven algorithm of Sn sweeps has not exploited potential advantages of multi-threading of multicore on shared memory. These advantages, however, as we shall demonstrate, could provide an elegant solution resolving problems in the previous MPI-only design. In this paper, we give a new design of data-driven parallel Sn sweeps using hybrid MPI and Pthread programming, namely Sweep-H, to exploit hierarchical parallelism of processes and threads. With special multi-threading techniques and vertex schedule policy, Sweep-H gets more efficient communication and better load balance. We further present an analytical performance model for Sweep-H to reveal why and when it is advantageous over former MPI counterpart. On a 64-node multicore cluster system with 12 cores per node, 768 cores in total, Sweep-H achieves nearly linear scalability for moderate problem sizes, and better absolute performance than the previous times speedup on 64 nodes）. MPI algorithm on more than 16 nodes （by up to two展开更多
E lement- partition- based methods for visualization of 3D unstructured grid data are presented. First, partition schemes for common elements, including curvilinear tetrahedra, pentahedra, hexahedra, etc., are given, ...E lement- partition- based methods for visualization of 3D unstructured grid data are presented. First, partition schemes for common elements, including curvilinear tetrahedra, pentahedra, hexahedra, etc., are given, so that complex elements can be divided into several rectilinear tetrahedra, and the visualization processes can be simplified.Then, a slice method for cloud map and an iso-surface method based on the partition schemes are described.展开更多
基金supported by the National Natural Science Foundation of China(41922027,4214200052)by the Macao Foundation+1 种基金by the Pre-research Project on Civil Aerospace Technologies No.D020308/D020303 funded by China National Space Administrationby the Macao Science and Technology Development Fund,grant No.0001/2019/A1。
文摘Global electromagnetic induction provides an efficient way to probe the electrical conductivity in the Earth’s deep interior.Owing to the increasing geomagnetic data especially from high-accuracy geomagnetic satellites,inverting the Earth’s three-dimensional conductivity distribution on a global scale becomes attainable.A key requirement in the global conductivity inversion is to have a forward solver with high-accuracy and efficiency.In this study,a finite volume method for global electromagnetic induction forward modeling is developed based on unstructured grids.Arbitrary polyhedral grids are supported in our algorithms to obtain high geometric adaptability.We employ a cell-centered collocated variable arrangement which allows convenient discretization for complex geometries and straightforward implementation of multigrid technique.To validate the method,we test our code with two synthetic models and compare our finite volume results with an analytical solution and a finite element numerical solution.Good agreements are observed between our solution and other results,indicating acceptable accuracy of the proposed method.
基金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.
基金supported by the National Natural Science Foundation of China (Grants No. 50909065 and 51109039)the National Basic Research Program of China (973 Program, Grant No. 2012CB417002)
文摘In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosity source term appears in the momentum equation. The numerical model of the shallow water equations with porosity is presented with the finite volume method on unstructured grids and the modified Roe-type approximate Riemann solver. The source terms of the bed slope and porosity are both decomposed in the characteristic direction so that the numerical scheme can exactly satisfy the conservative property. The present model was tested with a dam break with discontinuous porosity and a flash flood in the Toce River Valley. The results show that the model can simulate the influence of obstructions, and the numerical scheme can maintain the flux balance at the interface with high efficiency and resolution.
基金supported by the National Basic Research Program of China (2009CB724104)the National Natural Science Foundation of China (90716010)
文摘A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.
文摘This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.
基金Supported by National Natural Science Foundation of China (No. 10632050)
文摘A second-order mixing difference scheme with a limiting factor is deduced with the reconstruction gradient method and applied to discretizing the Navier-Stokes equation in an unstructured grid.The transform of nonorthogonal diffusion items generated by the scheme in discrete equations is provided.The Delaunay triangulation method is improved to generate the unstructured grid.The computing program based on the SIMPLE algorithm in an unstructured grid is compiled and used to solve the discrete equations of two types of incompressible viscous flow.The numerical simulation results of the laminar flow driven by lid in cavity and flow behind a cylinder are compared with the theoretical solution and experimental data respectively.In the former case,a good agreement is achieved in the main velocity and drag coefficient curve.In the latter case,the numerical structure and development of vortex under several Reynolds numbers match well with that of the experiment.It is indicated that the factor difference scheme is of higher accuracy,and feasible to be applied to Navier-Stokes equation.
基金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.
文摘A depth-integrated model for simulating wave-induced longshore current was developed with unstructured grids. Effects of surface roller and horizontal mixing under combined waves and currents were incorporated in the numerical model. Recommended values of model coefficients were also proposed based on sensitivity analysis. Field observations and three series of laboratory measurements including two cases conducted on the plane beach and one implemented on the ideal inlet were employed to examine the predictive capability of this model. For the field case and laboratory cases conducted on the plane beach, numerical results were compared favorably with the measured data. For the case with an ideal inlet, simulated circulation pattern is supposed to be reasonable although some deviations between numerical results and measured data still can be detected.
基金financially supported by the the National Natural Science Foundation of China(Grant No.51709054)the Public Science and Technology Research Funds Projects of Ocean(Grant Nos.201405025 and 201505019)
文摘By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.
基金financially supported by the Science and Technology Project of the Ministry of Transport (Grant No. 2011329224170)
文摘A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.41306078 and 41301414)the National Engineering Research Center for Inland Waterway Regulation and Key Laboratory of Hydraulic and Waterway Engineering of the Ministry of Education Program(Grant No.SLK2016B03)the Key Laboratory of the Inland Waterway Regulation of the Ministry of Transportation Program(Grant No.NHHD-201514)
文摘In most TVD schemes, the r-factors were proposed according to the cell-centered（CC） finite volume method（FVM） framework for the numerical approximation to the convective term. However, it is questionable whether those r-factors would be appropriate and effective for the vertex-centered（VC） FVM. In the paper, we collected five kinds of r-factor formulae and found out that only three of those, respectively by Bruner（1996）, Darwish and Moukalled（2003） and Cassuli and Zanolli（2005） can be formally extended to a context of the VC FVM. Numerical tests indicate that the TVD schemes and r-factors, after being extended and introduced to a context of the VC FVM, maintained their similar characteristics as in a context of the CC FVM. However, when the gradient-based r-factors and the SUPERBEE scheme were applied simultaneously, non-physical oscillations near the sharp step would appear. In the transient case, the oscillations were weaker in a context of the VC FVM than those in a context of the CC FVM, while the effect was reversed in the steady case. To eliminate disadvantages in the gradient-based r-factor formula, a new modification method by limiting values on the virtual node, namely Фu in the paper, was validated by the tests to effectively dissipate spurious oscillations.
基金The National Natural Science Foundation of China under contract No.42076214.
文摘Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important means to reduce storm surge-related losses.Storm surge numerical models are important for storm surge forecasting.To further improve the performance of the storm surge forecast models,we developed a numerical storm surge forecast model based on an unstructured spherical centroidal Voronoi tessellation(SCVT)grid.The model is based on shallow water equations in vector-invariant form,and is discretized by Arakawa C grid.The SCVT grid can not only better describe the coastline information but also avoid rigid transitions,and it has a better global consistency by generating high-resolution grids in the key areas through transition refinement.In addition,the simulation speed of the model is accelerated by using the openACC-based GPU acceleration technology to meet the timeliness requirements of operational ensemble forecast.It only takes 37 s to simulate a day in the coastal waters of China.The newly developed storm surge model was applied to simulate typhoon-induced storm surges in the coastal waters of China.The hindcast experiments on the selected representative typhoon-induced storm surge processes indicate that the model can reasonably simulate the distribution characteristics of storm surges.The simulated maximum storm surges and their occurrence times are consistent with the observed data at the representative tide gauge stations,and the mean absolute errors are 3.5 cm and 0.6 h respectively,showing high accuracy and application prospects.
文摘The three-dimensional Navier-Stokes equations were solved with the fractional step method where the hydrostatic pressure component was determined first, while the non-hydrostatic component of the pressure was computed from the pressure Poisson equation in which the coefficient matrix is positive definite and symmetric. The eddy viscosity was calculated from the efficient k-ε turbulence model. The resulting model is computationally efficient and unrestricted to the CFL condition. Computations with and without hydrostatic approximation were compared for the same cases to test the validity of the conventional hydrostatic pressure assumption. The model was verified against analytical solutions and experimental data, with excellent agreement.
基金the National Nature Science Foundation of China (Grant Nos. 50679019 and 50009001)the National Basic Research Program of China (973 Program, Grant No. 2008CB418202)the Social Technology Development Foundation of Jiangsu Province (Grant No.BS2006095)
文摘The modified QUICK scheme on unstructured grid was used to improve the advection flux approximation, and the depth-averaged κ-ε turbulence model with the scheme based on FVM by SIMPLE series algorithm was established and applied to spur-dike flow computation. In this model, the over-relaxed approach was adopted to estimate the diffusion flux in view of its advantages in reducing errors and sustaining numerical stability usually encountered in non-orthogonal meshes. Two spur-dike cases with different defection angles （90°and 135°） were analyzed to validate the model. Computed results show that the predicted velocities and recirculation lengths are in good agreement with the observed data. Moreover, the computations on structured and unstructured grids were compared in terms of the approximately equivalent grid numbers. It can be concluded that the precision with unstructured grids is higher than that with structured grids in spite that the CPU time required is slightly more with unstructured grids Thus, it is significant to apply the method to numerical simulation of practical hydraulic engineering.
文摘In this paper, an unstructured, collocated finite volume method for solvingthe Navier-Stokes equations was developed by virtue of auxiliary points. The derivatives weredetermined by the Gauss theorem. The proposed method could provide control volumes with arbitrarygeometry and preserve the second-order accuracy even if highly distorted grids are used. Althougharbitrary number of cell faces can be used, the hybrid quadrilateral/triangular grids are moredesirable for the simplicity of implementation and applications to engineering problems. Thepressure-velocity coupling was treated using a SIMPLE-like algorithm. The Generalized MinimumResidual (GMRES) method with the Incomplete LU (ILU) preconditioner was used to solve linearequations. Four test cases were studied for validating the proposed method. In using this method,grid quality is not important. Thus, engineers can pay mostly attention to physical mechanism ofproblems. Turbulence models can be simply integrated and the method can be straightforwardlyextended to treat three-dimensional problems.
基金Project supported by the China Postdoctoral Science Foundation (Grant No. 20100471016)the Fundamental Research Funds for Major Universities (Grant No. HEUCF 100307)
文摘A new method for sloshing simulation in a sway tank is present, in which the two phase interface is treated as a physical discontinuity, which can be captured by a well-designed high order scheme. Based on Normalized Variable Diagram （NVD）, a high order discretization scheme with unstructured grids is realized, together with a numerical method for free surface flow with a fixed grid. This method is implemented in an in-house code General Transport Equation Analyzer （ GTEA ） which is an unstructured grids finite volume solver. The present method is first validated by available analytical solutions. A simulation for a 2-D rectangular tank at different excitation frequencies of the sway is carried out. A comparison with experimental data in literature and results obtained by commercial software CFX shows that the sloshing load on the monitor points agrees well with the experimental data, with the same grids, and the present method gives better results on the secondary peak. It is shown that the present method can simulate the free surface overturning and breakup phenomena.
基金supported by the National Natural Science Foundation of China(Nos.91752114 and 11672160)
文摘In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.
基金the National Key Research&Development Program of China(2016YFB020071)the National Natural Science Foundation of China(Grants 11532016 and 11702315).
文摘The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the bow shock.To improve the accuracy of aero-heating predictions,three multi-dimensional approaches on unstructured grids are coupled in our Reynolds-averaged Navier-Stokes(RANS)solver,including multi-dimensional upwind flux reconstruction(MUP),multi-dimensional limiter(MLP-u2)and multi-dimensional gradient reconstruction(MLR).The coupled multi-dimensional RANS solver is validated by several typical verification and validation(V&V)cases,including hypersonic flows over a cylinder,a blunt biconic,and a double-ellipsoid,with commonly used prism/tetrahedral hybrid grids.Finally,the coupled multi-dimensional solver is applied to simulating the heat flux distribution over a 3D engineering configuration,i.e.a Hermes-like space shuttle model.The obtained numerical results are compared with experimental data.The predicted results demonstrate that the coupled multi-dimensional approach has a good prediction capability on aerodynamic heating over a wide range of complex engineering configurations.
基金supported by the National Natural Science Foundation of China under Grant Nos.60925009,61003062,61033009,61272134the National Basic 973 Program of China under Grant Nos.2011CB302502 and 2012CB316502
文摘In particle transport simulations, radiation effects are olden described by the discrete ordinates （Sn） form of Boltzmann equation. In each ordinate direction, the solution is computed by sweeping the radiation flux across the grid. Parallel Sn sweep on an unstructured grid can be explicitly modeled as topological traversal through an equivalent directed acyclic graph （DAG）, which is a data-driven algorithm. Its traditional design using MPI model results in irregular communication of massive short messages which cannot be efficiently handled by MPI runtime. Meanwhile, in high-end HPC cluster systems, multicore has become the standard processor configuration of a single node. The traditional data-driven algorithm of Sn sweeps has not exploited potential advantages of multi-threading of multicore on shared memory. These advantages, however, as we shall demonstrate, could provide an elegant solution resolving problems in the previous MPI-only design. In this paper, we give a new design of data-driven parallel Sn sweeps using hybrid MPI and Pthread programming, namely Sweep-H, to exploit hierarchical parallelism of processes and threads. With special multi-threading techniques and vertex schedule policy, Sweep-H gets more efficient communication and better load balance. We further present an analytical performance model for Sweep-H to reveal why and when it is advantageous over former MPI counterpart. On a 64-node multicore cluster system with 12 cores per node, 768 cores in total, Sweep-H achieves nearly linear scalability for moderate problem sizes, and better absolute performance than the previous times speedup on 64 nodes）. MPI algorithm on more than 16 nodes （by up to two
文摘E lement- partition- based methods for visualization of 3D unstructured grid data are presented. First, partition schemes for common elements, including curvilinear tetrahedra, pentahedra, hexahedra, etc., are given, so that complex elements can be divided into several rectilinear tetrahedra, and the visualization processes can be simplified.Then, a slice method for cloud map and an iso-surface method based on the partition schemes are described.
