$-----------------------------------------------------------------------------
$
$ Example provided by Peggy and Iñaki (LST-Ansys)
$
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA MPP R12.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (S.I)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*INCLUDE
mesh.k
*PARAMETER
$--- Fluid
$
R v_inlet 0.3638
Rrho_fluid 1035.0
R mu_fluid 0.0035
R dt_fluid 1.0e-5
R turb_I 0.07
$
$---- Rotor
$
R omega 366.52
$
$ ----------------- Control General --------------------------------------------
*ICFD_CONTROL_OUTPUT
$# msgl
4
*ICFD_CONTROL_TURBULENCE
3
*ICFD_CONTROL_MESH
1.2
*ICFD_CONTROL_TIME
$# ttm dt cfl lcidsf dtmin dtmax dtinit tdeath
0.34288 0 1.0 0 0.0 0.0 0.0 0.0
$ ----------------- Control Sliding --------------------------------------------
*ICFD_BOUNDARY_PERIODIC
4,3,5,
*ICFD_CONTROL_IMPOSED_MOVE
15
,,,100
*ICFD_CONTROL_IMPOSED_MOVE
5
,,,100
*ICFD_CONTROL_IMPOSED_MOVE
6
,,,100
$ ----------------- Parts and Boundary conditions-------------------------------
*ICFD_BOUNDARY_NONSLIP
$# pid
3
6
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 4 1 20 1.0 01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
2 2 1.01.00000E28 0.0
*ICFD_MAT
$# mid flg ro vis st
1 1&rho_fluid &mu_fluid 0.0
*ICFD_PART_TITLE
inlet
$# pid secid mid
1 1 1
*ICFD_PART_TITLE
outle
$# pid secid mid
2 1 1
*ICFD_PART_TITLE
rotor
$# pid secid mid
6 1 1
*ICFD_PART_TITLE
housing
$# pid secid mid
3 1 1
*ICFD_PART_TITLE
sliding_out
$# pid secid mid
4 1 1
*ICFD_PART_TITLE
sliding_in
$# pid secid mid
5 1 1
*ICFD_PART_VOL_title
fluid_volume_out
$# pid secid mid
15 1 1
5,6
*ICFD_PART_VOL_TITLE
fluid_volume_in
$# pid secid mid
16 1 1
1,2,3,4
*ICFD_SECTION
$# sid
1
$-----------------------------Output--------------------------------------------
*ICFD_CONTROL_TAVERAGE
$# dt
0.017143
*ICFD_DATABASE_FLUX
1
*ICFD_DATABASE_FLUX
2
*ICFD_DATABASE_DRAG
6
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
0.017143
$# ioopt
0
$-----------------------------MESH----------------------------------------------
*MESH_VOLUME
$# volid
15
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
5,6
*MESH_VOLUME
$# volid
16
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1,2,3,4
*MESH_BL
$# pid nelth blth blfe blst
3 4
*MESH_BL
$# pid nelth
6 4
*MESH_SIZE_SHAPE
$# sname
box
$# msize pminx pminy pminz pmaxx pmaxy pmaxz
0.00035 0.0 -0.03 0.004 0.024 -0.026 0.009
$ ----------------- Define Curves ----------------------------------------------
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 1.0
10000.0 1.0
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp lcint
2 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
10000.0 0.0
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp lcint
3 0 1.0 &turb_I 0.0 0.0 0 0
$# a1 o1
0.0 1.0
10000.0 1.0
*DEFINE_FUNCTION
20,x-velo
f(y,z)=-4.1977e14*(y**2+(z-0.294)**2)**3+5.5e12*sqrt(y**2+(z-0.294)**2)**5
-2.65e10*(y**2+(z-0.294)**2)**2+5.64e7*sqrt(y**2+(z-0.294)**2)**3
-5.86e4*(y**2+(z-0.294)**2)+2.365*sqrt(y**2+(z-0.294)**2)+1.157
*DEFINE_CURVE_TITLE
Rotation velocity
$# lcid sidr sfa sfo offa offo dattyp
100 &omega
$# a1 o1
0.0 0.0
0.0857125 1.0
10000.0 1.0
*END
