ICFD model for sloshing FSI in a container- initial velocity issue

- September 11, 2023 at 6:05 pm
  
  Md Zunayed Habib
  Subscriber
  
  Hello,
  I was trying to model a tank partially filled with water to impact. When I tried the ICFD fluid model separately, the models works fine, but when I try to couple them, the model doesn't coverge and after several iterations it says it has some intersecting faces, which should not happen. I tried to add initial velocity to the solid using *INITIAL_VELOCITY_GENERATION and to the fluid using *ICFD_INITIAL. When i remove the *ICFD_INITIAL card, the FSI model runs fine, but the fluid starts at a 0 velocity. I want the whole fluid volume to have an initial velocity same to the solid container.
  Can anyone please suggest how should I try to solve this issue?
  
  I am adding my keywords (with out the mesh) below:
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  Solid
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  *CONTROL_CONTACT
  $ slsfac rwpnal islchk shlthk penopt thkchg orien
  0.1 1.0 1 0 1 0 1 0
  $ usrstr usrfrc nsbcs interm xpene ssthk ecdt tiedprj
  0 0 0 0 4.0 1 1 0
  $# sfric dfric edc vfc th th_sf pen_sf ptscl
  0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0
  $# ignore frceng skiprwg outseg spotstp spotdel spothin
  0 0 0 0 0 0 0.0
  $# isym nserod rwgaps rwgdth rwksf icov swradf ithoff
  0 0 1 0.0 1.0 0 0.0 0
  $# shledg pstiff ithcnt tdcnof ftall unused shltrw igactc
  0 0 0 0 0 0.0 0
  *CONTROL_ENERGY
  $# hgen rwen slnten rylen irgen maten drlen disen
  2 2 2 1 2 1 1 1
  *CONTROL_IMPLICIT_AUTO
  $# iauto iteopt itewin dtmin dtmax dtexp kfail kcycle
  1 100 20 0.0 -700 0.0 0 0
  *CONTROL_IMPLICIT_DYNAMICS
  $# imass gamma beta tdybir tdydth tdybur irate alpha
  1 0.55 0.27563 0.01.00000E281.00000E28 1 0.0
  *CONTROL_IMPLICIT_GENERAL
  $# imflag dt0 imform nsbs igs cnstn form zero_v
  1 0.1 2 0 1 0 0 0
  *CONTROL_IMPLICIT_SOLUTION
  $# nsolvr ilimit maxref dctol ectol rctol lstol abstol
  12 11 15 0.001 0.01 0.0 0.01.0000E-20
  $# dnorm diverg istif nlprint nlnorm d3itctl cpchk
  2 1 1 3 2 10 0
  $# arcctl arcdir arclen arcmth arcdmp arcpsi arcalf arctim
  0 0 0.0 1 2 0.0 0.0 0.0
  $# lsmtd lsdir irad srad awgt sred
  5 2 0.0 0.0 0.0 0.0
  *CONTROL_OUTPUT
  $# npopt neecho nrefup iaccop opifs ipnint ikedit iflush
  1 3 0 0 0.0 0 100 5000
  *CONTROL_SHELL
  $# wrpang esort irnxx istupd theory bwc miter proj
  20.0 1 -1 0 2 1 1 1
  $# rotascl intgrd lamsht cstyp6 thshel
  1.0 0 0 1 0
  $# psstupd sidt4tu cntco itsflg irquad w-mode stretch icrq
  0 0 0 0 2 0.0 0.0 0
  $# nfail1 nfail4 psnfail keepcs delfr drcpsid drcprm intperr
  0 0 0 0 0 0 1.0 0
  *CONTROL_TERMINATION
  $# endtim endcyc dtmin endeng endmas nosol
  100.0 0 0.0 0.01.000000E8 0
  *DATABASE_GLSTAT
  $# dt binary lcur ioopt
  0.1 0 0 1
  *DATABASE_MATSUM
  $# dt binary lcur ioopt
  0.1 0 0 1
  *DATABASE_NODOUT
  $# dt binary lcur ioopt option1 option2
  0.01 0 0 1 0.0 0
  *DATABASE_RCFORC
  $# dt binary lcur ioopt
  0.1 0 0 1
  *DATABASE_BINARY_D3PLOT
  $# dt lcdt beam npltc psetid
  2.0 0 0 0 0
  *DATABASE_EXTENT_BINARY
  $# neiph neips maxint strflg sigflg epsflg rltflg engflg
  0 0 3 0 2 2 2 2
  $# cmpflg ieverp beamip dcomp shge stssz n3thdt ialemat
  0 1 0 1 1 1 2 1
  *MAT_PIECEWISE_LINEAR_PLASTICITY
  $# mid ro e pr sigy etan fail tdel
  1007.86000E-6 200.0 0.28 0.326 0.01.00000E21 0.0
  $# c p lcss lcsr vp
  40.0 5.0 0 0 1.0
  $ plastic stress strain curve
  $# eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
  0.0 0.0152 0.0226 0.0407 0.0691 0.0983 0.1345 0.7093
  $# es1 es2 es3 es4 es5 es6 es7 es8
  0.326 0.328 0.3788 0.4414 0.497 0.53 0.5557 0.7604
  *PART
  $# title
  Container shell
  $# pid secid mid eosid hgid grav adpopt tmid
  100 100 100 0 0 0 0 0
  *SECTION_SHELL
  $# secid elform shrf nip propt qr/irid icomp setyp
  100 2 1.0 5 1.0 0 0 1
  $ t1 t2 t3 t4 nloc
  $# t1 t2 t3 t4 nloc marea idof edgset
  3.0 3.0 3.0 3.0 0.0 0.0 0.0 0
  *INITIAL_VELOCITY_GENERATION
  $# id styp omega vx vy vz ivatn icid
  100 2 0.0 20.0 0.0 0.0 0 0
  $# xc yc zc nx ny nz phase irigid
  0.0 0.0 0.0 0.0 0.0 0.0 0 0
  *DEFINE_CURVE
  $# lcid sidr sfa sfo offa offo dattyp lcint
  700 0 1.0 1.0 0.0 0.0 0 0
  $# a1 o1
  0.0 0.1
  1.0 0.1
  100.0 0.1
  *RIGIDWALL_PLANAR_ID
  $# id title
  1
  $# nsid nsidex boxid offset birth death rwksf
  0 0 0 0.0 0.01.00000E20 1.0
  $# xt yt zt xh yh zh fric wvel
  1100.0 0.0 0.0 1000.0 100.0 0.0 0.0 0.0
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  Fluid and FSI
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  *ICFD_BOUNDARY_FSI
  $# pid
  1
  2
  *ICFD_CONTROL_FSI
  $# owc bt dt idc lcidsf xproj
  0 0.01.00000E28 0.25 0 0
  $# nsub
  0
  *ICFD_INITIAL
  $# pid vx vy vz t p dfunc
  0 20.0 0.0 0.0 0.0 0.0 0 0
  *LOAD_BODY_Z
  $# lcid sf lciddr xc yc zc cid
  2 1.0 0 0.0 0.0 0.0 0
  *DEFINE_CURVE_TITLE
  Velocity imposed
  $# lcid sidr sfa sfo offa offo dattyp lcint
  1 0 100.0 20.0 0.0 0.0 0 0
  $# a1 o1
  0.0 1.0
  10.0 1.0
  *DEFINE_CURVE_TITLE
  Gravity force
  $# lcid sidr sfa sfo offa offo dattyp lcint
  2 0 1.0 0.00981 0.0 0.0 0 0
  $# a1 o1
  0.0 1.0
  10000.0 1.0
  *INCLUDE
  ICFD_mesh.k
  *ICFD_BOUNDARY_NONSLIP
  $# pid
  1
  *ICFD_BOUNDARY_NONSLIP
  $# pid
  2
  *ICFD_CONTROL_TIME
  $# ttm dt cfl lcidsf dtmin dtmax dtinit tdeath
  100.0 0.0 1.0 01.00000E-91.00000E28 0.01.00000E28
  *ICFD_MAT
  $# mid flg ro vis st stsflcid ca
  1 11.00000E-61.00500E-9 0.0 0 0.0
  *ICFD_MAT
  $# mid flg ro vis st stsflcid ca
  2 0 0.0 0.0 0.0 0 0.0
  *ICFD_PART_TITLE
  Water
  $# pid secid mid
  1 1 1
  *ICFD_PART_TITLE
  Air
  $# pid secid mid
  2 1 2
  *ICFD_PART_TITLE
  Interface
  $# pid secid mid
  3 1 1
  *ICFD_PART_VOL
  $# pid secid mid
  10 1 1
  $# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
  1 3 0 0 0 0 0 0
  *ICFD_PART_VOL
  $# pid secid mid
  20 1 2
  $# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
  2 3 0 0 0 0 0 0
  *ICFD_SECTION
  $# sid
  1
  *MESH_INTERF
  $# volid
  10
  $# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
  3 0 0 0 0 0 0 0
  *MESH_VOLUME
  $# volid
  10
  $# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
  1 2 0 0 0 0 0 0
- September 12, 2023 at 4:12 pm
  
  Reno Genest
  Ansys Employee
  
  Hello,
  Have you checked the iCFD examples on sloshing?
  https://www.dynaexamples.com/icfd
  https://www.dynaexamples.com/icfd/basics-examples/global_imposed
  https://www.dynaexamples.com/icfd/intermediate-examples/slosh2
  
  There is also a course on iCFD on the Ansys Learning Hub (ALH):
  https://jam8.sapjam.com/groups/QyHAlhciRCZ0Xt6xKZ2MHi/overview_page/jPMouxqQ2AOP6jVZLBMZCl
  
  Reno.

Viewing 1 reply thread

You must be logged in to reply to this topic.

Ansys Innovation Space

Boost Ansys Fluent Simulations with AWS

Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.

Earth Rescue – An Ansys Online Series

The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.

Ansys Blog

Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.

Trending discussions