*KEYWORD *TITLE slosh3.k = ICFD Sloshing with global imposed move kg-m-s-K *INCLUDE geo1.k $=============================================================================== *PARAMETER R xminstr -0.5 R yminstr -0.5 R zminstr 0.0 R xmaxstr 0.5 R ymaxstr 0.5 R zmaxstr 1.0 $=============================================================================== R xminale -0.6 R yminale -0.6 R zminale -0.1 R xmaxale 0.6 R ymaxale 0.6 R zmaxale 1.1 $=============================================================================== R dl 0.02 $=============================================================================== *PARAMETER_EXPRESSION R nxale nint((xmaxale-xminale)/dl + 1) R nyale nint((ymaxale-yminale)/dl + 1) R nzale nint((zmaxale-zminale)/dl + 1) $=============================================================================== *PARAMETER R T_end 1.0 R dt_plot 0.02 $=============================================================================== R v_in 1. R v_back -1. Rrho_fluid 1000 R mu_fluid 0.001 R dt_fluid 0.020 R grav 9.81 $ $=============================================================================== $ [ ] EXECUTION CONTROLS $=============================================================================== *CONTROL_TERMINATION $ endtim endcyc dtmin endeng endmas unused &T_end *CONTROL_TIMESTEP $ dtinit tssfac isdo tslimt dt2ms lctm erode ms1st 0 0.6 0 0 0 0 0 0 $ dt2msf dt2mslc imscl unused 0 0 0 *CONTROL_ENERGY $ hgen rwen slnten rylen unused 2 2 2 2 $=============================================================================== $ [ ] OUTPUT CONTROLS $=============================================================================== *DATABASE_GLSTAT $ dt binary lcur ioopt dthff binhf unused1 1E-05 3 0 0 0 0 0 *DATABASE_MATSUM $ dt binary lcur ioopt dthff binhf unused1 1E-05 3 0 0 0 0 0 *DATABASE_ELOUT $ dt binary lcur ioopt dthff binhf unused1 1E-05 3 0 0 0 0 0 *DATABASE_EXTENT_BINARY $ neiph neips maxint strflg sigflg epsflg rltflg engflg 0 0 0 0 1 1 0 0 $ cmpflg ieverp beamip dcomp shge stssz n3thdt ialemat 0 0 0 0 0 0 0 0 $ nintsld pkp_sen sclp hydro msscl therm intout nodout 0 0 0 0 2 0 *DATABASE_BINARY_D3PLOT $ dt lcdt beam npltc psetid unused1 0.010 0 0 200 0 $ ioopt unused2 0 *DATABASE_BINARY_INTFOR $ dt lcdt beam npltc psetid unused 0.0005 0 0 0 0 $=============================================================================== $ [ ] |PART_ID|SECTION_ID|MAT_ID|EOS_ID|HG_ID|GRAV|ADPOPT|THERM_MID| DEFINITIONS $=============================================================================== *PART S1 = rigid can $ pid secid mid eosid hgid grav adpopt tmid 1 1 1 *SECTION_SHELL $ secid elform shrf nip propt qr/irid icomp setyp 1 2 0.000 1 $ t1 t2 t3 t4 nloc marea 0.008 0.008 0.008 0.008 *MAT_ELASTIC $ MID RHO E PR DA DB K 1 7850.0 2.1E+8 0.3 *MAT_RIGID $ mid ro e pr n couple m alias 11 7850.0000 210.00E9 0.300000 0.000 0.000 0.000 $ cmo con1 con2 1.000 5 7 $ lco or a1 a2 a3 v1 v2 v3 $=============================================================================== *PART H5 = AMMG1 = liq water $ pid secid mid eosid hgid grav adpopt tmid 5 5 5 5 *SECTION_SOLID $ secid elform 5 11 *MAT_NULL $ mid ro pc mu terod cerod ym pr 5 998.21 -8.7E6 0.8684E-3 0.0 0.0 0.0 0.0 *EOS_LINEAR_POLYNOMIAL $ eosid c0 c1 c2 c3 c4 c5 c6 5 101325.0 2.25E9 0.0 0.0 0.0 0.0 0.0 $ e0 v0 0.0 0.0 $------------------------------------------------------------------------------- *PART H6 = AMMG2 = air inside $ pid secid mid eosid hgid grav adpopt tmid 6 6 6 6 *SECTION_SOLID $ secid elform 6 11 *MAT_NULL $ MID RO PC MU TEROD CEROD YMBEAM PRBEAM 6 1.1845 0.0 1.8444E-5 0.0 0.0 0.0 0.0 *EOS_LINEAR_POLYNOMIAL $ EOSID C0 C1 C2 C3 C4 C5 C6 6 0.0 0 0 0 0.4 0.4 0 $ E0 V0 2.533125E5 1.0 $------------------------------------------------------------------------------- *PART H7 = AMMG3 = air outside $ pid secid mid eosid hgid grav adpopt tmid 7 6 6 6 $=============================================================================== $ [ ] EULERIAN & ALE & COUPLING CONSTRAINTS, ... $ SSTYP = SLAVE elm type: 0=PSID=part-set-id; 1=PID=part-id; 2=segm-set-id $=============================================================================== *ALE_STRUCTURED_MESH $ mshid dpid nbid ebid 3 4 40001 40001 $ nptx npty nptz nid0 lcsid 1001 1002 1003 1 101 $------------------------------------------------------------------------------- *ALE_STRUCTURED_MESH_CONTROL_POINTS $ CPID NONE NONE SFO NONE OFFO 1001 0.5 $ x1 x2 x3 x4 1 &xminale &nxale &xmaxale *ALE_STRUCTURED_MESH_CONTROL_POINTS $ CPID NONE NONE SFO NONE OFFO 1002 0.5 $ x1 x2 x3 x4 1 &yminale &nyale &ymaxale *ALE_STRUCTURED_MESH_CONTROL_POINTS $ CPID NONE NONE SFO NONE OFFO 1003 $ x1 x2 x3 x4 1 &zminale &nzale &zmaxale $---------------------------------- *DEFINE_COORDINATE_NODES $ CID NID1 NID1 NID1 FLAG DIR 101 1 2 52 0 $---------------------------------- *ALE_STRUCTURED_MESH_REFINE $ MSHID IFX IFY IFZ 1 3 3 3 $=============================================================================== *comment ALE_STRUCTURED_MESH_MOTION $ mshid command ammgset explim 3 FOLLOW_GC 12 1.1 $=============================================================================== *CONTROL_ALE $ dct nadv meth afac bfac cfac dfac efac -1 1 2 -1.0 0.0 0.0 0.0 0.0 $ start end aafac vfact prit ebc pref nsidebc 0.01.00000E20 0.0 0.0 0 0 101325.0 0 $ ncpl nbkt imascl checkr beamin mmgpref pdifmx dtmufac 1 50 0 0.0 0.0 0 0.0 0.0 $ optimpp ialedr bndflx minmas 0 0 0 $------------------------------------------------------------------------------- *ALE_MULTI-MATERIAL_GROUP $ pid typ 5 1 1h2o 6 1 2airin 7 1 3airout $---5----0----5----0----5----0----5----0----5----0----5----0----5----0----5----0 *ALE_STRUCTURED_MESH_VOLUME_FILLING $ MSHID notused AMMGTO notused NSAMPLE notused notused VID 3 1 4 $ GEOM IN/OUT E1 E2 E3 E4 E5 ALL 0 $---5----0----5----0----5----0----5----0----5----0----5----0----5----0----5----0 *ALE_STRUCTURED_MESH_VOLUME_FILLING $ MSHID notused AMMGTO notused NSAMPLE notused notused VID 3 2 4 $ sgs=GEOM IN/OUT nid=E1 offset=E2 E3 E4 E5 PLANE 0 20001 20002 $---5----0----5----0----5----0----5----0----5----0----5----0----5----0----5----0 *ALE_STRUCTURED_MESH_VOLUME_FILLING $ MSHID notused AMMGTO notused NSAMPLE notused notused VID 3 3 4 $ ell=GEOM IN/OUT SID=E1 offset=E2 RY=E3 RZ=E4 E5 PART 1 1 0.0 $---5----0----5----0----5----0----5----0----5----0----5----0----5----0----5----0 $ NODE ID= 5738 , xyz= 0, -0.5, 0.24 *NODE $ nid x y z tc rc 20001 0.0 -0.5 0.25 7 7 20002 0.0 -0.5 0.26 7 7 $------------------------------------------------------------------------------- *ALE_STRUCTURED_FSI $ SLAVE dpidMASTR SSTYP MSTYP notused notused notused MCOUP 1 4 1 1 -12 $ start end pfac notused notused norm notused notused 0.0 0.0 -12 $------------------------------------------------------------------------------- *SET_MULTI-MATERIAL_GROUP_LIST 12 1 2 *DEFINE_CURVE_TITLE Velocity imposed $# lcid sidr sfa sfo offa offo dattyp 12 $# a1 o1 0.0 0.0 0.020 10132.5 $------------------------------------------------------------------------------- *ALE_STRUCTURED_FSI $ SLAVE dpidMASTR SSTYP MSTYP notused notused notused MCOUP 1 4 1 1 -33 $ start end pfac notused notused norm notused notused 0.0 0.0 -12 1 $------------------------------------------------------------------------------- *SET_MULTI-MATERIAL_GROUP_LIST 33 3 $=============================================================================== $ [ ] BC's + IC's + BODY LOADS + FORCE FIELDS $=============================================================================== $ DOF:1=x;2=y;3=z;4=move along VID; VAD: 0=vel;1=accel;2=displacement $BOUNDARY_PRESCRIBED_MOTION_RIGID *BOUNDARY_PRESCRIBED_MOTION_SET $NID|NSID|PID DOF VAD LCID SF VID DEATH BIRTH 1 1 0 1 $ OFFSET1 OFFSET2 MRB NODE1 NODE2 DOF=9|10|11 or VAD=4 OPTCARD $ $------------------------------------------------------------------------------- *DEFINE_CURVE_TITLE Velocity imposed $# lcid sidr sfa sfo offa offo dattyp 1 $# a1 o1 0.0 0.0 1.0E-4 &v_in 0.5 &v_in 0.52 &v_back 0.80 &v_back 0.82 0 2.00 0 $------------------------------------------------------------------------------- *BOUNDARY_SPC_SET $ NID/NSID CID DOFX DOFY DOFZ DOFRX DOFRY DOFRZ 1 0 0 1 1 0 0 0 $------------------------------------------------------------------------------- *LOAD_BODY_Z $# lcid sf 2 1 *DEFINE_CURVE_TITLE Gravity force $# lcid sidr sfa sfo offa offo dattyp 2 &grav $# a1 o1 0.0 1.0 1.0 1.0 $=============================================================================== *END $===============================================================================