90 likes | 102 Views
Explore the capabilities of OpenFOAM, an open-source platform with over 50 standard solvers for fluid dynamics simulations. Benefit from access to source code, versatility in solver implementation, and applicability to various geometries. While lacking comprehensive documentation, OpenFOAM offers pre-processing tools like blockMesh and post-processing support through ParaView. The platform supports intricate simulations including incompressible and compressible flows, multiphase flows, combustion, heat transfer, and more. With linear solvers like PCG/PBICG and GAMG, OpenFOAM provides a powerful environment for conducting a wide range of computational analyses.
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