90 likes | 94 Vues
dr inż. Ziemowit Malecha TOWARDS A THERMAL MODEL FOR THE DIPOLE BY OpenFOAM. OpenFOAM. Open Source Field Operation and Manipulation C++ libraries ADVANTAGES: More than 50 standard solvers ( open source )
E N D
dr inż. Ziemowit MalechaTOWARDS A THERMAL MODEL FOR THE DIPOLE BY OpenFOAM
OpenFOAM • Open Source Field Operation and Manipulation C++ libraries • ADVANTAGES: • More than 50 standard solvers (open source) • Possibilityto implement new solver from the beginning or use standard solver as a frame • Access to source code • Applicable to any geometry DISADVANTAGE: • lack of full documentation
STANDARD SOLVERS • Incompressible flows RANS, LES, non Newtonian, PISO, SIMPLE • Compressible flows SIMPLE, PIMPLE, PISO, sonic/supersonic • Multiphase flows compressible, imcompressible, VOF (volume of fluid) • Combustion • Heat transfer and buoyancy-driven flows laminar/turbulent, compressible/incompressible • Molecular dynamics methods • Monte Carlo methods • Electromagnetics • Particle tracking flows • Stress analysis of solids • Finance
PREPROCESING • blockMesh • snappyHexMesh (hexahedral, stl geometry) • Transport mesh from commercial meshers • any mesh cell shape • Mesh, courtesy of: IT'IS Foundation for Research, Zürich, Switzerland
POSTPROCESING ParaView – Open Source Scientific Visualization
OpenFOAM syntax Navier – Stokes equation solve ( fvm::ddt(rho, U) + fvm::div(phi, U) - fvm::laplacian(mu, U) == - fvc::grad(p) ); • Linear solvers: • PCG/PBICG – preconditioned (bi-)conjugate gradient (DIC, FDIC, DILU, diagonal) • GAMG – generalised geometric-algebraic multi-grid
boundary condition fixedValueValue is specified fixedGradientNormal gradient is specified calculatedBoundary field derived from other fields mixedMixed fixedValue/ fixedGradient condition inletOutletSwitches and between fixedValue and zeroGradientdepending on direction of the flow pressureInletOutletVelocityCombination of pressureInletVelocity and inletOutlet Time and space varying
Two Fluid model ///PISO LOOP fvm::ddt(rho, U) + fvm::div(phi, U) == - fvc::grad(p) + fvc::grad(T) fvm::ddt(rho, Un) + fvm::div(phi, Un) == - fvc::grad(p) - fvc::grad(T) + laplacian(ni,Un) ///END PISO LOOP fvm::ddt(T) + fvm::div(phi, T) -laplacian(k, T) )
Fountain effect Superfluid normal fluid