Evaluation of Openfoam For Ship Hydrodynamics: Sung-Eun Kim

Download as pdf or txt
Download as pdf or txt
You are on page 1of 16

Evaluation of OpenFOAM for

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

3. Evaluation of free-surface capability “rasInterFoam”


ƒ “rasInterFoam” Based on VOF with URANS
ƒ The computation with a surface ship is underway.
4. Force/moment monitoring and reporting
ƒ “forceMoment” has been written and validated

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

You might also like