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

@@ -179,9 +179,9 @@ void Body::setNodes(std::vector<std::vector<double>> const &coords)
         throw std::runtime_error("sdpm::Problem: size of coordinates vector must be equal to the number of nodes!\n");
     // Set nodes coordinates
     for (size_t i = 0; i < _nodes.size(); ++i)
-        _nodes[i]->setCoords(sdpmVector3d(coords[i].data()));
+        _nodes[i]->setCoords(sdpmVector3d(coords[i][0], coords[i][1], coords[i][2]));
     // Update wake
-    double dx = _wake->getLagMap().at(_wake->getElements()[0]); // dx = dt of first cell
+    sdpmDouble dx = _wake->getLagMap().at(_wake->getElements()[0]); // dx = dt of first cell
     size_t nTe = _teNodes.size();
     std::vector<Node *> wkNodes = _wake->getNodes();
     for (size_t i = 1; i < wkNodes.size() / nTe; ++i)