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