foamedOver - A Library to Add a Dynamic Overset Grid Capability to ...
foamedOver - A Library to Add a Dynamic Overset Grid Capability to ...
foamedOver - A Library to Add a Dynamic Overset Grid Capability to ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
<strong>foamedOver</strong><br />
A <strong>Library</strong> <strong>to</strong> <strong>Add</strong> a <strong>Dynamic</strong> <strong>Overset</strong> <strong>Grid</strong> <strong>Capability</strong> <strong>to</strong> OpenFOAM<br />
David Boger Ralph Noack Eric Paterson<br />
Penn State University, Applied Research Lab<br />
State College, PA USA<br />
5th OpenFOAM Workshop<br />
21 - 24 June 2010<br />
Chalmers University of Technology<br />
Göteborg, Sweden<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 1 / 25
Outline<br />
1 Background<br />
<strong>Overset</strong> <strong>Grid</strong>s<br />
<strong>Overset</strong> <strong>Grid</strong> Approach<br />
Suggar++<br />
DiRTlib<br />
2 <strong>foamedOver</strong><br />
3 Results<br />
Standard Applications<br />
<strong>Dynamic</strong> Meshes<br />
4 Summary<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 2 / 25
Definition<br />
An overset grid assembly is<br />
comprised of a set of grids<br />
(structured, unstructured, or mixed)<br />
The grids overlap each other<br />
arbitrarily <strong>to</strong> cover the domain<br />
Component grids are constructed<br />
(mostly) independent of one another<br />
Field equations are connected by<br />
interpolation across the overlapping<br />
boundaries<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 3 / 25
Advantages<br />
<strong>Overset</strong> grids are a major enabling<br />
technology for<br />
◮ Complex geometries<br />
⋆ Component grids can be built<br />
independently<br />
⋆ Components can be easily<br />
modified, added, or removed<br />
◮ Bodies in relative motion<br />
◮ Tracking of important flow features<br />
⋆ <strong>Grid</strong> points are added only where<br />
they are needed<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 4 / 25
Advantages<br />
<strong>Overset</strong> grids are a major enabling<br />
technology for<br />
◮ Complex geometries<br />
⋆ Component grids can be built<br />
independently<br />
⋆ Components can be easily<br />
modified, added, or removed<br />
◮ Bodies in relative motion<br />
◮ Tracking of important flow features<br />
⋆ <strong>Grid</strong> points are added only where<br />
they are needed<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 4 / 25
Advantages<br />
<strong>Overset</strong> grids are a major enabling<br />
technology for<br />
◮ Complex geometries<br />
⋆ Component grids can be built<br />
independently<br />
⋆ Components can be easily<br />
modified, added, or removed<br />
◮ Bodies in relative motion<br />
◮ Tracking of important flow features<br />
⋆ <strong>Grid</strong> points are added only where<br />
they are needed<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 4 / 25
Advantages<br />
<strong>Overset</strong> grids are a major enabling<br />
technology for<br />
◮ Complex geometries<br />
⋆ Component grids can be built<br />
independently<br />
⋆ Components can be easily<br />
modified, added, or removed<br />
◮ Bodies in relative motion<br />
◮ Tracking of important flow features<br />
⋆ <strong>Grid</strong> points are added only where<br />
they are needed<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 4 / 25
<strong>Overset</strong> <strong>Grid</strong> Approach<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 5 / 25
<strong>Overset</strong> Hole Cutting<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 6 / 25
<strong>Overset</strong> Fringe Boundary Points<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 7 / 25
<strong>Overset</strong> Overlap Minimization<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 8 / 25
Requirements<br />
Obtain the Domain Connectivity Information (DCI)<br />
◮ Determine in/out/fringe status and interpolation stencils<br />
Read, access, and utilize the DCI<br />
◮ Perform interpolation of the dependent variables<br />
Modify the solver <strong>to</strong> exclude out cells<br />
◮ Equation solvers, residual measures, volume averaging, etc.<br />
Modify the field equations at fringe points<br />
◮ Replace governing equations by interpolation<br />
Move the grids<br />
Fix force and moment integrals on overlapping surfaces<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 9 / 25
Suggar++ 1<br />
<strong>Overset</strong> grid assembly software<br />
Performs hole-cutting, donor search,<br />
overlap minimization<br />
Static or dynamic (moving body)<br />
assemblies<br />
Structured and unstructured grids<br />
Node-centered and cell-centered flow<br />
solvers<br />
Stand-alone executable (static) or<br />
library calls for dynamic grids<br />
(libsuggar.so)<br />
1 Ralph W. Noack, David A. Boger, Robert F. Kunz, and Pablo M. Carrica, "Suggar++: An Improved General <strong>Overset</strong> <strong>Grid</strong><br />
Assembly <strong>Capability</strong>," AIAA Paper 2009-3992<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 10 / 25
DiRTlib 2<br />
DiRTlib is a solver-neutral library that simplifies the addition of an<br />
overset capability <strong>to</strong> a flow solver by encapsulating the required<br />
operations<br />
Acquire interpolation stencils via file<br />
I/O or direct communication with<br />
libSuggar++<br />
Provides higher-level methods <strong>to</strong><br />
transfer field data from donors <strong>to</strong><br />
recep<strong>to</strong>rs and interpolate<br />
Provides lower-level access <strong>to</strong> donor<br />
member indices and weights <strong>to</strong> help<br />
build implicit global matrix<br />
drt_get_dci()<br />
drt_generate_transmit_apply()<br />
drt_get_donor_members_donor_weights()<br />
2 Ralph W. Noack, “DiRTlib: A <strong>Library</strong> <strong>to</strong> <strong>Add</strong> an <strong>Overset</strong> <strong>Capability</strong> <strong>to</strong> Your Flow Solver” AIAA 2005-5116<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 11 / 25
<strong>foamedOver</strong><br />
FoamedOver is an interface between OpenFOAM and other<br />
specialized libraries<br />
DiRTlib – <strong>Library</strong> <strong>to</strong> facilitate<br />
the addition of overset<br />
capability <strong>to</strong> any flow solver<br />
Suggar++ – <strong>Overset</strong> grid<br />
assembly software<br />
PETSc – <strong>Library</strong> of data<br />
structures and routines for the<br />
parallel solution of large<br />
systems of linear and<br />
non-linear equations<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 12 / 25
<strong>foamedOver</strong><br />
FoamedOver is a collection of cus<strong>to</strong>m classes, solvers, and<br />
applications.<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 13 / 25
<strong>foamedOver</strong><br />
FoamedOver is a stand-alone library that provides a dynamic overset<br />
grid capability <strong>to</strong> any OpenFOAM solver<br />
No changes are required <strong>to</strong><br />
the OpenFOAM library itself<br />
Any OpenFOAM solver is<br />
made overset-capable by the<br />
insertion of a half-dozen lines<br />
of code . . .<br />
. . . and the use of cus<strong>to</strong>m<br />
run-time selectable objects<br />
# include “oversetObject.H”<br />
# include “create<strong>Overset</strong>.H”<br />
overset.updateFringeValues(U);<br />
dynamicFvMesh oversetFvMesh;<br />
mesh.update()<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 14 / 25
Run–Time Selectable Objects<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 15 / 25
Potential Foam<br />
cylinder<br />
steady<br />
potential flow<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 16 / 25
interFoam<br />
submerged hydrofoil a<br />
steady<br />
incompressible multiphase<br />
a J. Duncan, "The Breaking and Non-Breaking Wave<br />
Resistance of a Two-Dimensional Hydrofoil," J. Fluid Mech,<br />
126:507–520, 1983.<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 17 / 25
interFoam<br />
extension of damBreak tu<strong>to</strong>rial<br />
unsteady<br />
incompressible multiphase<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 18 / 25
compressibleInterFoam<br />
extension of depthCharge<br />
tu<strong>to</strong>rial<br />
unsteady<br />
compressible multiphase<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 19 / 25
Unsteady potential flow<br />
cylinder<br />
prescribed mesh motion<br />
unsteady<br />
unsteady potentail flow is a<br />
series of steady solutions<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 20 / 25
icoDyMFoam<br />
pitching foil<br />
prescribed mesh motion<br />
incompressible, laminar flow<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 21 / 25
interDyMFoam<br />
falling cylinder<br />
6DOF mesh motion<br />
incompressible, laminar flow<br />
Compare <strong>to</strong> simplified<br />
analytical solution for<br />
free-falling body with constant<br />
CD<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 22 / 25
Two Paddles<br />
interlacing rotating paddles<br />
motion would be very difficult<br />
<strong>to</strong> resolve with GGI, RBF, or<br />
Laplacian dynamicFvMesh<br />
methods.<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 23 / 25
Water Entry<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 24 / 25
Summary<br />
<strong>foamedOver</strong> is a collection of cus<strong>to</strong>m classes, solvers, and<br />
applications<br />
It is a stand-alone library that provides a dynamic overset grid<br />
capability <strong>to</strong> any OpenFOAM solver<br />
motionObjects have been developed which permit prescribed,<br />
analytical, and 6DOF motions<br />
For now, it is controlled by ITAR 3 , and not available outside the<br />
USA<br />
Future work<br />
◮ Develop interface with native OpenFOAM algebraic solvers.<br />
◮ Parallel efficiency<br />
3 International Traffic in Arms Regulations<br />
Boger, Noack, and Paterson (PSU/ARL) <strong>foamedOver</strong> OFW5 25 / 25