|
|
(6 intermediate revisions by the same user not shown) |
Line 1: |
Line 1: |
− | == Two Dimensional Equations for fixed bodies in the time domain ==
| |
| | | |
− | We consider a two-dimensional fluid domain of constant depth, which
| + | {{coordinate definition in two dimension}} |
− | contains a finite number of fixed bodies of arbitrary geometry. We
| + | |
− | denote the fluid domain by <math>\Omega</math>, the boundary of the fluid domain
| + | {{equations of motion time domain without body condition}} |
− | which touches the fixed bodies by <math>\partial\Omega</math>, and the free
| + | |
− | surface by <math>F.</math> The <math>x</math> and <math>z</math> coordinates are such that <math>x</math> is
| + | The body boundary condition for a fixed body is |
− | pointing in the horizontal direction and <math>z</math> is pointing in the
| |
− | vertical upwards direction (we denote <math>\mathbf{x}=\left( x,z\right) ).</math> The
| |
− | free surface is at <math>z=0</math> and the sea floor is at <math>z=-h</math> (the equations
| |
− | would be almost identical if the sea floor depth varied). The equations
| |
− | of motion in the time domain are | |
− | <center><math>
| |
− | \Delta\Phi\left( \mathbf{x,}t\right) =0,\ \ \mathbf{x}\in\Omega,
| |
− | </math></center>
| |
| <center><math> | | <center><math> |
| \partial_{n}\Phi=0,\ \ \mathbf{x}\in\partial\Omega, | | \partial_{n}\Phi=0,\ \ \mathbf{x}\in\partial\Omega, |
| </math></center> | | </math></center> |
− | <center><math>
| + | |
− | \partial_{n}\Phi=0,\ \ z=-h,
| + | |
− | </math></center>
| + | The initial conditions are |
− | where <math>\Phi</math> is the velocity potential for the fluid. At the free
| + | {{initial free surface time domain}} |
− | surface we have the kinematic condition
| |
− | <center><math>
| |
− | \partial_{t}\zeta=\partial_{n}\Phi,\ \ z=0,\ x\in F,
| |
− | </math></center>
| |
− | and the dynamic condition (the linearized Bernoulli equation)
| |
− | <center><math>
| |
− | \zeta = -\partial_{t}\Phi,\ \ z=0,\ x\in F,
| |
− | </math></center>
| |
− | where <math>\zeta</math> is the free-surface elevation. These equations are in non-dimensional
| |
− | form (so that the fluid density and gravity are both unity). They are
| |
− | also subject to initial conditions
| |
− | <center><math>
| |
− | \left.\zeta\right|_{t=0} = \zeta_0(x)\,\,\,\mathrm{and}\,\,\,
| |
− | \left.\partial_t\zeta\right|_{t=0} = \partial_t\zeta_0(x).
| |
− | </math></center>
| |
We consider a two-dimensional fluid domain of constant depth, which
contains a finite number of fixed bodies of arbitrary geometry. We
denote the fluid domain by [math]\displaystyle{ \Omega }[/math], the boundary of the fluid domain
which touches the fixed bodies by [math]\displaystyle{ \partial\Omega }[/math], and the free
surface by [math]\displaystyle{ F. }[/math] The [math]\displaystyle{ x }[/math] and [math]\displaystyle{ z }[/math] coordinates are such that [math]\displaystyle{ x }[/math] is
pointing in the horizontal direction and [math]\displaystyle{ z }[/math] is pointing in the
vertical upwards direction (we denote [math]\displaystyle{ \mathbf{x}=\left( x,z\right) ). }[/math] The
free surface is at [math]\displaystyle{ z=0 }[/math] and the sea floor is at [math]\displaystyle{ z=-h }[/math]. The
fluid motion is described by a velocity potential [math]\displaystyle{ \Phi }[/math] and free surface by
[math]\displaystyle{ \zeta }[/math].
The equations
of motion in the time domain are
Laplace's equation through out the fluid
[math]\displaystyle{
\Delta\Phi\left( \mathbf{x,}t\right) =0,\ \ \mathbf{x}\in\Omega.
}[/math]
At the bottom surface we have no flow
[math]\displaystyle{
\partial_{n}\Phi=0,\ \ z=-h.
}[/math]
At the free
surface we have the kinematic condition
[math]\displaystyle{
\partial_{t}\zeta=\partial_{n}\Phi,\ \ z=0,\ x\in \partial\Omega_{F},
}[/math]
and the dynamic condition (the linearized Bernoulli equation)
[math]\displaystyle{
\partial_{t}\Phi = -g\zeta ,\ \ z=0,\ x\in \partial\Omega_{F}.
}[/math]
The body boundary condition for a fixed body is
[math]\displaystyle{
\partial_{n}\Phi=0,\ \ \mathbf{x}\in\partial\Omega,
}[/math]
The initial conditions are
[math]\displaystyle{
\left.\zeta\right|_{t=0} = \zeta_0(x)\,\,\,\mathrm{and}\,\,\,
\left.\partial_t\zeta\right|_{t=0} = \partial_t\zeta_0(x).
}[/math]