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/sdpmBody.cpp

@@ -157,6 +157,8 @@ Body::Body(Mesh &msh, std::string const &name, std::string const &teName,
     _nLoads.resize(_nodes.size(), sdpmVector3d::Zero());
     _nLoads1Re.resize(nM, std::vector<std::vector<sdpmVector3d>>(nF, std::vector<sdpmVector3d>(_nodes.size(), sdpmVector3d::Zero())));
     _nLoads1Im.resize(nM, std::vector<std::vector<sdpmVector3d>>(nF, std::vector<sdpmVector3d>(_nodes.size(), sdpmVector3d::Zero())));
+    _iVelocityS.resize(nM, std::vector<sdpmVector3d>(getElements().size(), sdpmVector3d::Zero()));
+    _iVelocityH.resize(nM, std::vector<sdpmVector3d>(getElements().size(), sdpmVector3d::Zero()));
 }
 
 Body::~Body()
@@ -233,8 +235,8 @@ void Body::computeVelocity(size_t m, sdpmVector3d const &ui)
         ucS[i] = _motions[m]->computeSteady(*elems[i], ui);
         ucH[i] = _motions[m]->computeHarmonic(*elems[i]);
     }
-    _iVelocityS.push_back(ucS);
-    _iVelocityH.push_back(ucH);
+    _iVelocityS[m] = ucS;
+    _iVelocityH[m] = ucH;
 }
 
 /**