Evaluation of Openfoam For Ship Hydrodynamics: Sung-Eun Kim
Evaluation of Openfoam For Ship Hydrodynamics: Sung-Eun Kim
Evaluation of Openfoam For Ship Hydrodynamics: Sung-Eun Kim
Ship Hydrodynamics
Sung-Eun Kim
CFD Group, NSWCCD, U.S.A.
June 8, 2007
Outline
• Motivation
• Pilot Studies Using OpenFOAM
• Summary & Conclusions
2
Motivation
• Ideal framework for in-house CFD capability
development
• Collaboration and knowledge sharing
• Rapid prototyping and technology transfer
• Hands-on learning in modern software
engineering applied to CFD
3
Outline
• Motivation
• Pilot Studies Using OpenFOAM
• Summary & Conclusions
4
Pilot Studies - Objectives
• Usability
Learning curve
User-friendliness
• Programmability
Impact of language (C++) language barrier
Tasks for different levels of programming skills
• Extensibility – code architecture & design
Solver applications
Classes and libraries
• Solver performance – benchmarking against
commercial and in-house codes
Speed
Accuracy
Robustness
5
Pilot Studies Using OpenFOAM
Tasks
1. Implementation of a Wilcox’ 1998 k-ω turbulence
model
A new derived class with runtime selection
Wall boundary conditions for k and ω (2-layer wall functions)
User-interface for model constants
Validation for a body of revolution at incidence
2. RANS solver in a rotating frame of reference
1. Source terms
2. Transformation of boundary conditions in rotating coordinates
3. Validation for an open-water propeller
6
Pilot Development Projects
7
Task1 1 - Validation with a BOR
• Body of revolution
ReL = 11.7 x 106, L/D = 7.3
• Lift and moment measured for -30° < α < 30°
Earlier RANS computations
• Rhee and Sung (2004)
• Wilson et al. (2004)
1.3M-cell structured mesh – the same mesh as the one
used with FLUENT
• Y+ < 0.7 at wall-adjacent cells
• Minimum 15 grid points in BL
QUICK scheme was used for for convection discretization.
8
BOR - Pressure Distribution (ParaFOAM)
9
BOR – Pathlines (FLUENT)
10
BOR - Normal Force Coefficient
vs. Angle of Attack
11
BOR -Normal Force Coeff.
at α = 18°
12
BOR - Moment Coefficient
vs. Angle of Attack
13
Task 2 – Validation for a Marine
Propeller
• Computational mesh
450K-cell tet. mesh
Coarse near-wall (wall function) mesh
• FLUENT runs are complete with realizable k-ε and
Wilcox’ k-ω model
• Solved in non-inertial (rotating) frame of reference
• OpenFOAM computations are in progress
Run into convergence difficulty
• Possibly due to strong body-force terms
The solver (sransFoam) and the boundary/initial conditions
14
Task 3 – rasInterFoam Validation
A Surface Ship at Fr = 0.28
• Computational mesh
1.8M-cell structured mesh
• FLUENT has been run on the
same case.
VOF with modified HRIC
• OpenFOAM computations
Mesh translation and problem-
setup completed
The computation in progress
15
Conclusions
• Usability
User-friendliness – satisfactory
• Programmability
A wide range of development (coding) tasks can be
accomplished without deep knowledge in C++
Learning curve – fairly reasonable with the precompiled solver
applications
C++ knowledge is obviously a great plus and opens up many
possibilities.
Greatly aided by “learn-by-examples”
• Extensibility
Great!
• Solver performance
Slightly less robust and less efficient than commercial codes
Comparable accuracy
16