A staggered control volume scheme for unstructured triangular grids

A Brinkman penalization method for compressible flows in complex geometries

Large scale unsteady motions in many practical engineering flows play a very important role and it is very unlikely that these unsteady flow features can be captured within the framework of Reynolds averaged Navier–Stokes approach. Large‐eddy simulation (LES) has become, arguably, the only practical numerical tool for predicting those flows more accurately since it is still not realistic to apply DNS to practical engineering flows with the current and near future available computing power.Numerical methods for the LES of turbulent flows in complex geometry have been developed and applied to predict practical engineering flows successfully. The method is based on body‐fitted curvilinear co‐ordinates with the contravariant velocity components of the general Navier–Stokes equations discretized on a staggered orthogonal mesh. For incompressible flow simulations the main source of computational expense is due to the solution of a Poisson equation for pressure. This is especially true for flows in complex geometry. A multigrid 3D pressure solver is developed to speed up the solution. In addition, the Poisson equation for pressure takes a simpler form with no cross‐derivatives when orthogonal mesh is used and hence resulting in increased convergence rate and producing more accurate solutions. Copyright © 2004 John Wiley & Sons, Ltd.

The development of better laser-based experimental methods and the fast rise in computer power has created an unprecedented shift in turbulent combustion research. The range of species and quantities measured and the advent of kHz-level planar diagnostics are now providing great insights in important phenomena and applications such as local and global extinction, pollutants, and spray combustion that were hitherto unavailable. In simulations, the shift to LES allows better representation of the turbulent flow in complex geometries, but despite the fact that the grid size is smaller than in RANS, the push towards realistic conditions and the need to include more detailed chemistry that includes very fast species and thin reaction zones emphasize the necessity of a sub-grid turbulent combustion model. The paper discusses examples from current research with experiments and modelling that focus on flame transients (self-excited oscillations, local extinction), sprays, soot emissions, and on practical applications. These demonstrate how current models are being validated by experimental data and the concerted efforts the community is taking to promote the modelling tools to industry. In addition, the various coordinated International Workshops on non-premixed, premixed, and spray flames, and on soot are discussed and some of their target flames are explored. These comprise flames that are relatively simple to describe from a fluid mechanics perspective but contain difficult-to-model combustion problems such as extinction, pollutants and multi-mode reaction zones. Recently, swirl spray flames, which are more representative of industrial devices, have been added to the target flames. Typically, good agreement is found with LES and some combustion models such as the progress variable - mixture fraction flamelet model, the Conditional Moment Closure, and the Transported PDF method, but predicting soot emissions and the condition of complete extinction in complex geometries is still elusive.

The development of better laser-based experimental methods and the fast rise in computer power has created an unprecedented shift in turbulent combustion research. The range of species and quantities measured and the advent of kHz-level planar diagnostics are now providing great insights in important phenomena and applications such as local and global extinction, pollutants, and spray combustion that were hitherto unavailable. In simulations, the shift to LES allows better representation of the turbulent flow in complex geometries, but despite the fact that the grid size is smaller than in RANS, the push towards realistic conditions and the need to include more detailed chemistry that includes very fast species and thin reaction zones emphasize the necessity of a sub-grid turbulent combustion model. The paper discusses examples from current research with experiments and modelling that focus on flame transients (self-excited oscillations, local extinction), sprays, soot emissions, and on practical applications. These demonstrate how current models are being validated by experimental data and the concerted efforts the community is taking to promote the modelling tools to industry. In addition, the various coordinated International Workshops on non-premixed, premixed, and spray flames, and on soot are discussed and some of their target flames are explored. These comprise flames that are relatively simple to describe from a fluid mechanics perspective but contain difficult-to-model combustion problems such as extinction, pollutants and multi-mode reaction zones. Recently, swirl spray flames, which are more representative of industrial devices, have been added to the target flames. Typically, good agreement is found with LES and some combustion models such as the progress variable - mixture fraction flamelet model, the Conditional Moment Closure, and the Transported PDF method, but predicting soot emissions and the condition of complete extinction in complex geometries is still elusive.

A control volume approach for integral viscoelastic models and its application to contraction flow of polymer melts

A lattice Boltzmann (LB) framework is developed for simulation of three-phase viscoelastic fluid flows in complex geometries. This model is based on a Rothman-Keller type model for immiscible multiphase flows which ensures mass conservation of each component in porous media even for a high density ratio. To account for the viscoelastic effects, the Maxwell constitutive relation is correctly introduced into the momentum equation, which leads to a modified lattice Boltzmann evolution equation for Maxwell fluids by removing the normal but excess viscous term. Our simulation tests indicate that this excess viscous term may induce significant errors. After three benchmark cases, the displacement processes of oil by dispersed polymer are studied as a typical example of three-phase viscoelastic fluid flow. The results show that increasing either the polymer intrinsic viscosity or the elastic modulus will enhance the oil recovery.

A fully automatic algorithm for three‐dimensional mesh generation is presented. The algorithm preserves the robustness and time efficiency of the finite octree technique; replacing octrees by more general 2–4–8‐ trees, it is capable of generating oriented meshes. In a post‐processing step, meshes are transformed in order to satisfy the Delaunay criterion, and so that non‐overlapping, closed control volumes can be defined for each mesh point through edge bisectors. The method is shown to be appropriate for modelling semiconductor devices, where the control volume scheme is the method of choice due to the peculiarities of the partial differential equations involved, and where locally refined and oriented meshes are essential to describe the relevant internal physical quantities adequately while keeping the number of mesh points practical.

HemeLB: A high performance parallel lattice-Boltzmann code for large scale fluid flow in complex geometries

An immersed boundary method for direct and large eddy simulation of stratified flows in complex geometry

A strongly conservative finite-volume procedure is presented for flows in complex geometries. The technique is based on a complete transformation of the governing equations, and physical velocity components, rather than the traditionally used Cartesian velocity components, are used as primitive variables. It was found that projecting the discretized vector transport equation in the direction of the covariant base vectors eliminated two substantial difficulties associated with flows in complex geometries. These difficulties stem from the presence of cross-pressure gradient terms and the need for a transformation between the different types of curvilinear velocity components in the mass conservation equation. It is shown that the present formulation ensures that the computational scheme is diagonally dominant. It was found that partially implicit treatment of nonorthogonal diffusion terms improved the convergence rate primarily for high-cell-Reynolds-number values. For nonstaggered grids, a new sol...

Large Eddy Simulations of a Stirred Tank Using the Lattice Boltzmann Method on a Nonuniform Grid

Numerical simulation of a square cavity was conducted to validate an implementation of an Immersed Boundary Method (IBM). The cavity consists of two differentially heated side walls and adiabatic top and bottom walls. A Cartesian grid is used with a finite volume, fractional step pressure correction method. Simulations use Dirichlet boundaries for vertical walls and Neumann boundaries for horizontal walls. The Immersed Boundary Method involves modifying the Navier--Stokes equations to include a forcing function in the momentum and energy equations that creates a virtual boundary. This method is useful because the boundary does not necessarily have to coincide with grid points; however, it is much less computationally expensive than other similar methods such as the cut cell method. The IBM is commonly used in simulations involving complex objects and can also accommodate moving boundaries. A standard numerical simulation with grid aligned with the boundary is first compared with previous results. The same geometry is then simulated by tilting the grid at various angles and using the IBM for each of the walls, and comparing these results with those initially obtained. We detail the implementation method and common problems associated with this. Velocity and temperature profiles are presented and the IBM is shown to maintain second order spatial accuracy. References E. A. Fadlun, R. Verzicco, P. Orlandi, and J. Mohd-Yusof. Combined immersed-boundary finite-difference methods for three-dimensional complex flow simulations. Journal of Computational Physics, 161(1):35--60, 2000. doi:10.1006/jcph.2000.6484 T. Gao, Y. H. Tseng, and X. Y. Lu. An improved hybrid cartesian/immersed boundary method for fluid-solid flows. International Journal for Numerical Methods in Fluids, 55(12):1189--1211, 2007. doi:10.1002/fld.1522 M. P. Kirkpatrick and S. W. Armfield. Experimental and large eddy simulation results for the purging of salt water from a cavity by an overflow of fresh water. International Journal of Heat and Mass Transfer, 48(2):341--359, 2005. doi:10.1016/j.ijheatmasstransfer.2004.08.016 M. P. Kirkpatrick, S. W. Armfield, and J. H. Kent. A representation of curved boundaries for the solution of the navier-stokes equations on a staggered three-dimensional cartesian grid. Journal of Computational Physics, 184(1):1--36, 2003. doi:10.1016/S0021-9991(02)00013-X B. P. Leonard and Simin Mokhtari. Beyond first-order upwinding: The ultra-sharp alternative for non-oscillatory steady-state simulation of convection. International Journal for Numerical Methods in Engineering, 30(4):729--766, 1990. doi:10.1002/nme.1620300412 S. H. Peng and L. Davidson. Large eddy simulation for turbulent buoyant flow in a confined cavity. volume 22, pages 323--331, 2001. doi:10.1016/S0142-727X(01)00095-9 C. S. Peskin. Numerical-analysis of blood-flow in heart. Journal of Computational Physics, 25(3):220--252, 1977. doi:10.1016/0021-9991(77)90100-0 J. Salat, S. Xin, P. Joubert, A. Sergent, F. Penot, and P. Le Quere. Experimental and numerical investigation of turbulent natural convection in a large air-filled cavity. International Journal of Heat and Fluid Flow, 25:824--832, 2004. doi:10.1016/j.ijheatfluidflow.2004.04.003 Y. H. Tseng and J. H. Ferziger. A ghost-cell immersed boundary method for flow in complex geometry. Journal of Computational Physics, 192(2):593--623, 2003. doi:10.1016/j/jcp.2003.07.024 N. Zhang and Z. C. Zheng. An improved direct-forcing immersed-boundary method for finite difference applications. Journal of Computational Physics, 221(1):250--268, 2007. doi:10.1016/j.jcp.2006.06.012 N. Zhang, Z. C. Zheng, and S. Eckels. Study of heat-transfer on the surface of a circular cylinder in flow using an immersed-boundary method. International Journal of Heat and Fluid Flow, 29(6):1558--1566, 2008. doi:10.1016/j.ijheatfluidflow.2008.08.009

Horizontal wells can have very complex flow geometries, in part due to interaction between the main flow stream and the influxes along the wellbore, and also due to completion type. In this study, the flow behavior in horizontal wells with a single perforation and with multiple perforations of perforation densities equivalent to 1, 2 and 4 shots per foot were investigated. A new test facility was designed and constructed to simulate single phase liquid flow in horizontal wells. Experiments were conducted with Reynolds numbers ranging from 5,000 to 60,000 and influx to main flow rate ratios ranging from 1/5 to 1/100 for the single injection case and from 1/100 to 1/2000 for the multiple injection case and for the no influx case. A general horizontal well friction factor expression was developed using the principles of conservation of mass and momentum. Horizontal well friction factor correlations were developed by applying experimental data to the general friction factor expressions. It was observed that the friction factor for a perforated pipe with fluid injection can be either smaller or greater than that for a smooth pipe, depending on influx to main flow rate ratios. The proposed friction factor correlation can be used in any horizontal well model which considers pressure variation along the wellbore. Introduction The most commonly used assumptions in studying horizontal well production behavior are: infinite conductivity, and uniform influx. Infinite conductivity assumes no pressure drop along a horizontal well, and uniform influx assumes that the influx from the reservoir is constant along a horizontal well. It has been argued in the literature that the infinite conductivity wellbore assumption is adequate for describing flow behavior in horizontal wells. Although this may be a good assumption in situations where the pressure drop along the horizontal section of the wellbore is negligible compared to that in the reservoir, it is also reasonable to expect the friction and acceleration effects to cause noticeable pressure drops in long horizontal wellbores. In a horizontal well, depending upon the completion method, fluid may enter the wellbore at various locations along the well length. The distance between perforations may not be sufficient to achieve a stabilized velocity profile, and this may lead to different pressure behavior than for fully developed flow. The pressure distribution in a horizontal well can influence the well completion and well profile design, as well as having an impact on the production behavior of the well. Therefore, both the pressure drop vs. flow behavior along the well and the relationship between the pressure drop along the well and the influx from the reservoir need to be understood. Flow behavior in porous wall pipes or plates has been studied by several investigators in other disciplines. Although they provide useful information to the study of horizontal well flow behavior, the results can not readily be applied in the petroleum industry. The petroleum industry started to investigate horizontal wellbore hydraulics in the late 1980's. Investigators conducted analytical or experimental studies to investigate different aspects of horizontal well flow behavior. A friction factor correlation for horizontal wellbores was proposed by Asheim et al. which includes accelerational pressure losses due to continuous fluid influx along the wellbore. However, the literature survey reveals that experimental data are limited and theoretical studies are inconclusive for predicting frictional pressure loss in horizontal wells. In this study, an experimental and theoretical investigation of the flow behavior in a horizontal pipe with fluid injection from a single injection point in the pipe wall and from multiple injection points in the pipe wall with perforation densities equivalent to 1, 2 and 4 shots per foot has been conducted. New correlations based on acquired experimental data were developed for the prediction of friction factors in a horizontal wellbore. Model Development A general model was developed to predict horizontal well friction factors. P. 603

Meshfree-enriched simplex elements with strain smoothing for the finite element analysis of compressible and nearly incompressible solids

In the present study, the dynamic Vreman model based on the global equilibrium between the subgrid-scale dissipation and the viscous dissipation with a global model coefficient are applied to large eddy simulation of turbulent flow in complex geometries. Distributions of pressure, velocity and vorticity as well as some flow structure are gained, which is helpful to examine the performance of SGS model and understand the flow characters in complex geometries.