Solution of the implicitly discretised fluid flow equations by operator-splitting

Issa, R. I.

doi:10.1016/0021-9991(86)90099-9

Cited by 4,474 publications

(1,833 citation statements)

References 8 publications

Supporting

Mentioning

1,782

Contrasting

Unclassified

Order By: Relevance

“…The approach features a variable predictor-corrector variant (with splitting error control) of the PISO [81] algorithm for the velocity-pressure coupling with a TVD scheme [82] used to reduce numerical diffusion. A more detailed description of the numerical procedure is available elsewhere [31].…”

Section: Numerical Proceduresmentioning

confidence: 99%

The impact of dilatation, scrambling, and pressure transport in turbulent premixed flames

Tian

Lindstedt

2017

Combustion Theory and Modelling

View full text Add to dashboard Cite

Premixed turbulent flames feature strong interactions between chemical reactions and turbulence that affect scalar and turbulence statistics. The focus of the present work is on clarifying the impact of pressure dilatation/flamelet scrambling effects with a comprehensive second-moment closure used for evaluation purposes. Model extensions that take into account flamelet orientation and molecular diffusion are derived. Isothermal pressure transport is included with an additional variable density contribution derived for the flamelet regime of combustion. The full closure is assessed by comparisons with direct numerical simulations (DNS) of statistically "steady" fully developed premixed turbulent planar flames at different expansion ratios. Subsequently, the prediction of lean premixed turbulent methane-air flames featuring fractal grid generated turbulence in an opposed jet geometry is considered. The overall agreement shows that "dilatation" effects contribute to counter-gradient transport and can also increase the turbulent kinetic energy significantly. Levels of anisotropy are broadly consistent with the DNS data and key aspects of opposed jet flames are well predicted. However, it is also shown that complications arise due to interactions between the imposed pressure gradient and combustion and that redistribution is affected along with the scalar flux at the leading edge. The latter is strongly affected by the reaction rate closure and, potentially, by pressure transport. Overall, the derived models offer significant improvements and can readily be applied to the modelling of premixed turbulent flames at practical rates of heat release.

show abstract

Section: Numerical Proceduresmentioning

confidence: 99%

The impact of dilatation, scrambling, and pressure transport in turbulent premixed flames

Tian

Lindstedt

2017

Combustion Theory and Modelling

View full text Add to dashboard Cite

show abstract

“…It is therefore not necessary to update the velocity more than once per time step. This is essentially the Pressure-Implicit with Splitting of Operators (PISO) algorithm [32] and is chosen in OpenFOAM by setting the nOuterCorrectors to 1, and performing several pressure corrections, typically 3. PISO is efficient, and the accuracy is good with the strict Courant number limit.…”

Section: Openfoam Solver Interdymfoammentioning

confidence: 99%

Numerical Simulations of Sloshing and the Thermodynamic Response Due to Mixing

Grotle¹,

Æsøy²

2017

Energies

View full text Add to dashboard Cite

Abstract:In this paper, we apply computational fluid dynamics (CFD) to study the thermodynamic response enhanced by sloshing inside liquefied natural gas (LNG) fuel tanks. An existing numerical solver provided by OpenFOAM is used to simulate sloshing in a model scaled tank of similar form to an LNG fuel tank. The interface area has been estimated for different sloshing regimes on three different numerical grids representing the tank in 3D. Estimating the interface area is done by performing a grid-independence study. In the most severe sloshing conditions, convergence is not achieved. By combining the results from experiments and CFD, it is found that the interface area and the condensation mass flow rate are in phase for the most severe sloshing condition. The existing CFD solver is modified to determine the pressure drop. The simulation results are compared to the experimental data, and the results are acceptable and thereby show a potential in applying CFD to predict the thermodynamic response due to sloshing. By plotting the temperature contours, indications are found that the exchange of cold bulk and saturated liquid due to sloshing has a significant influence on the thermodynamic response.

show abstract

“…The solution procedure is based on the well known PISO procedure [27] and consists of the following steps: The efficiency of the above described solution procedure can be substantially increased if only interfacial mesh points are moved in step (b) instead of moving all mesh points. In that case, the entire mesh is deformed before the start of the outer iteration loop, according to the total interface displacement obtained in the previous time step.…”

Section: Solution Proceduresmentioning

confidence: 99%

“…The incompressible fluid flow model discretised on moving collocated mesh is solved using the iterative PISO algorithm [27] where the Rhie-Chow momentum interpolation method [28] is employed to evaluate cell face velocity in order to eliminate the checkerboard pressure field effect. Here we came across a limitation of the original Rhie and Chow interpolation connected with the application of very small time steps (due to semi-implicit treatment of surface tension forces).…”

Section: Introductionmentioning

confidence: 99%

A moving mesh finite volume interface tracking method for surface tension dominated interfacial fluid flow

2012

View full text Add to dashboard Cite

a b s t r a c tThis paper describes a moving mesh interface tracking method implemented in OpenFOAM for simulating three-dimensional (3-D) incompressible and immiscible two-phase interfacial fluid flows with dominant surface tension forces. Collocated finite volume (FV) method is used for spatial discretisation of Navier-Stokes equations on moving polyhedral mesh. The mesh consists of two parts separated on interface. Fluid flow is solved on each mesh separately and coupling is accomplished in an iterative manner by enforcing the kinematic and dynamic condition at the interface. Surface tension force is calculated on arbitrary polygonal surface mesh with second order accuracy using a ''force-conservative'' approach. Arbitrary polyhedral mesh adapts to the time-varying shape of the interface using vertex-based automatic mesh motion solver which calculates the motion of internal points based on the prescribed motion of interface points by solving the variable diffusivity Laplace equation discretised using the finite element method. The overall solution procedure based on iterative PISO algorithm with modified Rhie-Chow interpolation is second-order accurate in space and time, as is confirmed by numerical experiments on small amplitude sloshing in a two-dimensional (2-D) tank, 3-D droplet oscillation and buoyant rise of a 3-D air bubble in water. Numerical results are found to be in excellent agreement with available theoretical and experimental results.

show abstract

Solution of the implicitly discretised fluid flow equations by operator-splitting

Cited by 4,474 publications

References 8 publications

The impact of dilatation, scrambling, and pressure transport in turbulent premixed flames

The impact of dilatation, scrambling, and pressure transport in turbulent premixed flames

Numerical Simulations of Sloshing and the Thermodynamic Response Due to Mixing

A moving mesh finite volume interface tracking method for surface tension dominated interfacial fluid flow

Contact Info

Product

Resources

About