OpenFOAM Tutorials: Basic Session

Hrvoje Jasak
[email protected]

Wikki Ltd, United Kingdom and FSB, University of Zagreb, Croatia 7-9th June 2007

Open FOAM

OpenFOAM Tutorials: Basic Session p.1/19

Open Source CFD Platform

OpenFOAM: Open Source Computational Continuum Mechanics Commercial CFD vendors do not provide exibility for customisation and add-on developments to answer the needs of commercial and research users Proprietary approach to solution methodology is the limiting factor: closed software architecture; no customer-specic proprietary models etc. Complete CFD methodology is already in the public domain (research papers, model formulation, numerical schemes, linear equation solvers etc.) The way forward involves educated users and customised simulations Objective: open source implementation of existing knowledge and an object-oriented platform for easy future development 1. Completely open software platform using object-oriented design 2. Extensive modelling capabilities in library form: component re-use 3. Fast, robust and accurate numerical solver 4. State of the art complex geometry handling 5. On-order targeted and customer-driven model development . . . but the mode of operation changes considerably

Open FOAM

OpenFOAM Tutorials: Basic Session p.2/19

Implementing Continuum Models

How to Handle Complex Continuum Models in Software? Natural language of continuum mechanics: partial differential equations Example: turbulence kinetic energy equation k + (uk) [( + t )k] = t t 1 (u + u T ) 2 2 o k ko

Objective: represent differential equations in their natural language

solve ( fvm::ddt(k) + fvm::div(phi, k) - fvm::laplacian(nu() + nut, k) == nut*magSqr(symm(fvc::grad(U))) - fvm::Sp(epsilon/k, k) );

Correspondence between the implementation and the original equation is clear

Open FOAM

OpenFOAM Tutorials: Basic Session p.3/19

Object Orientation
Object-Oriented Software Analysis of numerical simulation software through object orientation: Recognise main objects from the numerical modelling viewpoint Objects consist of data they encapsulate and functions which operate on the data Example: Sparse Matrix Class Data members Sparse addressing pattern (CR format, arrow format) Diagonal coefcients, off-diagonal coefcients Operations on matrices or data members Matrix algebra operations: +, , , /, Matrix-vector product, transpose, triple product, under-relaxation Actual data layout and functionality is important only internally: efciency Example: Linear Equation Solver Operate on a system of linear equations Ax = b to obtain x It is irrelevant how the matrix was assembled or what shall be done with solution Ultimately, even the solver algorithm is not of interest: all we want is new x! Gauss-Seidel, AMG, direct solver: all answer to the same interface

Open FOAM

OpenFOAM Tutorials: Basic Session p.4/19

Object Orientation
Basic Components Scalars, vectors and tensors with algebra Computational mesh; mesh motion, adaptive renement, topological changes Fields (scalar, vector, tensor) and boundary conditions: Dirichlet, Neumann etc. Sparse matrix support with linear solver technology Discretisation Classes Implemented as interpolation, differentiation and discretisation operators All discretisation methods use identical basic components, e.g. common mesh and matrix support. Better testing and more compact software implementation Physical Modelling Libraries and Top-Level Solvers Libraries encapsulate interchangeable models with common interfaces New models provide functionality adhering to common interface Custom-written and optimised top-level solvers for class of physics Utilities Common functionality is needed for most simulation methods Example: problem setup, mesh manipulation, data analysis

Open FOAM

OpenFOAM Tutorials: Basic Session p.5/19

Tutorials
Running OpenFOAM 1. Laminar uid ow in a lid-driven cavity Pre-processing with FoamX Case structure and dictionary format Running OpenFOAM Basic post-processing Create a graded mesh, map a solution and repeat a simulation 2. Free surface ow: dam break problem Pre-processing using le editors and command line Multi-block mesh generation with blockMesh Setting non-uniform initial eld Post-processing: create an animation Running OpenFOAM on multiple processors

Open FOAM

OpenFOAM Tutorials: Basic Session p.6/19

Lid-Driven Cavity
Flow in a Lid-Driven Cavity: Case Setup

Ux = 1 m/s

d = 0.1 m y x
Open FOAM

OpenFOAM Tutorials: Basic Session p.7/19

Lid-Driven Cavity
Mesh Generation: blockMesh

3 7

2 6

y
0

x z
4

1 5
OpenFOAM Tutorials: Basic Session p.8/19

Open FOAM

Lid-Driven Cavity
Dimensions and Properties Domain size: 0.1 0.1 m Lid velocity 1 m/s Kinematic viscosity: 0.01 m2 /s Simulation Settings Start time = 0 End time = 0.5 s Time step: t = 0.005 s Sampling Lines (0.0500.0005) to (0.050.10.0005) (0.0500.0005) to (0.050.10.0005)

Open FOAM

OpenFOAM Tutorials: Basic Session p.9/19

Lid-Driven Cavity
Phases of Simulation Mesh generation Boundary conditions Discretisation control Running the simulation Graphical post-processing Velocity magnitude Sampling lines Mesh Renement Mesh generation Solution mapping Running the simulation Graphical post-processing

Open FOAM

OpenFOAM Tutorials: Basic Session p.10/19

Case Data Format and Organisation

Structure of an OpenFOAM Case Unlike other applications, mesh, case setup and results data in OpenFOAM is split into a number of les The case is stored in a directory, with a prescribed structure. A case consists of: Mesh or series of meshes (deformation, topological changes) Material properties Choice of models Solver settings: discretisation choice, convergence tolerances, etc. Results of simulation, organised in time or iteration directories Idea: keep similar cases together: root directory Case data separated into units: case directory Typically, 1 high-level object lives in 1 le Mesh separated into several les: point motion/topological change

Open FOAM

OpenFOAM Tutorials: Basic Session p.11/19

Case Data Format and Organisation

Data File Format All OpenFOAM les annotated by a version-control and context header Header contains File version, format and location (root/cas/instance) Class type and object name of the class the data belongs to C++ comment stype applies to all parsed data
FoamFile { version format root case instance local class object

2.0; ascii; ""; ""; ""; ""; dictionary; transportProperties;

Open FOAM

OpenFOAM Tutorials: Basic Session p.12/19

Case Data Format and Organisation

Dictionary Format Main input system: dictionary List of keyword value pairs delimited by semicolon Order of entries is irrelevant Embedded dictionary levels allowed Field data also in dictionary format: uniform or expanded data, using supertokens for efcient reading of expanded data Embedded dictionaries use curly brackets delimitation Main dictionary, controlling execution: system/controlDict All other dictionaries tailored to applications Example: transport properties for icoFoam
nu nu [0 2 -1 0 0 0 0] 0.01;

Open FOAM

OpenFOAM Tutorials: Basic Session p.13/19

Case Data Format and Organisation

Dictionary Format: Field Data Example
dimensions [0 1 -1 internalField uniform boundaryField { movingWall { type value } fixedWalls { type value } frontAndBack { type } } 0 0 0 0]; (0 0 0);

fixedValue; uniform (1 0 0);

fixedValue; uniform (0 0 0);

empty;

Open FOAM

OpenFOAM Tutorials: Basic Session p.14/19

Case Data Format and Organisation

OpenFOAM Case Structure

<case>
system controlDict fvSchemes fvSolution constant . . . Properties polyMesh points cells faces boundary time directories

Open FOAM

OpenFOAM Tutorials: Basic Session p.15/19

Dam Break Free Surface Flow

0.584 m

0.292 m

1111 0000 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111 0000 1111
0.146 m

water column

0.584 m

0.048 m 0.024 m
OpenFOAM Tutorials: Basic Session p.16/19

Open FOAM

Dam Break Free Surface Flow

Material Properties 1 = 1000 kg/m3 2 = 1 kg/m3 1 = 1 106 Pa/s 2 = 1.48 106 Pa/s Gravity vector (0, 9.81, 0) m/s2 Initial eld Column of water, 0.292 m high and 0.1461 m wide Both uids at rest at t = 0 s

Open FOAM

OpenFOAM Tutorials: Basic Session p.17/19

OpenFOAM: Implemented Capabilities

Discretisation Methods and Support Second and fourth-order Finite Volume with mesh motion and topological changes Polyhedral Finite Element solver (mesh motion) Lagrangian particle tracking (discrete particle model); Diesel spray model Finite Area Method: FVM on a curved surface in 3-D A-posteriori error estimation Dynamic mesh handling and topology changes; automatic mesh motion Physical Modelling Libraries Thermo-physical models (liquids and gasses) Chemical reaction library interface (Chemkin) Non-Newtonian viscosity models Turbulence models (RANS and LES, compressible and incompressible); DNS Diesel spray (atomisation, dispersion, heat transfer, evaporation, spray-wall etc.) High Performance Computing Support Massively parallel computing: domain decomposition approach Next-generation of linear equation solver technology: up to 3 times faster!

Open FOAM

OpenFOAM Tutorials: Basic Session p.18/19

OpenFOAM: Implemented Capabilities

Standard Top-Level Solvers Basic: Laplace, potential ow, transport Incompressible ow, compressible ow Heat transfer: buoyancy-driven ows Multiphase: Euler-Euler, surface capturing and surface tracking DNS and LES turbulent ows, aero-acoustics Pre-mixed and Diesel combustion, spray and in-cylinder ows Stress analysis, uid-structure interaction, electromagnetics, MHD, etc. Utilities Pre-processing, data manipulation Mesh import and export, mesh generation and manipulation Parallel processing tools: decomposition and reconstruction Post processor hook-up (reader module) and data export This is just a standard set: users write their own applications using the library

Open FOAM

OpenFOAM Tutorials: Basic Session p.19/19