2D Pole Deformation
This is a 2D fluid-structure interaction (FSI) problem. In this test, a structure pole (shell element) stands in the middle of a flow field and its bottom end is fixed. When a strong wind flows from left to right, it pushes the pole and forces it to deform. It is simulated by the dual CESE FSI-mmm solver (this FSI solver is based on the moving mesh method (mmm), and can be more accurate than the immersed boundary method (ibm) in capturing small deformation problems).
*DUALCESE_MODEL
*DUALCESE_CONTROL_SOLVER
*DUALCESE_CONTROL_TIMESTEP
*DUALCESE_CONTROL_LIMITER
*DUALCESE_INCLUDE_MODEL
*DUALCESE_BOUNDARY_PRESCRIBED_SEGMENT_SET
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
*DUALCESE_INITIAL
*DUALCESE_PART
*DUALCESE_CONTROL_MESH_MOV
*DUALCESE_MAT_GAS
*DUALCESE_EOS_IDEAL_GAS
*DUALCESE_D3PLOT
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3DUMP
*INCLUDE
*PART
*SECTION_SHELL
*MAT_ELASTIC
*CONTROL_CONTACT
*BOUNDARY_SPC_SET
*CONTROL_STRUCTURED
*KEYWORD
*TITLE
shell_swing_in_flow
$===========================================================================
$ Units:
$ rho (density): Kg/m^3
$ V (velocity): m/s
$ P (pressure): Pa
$ T (temperature): K
$============================================================================
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
0.015 20000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state and restart data
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
1.0e-4 0
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
10000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
struct_setup.k
$
*DUALCESE_MODEL
fluid_setup.k
*END
$===========================================================================
*KEYWORD
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_CONTROL_SOLVER
$ ieqns igeom iframeref
Euler 2D fixed
$ NS 3D fixed
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2 0.9 0.1e-5
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
0 2.0 1.0 0.2
*DUALCESE_INCLUDE_MODEL
fluid_mesh1.k
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup the boundary conditions for fluid
$ Prescribed inlet BC (left)
*DUALCESE_BOUNDARY_PRESCRIBED_SEGMENT_SET
$ ssid
1
$ lcid_u lcid_v lcid_w lcid_d lcid_p lcid_t
-1
$ sf_u sf_v sf_w sf_d sf_p sf_t
1.0e+2 0.0 0.0 1.5 2.0e5
$
$ open BCs
$
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
$ ssid
2
$ 5
$
$ bottom solid BC_7
$
$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
$ ssid
3
5
*DUALCESE_BOUNDARY_FSI_SEGMENT_SET
$ ssid
11
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup the initial conditions for fluid
*DUALCESE_INITIAL
$ uic vic wic rhoic pic tic hic
0.0 0.0 0.0 1.225 1.01325e5
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup fluid properties
*DUALCESE_PART
$ pid mid eosid FSI-type MOVMSHALG
1 4 3 MMM 21
*DUALCESE_CONTROL_MESH_MOV
$ ID IALG NITER RELERRTOL MXDISP_R
21 9 1.5e-4
*DUALCESE_MAT_GAS
$ mid c1 c2 prnd
4 1.458e-6 110.4 0.72
*DUALCESE_EOS_IDEAL_GAS
$ eosid cv cp t0 v0
3 713.5 1001.5
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state and restart data
*DUALCESE_D3PLOT
density
pressure
velocity
vorticity
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END