Version 1.2

Version 1.2 - Transonic correction and analytical gradients

This version mainly features a transonic solver and analytical pressure derivatives.

Changes

• Add transonic correction
• Add analytical gradients
• Add feature to deform the body
• Fix calculation of motion-induced velocity and unsteady pressure
• Refactor APIs and utilities

Notes

Transonic correction
The unsteady pressure solution can now be corrected by higher-fidelity data (typically obtained from CFD) to include transonic and viscous effects.
Analytical gradients
Although CoDiPack works well, its memory footprint is too large for practical unsteady cases (see #6). Analytical pressure and aerodynamic forces gradients have therefore been implemented so that flutter constrained optimization can be performed within a reasonable computational time.
Body deformation
The body nodes can now be set to new positions for aeroelastic and optimization purposes.
Formulation fixes
The calculation of the surface velocity induced by a rigid pitching motion has been corrected to use original instead of compressible coordinates. Additionally, a mistake in the calculation of the unsteady pressure coefficient has been fixed.
Refactoring
The default values in the API have been refactored. The slicing utility now also outputs z-coordinates.

Tests

• Ubuntu 22.04 (gcc 11, python 3.10)

sdpm/src/sdpmWake.h

@@ -31,6 +31,7 @@ namespace sdpm
  */
 class SDPM_API Wake : public Group
 {
+    std::vector<Node *> _nodes;                                  ///< nodes of the wake
     std::map<Element *, std::pair<Element *, Element *>> _wwMap; ///< wing-wake map
     std::map<Element *, sdpmDouble> _lagMap;                     ///< time lag map
 
@@ -39,10 +40,11 @@ public:
     ~Wake() {}
 
 #ifndef SWIG
-    std::map<Element *, std::pair<Element *, Element *>> const &getElMap() const
+    std::vector<Node *> const &getNodes() const
     {
-        return _wwMap;
+        return _nodes;
     }
+    std::map<Element *, std::pair<Element *, Element *>> const &getElMap() const { return _wwMap; }
     std::map<Element *, sdpmDouble> const &getLagMap() const { return _lagMap; }
     virtual void write(std::ostream &out) const override;
 #endif