Vacuum Tube
In this FSI test case, three pieces of flimsy placed on a table are sucked by a vacuuming tube. The sucking process and the interaction between the flimsy and the fluid is simulated including the deformation pattern of the flimsy.
*DUALCESE_MODEL
*DUALCESE_CONTROL_SOLVER
*DUALCESE_CONTROL_TIMESTEP
*DUALCESE_CONTROL_LIMITER
*DUALCESE_INCLUDE_MODEL
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET
*DUALCESE_INITIAL
*DUALCESE_PART
*DUALCESE_MAT_GAS
*DUALCESE_EOS_IDEAL_GAS
*DUALCESE_D3PLOT
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3DUMP
*INCLUDE
*DEFINE_CURVE
*PART
*SECTION_SHELL
*MAT_ELASTIC
*CONTACT_AUTOMATIC_GENERAL_MPP
*CONTROL_CONTACT
*SET_PART_LIST
*LOAD_BODY_Z
*LOAD_BODY_PARTS
*BOUNDARY_SPC_SET
*KEYWORD
*TITLE
3_Flimsy_parallel_1
$===========================================================================
$ Units:
$ rho (density): Kg/m^3
$ V (velocity): m/s
$ P (pressure): Pa
$ T (temperature): K
$============================================================================
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
0.07 60000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Define the load curves
*DEFINE_CURVE
$ lcid sidr sfa sfo offa offo dattyp
1
$ a1 o1
0.0 1.0
$ 1000000.0 1.0
10.0 1.0
*DEFINE_CURVE
$ lcid sidr sfa sfo offa offo dattyp
2
$ a1 o1
0.0 0.0
$ 1000000.0 0.0
10.0 0.0
*DEFINE_CURVE
$ lcid sidr sfa sfo offa offo dattyp
101
$ a1 o1
0.0 1.0
1000000.0 1.0
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state, time-history, and restart data
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
0.1e-2 0
$ 0.4e-3 0
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
30000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
struct_setup.k
*DUALCESE_MODEL
dcese_setup.k
*END
$===========================================================================
*KEYWORD
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup the CE/SE solver
*DUALCESE_CONTROL_SOLVER
$ ieqns igeom iframeref
Euler 3D fixed
$ NS 3D fixed
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2 0.9 0.0001
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
0 1.0 1.0 0.01
*DUALCESE_INCLUDE_MODEL
fluid_m5h.k
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_FSI_EXCLUDE
1005
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup the boundary conditions
$ Prescribed inlet BC (left)
*DUALCESE_BOUNDARY_PRESCRIBED_SEGMENT_SET
$ ssid
10
$ lcid_u lcid_v lcid_w lcid_d lcid_p lcid_t
-1 -1 -1 -1 -1
$ sf_u sf_v sf_w sf_d sf_p sf_t
$ 0.98325e+5
0.98325e+5
$
$ open BCs
$
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
$ ssid
1
2
3
5
12
$
$ bottom solid BC_7
$
*DUALCESE_BOUNDARY_SOLID_WALL_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 719.0 1006.0
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state, time-history, and restart data
*DUALCESE_D3PLOT
density
pressure
velocity
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END
$============================================================================
*KEYWORD
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*MAT_ELASTIC
5,2.4e+2,5.15e+8,0.1,0.0,0.0,0.0
*MAT_ELASTIC
9,5.0e+5,1.0e+9,0.1,0.0,0.0,0.0
*SECTION_SHELL
6,2,0.0,0.0,0.0
0.00035,0.00035,0.00035,0.00035
*SECTION_SHELL
7,16,0.0,0.0,0.0
0.002,0.002,0.002,0.002
*PART
flimsy
3,6,5,0,0,0,0,0
*PART
flimsy
4,6,5,0,0,0,0,0
*PART
flimsy
5,6,5,0,0,0,0,0
*PART
nozzle
501,7,9,0,0,0,0,0
*PART
table
1005,7,9,0,0,0,0,0
$===============================================================================
*CONTACT_AUTOMATIC_GENERAL_MPP
$# cid title
$# iptrack bsortfq unused track initier parmax unused beampart
0 1 0 10 0 1.100000
$# ssid msid sstyp mstyp sboxid mboxid spr mpr
910 0 2 2 1 1
$# fs fd dc vc vdc penchk bt dt
0.0 60.0 0 0.0 0.0
$# sfs sfm sst mst sfst sfmt fsf vsf
1E-6 1E-6
1
$===============================================================================
*CONTROL_CONTACT
0.10000 1 1
1 1
2
$===============================================================================
*SET_PART_LIST
910,0.0,0.0,0.0,0.0
3,4,5,501,1005
*SET_PART_LIST
922,0.0,0.0,0.0,0.0
3,4,5
*LOAD_BODY_Z
$ lcid sf lciddr xc yc zc
101,9.81
*LOAD_BODY_PARTS
922
*BOUNDARY_SPC_SET
79,0,1,1,1,1,1,1
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
struct_m1h.k
*END
$============================================================================