add icfd model:mesh morphing

Author: zhangzhanqun

examples/ICFD/icfd_mesh_morphing.py

@@ -0,0 +1,89 @@
+"""
+Mesh Morphing
+=============
+
+This example shows a simple ICFD problem with mesh morphing.
+"""
+
+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
+)
+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.mesh_morphing + os.sep
+fns.append(path + "mesh_morphing.k")
+solution.open_files(fns)
+solution.set_termination(termination_time=40)
+icfd = DynaICFD()
+solution.add(icfd)
+
+icfdanalysis = ICFDAnalysis()
+icfdanalysis.set_timestep(0.05)
+icfd.add(icfdanalysis)
+
+# define model
+mat = MatICFD(flow_density=1.0, dynamic_viscosity=0.005)
+
+part_inflow = ICFDPart(1)
+part_inflow.set_material(mat)
+part_inflow.set_prescribed_velocity(dof=ICFDDOF.X, motion=Curve(x=[0, 5, 6, 10000], y=[0, 0, 1, 1]))
+part_inflow.set_prescribed_velocity(dof=ICFDDOF.Y, motion=Curve(x=[0, 10000], y=[0, 0]))
+icfd.parts.add(part_inflow)
+
+part_outflow = ICFDPart(2)
+part_outflow.set_material(mat)
+part_outflow.set_prescribed_pressure(pressure=Curve(x=[0, 10000], y=[0, 0]))
+icfd.parts.add(part_outflow)
+
+part_symmetric = ICFDPart(3)
+part_symmetric.set_material(mat)
+part_symmetric.set_free_slip()
+icfd.parts.add(part_symmetric)
+
+part_wall = ICFDPart(4)
+part_wall.set_material(mat)
+part_wall.set_non_slip()
+part_wall.compute_drag_force()
+part_wall.set_boundary_layer(number=2)
+part_wall.set_imposed_move(vy=Curve(func="2*3.14/10*sin(2*3.14/10*TIME+3.14/2)"))
+icfd.parts.add(part_wall)
+
+part_meshmorph = ICFDPart(5)
+part_meshmorph.set_material(mat)
+icfd.parts.add(part_meshmorph)
+
+partvol1 = ICFDVolumePart(surfaces=[1, 2, 3, 5])
+partvol1.set_material(mat)
+icfd.parts.add(partvol1)
+partvol2 = ICFDVolumePart(surfaces=[5, 4])
+partvol2.set_material(mat)
+partvol2.set_imposed_move(vy=Curve(func="2*3.14/10*sin(2*3.14/10*TIME+3.14/2)"))
+icfd.parts.add(partvol2)
+# define the volume space that will be meshed,The boundaries
+# of the volume are the surfaces "spids"
+meshvol = MeshedVolume(surfaces=[1, 2, 3, 4])
+meshvol.set_fluid_interfaces([5])
+icfd.add(meshvol)
+
+solution.create_database_binary(dt=0.5)
+solution.save_file()

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

@@ -619,15 +619,37 @@ def __init__(self, surfaces):
         self.secid = 1
         self.mid = 0
         self.surfaces = surfaces
+        self.defined_imposed_move = False
 
     def set_material(self, mat):
         """Set material."""
         self.mid = mat.material_id
 
+    def set_imposed_move(self, vx=None, vy=None, vz=None):
+        """Impose a velocity on specific ICFD part."""
+        self.defined_imposed_move = True
+        self.vx = vx
+        self.vy = vy
+        self.vz = vz
+
     def create(self):
         """Create ICFD volume part."""
         ret = self.stub.ICFDCreatePartVol(ICFDPartVolRequest(secid=1, mid=self.mid, spids=self.surfaces))
         self.id = ret.id
+        if self.defined_imposed_move:
+            lcvx, lcvy, lcvz = 0, 0, 0
+            if self.vx != None:
+                self.vx.create(self.stub)
+                lcvx = self.vx.id
+            if self.vy != None:
+                self.vy.create(self.stub)
+                lcvy = self.vy.id
+            if self.vz != None:
+                self.vz.create(self.stub)
+                lcvz = self.vz.id
+            ret = self.stub.ICFDCreateControlImposedMove(
+                ICFDControlImposedMoveRequest(pid=self.id, lcvx=lcvx, lcvy=lcvy, lcvz=lcvz)
+            )
         logging.info(f"ICFD part volume {self.id} Created...")
         return ret
 

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

@@ -27,6 +27,7 @@
 strong_fsi = os.path.join(_module_path, "icfd", "strong_fsi")
 imposed_move = os.path.join(_module_path, "icfd", "imposed_move")
 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")
 
 iga_sample = os.path.join(_module_path, "iga", "iga_sample")