add icfd model:sloshing

Author: zhangzhanqun

examples/ICFD/icfd_sloshing.py

@@ -0,0 +1,85 @@
+"""
+Sloshing
+========
+
+This example shows how to directly impose a displacement of the entire volume mesh through the use of the keyword ICFD_CONTROL_IMPOSED_MOVE.
+"""
+
+import os
+import sys
+
+
+from ansys.dyna.core.pre.dynasolution import DynaSolution
+from ansys.dyna.core.pre.dynaicfd import (
+    DynaICFD,
+    MatICFD,
+    ICFDPart,
+    ICFDDOF,
+    Curve,
+    ICFDVolumePart,
+    MeshedVolume,
+    ICFDAnalysis,
+    Curve,
+    Gravity,
+    GravityOption,
+    Compressible
+    )
+from ansys.dyna.core.pre import examples
+
+
+hostname = "localhost"
+if len(sys.argv) > 1:
+    hostname = sys.argv[1]
+
+solution = DynaSolution(hostname)
+# Import the initial mesh data(nodes and elements)
+fns = []
+path = examples.sloshing + os.sep
+fns.append(path + "sloshing.k")
+solution.open_files(fns)
+solution.set_termination(termination_time=1)
+icfd = DynaICFD()
+solution.add(icfd)
+
+icfdanalysis = ICFDAnalysis()
+icfdanalysis.set_timestep(0.02)
+icfd.add(icfdanalysis)
+
+# define model
+mat1 = MatICFD(flow_density=1000, dynamic_viscosity=0.001)
+mat2 = MatICFD(flag=Compressible.VACUUM)
+
+part_bottom = ICFDPart(1)
+part_bottom.set_material(mat1)
+part_bottom.set_non_slip()
+icfd.parts.add(part_bottom)
+
+part_top = ICFDPart(2)
+part_top.set_material(mat2)
+part_top.set_non_slip()
+icfd.parts.add(part_top)
+
+part_mid = ICFDPart(3)
+part_mid.set_material(mat1)
+icfd.parts.add(part_mid)
+
+g = Gravity(dir=GravityOption.DIR_Z, load=Curve(x=[0, 10000], y=[1, 1]))
+icfd.add(g)
+
+icfd.set_imposed_move(vx=Curve(x=[0, 0.5, 0.52, 0.8, 0.82, 2.0], y=[1, 1, -1, -1, 0, 0]))
+
+partvol_bottom = ICFDVolumePart(surfaces=[1, 3])
+partvol_bottom.set_material(mat1)
+icfd.parts.add(partvol_bottom)
+
+partvol_top = ICFDVolumePart(surfaces=[2, 3])
+partvol_top.set_material(mat2)
+icfd.parts.add(partvol_top)
+# define the volume space that will be meshed,The boundaries
+# of the volume are the surfaces "spids"
+meshvol = MeshedVolume(surfaces=[1, 2])
+meshvol.set_fluid_interfaces([3])
+icfd.add(meshvol)
+
+solution.create_database_binary(dt=0.02)
+solution.save_file()

src/ansys/dyna/core/pre/dynaicfd.py

@@ -186,6 +186,23 @@ def set_initial(self, velocity=Velocity(0, 0, 0), temperature=0, pressure=0):
         logging.info("ICFD_INIT Created...")
         return ret
 
+    def set_imposed_move(self, vx=None, vy=None, vz=None):
+        """Impose a velocity on the whole volume mesh."""
+        lcvx, lcvy, lcvz = 0, 0, 0
+        if vx != None:
+            vx.create(self.stub)
+            lcvx = vx.id
+        if vy != None:
+            vy.create(self.stub)
+            lcvy = vy.id
+        if vz != None:
+            vz.create(self.stub)
+            lcvz = vz.id
+        ret = self.stub.ICFDCreateControlImposedMove(
+            ICFDControlImposedMoveRequest(pid=0, lcvx=lcvx, lcvy=lcvy, lcvz=lcvz)
+        )
+        return ret
+
     def save_file(self):
         """Save keyword files."""
         self.create_section_icfd(1)

src/ansys/dyna/core/pre/examples/examples.py

@@ -29,6 +29,7 @@
 mesh_adaptivity = os.path.join(_module_path, "icfd", "mesh_adaptivity")
 mesh_morphing = os.path.join(_module_path, "icfd", "mesh_morphing")
 dem_coupling = os.path.join(_module_path, "icfd", "dem_coupling")
+sloshing = os.path.join(_module_path, "icfd", "sloshing")
 
 iga_sample = os.path.join(_module_path, "iga", "iga_sample")
 camry_rc = os.path.join(_module_path, "implicit", "camry_rc")

src/ansys/dyna/core/pre/examples/icfd/cylinder_flow/lspost.cfile

@@ -0,0 +1,4 @@
+$# LS-PrePost command file created by LS-PrePost(R) V4.10 (Beta)-30Dec2022 -64bit-Window
+$# Created on Jan-9-2023 (09:30:19)
+open keyword "D:\pyDynavenv\pyDyna\src\ansys\dyna\core\pre\examples\icfd\cylinder_flow\cylinder_flow.k"
+open keyword "D:\!WORK\ICFD\ICFDPRE_EXAMPLE\Basics\7_User_defined_mesh\i.k"

src/ansys/dyna/core/pre/examples/icfd/cylinder_flow/lspost.msg

@@ -0,0 +1,10 @@
+open keyword "D:\pyDynavenv\pyDyna\src\ansys\dyna\core\pre\examples\icfd\cylinder_flow\cylinder_flow.k"
+ Time to read keyword file =  0.141 seconds
+ Time to process element data =  0.109 seconds
+Finished reading model
+open keyword "D:\!WORK\ICFD\ICFDPRE_EXAMPLE\Basics\7_User_defined_mesh\i.k"
+LS-PrePost encounters an error when open D:\pyDynavenv\pyDyna\src\ansys\dyna\core\pre\examples\icfd\cylinder_flow\mesh.k
+because the path is not found
+Error - Include File mesh.k Not open
+  Error occurs in file D:\!WORK\ICFD\ICFDPRE_EXAMPLE\Basics\7_User_defined_mesh\i.k line # 31
+Save Config. parameters to file C:\Users\dell\AppData\Roaming\LSTC\LS-PrePost4.10\lsppconf

src/ansys/dyna/core/pre/examples/icfd/sloshing/sloshing.k

@@ -827,4 +827,61 @@ def test_icfd_dem_coupling(resolve_icfd_path, resolve_server_path, resolve_stand
     solution.save_file()
     outputfile = os.path.join(resolve_server_path, "output", "test_dem_coupling.k")
     standardfile = os.path.join(resolve_standard_path, "icfd", "dem_coupling.k")
+    assert comparefile(outputfile, standardfile)
+
+def test_icfd_sloshing(resolve_icfd_path, resolve_server_path, resolve_standard_path):
+    solution = DynaSolution("localhost")
+    icfd_initialfile = os.path.join(resolve_icfd_path, "test_sloshing.k")
+    fns = []
+    fns.append(icfd_initialfile)
+    solution.open_files(fns)
+    solution.set_termination(termination_time=1)
+    icfd = DynaICFD()
+    solution.add(icfd)
+
+    icfdanalysis = ICFDAnalysis()
+    icfdanalysis.set_timestep(0.02)
+    icfd.add(icfdanalysis)
+
+    # define model
+    mat1 = MatICFD(flow_density=1000, dynamic_viscosity=0.001)
+    mat2 = MatICFD(flag=Compressible.VACUUM)
+
+    part_bottom = ICFDPart(1)
+    part_bottom.set_material(mat1)
+    part_bottom.set_non_slip()
+    icfd.parts.add(part_bottom)
+
+    part_top = ICFDPart(2)
+    part_top.set_material(mat2)
+    part_top.set_non_slip()
+    icfd.parts.add(part_top)
+
+    part_mid = ICFDPart(3)
+    part_mid.set_material(mat1)
+    icfd.parts.add(part_mid)
+
+    g = Gravity(dir=GravityOption.DIR_Z, load=Curve(x=[0, 10000], y=[1, 1]))
+    icfd.add(g)
+
+    icfd.set_imposed_move(vx=Curve(x=[0, 0.5, 0.52, 0.8, 0.82, 2.0], y=[1, 1, -1, -1, 0, 0]))
+
+    partvol_bottom = ICFDVolumePart(surfaces=[1, 3])
+    partvol_bottom.set_material(mat1)
+    icfd.parts.add(partvol_bottom)
+
+    partvol_top = ICFDVolumePart(surfaces=[2, 3])
+    partvol_top.set_material(mat2)
+    icfd.parts.add(partvol_top)
+    # define the volume space that will be meshed,The boundaries
+    # of the volume are the surfaces "spids"
+    meshvol = MeshedVolume(surfaces=[1, 2])
+    meshvol.set_fluid_interfaces([3])
+    icfd.add(meshvol)
+
+    solution.create_database_binary(dt=0.02)
+    solution.save_file()
+
+    outputfile = os.path.join(resolve_server_path, "output", "test_sloshing.k")
+    standardfile = os.path.join(resolve_standard_path, "icfd", "sloshing.k")
     assert comparefile(outputfile, standardfile)