Advanced : Irregular Waves
Fluid Velocity fringes
*KEYWORD *TITLE *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_DYNAMICS *CONTROL_IMPLICIT_SOLUTION *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DEFINE_CURVE *ELEMENT_SOLID
*ICFD_CONTROL_FSI
*ICFD_CONTROL_OUTPUT
*ICFD_CONTROL_MESH
*ICFD_CONTROL_TIME *ICFD_BOUNDARY_FREESLIP *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_FSWAVE *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME
*ICFD_INITIAL_LEVELSET *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION
*ICFD_SOLVER_TOL_LSET *ICFD_SOLVER_TOL_MOM *ICFD_SOLVER_TOL_PRE *INCLUDE *LOAD_BODY *MAT_ELASTIC *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *PART *SECTION_SOLID
*END
*KEYWORD
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_TIME
$# ttm dt
&T_end &dt_fluid
*ICFD_CONTROL_FSI
$# owc bt
0 &fsi_bt
*ICFD_CONTROL_OUTPUT
$# msglv
3
*ICFD_CONTROL_MESH
1.7
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD SOLVER TOLERANCES $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_SOLVER_TOL_MOM
1.e-10,1.e-10,1e3,1000
*ICFD_SOLVER_TOL_LSET
1.e-10,1.e-10,1e3,1000
*ICFD_SOLVER_TOL_PRE
1.e-10,1.e-10,1e3,1000
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_PART
$# pid secid mid
1 1 1
*ICFD_PART
$# pid secid mid
2 1 1
*ICFD_PART
$# pid secid mid
3 1 2
*ICFD_PART
$# pid secid mid
4 1 1
*ICFD_PART
$# pid secid mid
5 1 1
*ICFD_PART
$# pid secid mid
6 1 1
*ICFD_PART_VOL
10 1 1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1 2 3 4 5 6
*ICFD_MAT
$# mid flg ro vis
1 0
*ICFD_MAT
$# mid flg ro vis
2 1&rho_fluid &mu_fluid
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_INITIAL_LEVELSET
1 0 0 1 0 0 2
*ICFD_BOUNDARY_FREESLIP
1
*ICFD_BOUNDARY_FSWAVE
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ pid wtype h0 Hs wmin wmax Nmodes betadir
2 5 &H_0 &H_S &w_min &w_max &Nmodes &betadir
*ICFD_BOUNDARY_FREESLIP
3
*ICFD_BOUNDARY_NONSLIP
4
*ICFD_BOUNDARY_FSI
4
*ICFD_BOUNDARY_FREESLIP
5
*ICFD_BOUNDARY_PRESCRIBED_PRE
6,2
*DEFINE_CURVE
2 , 0, 1. , 1. , 0. , 0.
0.,0.
1e3,0.
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
20
1 2 3 4 5 6
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_DRAG
4
$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END
This example shows irregular (wind induced) ocean waves impacting a pillar using the wave B.C Generator. The irregular waves use the JONSWAP spectra. Freesurface is initialized directly via keyword rather than the classic way of building the initial interface mesh.