dart/api/core.py

@@ -203,21 +203,30 @@ def init_dart(cfg, scenario='aerodynamic', task='analysis', viscous=False):
                 _pbl.add(dart.Wake(_msh, [cfg['Wakes'][i], cfg['Wakes'][i]+'_', cfg['Fluid']])) # 2.5D
             _pbl.add(dart.Kutta(_msh, [cfg['Tes'][i], cfg['Wakes'][i]+'_', cfg['Wings'][i], cfg['Fluid']]))
     # add transpiration (blowing) boundary conditions
+    _blw = {}
     if viscous:
         if _dim == 2:
-            _blwb = dart.Blowing(_msh, cfg['Wing'])
-            _blww = dart.Blowing(_msh, cfg['Wake'])
+            _blwb = dart.Blowing(_msh, cfg["Wing"])
+            _blww = dart.Blowing(_msh, cfg["Wake"])
             _pbl.add(_blwb)
             _pbl.add(_blww)
-            _blw = [[_blwb, _blww]]
+            print(f'Adding blowing boundary {cfg["Wing"]}')
+            _blw[cfg['Wing']] = {'wing':_blwb,
+                                 'wake':_blww}
         else:
-            _blw = []
             for ibody in range(len(cfg['Wings'])):
-                _blwb = dart.Blowing(_msh, cfg['Wings'][ibody])
-                _blww = dart.Blowing(_msh, cfg['Wakes'][ibody])
+                print(f'Adding blowing boundary {cfg["Wings"][ibody]}')
+                _blwb = dart.Blowing(_msh, cfg["Wings"][ibody])
+                _blww = dart.Blowing(_msh, cfg["Wakes"][ibody])
                 _pbl.add(_blwb)
                 _pbl.add(_blww)
-                _blw.append([_blwb, _blww])
+                _blw[cfg["Wings"][ibody]] = {'wing':_blwb,
+                                    'wake':_blww}
+            if 'Fuselage' in cfg:
+                print(f'Adding blowing boundary {cfg["Fuselage"]}')
+                _blwf = dart.Blowing(_msh, cfg["Fuselage"])
+                _pbl.add(_blwf)
+                _blw[cfg["Fuselage"]] = {'fuselage':_blwf}
     else:
         _blw = None
 

dart/src/wSolver.h

@@ -81,6 +81,8 @@ public:
     void save(tbox::MshExport *mshWriter, std::string const &suffix = "");
 
     int getRow(size_t i) const { return rows[i]; };
+    void setRelTol(double const rtol) { relTol = rtol; }
+    void setAbsTol(double const atol) { absTol = atol; }
 
 protected:
     void computeFlow();