The finite volume method has been successfully applied in several engineering fields and has shown outstanding performance in fluid dynamics simulation. In this paper, the general framework for the simulation ofnear-w...The finite volume method has been successfully applied in several engineering fields and has shown outstanding performance in fluid dynamics simulation. In this paper, the general framework for the simulation ofnear-wellbore systems using the finite volume method is described. The mathematical model and the numerical model developed by the authors are presented and discussed. A radial geometry in the vertical plane was implemented so as to thoroughly describe near-wellbore phenomena. The model was then used to simulate injection tests in an oil reservoir through a horizontal well and proved very powerful to correctly reproduce the transient pressure behavior. The reason for this is the robustness of the method, which is independent of the gridding options because the discretization is performed in the physical space. The model is able to describe the phenomena taking place in the reservoir even in complex situations, i.e. in the presence of heterogeneities and permeability barriers, demonstrating the flexibility of the finite volume method when simulating non-conventional tests. The results are presented in comparison with those obtained with the finite difference numerical approach and with analytical methods, if possible.展开更多
The examination of an unstructured finite volume method for structural dynamics is assessed for simulations of systematic impact dynamics. A robust display dual-time stepping method is utilized to obtain time accurate...The examination of an unstructured finite volume method for structural dynamics is assessed for simulations of systematic impact dynamics. A robust display dual-time stepping method is utilized to obtain time accurate solutions. The study of impact dynamics is a complex problem that should consider strength models and state equations to describe the mechanical behavior of materials. The current method has several features, l） Discrete equations of unstructured finite volume method naturally follow the conservation law. 2） Display dual-time stepping method is suitable for the analysis of impact dynamic problems of time accurate solutions. 3） The method did not produce grid distortion when large deformation appeared. The method is validated by the problem of impact dynamics of an elastic plate with initial conditions and material properties. The results validate the finite element numerical data展开更多
Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mec...Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mechanical loadings with a tensile strain and a compressive strain on the hysteresis and butterfly loops is studied numerically. Different from the traditional finite element and finite difference methods, the finite volume method is applicable to simulate the ferroelectric phase transitions and properties of ferroelectric materials even for more realistic and physical problems.展开更多
Underground coal mining frequently uses longwalls.The occurrence of a potentially explosive mixture of methane and air is one of the most serious hazards.A large number of papers have applied numerical modeling of met...Underground coal mining frequently uses longwalls.The occurrence of a potentially explosive mixture of methane and air is one of the most serious hazards.A large number of papers have applied numerical modeling of methane propagation in research aimed at this problem.To date,none of the CFD simulations has considered the movement of the shearer in the analyses.This paper proposes an adaptation of a method used for the description of the movement of trains in tunnels to a specific geometry of a longwall district.The flow of the air-methane mixture was calculated using the finite volume method,in particular the k-w SST and SAS turbulence models.Due to the movement of the shearer,moving and deforming meshes were used for simulation of unsteady flows.Examples of solutions for two hypothetical cases are presented.Finally,the drawbacks and advantages of presented methods are discussed.Further development with the application of either local mesh variability or overset meshes is outlined.展开更多
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.展开更多
This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, where...This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time （SRT） lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.展开更多
This article proposes a high-order numerical method for a space distributed-order time-fractional diffusion equation.First,we use the mid-point quadrature rule to transform the space distributed-order term into multi-...This article proposes a high-order numerical method for a space distributed-order time-fractional diffusion equation.First,we use the mid-point quadrature rule to transform the space distributed-order term into multi-term fractional derivatives.Second,based on the piecewise-quadratic polynomials,we construct the nodal basis functions,and then discretize the multi-term fractional equation by the finite volume method.For the time-fractional derivative,the finite difference method is used.Finally,the iterative scheme is proved to be unconditionally stable and convergent with the accuracy O(σ^(2)+τ^(2-β)+h^(3)),whereτand h are the time step size and the space step size,respectively.A numerical example is presented to verify the effectiveness of the proposed method.展开更多
An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D...An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.展开更多
A flux difference splitting numerical scheme based on the finite volume method is applied to study ideal/resistive magnetohydrodynamics. The ideal/resistive MHD equations are cast as a set of hyperbolic conservation l...A flux difference splitting numerical scheme based on the finite volume method is applied to study ideal/resistive magnetohydrodynamics. The ideal/resistive MHD equations are cast as a set of hyperbolic conservation laws, and we develop a numerical capability to solve the weak solutions of these hyperbolic conservation laws by combining a multi-state Harten-Lax-Van Leer approximate Riemann solver with the hyperbolic divergence cleaning technique, high order shock-capturing reconstruction schemes, and a third order total variance diminishing Runge-Kutta time evolving scheme. The developed simulation code is applied to study the long time nonlinear evolution of the coalescence instability. It is verified that small structures in the instability oscillate with time and then merge into medium structures in a coherent manner. The medium structures then evolve and merge into large structures, and The physics of this interesting nonlinear dynamics this trend continues through all scale-lengths is numerically analyzed.展开更多
A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the gov...A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.展开更多
In this paper, a reactive dynamic user equilibrium model is extended to simulate two groups of pedestrians traveling on crossing paths in a continuous walking facility. Each group makes path choices to minimize the tr...In this paper, a reactive dynamic user equilibrium model is extended to simulate two groups of pedestrians traveling on crossing paths in a continuous walking facility. Each group makes path choices to minimize the travel cost to its destination in a reactive manner based on instantaneous information. The model consists of a conservation law equation coupled with an Eikonal-type equation for each group. The velocity-density relationship of pedestrian movement is obtained via an experimental method. The model is solved using a finite volume method for the conservation law equation and a fast-marching method for the Eikonal-type equation on unstructured grids. The numerical results verify the rationality of the model and the validity of the numerical method. Based on this continuum model, a number of results, e.g., the formation of strips or moving clusters composed of pedestrians walking to the same destination, are also observed.展开更多
To deal with the effect of compressible fluids on the supercavitating flow over the subsonic disk cavitator of a projectile, a finite volume method is formulated based on the ideal compressible potential theory. By us...To deal with the effect of compressible fluids on the supercavitating flow over the subsonic disk cavitator of a projectile, a finite volume method is formulated based on the ideal compressible potential theory. By using the continuity equation and Tait state equation as well as Riabouchinsky closure model, an“inverse problem”solution is presented for the supercavitating flow. According to the impenetrable condition on the surface of supercavity, a new iterative method for the supercavity shape is designed to deal with the effect of compressibility on the supercavity shape, pressure drag coefficient and density field. By this method, the very low cavitation number can be computed. The calculated results agree well with the experimental data and empirical formula. At the subsonic condition, the fluid compressibility will make supercavity length and radius increase. The supercavity expands, but remains spheroid. The effect on the first 1/3 part of supercavity is not obvious. The drag coefficient of projectile increases as the cavitation number or Mach number increases. With Mach number increasing, the compressibility is more and more significant. The compressibility must be considered as far as the accurate calculation of supercavitating flow is concerned.展开更多
In the present study, a three-dimensional computational fluid dynamics simulation together with experimental field measurements was applied to optimize the performance of an industrial hydrocyclone at Sarcheshmeh copp...In the present study, a three-dimensional computational fluid dynamics simulation together with experimental field measurements was applied to optimize the performance of an industrial hydrocyclone at Sarcheshmeh copper complex. In the simulation, the Eulerian–Eulerian approach was used for solid and liquid phases, the latter being water. In this approach, nine continuous phases were considered for the solid particles with different sizes and one continuous phase for water. The continuity and momentum equations with inclusion of buoyancy and drag forces were solved by the finite volume method. The k–e RNG turbulence model was used for modeling of turbulency. There was a good agreement between the simulation results and the experimental data. After validation of the model accuracy, the effect of inlet solid percentage, pulp inlet velocity, rod inserting in the middle of the hydrocyclone and apex diameter on hydrocyclone performance was investigated. The results showed that by decreasing the inlet solid percentage and increasing the pulp inlet velocity, the efficiency of hydrocyclone increased. Decreasing the apex diameter caused an increase in the hydrocyclone efficiency.展开更多
Offshore observation platforms are required to have great ability to resist waves when they are operating at sea. Investigation on the motion characteristics of the platforms in the sea can provide significant referen...Offshore observation platforms are required to have great ability to resist waves when they are operating at sea. Investigation on the motion characteristics of the platforms in the sea can provide significant reference values during the platform design procedure. In this paper, a series of numerical simulation on the interaction of a triple-hulled offshore observation platform with different incident waves is carried out. All of the simulations are implemented utilizing our own solver naoe-FOAM-SJTU, which is based and developed on the open source tools of OpenFOAM. Duration curves of motion characteristics and loads acting on the platform are obtained, and a comparison between the results of the amplitude in different incident waves is presented. The results show that the solver is competent in the simulation of motion response of platforms in waves.展开更多
In this paper, we report our study on a numerical fluid-structure interaction problem originally presented by Mok et al.(2001) in two dimensions and later studied in three dimensions by Valdés Vazquez(2007), Lomb...In this paper, we report our study on a numerical fluid-structure interaction problem originally presented by Mok et al.(2001) in two dimensions and later studied in three dimensions by Valdés Vazquez(2007), Lombardi(2012), and Trimarchi(2012). We focus on a 3D test case in which we evaluated the sensitivity of several input parameters on the fluid and structural results. In particular, this analysis provides a starting point from which we can look deeper into specific aspects of these simulations and analyze more realistic cases, e.g., in sails design. In this study, using the commercial software ADINATM, we addressed a well-known unsteadiness problem comprising a square box representing the fluid domain with a flexible bottom modeled with structural shell elements. We compared data from previously published work whose authors used the same numerical approach, i.e., a partitioned approach coupling a finite volume solver(for the fluid domain) and a finite element solver(for the solid domain). Specifically, we established several benchmarks and made comparisons with respect to fluid and solid meshes, structural element types, and structural damping, as well as solution algorithms. Moreover, we compared our method with a monolithic finite element solution method. Our comparisons of new and old results provide an outline of best practices for such simulations.展开更多
This paper proposes and analyzes a novel heating coil bundle with the tubes arranged in a multi-level manner.The bundle generates a heated cargo large-scale circulation that enables a superposition of the circulation-...This paper proposes and analyzes a novel heating coil bundle with the tubes arranged in a multi-level manner.The bundle generates a heated cargo large-scale circulation that enables a superposition of the circulation-driven forced convection on the buoyancy-driven natural convection,providing a more efficient mixed convection heat transfer mechanism.A simulationbased comparison of the proposed design and the conventional design is provided.The test case comprising an actual tank heating of an RMH 45 residual fuel oil by an 8-bar steam is simulated by a finite volume method and an OpenFOAM computational fluid dynamics software.The simulation results reveal that a 47.1%higher average heat transfer coefficient may be achieved,allowing a 32.0%reduction of the required heating coil area.展开更多
The exhaust smoke dispersion for a generic frigate is investigated numerically through the numerical solution of the governing fluid flow, energy, species and turbulence equations. The main objective of this work is t...The exhaust smoke dispersion for a generic frigate is investigated numerically through the numerical solution of the governing fluid flow, energy, species and turbulence equations. The main objective of this work is to obtain the effects of the yaw angle velocity ratio and buoyancy on the dispersion of the exhaust smoke The numerical method is based on the fully conserved control-volume representation of the fully elliptic Navier-Stokes equations. Turbulence is modeled using a two-equation （k-ε） model The flow visualization tests using a 1/100 scale model of the frigate in the wind tunnel were also carried out to determine the exhaust plume path and to validatc the computational results. The results show that down wash phenomena occurs for the yaw angles between ψ=10° and 20°. The results with different exhaust gas temperatures show that the buoyancy effect increases with the increasing of the exhaust gas temperature. However, its effect on the plume rise is less significant in comparison with its momentum. A good agreement between the predictions and experiment results is obtained.展开更多
Personal desktop platform with teraflops peak performance of thousands of cores is realized at the price of conventional workstations using the programmable graphics processing units(GPUs).A GPU-based parallel Euler/N...Personal desktop platform with teraflops peak performance of thousands of cores is realized at the price of conventional workstations using the programmable graphics processing units(GPUs).A GPU-based parallel Euler/Navier-Stokes solver is developed for 2-D compressible flows by using NVIDIA′s Compute Unified Device Architecture(CUDA)programming model in CUDA Fortran programming language.The techniques of implementation of CUDA kernels,double-layered thread hierarchy and variety memory hierarchy are presented to form the GPU-based algorithm of Euler/Navier-Stokes equations.The resulting parallel solver is validated by a set of typical test flow cases.The numerical results show that dozens of times speedup relative to a serial CPU implementation can be achieved using a single GPU desktop platform,which demonstrates that a GPU desktop can serve as a costeffective parallel computing platform to accelerate computational fluid dynamics(CFD)simulations substantially.展开更多
In order to improve lubricating characteristics of slippers in an axial piston pump, the combining center cavity slipper approach was proposed based on slipper shape and moving characteristic. The cylindrical coordina...In order to improve lubricating characteristics of slippers in an axial piston pump, the combining center cavity slipper approach was proposed based on slipper shape and moving characteristic. The cylindrical coordinate was used in the lubricant area and mesh was made. The blockweight approach was implemented to deal with non-coincidence of mesh and shallow recess border in numerical method. The finite control volume method was applied in calculating pressure distribution. The flow conservation equation and film thickness model were resolved through Gauss-Siedel relaxation iteration. The calculation and analysis results indicate that compared to the slipper （1） slip- per pressure distribution is improved; （2） hydrodynamic pressure of the combining slipper is greatly increased; （3） inclining degree is greatly reduced; （4） negative pressure in lubricant film disappear. So the combining center cavity slipper is lubricated better.展开更多
文摘The finite volume method has been successfully applied in several engineering fields and has shown outstanding performance in fluid dynamics simulation. In this paper, the general framework for the simulation ofnear-wellbore systems using the finite volume method is described. The mathematical model and the numerical model developed by the authors are presented and discussed. A radial geometry in the vertical plane was implemented so as to thoroughly describe near-wellbore phenomena. The model was then used to simulate injection tests in an oil reservoir through a horizontal well and proved very powerful to correctly reproduce the transient pressure behavior. The reason for this is the robustness of the method, which is independent of the gridding options because the discretization is performed in the physical space. The model is able to describe the phenomena taking place in the reservoir even in complex situations, i.e. in the presence of heterogeneities and permeability barriers, demonstrating the flexibility of the finite volume method when simulating non-conventional tests. The results are presented in comparison with those obtained with the finite difference numerical approach and with analytical methods, if possible.
文摘The examination of an unstructured finite volume method for structural dynamics is assessed for simulations of systematic impact dynamics. A robust display dual-time stepping method is utilized to obtain time accurate solutions. The study of impact dynamics is a complex problem that should consider strength models and state equations to describe the mechanical behavior of materials. The current method has several features, l） Discrete equations of unstructured finite volume method naturally follow the conservation law. 2） Display dual-time stepping method is suitable for the analysis of impact dynamic problems of time accurate solutions. 3） The method did not produce grid distortion when large deformation appeared. The method is validated by the problem of impact dynamics of an elastic plate with initial conditions and material properties. The results validate the finite element numerical data
基金Supported by the Research Starting Funds for Imported Talents of Ningxia University under Grant No BQD2012011
文摘Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mechanical loadings with a tensile strain and a compressive strain on the hysteresis and butterfly loops is studied numerically. Different from the traditional finite element and finite difference methods, the finite volume method is applicable to simulate the ferroelectric phase transitions and properties of ferroelectric materials even for more realistic and physical problems.
基金The results presented in this paper are the result of the PICTO research project titled"Production Face Environmental Risk Minimization in Coal and Lignite Mines”,No.800711,financed by the Research Programme of the Research Fund for Coal and Steel(RFCS)and Polish MNiSW No W93/FBWiS/2018 and the statutory research funds of the Institute.
文摘Underground coal mining frequently uses longwalls.The occurrence of a potentially explosive mixture of methane and air is one of the most serious hazards.A large number of papers have applied numerical modeling of methane propagation in research aimed at this problem.To date,none of the CFD simulations has considered the movement of the shearer in the analyses.This paper proposes an adaptation of a method used for the description of the movement of trains in tunnels to a specific geometry of a longwall district.The flow of the air-methane mixture was calculated using the finite volume method,in particular the k-w SST and SAS turbulence models.Due to the movement of the shearer,moving and deforming meshes were used for simulation of unsteady flows.Examples of solutions for two hypothetical cases are presented.Finally,the drawbacks and advantages of presented methods are discussed.Further development with the application of either local mesh variability or overset meshes is outlined.
文摘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 the Innovation Fund for Aerospace Science and Technology of China(Grant No.2009200066)the Aeronautical Science Fund of China(Grant No.20111453012)
文摘This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time （SRT） lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.
基金supported by the Natural and Science Foundation Council of China(11771059)Hunan Provincial Natural Science Foundation of China(2018JJ3519)Scientific Research Project of Hunan Provincial office of Education(20A022)。
文摘This article proposes a high-order numerical method for a space distributed-order time-fractional diffusion equation.First,we use the mid-point quadrature rule to transform the space distributed-order term into multi-term fractional derivatives.Second,based on the piecewise-quadratic polynomials,we construct the nodal basis functions,and then discretize the multi-term fractional equation by the finite volume method.For the time-fractional derivative,the finite difference method is used.Finally,the iterative scheme is proved to be unconditionally stable and convergent with the accuracy O(σ^(2)+τ^(2-β)+h^(3)),whereτand h are the time step size and the space step size,respectively.A numerical example is presented to verify the effectiveness of the proposed method.
基金the Youth Program of the National Natural Science Foundation of China(Grant Nos.11972272,12072246,and 12202331)the National Key Project,China(Grant No.GJXM92579)the Natural Science Basic Research Program of Shaanxi Province,China(Program No.2022JQ-028)。
文摘An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.
基金supported by the National Magnetic Confinement Fusion Science Program of China(Nos.2013GB111002,2013GB105003,2013GB111000,2014GB124005,2015GB111003)National Natural Science Foundation of China(Nos.11305171,11405208)+2 种基金JSPSNRF-NSFC A3 Foresight Program in the field of Plasma Physics(NSFC-11261140328)the Science Foundation of the Institute of Plasma Physics,Chinese Academy of Sciences(DSJJ-15-JC02)the CAS Program for the Interdisciplinary Collaboration Team
文摘A flux difference splitting numerical scheme based on the finite volume method is applied to study ideal/resistive magnetohydrodynamics. The ideal/resistive MHD equations are cast as a set of hyperbolic conservation laws, and we develop a numerical capability to solve the weak solutions of these hyperbolic conservation laws by combining a multi-state Harten-Lax-Van Leer approximate Riemann solver with the hyperbolic divergence cleaning technique, high order shock-capturing reconstruction schemes, and a third order total variance diminishing Runge-Kutta time evolving scheme. The developed simulation code is applied to study the long time nonlinear evolution of the coalescence instability. It is verified that small structures in the instability oscillate with time and then merge into medium structures in a coherent manner. The medium structures then evolve and merge into large structures, and The physics of this interesting nonlinear dynamics this trend continues through all scale-lengths is numerically analyzed.
基金Supported by the National Natural Science Foundation of China(11272153)
文摘A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.
基金supported by grants from the Research Grants Council of the Hong Kong Special Administrative Region, China (HKU 7183/06E)the University of Hong Kong (10207394)the National Natural Science Foundation of China (70629001 and 10771134)
文摘In this paper, a reactive dynamic user equilibrium model is extended to simulate two groups of pedestrians traveling on crossing paths in a continuous walking facility. Each group makes path choices to minimize the travel cost to its destination in a reactive manner based on instantaneous information. The model consists of a conservation law equation coupled with an Eikonal-type equation for each group. The velocity-density relationship of pedestrian movement is obtained via an experimental method. The model is solved using a finite volume method for the conservation law equation and a fast-marching method for the Eikonal-type equation on unstructured grids. The numerical results verify the rationality of the model and the validity of the numerical method. Based on this continuum model, a number of results, e.g., the formation of strips or moving clusters composed of pedestrians walking to the same destination, are also observed.
基金Foundation item： Supported by the National Natural Science Foundation of China （Grant No. 51309230）, and China Postdoctoral Science Foundation （Nos. 2014T70992 and 2013 M542531）We would like to thank Dr. Tao Miao for closely following our work and making several useful suggestions.
文摘To deal with the effect of compressible fluids on the supercavitating flow over the subsonic disk cavitator of a projectile, a finite volume method is formulated based on the ideal compressible potential theory. By using the continuity equation and Tait state equation as well as Riabouchinsky closure model, an“inverse problem”solution is presented for the supercavitating flow. According to the impenetrable condition on the surface of supercavity, a new iterative method for the supercavity shape is designed to deal with the effect of compressibility on the supercavity shape, pressure drag coefficient and density field. By this method, the very low cavitation number can be computed. The calculated results agree well with the experimental data and empirical formula. At the subsonic condition, the fluid compressibility will make supercavity length and radius increase. The supercavity expands, but remains spheroid. The effect on the first 1/3 part of supercavity is not obvious. The drag coefficient of projectile increases as the cavitation number or Mach number increases. With Mach number increasing, the compressibility is more and more significant. The compressibility must be considered as far as the accurate calculation of supercavitating flow is concerned.
文摘In the present study, a three-dimensional computational fluid dynamics simulation together with experimental field measurements was applied to optimize the performance of an industrial hydrocyclone at Sarcheshmeh copper complex. In the simulation, the Eulerian–Eulerian approach was used for solid and liquid phases, the latter being water. In this approach, nine continuous phases were considered for the solid particles with different sizes and one continuous phase for water. The continuity and momentum equations with inclusion of buoyancy and drag forces were solved by the finite volume method. The k–e RNG turbulence model was used for modeling of turbulency. There was a good agreement between the simulation results and the experimental data. After validation of the model accuracy, the effect of inlet solid percentage, pulp inlet velocity, rod inserting in the middle of the hydrocyclone and apex diameter on hydrocyclone performance was investigated. The results showed that by decreasing the inlet solid percentage and increasing the pulp inlet velocity, the efficiency of hydrocyclone increased. Decreasing the apex diameter caused an increase in the hydrocyclone efficiency.
基金Supported by the National Natural Science Foundation of China （Grant No. 50739004 and 11072154） Foundation of State Key Laboratory of Ocean Engineering of China （GKZD010059）+1 种基金 the Program for Professor of Special Appointment （Eastern Scholar） at Shanghai Institutions of Higher Learning （2008007） The Lloyd＇s Register Educational Trust （The LRET）
文摘Offshore observation platforms are required to have great ability to resist waves when they are operating at sea. Investigation on the motion characteristics of the platforms in the sea can provide significant reference values during the platform design procedure. In this paper, a series of numerical simulation on the interaction of a triple-hulled offshore observation platform with different incident waves is carried out. All of the simulations are implemented utilizing our own solver naoe-FOAM-SJTU, which is based and developed on the open source tools of OpenFOAM. Duration curves of motion characteristics and loads acting on the platform are obtained, and a comparison between the results of the amplitude in different incident waves is presented. The results show that the solver is competent in the simulation of motion response of platforms in waves.
文摘In this paper, we report our study on a numerical fluid-structure interaction problem originally presented by Mok et al.(2001) in two dimensions and later studied in three dimensions by Valdés Vazquez(2007), Lombardi(2012), and Trimarchi(2012). We focus on a 3D test case in which we evaluated the sensitivity of several input parameters on the fluid and structural results. In particular, this analysis provides a starting point from which we can look deeper into specific aspects of these simulations and analyze more realistic cases, e.g., in sails design. In this study, using the commercial software ADINATM, we addressed a well-known unsteadiness problem comprising a square box representing the fluid domain with a flexible bottom modeled with structural shell elements. We compared data from previously published work whose authors used the same numerical approach, i.e., a partitioned approach coupling a finite volume solver(for the fluid domain) and a finite element solver(for the solid domain). Specifically, we established several benchmarks and made comparisons with respect to fluid and solid meshes, structural element types, and structural damping, as well as solution algorithms. Moreover, we compared our method with a monolithic finite element solution method. Our comparisons of new and old results provide an outline of best practices for such simulations.
文摘This paper proposes and analyzes a novel heating coil bundle with the tubes arranged in a multi-level manner.The bundle generates a heated cargo large-scale circulation that enables a superposition of the circulation-driven forced convection on the buoyancy-driven natural convection,providing a more efficient mixed convection heat transfer mechanism.A simulationbased comparison of the proposed design and the conventional design is provided.The test case comprising an actual tank heating of an RMH 45 residual fuel oil by an 8-bar steam is simulated by a finite volume method and an OpenFOAM computational fluid dynamics software.The simulation results reveal that a 47.1%higher average heat transfer coefficient may be achieved,allowing a 32.0%reduction of the required heating coil area.
基金the financial support received from the Scientific Research Projects Office of Istanbul Technical University,stanbul, Turkey for the project titled as ‘Numerical and Experimental Investigation of Exhaust Smoke Dispersion for Naval Ships’
文摘The exhaust smoke dispersion for a generic frigate is investigated numerically through the numerical solution of the governing fluid flow, energy, species and turbulence equations. The main objective of this work is to obtain the effects of the yaw angle velocity ratio and buoyancy on the dispersion of the exhaust smoke The numerical method is based on the fully conserved control-volume representation of the fully elliptic Navier-Stokes equations. Turbulence is modeled using a two-equation （k-ε） model The flow visualization tests using a 1/100 scale model of the frigate in the wind tunnel were also carried out to determine the exhaust plume path and to validatc the computational results. The results show that down wash phenomena occurs for the yaw angles between ψ=10° and 20°. The results with different exhaust gas temperatures show that the buoyancy effect increases with the increasing of the exhaust gas temperature. However, its effect on the plume rise is less significant in comparison with its momentum. A good agreement between the predictions and experiment results is obtained.
基金supported by the National Natural Science Foundation of China (No.11172134)the Funding of Jiangsu Innovation Program for Graduate Education (No.CXLX13_132)
文摘Personal desktop platform with teraflops peak performance of thousands of cores is realized at the price of conventional workstations using the programmable graphics processing units(GPUs).A GPU-based parallel Euler/Navier-Stokes solver is developed for 2-D compressible flows by using NVIDIA′s Compute Unified Device Architecture(CUDA)programming model in CUDA Fortran programming language.The techniques of implementation of CUDA kernels,double-layered thread hierarchy and variety memory hierarchy are presented to form the GPU-based algorithm of Euler/Navier-Stokes equations.The resulting parallel solver is validated by a set of typical test flow cases.The numerical results show that dozens of times speedup relative to a serial CPU implementation can be achieved using a single GPU desktop platform,which demonstrates that a GPU desktop can serve as a costeffective parallel computing platform to accelerate computational fluid dynamics(CFD)simulations substantially.
基金Supported by the National Key Laboratory Foundation Project(9140C3403010903)
文摘In order to improve lubricating characteristics of slippers in an axial piston pump, the combining center cavity slipper approach was proposed based on slipper shape and moving characteristic. The cylindrical coordinate was used in the lubricant area and mesh was made. The blockweight approach was implemented to deal with non-coincidence of mesh and shallow recess border in numerical method. The finite control volume method was applied in calculating pressure distribution. The flow conservation equation and film thickness model were resolved through Gauss-Siedel relaxation iteration. The calculation and analysis results indicate that compared to the slipper （1） slip- per pressure distribution is improved; （2） hydrodynamic pressure of the combining slipper is greatly increased; （3） inclining degree is greatly reduced; （4） negative pressure in lubricant film disappear. So the combining center cavity slipper is lubricated better.
