Embedded FSI example (DualCESE version)
This examples features a simple immersed FSI example. The active structure needs to be entirely embedded in the fluid for the immersed method. Interfaces are tracked automatically. It is advised to use a finer mesh for the fluid than for the solid. For FSI problems, the two solvers will use a synchronized timestep i.e the smallest value computed from their respective domains (CFL condition). No leakage can occur.
*DUALCESE_MODEL
*DUALCESE_CONTROL_SOLVER
*DUALCESE_CONTROL_TIMESTEP
*DUALCESE_CONTROL_LIMITER
*DUALCESE_INCLUDE_MODEL
*DUALCESE_FSI_EXCLUDE
*DUALCESE_BOUNDARY_PRESCRIBED_SEGMENT_SET
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
*DUALCESE_INITIAL
*DUALCESE_PART
*DUALCESE_MAT_GAS
*DUALCESE_EOS_IDEAL_GAS
*DUALCESE_D3PLOT
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3DUMP
*INCLUDE
*PART
*SECTION_SOLID
*SECTION_SHELL
*MAT_ELASTIC
*HOURGLASS
*CONTACT_AUTOMATIC_GENERAL_INTERIOR_MPP
*CONTROL_CONTACT
*SET_PART_LIST
*LOAD_BODY_Z
*LOAD_BODY_PARTS
*BOUNDARY_SPC_SET
*DEFINE_CURVE
*KEYWORD
*TITLE
fluid_solid_shell_interaction
$===========================================================================
$ Units:
$ rho (density): Kg/m^3
$ V (velocity): m/s
$ P (pressure): Pa
$ T (temperature): K
$============================================================================
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
0.03 20000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state and restart data
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
5.0e-4 0
$ 1.0e-7 0
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
1000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
struct_setup.k
*DUALCESE_MODEL
dcese_setup.k
*END
$===========================================================================
*KEYWORD
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_CONTROL_SOLVER
$ ieqns igeom iframeref
Euler 3D fixed
$ NS 3D fixed
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2 .99 0.1e-5
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
0 1.0 1.0 .1
*DUALCESE_INCLUDE_MODEL
dcese_m1h.k
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_FSI_EXCLUDE
10005
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ 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
- 0.0 0.0 1.22 1.51988e5
$ open BCs
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
$ ssid
2
3
5
10
$ bottom solid BC_7
$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
$ ssid
7
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup the initial conditions for fluid
*DUALCESE_INITIAL
$ uic vic wic rhoic pic tic hic
0.0 0.0 0.0 1.22 1.01325e5
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup fluid properties
*DUALCESE_PART
$ pid mid eosid FSI-type MOVMSHALG
1 4 3 IBM
*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
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END
$============================================================================
*KEYWORD
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*PART
solid
$ PID SECID MID EOSID HGID GRAV ADPOPT TMID
3 6 5 0 1
*SECTION_SOLID
6 1
*MAT_ELASTIC
5,1.4e+3,1.15e+10,0.1,0.0,0.0,0.0
*HOURGLASS
1 4
$===============================================================================
*PART
flimsy
15,17,19,0,0,0,0,0
*SECTION_SHELL
17,2,0.0,0.0,0.0
0.003,0.003,0.003,0.003
*MAT_ELASTIC
19,6.4e+2,5.15e+8,0.1,0.0,0.0,0.0
$===============================================================================
*PART
bottom_table
10005,7,9,0,0,0,0,0
*SECTION_SHELL
7,16,0.0,0.0,0.0
0.002,0.002,0.002,0.002
*MAT_ELASTIC
9,5.0e+3,1.0e+11,0.1,0.0,0.0,0.0
$===============================================================================
*CONTACT_AUTOMATIC_GENERAL_INTERIOR_MPP
$*CONTACT_AUTOMATIC_GENERAL_INTERIOR
$*CONTACT_AUTOMATIC_SINGLE_SURFACE
,10
910,0,2
0.0,0.0,0.0,0.0,50.0,2,0.0,0.0
1
$===============================================================================
*CONTROL_CONTACT
0.10000
2
$===============================================================================
*SET_PART_LIST
910,0.0,0.0,0.0,0.0
3,15,10005
*SET_PART_LIST
922,0.0,0.0,0.0,0.0
3,15
$===============================================================================
*LOAD_BODY_Z
$ lcid sf lciddr xc yc zc
101,9.81e-0
*LOAD_BODY_PARTS
922
*BOUNDARY_SPC_SET
79,0,1,1,1,1,1,1
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Define the load curves
*DEFINE_CURVE
$ lcid sidr sfa sfo offa offo dattyp
101
$ a1 o1
0.0 1.0
1000000.0 1.0
*DEFINE_CURVE
$ lcid sidr sfa sfo offa offo dattyp
11
$ a1 o1
0.0 0.98325e5
1000000.0 0.98325e5
*DEFINE_CURVE
$ lcid sidr sfa sfo offa offo dattyp
21
$ a1 o1
0.0 1.51988e5
1000000.0 1.51988e5
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
struct_m1h.k
*END
$============================================================================