25.10.2014 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!