2D Pitching Fixed Wing

amr-wind/wind_energy/actuator/wing/ActuatorWing.H

@@ -37,6 +37,15 @@ struct WingBaseData
     //! Pitch angle for the wing
     amrex::Real pitch{0.0};
 
+    //! File name for pitch actuation table
+    std::string pitch_timetable_file;
+
+    //! Pitch actuation: time table
+    amrex::Vector<amrex::Real> time_table;
+
+    //! Pitch actuation: pitch angle table
+    amrex::Vector<amrex::Real> pitch_table;
+
     //! Relative velocity at the actuator node
     VecList vel_rel;
 

amr-wind/wind_energy/actuator/wing/fixed_wing_ops.H

@@ -46,6 +46,24 @@ struct ReadInputsOp<FixedWing, ActSrcLine>
         bool use_fllc = false;
         pp.query("fllc", use_fllc);
 
+        // Initialize tables for pitch actuation (in degrees)
+        pp.query("pitch_timetable", wdata.pitch_timetable_file);
+        if (!wdata.pitch_timetable_file.empty()) {
+            std::ifstream ifh(wdata.pitch_timetable_file, std::ios::in);
+            if (!ifh.good()) {
+                amrex::Abort(
+                    "Cannot find input file: " + wdata.pitch_timetable_file);
+            }
+            amrex::Real data_time;
+            amrex::Real data_pitch_deg;
+            ifh.ignore(std::numeric_limits<std::streamsize>::max(), '\n');
+            while (ifh >> data_time) {
+                ifh >> data_pitch_deg;
+                wdata.time_table.push_back(data_time);
+                wdata.pitch_table.push_back(data_pitch_deg);
+            }
+        }
+
         pp.query("motion_type", wdata.motion_type);
         {
             // do nothing for default, "none"

amr-wind/wind_energy/actuator/wing/wing_ops.H

@@ -7,6 +7,7 @@
 #include "amr-wind/wind_energy/actuator/aero/AirfoilTable.H"
 #include "amr-wind/wind_energy/actuator/ActParser.H"
 #include "amr-wind/wind_energy/actuator/FLLCOp.H"
+#include "amr-wind/utilities/linear_interpolation.H"
 
 namespace amr_wind::actuator {
 namespace wing {
@@ -164,8 +165,19 @@ struct ComputeForceOp<
             windvector[2] = (grid.vel[ip] - wdata.vel_tr) & blade_z;
 
             const auto vmag = vs::mag(windvector);
+
+            // Calculate the aoa based on the pitch actuation table if supplied
+            // (interpolate pitch angle at step n) or the initial pitch angle
+            amrex::Real aoa;
+            if (!wdata.pitch_timetable_file.empty()) {
+                wdata.pitch = amr_wind::interp::linear(
+                    wdata.time_table, wdata.pitch_table, time.current_time());
+                amrex::Print() << "Wing Ops: Current wing pitch is "
+                               << wdata.pitch << " degrees.\n";
+
             const auto aoa = std::atan2(windvector[2], windvector[0]) +
-                             amr_wind::utils::radians(wdata.pitch);
+                  amr_wind::utils::radians(wdata.pitch);
+
 
             // Make up some Cl, Cd values
             amrex::Real cl, cd;