Merge pull request #32 from ulgltas/lacroix

Adaptive Time Step Support

Author: rboman

README.md

@@ -1,25 +1,29 @@
 # CUPyDO
+
 FSI tools for partinioned coupling between generic solid and fluid solvers.
 
 [![Apache License Version 2.0](https://img.shields.io/badge/license-Apache_2.0-green.svg)](LICENSE)
 
 ## Solvers
+
 CUPyDO currently features interfaces for the following solvers:
 - **Solid**
-  - Metafor **[v3514]** (http://metafor.ltas.ulg.ac.be/dokuwiki/start)
+  - Metafor **[v3516]** (http://metafor.ltas.ulg.ac.be/dokuwiki/start)
   - A Nonlinear Finite Element solid solver developed at the University of Liège.
   ---
-  - RBMI **[1.1]** (https://github.com/ulgltas/NativeSolid)
+  - NativeSolid (a.k.a. RBM) **[v1.2]** (https://github.com/ulgltas/NativeSolid)
   - A dynamic 2dof pitch/plunge solid solver developed at the University of Liège.
   ---
-  - SU2 **[X]** (https://su2code.github.io/)
+  - SU2 **[-]** (https://su2code.github.io/)
   - Open-source CFD code developed at Stanford University.
+  - ⚠️ **This interface is currently broken**
   ---
   - Modali **[v2.0]** (https://github.com/ulgltas/modali)
   - A static/dynamic modal solver developed at the University of Liège.
   ---
-  - GetDP **[X]** (http://getdp.info/)
+  - GetDP **[-]** (http://getdp.info/)
   - A free finite element software and a general environment for the treatment of discrete problems, developed at University of Liège.
+   - ⚠️ **This interface is currently broken**
   ---
   - pyBeam **[master]** (https://github.com/pyBeam/pyBeam)
   - A nonlinear beam finite element solver developed for aeronautical design applications.
@@ -31,7 +35,7 @@ CUPyDO currently features interfaces for the following solvers:
   - SU2 **[fix_wrap_strong]** (https://github.com/ulgltas/SU2/tree/fix_wrap_strong)
   - Open-source CFD code developed at Stanford University.
   ---
-  - DART **[v1.2.0]** (https://gitlab.uliege.be/am-dept/dartflo)
+  - DART **[v1.2.1]** (https://gitlab.uliege.be/am-dept/dartflo)
   - Open-source transonic full potential finite element fluid solver, developed at the University of Liège.
   ---
   - VLM **[v2.0]** (https://github.com/ulgltas/VLM)
@@ -55,14 +59,17 @@ Finally, CUPyDO features two main couplers:
   - Multi-Vector Jacobian (IQN_MVJ)
 
 ## Compilation
+
 Detailed build instructions can be found in the [wiki](https://github.com/ulgltas/CUPyDO/wiki/Installation).
 
 ## Examples
+
 Examples of simulations are available in [CUPyDO/tests](https://github.com/ulgltas/CUPyDO/tree/master/tests) and [CUPyDO/cases](https://github.com/ulgltas/CUPyDO/tree/master/cases).
 
 ![Screenshot](/tests/fsi_examples.png)
 
-## Publications:
+## Publications
+
 Cerquaglia M.L., Thomas D., Boman R., Terrapon V.E., Ponthot J.-P., [A fully partitioned Lagrangian framework for FSI problems characterized by free surfaces, large solid deformations and displacements, and strong added-mass effects](https://doi.org/10.1016/j.cma.2019.01.021), Computer Methods in Applied Mechanics and Engineering, in press (2019)
 
 Thomas D., Cerquaglia M.L., Boman R., Economon T.D., Alonso J.J., Dimitriadis G., Terrapon V.E., [CUPyDO - An integrated Python environment for coupled multi-physics simulations](https://doi.org/10.1016/j.advengsoft.2018.05.007), Advances in Engineering Software 128:69-85 (2019)

cases/PFEM3D_Metafor/birdAxisym/input_fsi.py

@@ -38,9 +38,9 @@ def getFsiP():
 
     # FSI objects
 
-    p['criterion'] = 'Displacements'
-    p['interpolator'] = 'Matching'
-    p['algorithm'] = 'AitkenBGS'
+    p['criterion'] = 'displacement'
+    p['interpolator'] = 'matching'
+    p['algorithm'] = 'aitkenBGS'
 
     # FSI parameters
 
@@ -66,7 +66,7 @@ def main():
     cupydo = cupy.CUPyDO(param)
     cupydo.run()
 
-    cupydo.algorithm.FluidSolver.save(cupydo.algorithm.timeIter)
+    cupydo.algorithm.FluidSolver.save(cupydo.algorithm.step.timeIter)
     test(cupydo.algorithm.getMeanNbOfFSIIt())
 
 if __name__=='__main__':

cases/PFEM3D_Metafor/birdAxisym/input_meta.py

@@ -7,7 +7,6 @@
 def params(input):
 
     input['bndno'] = 12
-    input['saveAllFacs'] = False
     input['bctype'] = 'pydeadloads'
     return input
 

cases/PFEM_Metafor/StaticCylinder_fluid_Pfem.py

@@ -86,9 +86,7 @@ def getPfem():
     w.Boundary(msh, 7, 2, 0.0)
     w.Boundary(msh, 9, 1, 0.0)
     w.Boundary(msh, 9, 2, 0.0)
-    
-    scheme.savefreq=1
-    scheme.nthreads=3
+
     scheme.gamma = 0.6
     scheme.omega = 0.5
     scheme.addRemoveNodesOption = True

cases/PFEM_Metafor/VIV_cantileverBeam_air_Pfem.py

@@ -82,8 +82,6 @@ def getPfem():
     n3.isOnAxisOfSymmetry = True
     # ---
 
-    scheme.savefreq=100
-    scheme.nthreads=3
     scheme.gamma = 0.7
     scheme.omega = 0.7
     scheme.addRemoveNodesOption = True

cases/PFEM_Metafor/VIV_cantileverBeam_freeSlip_air_Pfem.py

@@ -85,9 +85,7 @@ def getPfem():
     n2.isOnFreeSlipBoundaryY = True
     n3.isOnFreeSlipBoundaryY = True
     # ---
-    
-    scheme.savefreq=100
-    scheme.nthreads=4
+
     scheme.gamma = 0.8
     scheme.omega = 0.7
     scheme.addRemoveNodesOption = True

cases/PFEM_Metafor/birdImpact_deformable_panel_bird_Pfem.py

@@ -66,9 +66,7 @@ def getPfem():
     bird = w.Group(msh, 16)
     loadingset = w.LoadingSet(msh)
     loadingset.add(1,w.InitialVelocity(msh,bird,0.,-U0,0.))
-    
-    scheme.savefreq=1
-    scheme.nthreads=3
+
     scheme.gamma = 0.5
     scheme.omega = 0.5
     scheme.addRemoveNodesOption = True

cases/PFEM_Metafor/birdImpact_deformable_panel_bird_Pfem_Axisym.py

@@ -72,8 +72,6 @@ def getPfem():
     loadingset = w.LoadingSet(msh)
     loadingset.add(1,w.InitialVelocity(msh,bird,0.,-U0,0.))
 
-    scheme.savefreq=1
-    scheme.nthreads=3
     scheme.gamma = 0.5
     scheme.omega = 0.5
     scheme.addRemoveNodesOption = True

cases/PFEM_Metafor/birdImpact_deformable_panel_bird_Pfem_Smojver.py

@@ -66,9 +66,7 @@ def getPfem():
     bird = w.Group(msh, 16)
     loadingset = w.LoadingSet(msh)
     loadingset.add(1,w.InitialVelocity(msh,bird,0.,-U0,0.))
-    
-    scheme.savefreq=1
-    scheme.nthreads=3
+
     scheme.gamma = 0.5
     scheme.omega = 0.5
     scheme.addRemoveNodesOption = True

cases/PFEM_Metafor/birdImpact_deformable_panel_bird_Pfem_Smojver_Axisym.py

@@ -72,8 +72,6 @@ def getPfem():
     loadingset = w.LoadingSet(msh)
     loadingset.add(1,w.InitialVelocity(msh,bird,0.,-U0,0.))
 
-    scheme.savefreq=1
-    scheme.nthreads=3
     scheme.gamma = 0.5
     scheme.omega = 0.5
     scheme.addRemoveNodesOption = True