This page gives a short introduction of panMARE's theoretical background. Thereby, panMARE is based on the principles and equations of the potential method described in detail by Katz and Plotkin in Low-Speed Aerodynamics. More information can also be found in the Institute for Fluid Dynamics and Ship Theory's various panMARE-related publications.
Potential Theory
The panel method panMARE uses potential theory to calculate the flow.
A shematic representation of the problem to be solved is visualised in Figure 1, including
the boundaries and used nomenclature. One or more bodies with the surface
are located in the fluid , which is enclosed by a domain boundary . It is
possible to further include boundaries such as walls, a free surface ,
or the seabed . These boundaries are counted among the fluid domain's
surrounding outer boundaries. The surface normal points into the fluid domain .
The body (or bodies) move(s) through the fluid domain, whereby translations, rotations and
deformations are possible in arbitrary combinations. The velocities resulting from the body-movement
are known at every point on the surface as . For lift-generating bodies an additional
wake surface is introduced, which prevents the flow around the hydrofoil's trailing edge, thus
satisfying the Kutta boundary condition.
Figure 1: Boundaries in panMARE
Since the flow is assumed to be inviscid, irrotational and incompressible in potential theory,
the conservation equations can be simplified resulting in the incompressible and irrotational
equation of continuity (LaPlace equation) formulated for the velocity potential :
and the Bernoulli equation for the pressure :
Hereby, is the potential, whose spatial gradient corresponds to the flow velocity,
whereas is the gravitational acceleration, the observation point's position and
the flow velocity. Assuming that there is no flow through the body (boundaries and )
and that the influence of the body vanishes at infinity (boundary ), a
Green's identity can be used as a solution of the LaPlace equation:
The fluid domain's boundary consists of all boundaries , and of the fluid .
The identity equation above can be rearranged and the following subsitution of a doublet strength and a source strength applied:
Hereby, the doublet strength corresponds to the jump of the velocity potential from
within the body to outside the body in normal direction , and the source strength
corresponds to the jump of the velocity potential's gradient, both at the boundary . The normal direction is positive pointing away from the body, i.e. into the fluid.
Based on this, the velocity potential at a point in the fluid can be described by the following equation:
Hereby, denotes the distance between the point of observation and the surface. For the wake surface
there is no jump in the velocity potential's gradient in normal direction, since the wake surface is not a
body. The fluid's perturbation caused by the body is represented by
the potential. The undisturbed potential describes the condition of the fluid without the
body's presence and includes possible currents as well as the influence of the distant boundary .
Since the distant boundary has no contribution and its influence vanishes at infinity, the undisturbed
potential is equal to the velocity potential of the underlying currents,
which is defined continuously throughout the fluid domain. Assuming that:
and substituting with in the Green's identity, the velocity potential at a point
in the fluid can be written as follows:
Thereby, the influence of the body is represented by the so-called induced potential .
The spatial gradient of the induced potential corresponds to the induced velocity of the flow.
Boundary Conditions
To determine the doublet strength and the source strength boundary conditions are required.
Since the induced potential vanishes inside the body, the inner potential is defined as:
This leads to the following formulation of the Dirichlet problem:
This equation is evaluated at the inner side of the body's boundary . If the equation is evaluated at
the distant boundary , the terms in the integrals vanish, so that the boundary condition is also fulfilled there.
Another boundary condition is that there must be no flow through the body (Neumann boundary condition).
To satisfy the Neumann boundary condition the velocities in normal direction must vanish.
The sum of the induced and the external potential's gradients and the body's motion velocity is calculated and multiplied with the normal vector. The result is set to zero.
This leads to the following relation:
The contribution of the motion velocity is included as negative, since the body is moved throught the resting fluid. As explained above, the source strength is defined as the jump of the potential's gradient. Since both the gradient as well as the induced potential vanish inside the body, the source strength is: . Thus the source strength can be written as follows:
Hydrofoil Wake Surface
For the calculation of a hydrofoil a wake surface is is inserted at the trailing edge of the hydrofoil, see Figure 1, in order
to prevent the fluid from flowing around the trailing edge. This wake surface can be viewed as the free vortex sheet. The Kutta
boundary condition demands that the pressure difference between suction and pressure side vanishes at the trailing edge. Linearisation
leads to the following Morino-Kutta condition:
This equation establishes a connection between the doublet strengths on the upper and lower side of the trailing edge and the wake surface.
As long as the flow is mainly perpendicular to the trailing edge, the linearisation does not lead to significant mistakes. The orientation
of the free vortex sheet is parallel to the respective flow lines, since there is no flow through this surface.