Add new solver and options for torque tube as a spring

.github/workflows/test_pvade.yaml

@@ -33,6 +33,8 @@ jobs:
         shell: bash -l {0}
         run: PYTHONPATH=. pytest -sv pvade/tests/
       - name: test all inputs
+        env:
+          OMPI_MCA_regx: "naive"
         shell: bash -l {0}
         run: pytest -sv  test_all_inputs.py 
 

examples/cylinderflow.py

@@ -45,7 +45,6 @@
 from dolfinx.io import VTXWriter, gmshio, XDMFFile
 from dolfinx.mesh import locate_entities_boundary
 
-
 ####################################################
 #                                                  #
 #           MESH PARAMETERS                        #

examples/poissoneq.py

@@ -10,7 +10,6 @@
 from timeit import default_timer as timer
 from dolfinx.common import TimingType, list_timings
 
-
 start = timer()
 elems = int(sys.argv[1])  # nelems
 # Create mesh and define function space

examples/synthetic_turbulent_inflow/generate_turbulent_inflow_h5_file.ipynb

@@ -24,12 +24,18 @@
     "import pandas as pd  # need this to load our data from the csv files\n",
     "\n",
     "import sys\n",
+    "\n",
     "# NOTE: to run this, you will need to replace the path below #\n",
-    "sys.path.append('/Users/bstanisl/repos/pyconturb/pyconturb') # path to the pyconturb repository location on your local machine\n",
+    "sys.path.append(\n",
+    "    \"/Users/bstanisl/repos/pyconturb/pyconturb\"\n",
+    ")  # path to the pyconturb repository location on your local machine\n",
     "from pyconturb import gen_turb, gen_spat_grid  # generate turbulence, useful helper\n",
     "from pyconturb.sig_models import iec_sig  # IEC 61400-1 turbulence std dev\n",
     "from pyconturb.spectral_models import kaimal_spectrum  # Kaimal spectrum\n",
-    "from pyconturb.wind_profiles import constant_profile, power_profile  # wind-speed profile functions\n",
+    "from pyconturb.wind_profiles import (\n",
+    "    constant_profile,\n",
+    "    power_profile,\n",
+    ")  # wind-speed profile functions\n",
     "\n",
     "from _nb_utils import plot_slice\n",
     "import h5py"
@@ -217,11 +223,13 @@
    "source": [
     "# generate spatial gridpoints dataframe\n",
     "ny = 80\n",
-    "nz = 80 #40\n",
+    "nz = 80  # 40\n",
     "y = np.linspace(-10.0, 10.0, ny)\n",
-    "z = np.linspace(0.00001, 20.0, nz) \n",
+    "z = np.linspace(0.00001, 20.0, nz)\n",
     "\n",
-    "spat_df = gen_spat_grid(y, z)  # if `comps` not passed in, assumes all 3 components are wanted\n",
+    "spat_df = gen_spat_grid(\n",
+    "    y, z\n",
+    ")  # if `comps` not passed in, assumes all 3 components are wanted\n",
     "spat_df.head()  # look at the first few rows"
    ]
   },
@@ -233,9 +241,9 @@
    "outputs": [],
    "source": [
     "# generate time vector\n",
-    "t_final = 1.0 #10.0 [s]\n",
-    "dt = 0.01 # [s]\n",
-    "t_steps = int(t_final/dt)\n",
+    "t_final = 1.0  # 10.0 [s]\n",
+    "dt = 0.01  # [s]\n",
+    "t_steps = int(t_final / dt)\n",
     "time = np.linspace(0, t_final, t_steps)"
    ]
   },
@@ -483,16 +491,30 @@
    "source": [
     "# reshape to 3D array and visualize first timestep\n",
     "data = {}\n",
-    "data['u'] = turb_df.filter(regex='u').values.reshape(len(turb_df),y.size,z.size).transpose((0, 2, 1))\n",
-    "data['v'] = turb_df.filter(regex='v').values.reshape(len(turb_df),y.size,z.size).transpose((0, 2, 1))\n",
-    "data['w'] = turb_df.filter(regex='w').values.reshape(len(turb_df),y.size,z.size).transpose((0, 2, 1))\n",
+    "data[\"u\"] = (\n",
+    "    turb_df.filter(regex=\"u\")\n",
+    "    .values.reshape(len(turb_df), y.size, z.size)\n",
+    "    .transpose((0, 2, 1))\n",
+    ")\n",
+    "data[\"v\"] = (\n",
+    "    turb_df.filter(regex=\"v\")\n",
+    "    .values.reshape(len(turb_df), y.size, z.size)\n",
+    "    .transpose((0, 2, 1))\n",
+    ")\n",
+    "data[\"w\"] = (\n",
+    "    turb_df.filter(regex=\"w\")\n",
+    "    .values.reshape(len(turb_df), y.size, z.size)\n",
+    "    .transpose((0, 2, 1))\n",
+    ")\n",
     "\n",
     "fig, ax = plt.subplots()\n",
-    "plt.imshow(data['u'][0,:,:],  # imshow requires nz-ny slice\n",
-    "            origin='lower',  # smallest y-z in lower left, not upper left\n",
-    "            extent=[y[0], y[-1], z[0], z[-1]],  # lateral and vertical limits\n",
-    "            interpolation='bilinear',\n",
-    "            cmap='coolwarm')  # image smoothing\n",
+    "plt.imshow(\n",
+    "    data[\"u\"][0, :, :],  # imshow requires nz-ny slice\n",
+    "    origin=\"lower\",  # smallest y-z in lower left, not upper left\n",
+    "    extent=[y[0], y[-1], z[0], z[-1]],  # lateral and vertical limits\n",
+    "    interpolation=\"bilinear\",\n",
+    "    cmap=\"coolwarm\",\n",
+    ")  # image smoothing\n",
     "plt.colorbar()"
    ]
   },
@@ -507,24 +529,24 @@
     "h5_filename = \"pct_turb_ny{}_nz{}_unconstrained_{}s_dt{}_uref{}.h5\".format(\n",
     "    ny, nz, t_final, dt, u_ref\n",
     ")\n",
-    "with h5py.File(\"../../input/\"+h5_filename, \"w\") as fp:\n",
+    "with h5py.File(\"../../input/\" + h5_filename, \"w\") as fp:\n",
     "    fp.create_dataset(\"time_index\", shape=(t_steps,))\n",
     "    fp[\"time_index\"][:] = time\n",
-    "    \n",
+    "\n",
     "    fp.create_dataset(\"y_coordinates\", shape=(ny,))\n",
     "    fp[\"y_coordinates\"][:] = y\n",
-    "        \n",
+    "\n",
     "    fp.create_dataset(\"z_coordinates\", shape=(nz,))\n",
     "    fp[\"z_coordinates\"][:] = z\n",
-    "        \n",
+    "\n",
     "    fp.create_dataset(\"u\", shape=(t_steps, nz, ny))\n",
-    "    fp[\"u\"][:] = data['u'][:]\n",
-    "    \n",
+    "    fp[\"u\"][:] = data[\"u\"][:]\n",
+    "\n",
     "    fp.create_dataset(\"v\", shape=(t_steps, nz, ny))\n",
-    "    fp[\"v\"][:] = data['v'][:]\n",
-    "    \n",
+    "    fp[\"v\"][:] = data[\"v\"][:]\n",
+    "\n",
     "    fp.create_dataset(\"w\", shape=(t_steps, nz, ny))\n",
-    "    fp[\"w\"][:] = data['w'][:]"
+    "    fp[\"w\"][:] = data[\"w\"][:]"
    ]
   },
   {
@@ -551,7 +573,7 @@
    "source": [
     "# read h5 file\n",
     "# Note: Replace the path below with the path on your local machine to the .h5 file that you want to read\n",
-    "fname = '/Users/bstanisl/Documents/repos/PVade/input/pct_turb_ny80_nz80_unconstrained_1.0s_dt0.01_uref20.h5'\n",
+    "fname = \"/Users/bstanisl/Documents/repos/PVade/input/pct_turb_ny80_nz80_unconstrained_1.0s_dt0.01_uref20.h5\"\n",
     "\n",
     "# Open the file (read-only)\n",
     "with h5py.File(fname, \"r\") as f:\n",

input/duramat_case_study.yaml

@@ -23,6 +23,13 @@ pv_array:
   elevation: 2.1
   tracker_angle: 0.0
   span_fixation_pts: [13.2]
+  # torque_tube_separation: 0.2 # gap between panel and tube center
+  # torque_tube_outer_radius: 0.1 # radius of the torque tube
+  # torque_tube_inner_radius: 0.09 # radius of the torque tube
+  modules_per_span: 1
+  fixed_location: 1   # fixed location of the panel along the span (from 0 (fixed at left) to modules_per_span (fixed at right))
+  block_chord_div_by_panel_chord: 0.02
+
 solver:
   dt: 0.01
   t_final: 20.0
@@ -58,4 +65,7 @@ structure:
   bc_list: []
   motor_connection: true
   tube_connection: true
-  beta_relaxation: 0.5
+  beta_relaxation: 0.5
+  elasticity_modulus_tube: 2.0e+11
+  poissons_ratio_tube: 0.3
+  rho_tube: 7800.0

input/flag2d.yaml

@@ -52,7 +52,7 @@ structure:
   body_force_x: 0
   body_force_y: 0
   body_force_z: 0 #100
-  bc_list: ["left"]
+  bc_list: ["panel_left"]
   motor_connection: False
   tube_connection: False
   beta_relaxation: 0.005

input/heated_panels2d.yaml

@@ -59,7 +59,7 @@ structure:
     body_force_x: 0
     body_force_y: -1
     body_force_z: 0 #100
-    bc_list: ["left"]
+    bc_list: ["panel_left"]
     motor_connection: False
     tube_connection: False
     beta_relaxation: 0.005

input/inflow_input.yaml

@@ -56,7 +56,7 @@ structure:
   body_force_x: 0
   body_force_y: 0
   body_force_z: 0 #100
-  bc_list: ["left"]
+  bc_list: ["panel_left"]
   motor_connection: False
   tube_connection: False
   beta_relaxation: 0.005

input/panels2d.yaml

@@ -46,7 +46,7 @@ structure:
     body_force_x: 0
     body_force_y: -1
     body_force_z: 0 #100
-    bc_list: ["top"]
+    bc_list: [panel_top]
     motor_connection: False
     tube_connection: False
     beta_relaxation: 0.005

input/panels3d.yaml

@@ -56,5 +56,5 @@ structure:
   body_force_x: 0
   body_force_y: 0
   body_force_z: -1 #100
-  bc_list: [left ]
+  bc_list: [panel_left]
   tube_connection: False

pvade/IO/input_schema.yaml

@@ -136,7 +136,7 @@ properties:
         minimum: 0.0
         maximum: 10.0
         type: "number"
-        description: "The vertical distance between the center of the panel and the ground."
+        description: "The vertical distance between the center of the torque tube and the ground."
         units: "meter"
       stream_spacing:
         default: 7.0
@@ -154,7 +154,7 @@ properties:
         units: "meter"
       span_fixation_pts:
         default: 3.5
-        minimum: 1.0
+        minimum: 0.2
         # maximum: 32.0
         type:
           - "number"
@@ -198,6 +198,39 @@ properties:
         be unrolled such that tracking_angle_0 is applied to the panel at (x_min, y_min), tracking_angle_1 is applied to the panel at
         (x_min + x_spacing, y_min), etc., i.e., x-major direction. If argument is a single value, that tracking angle will be applied to the every panel."
         units: "degree"
+      torque_tube_separation:
+        default: 0.0
+        minimum: 0.0
+        maximum: 1.0
+        type: "number"
+        description: "The distance between the bottom of the panel structure and the centerline of the torque tube, in meters."
+      torque_tube_outer_radius:
+        default: 0.0
+        minimum: 0.0
+        maximum: 1.0
+        type: "number"
+        description: "The radius of the torque tube, in meters."
+      torque_tube_inner_radius:
+        default: 0.0
+        minimum: 0.0
+        maximum: 1.0
+        type: "number"
+        description: "The radius of the torque tube, in meters."
+      modules_per_span:
+        default: 1
+        minimum: 1
+        type: "integer"
+        description: "The number of modules/panels per every panel_span distance. This can be thought of as the number of modules/panels per each table, counted in the spanwise direction only (i.e., 2-in-portait systems should report only the number in the spanwise direction, vs 2x that value)."
+      fixed_location:
+        default: 5
+        type: "integer"
+        description: "If an integer between 0 and modules_per_span, indicates the fixed location"
+      block_chord_div_by_panel_chord:
+        default: 0.1
+        minimum: 0.01
+        type: "number"
+        description: "block length or width divided by panel chord"
+
   solver:
     additionalProperties: false
     type: "object"
@@ -617,7 +650,7 @@ properties:
         type: "boolean"
         description: "Controls whether or not a boundary condition is applied along the lines defined by the spanwise fixation points which run in the streamwise direction and divide the panel into multiple rectangles in the spanwise direction."
       bc_list:
-        default: [front,back]
+        default: []#[panel_front,panel_back]
         # type: "list"
         description: "location for clamped boundary conditions"
       dt:
@@ -648,6 +681,27 @@ properties:
         type: "number"
         description: "poissons ratio of the panel structure."
         units: "None"
+      elasticity_modulus_tube:
+        default: 200.0e+9
+        minimum: 1.0e+2
+        maximum: 500.0e+9
+        type: "number"
+        description: "The effective Young's modulus of the torque tube."
+        units: "Pascal"  
+      poissons_ratio_tube: 
+        default: 0.3
+        minimum: 0.1
+        maximum: 0.9
+        type: "number"
+        description: "poissons ratio of the torque tube."
+        units: "None"
+      rho_tube:
+        default: 1.0
+        # minimum: 0.001
+        # maximum: 1000.0
+        type: "number"
+        description: "The density of the structure."
+        units: "kg/meter^3"
       body_force_x: 
         default:  100
         type: "number"