This is a 3D fluid-structure interaction (FSI) problem solved with the dual CESE FSI solver using the immersed boundary method (ibm). In this test, a folded bag is deployed by a high-pressure jet from the entry tube.
*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_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
*PART
*SECTION_SHELL
*MAT_FABRIC
*MAT_ELASTIC
*CONTACT_AIRBAG_SINGLE_SURFACE
*CONTROL_CONTACT
*SET_PART_LIST
*LOAD_BODY_Z
*LOAD_BODY_PARTS
*BOUNDARY_SPC_SET
*DEFINE_CURVE
*KEYWORD
*TITLE
fold_bag_deployment_test4
$
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
0.05 20000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state and restart data
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
1.0e-3 0
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
5000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
bag_setup4.k
*DUALCESE_MODEL
fluid_setup4.k
*END
$=================================================================================
*KEYWORD
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_CONTROL_SOLVER
$ ieqns igeom iframeref
Euler 3D fixed
$ NS 3D fixed
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2 0.99 1.e-4
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
0 2.0 1.0 0.05
*DUALCESE_INCLUDE_MODEL
fluid_mesh4.k
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Setup the boundary conditions for fluid
$ Prescribed inlet BC (bottom hole)
*DUALCESE_BOUNDARY_PRESCRIBED_SEGMENT_SET
$ ssid
10
$ 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
0.0 0.0 2.0e+2 2.5 3.03975e5
$
$ surrounding open BCs
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
$ ssid
1
2
3
5
7
12
*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET
$ ssid
15
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ 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 modeleqs FSI-type MOVMSHALG
$ 1 7 5
$ pid mid eosid FSI-type MOVMSHALG
1 7 5 IBM
*DUALCESE_MAT_GAS
$ mid c1 c2 prnd
7 1.458e-6 110.4 0.72
*DUALCESE_EOS_IDEAL_GAS
$ eosid cv cp t0 v0
5 717.5 1004.5
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ Handle output of state and restart data
*DUALCESE_D3PLOT
density
pressure
velocity
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END
$================================================================================
*KEYWORD
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*PART
bag
503,6,5,0,0,0,0,0
*SECTION_SHELL
6,2,0.0,0.0,0.0
0.00035,0.00035,0.00035,0.00035
*MAT_FABRIC
$# mid ro ea eb ec prba prca prcb
5 647.00000 2.8100E+8 2.8100E+8 2.8100E+8 0.300000 0.300000 0.300000
$# gab gbc gca cse el prl lratio damp
0.000 0.000 0.000 1.000000 0.000 0.000 0.100000 0.050000
$# aopt flc fac ela lnrc form fvopt tsrfac
$ 3.000000 0.000 1.000000 0.000 1.000000 0 0.000 0.000
3.000000 0.000 0.000 0.000 0.000000 0 0.000 0.000
$# xp yp zp a1 a2 a3
0.000 0.000 0.000 0.000 0.000 0.000
$# v1 v2 v3 d1 d2 d3 beta
1.000000 1.000000 1.000000 0.000 0.000 0.000 0.000
$——————————————————————————-
*PART
dashpanel
105,7,9,0,0,0,0,0
*SECTION_SHELL
7,2,0.0,0.0,0.0
0.00035,0.00035,0.00035,0.00035
*MAT_ELASTIC
9,1.2e+3,1.0e+9,0.1,0.0,0.0,0.0
$——————————————————————————-
*CONTACT_AIRBAG_SINGLE_SURFACE
$# cid title
$ SSID MSID SSTYP MSTYP SBOXID MBOXID SPR MPR
$# ssid msid sstyp mstyp sboxid mboxid spr mpr
910 0 2
$ FS FD DC VC VDC PENCHK BT DT
$# fs fd dc vc vdc penchk bt dt
0.100000 0.0 0.0 0.0 0.0 0 0.0 0.0
$ SFS SFM SST MST SFST SFMT FSF VSF
$# sfs sfm sst mst sfst sfmt fsf vsf
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$ SOFT SOFSCL LCIDAB MAXPAR EDGE DEPTH BSORT FRCFRQ
$# soft sofscl lcidab maxpar sbopt depth bsort frcfrq
2 0.0 0 0.0 0.0
$——————————————————————————-
*CONTROL_CONTACT
0.10000
2
$——————————————————————————-
*SET_PART_LIST
910,0.0,0.0,0.0,0.0
503,105
$——————————————————————————-
*SET_PART_LIST
922,0.0,0.0,0.0,0.0
503,105
$——————————————————————————-
*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
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ 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
1
$ a1 o1
0.0 3.03975e5
1000000.0 3.03975e5
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
bag_mesh5.k
*END