BLASTER v1.0

The viscous-inviscid interaction implementation was removed from DARTFlo and will be in BLASTER

Features

• Quasi-2D viscous solver

  • Transonic flows
  • Laminar and turbulent flows
  • Laminar to turbulent transition
  • Attached and mildly separated flows

• 2D and 3D viscous inviscid interface (Python)

  • 2D geometries with matching meshes
  • 3D geometries with Radial Basis Function interpolation
  • Solvers coupling
  • Quasi-simultaneous interaction law

• Interface with DART

• Interface with SU2 - Not tested to avoid cloning SU2

• 2D coupled adjoint methodology (developed with DART)

• Interface with OpenMDAO

  • NACA parametrisation
  • Free-form deformation with pyGeo

• Basic features (api, attached compressible flows, 2D adjoint, 3D transonic toy problem) are tested with scripts in the 'tests/' folder

• Examples for advanced cases (separated, 2D transonic, 3D transonic) are present in the 'validation/' folder

• 2D optimization cases (NACA subsonic and transonic) are present in the 'mdao/' folder

Notes

• 3D cases require a dedicated surface transfinite mesh to work. Examples are given in 'blast/models/' folder. This will be removed in the future
• SU2 interface is not tested for now and will be updated to work with the new v8 version (with reworked python api)

Future works

Future works will be focused on 3D improvements and the development of the 3D adjoint methodology. Additional features will also be included.

• Modification of the section generation in 3D. It will be automatic.

• New RBF interpolator with gradient computation for 3D adjoint

• 3D adjoint solver

• 3D optimization examples

• New SU2 v8 interface

• Interface with the panel method HSPM for fast 2D computation (Master thesis of Corentin Thomée)

• New semi-inverse boundary layer solver for increased robustness in separated cases (Master thesis of Massimo Volders)

blast/api/mda_api.py

@@ -24,7 +24,6 @@ from fwk.coloring import ccolors
 import openmdao.api as om
 import blast
 import blast.blUtils as vutils
-from pygeo import DVGeometryCST
 from scipy.optimize import least_squares
 import fwk
 
@@ -79,7 +78,7 @@ class BlasterSolver(om.ExplicitComponent):
         self.coupler.run(write=False)
         self.__write(sfx='baseline')
         self.coupler.reset()
-        print('object created')
+        print('BLASTER object created')
     
     def setup(self):
         self.add_input('aoa', val=0.0, desc='Flow angle of attack')
@@ -205,6 +204,7 @@ class BlasterSolver(om.ExplicitComponent):
     def getBoundaryFFDGeometry(self):
         """Airfoil FFD geometry
         """
+        from pygeo import DVGeometryCST
         with open('surface_DVGeo.dat', 'w') as f:
             for i, node_coord in enumerate(self.isol.vBnd[0][0].nodesCoord):
                 f.write(f"{node_coord[0]} {node_coord[1]}\n")