diff --git a/amr-wind/utilities/sampling/PlaneSampler.cpp b/amr-wind/utilities/sampling/PlaneSampler.cpp index 7c16e53ace..f7fb48697f 100644 --- a/amr-wind/utilities/sampling/PlaneSampler.cpp +++ b/amr-wind/utilities/sampling/PlaneSampler.cpp @@ -65,7 +65,7 @@ void PlaneSampler::check_bounds() const amrex::Vector points = { point, point + m_axis1[d], point + m_axis2[d]}; for (const auto& pt : points) { - if ((pt <= prob_lo[d]) || (pt >= prob_hi[d])) { + if ((pt < prob_lo[d]) || (pt >= prob_hi[d])) { amrex::Abort( "PlaneSampler: Point out of domain. Redefine your " "planes so they are completely inside the domain."); diff --git a/docs/sphinx/walkthrough/precursor_inp.txt b/docs/sphinx/walkthrough/precursor_inp.txt index 65c8dbb4b4..072d5f5d92 100644 --- a/docs/sphinx/walkthrough/precursor_inp.txt +++ b/docs/sphinx/walkthrough/precursor_inp.txt @@ -30,7 +30,7 @@ transport.turbulent_prandtl = 0.3333 geometry.prob_lo = 0. 0. 0. # Lo corner coordinates geometry.prob_hi = 2560. 2560. 1280. # Hi corner coordinates amr.n_cell = 128 128 64 # Grid cells at coarsest AMRlevel -amr.max_level = 0 # Max AMR level in hierarchy +amr.max_level = 0 # Max AMR level in hierarchy geometry.is_periodic = 1 1 0 # Periodicity x y z (0/1) zlo.type = wall_model @@ -74,21 +74,21 @@ sampling.output_frequency = 96 sampling.fields = velocity temperature #---- sample defs ---- -sampling.labels = xy-domain xz-domain +sampling.labels = xy-domain xz-domain -sampling.xy-domain.type = PlaneSampler -sampling.xy-domain.num_points = 256 256 -sampling.xy-domain.origin = 0.0 0.0 86.5 -sampling.xy-domain.axis1 = 2550.0 0.0 0.0 -sampling.xy-domain.axis2 = 0.0 2550.0 0.0 -sampling.xy-domain.offset_vector = 0.0 0.0 1.0 -sampling.xy-domain.offsets = -63.45 0.0 63.45 +sampling.xy-domain.type = PlaneSampler +sampling.xy-domain.num_points = 256 256 +sampling.xy-domain.origin = 0.0 0.0 86.5 +sampling.xy-domain.axis1 = 2550.0 0.0 0.0 +sampling.xy-domain.axis2 = 0.0 2550.0 0.0 +sampling.xy-domain.offset_vector = 0.0 0.0 1.0 +sampling.xy-domain.offsets = -63.45 0.0 63.45 -sampling.xz-domain.type = PlaneSampler -sampling.xz-domain.num_points = 256 128 -sampling.xz-domain.origin = 0.0 1280.0 0.0 -sampling.xz-domain.axis1 = 2550.0 0.0 0.0 -sampling.xz-domain.axis2 = 0.0 0.0 1270.0 +sampling.xz-domain.type = PlaneSampler +sampling.xz-domain.num_points = 256 128 +sampling.xz-domain.origin = 0.0 1280.0 0.0 +sampling.xz-domain.axis1 = 2550.0 0.0 0.0 +sampling.xz-domain.axis2 = 0.0 0.0 1270.0 #¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# # AVERAGING # diff --git a/docs/sphinx/walkthrough/spinup_inp.txt b/docs/sphinx/walkthrough/spinup_inp.txt index bf42e66e49..358a3f3808 100644 --- a/docs/sphinx/walkthrough/spinup_inp.txt +++ b/docs/sphinx/walkthrough/spinup_inp.txt @@ -11,7 +11,7 @@ time.checkpoint_interval = 1800 # Steps between checkpoi incflo.physics = ABL turbulence.model = OneEqKsgsM84 # For neutral ABL, use "Smagorinsky" -TKE.source_terms = KsgsM84Src +TKE.source_terms = KsgsM84Src incflo.gravity = 0. 0. -9.81 # Gravitational force (3D) incflo.density = 1.225 # Reference density; make sure this agrees with OpenFAST values transport.viscosity = 1.0e-5 # Dynamic viscosity [N-s/m^2] @@ -24,7 +24,7 @@ transport.turbulent_prandtl = 0.3333 geometry.prob_lo = 0. 0. 0. # Lo corner coordinates geometry.prob_hi = 2560. 2560. 1280. # Hi corner coordinates amr.n_cell = 128 128 64 # Grid cells at coarsest AMRlevel -amr.max_level = 0 # Max AMR level in hierarchy +amr.max_level = 0 # Max AMR level in hierarchy geometry.is_periodic = 1 1 0 # Periodicity x y z (0/1) zlo.type = wall_model @@ -63,22 +63,22 @@ ABL.surface_temp_flux = 0.05 # Surface temperature flux [K-m incflo.post_processing = sampling # --- Sampling parameters --- -sampling.output_frequency = 900 +sampling.output_frequency = 900 sampling.fields = velocity temperature #---- sample defs ---- -sampling.labels = xy-domain xz-domain +sampling.labels = xy-domain xz-domain -sampling.xy-domain.type = PlaneSampler -sampling.xy-domain.num_points = 256 256 -sampling.xy-domain.origin = 0.0 0.0 86.5 -sampling.xy-domain.axis1 = 2550.0 0.0 0.0 -sampling.xy-domain.axis2 = 0.0 2550.0 0.0 -sampling.xy-domain.offset_vector = 0.0 0.0 1.0 -sampling.xy-domain.offsets = -63.45 0.0 63.45 +sampling.xy-domain.type = PlaneSampler +sampling.xy-domain.num_points = 256 256 +sampling.xy-domain.origin = 0.0 0.0 86.5 +sampling.xy-domain.axis1 = 2550.0 0.0 0.0 +sampling.xy-domain.axis2 = 0.0 2550.0 0.0 +sampling.xy-domain.offset_vector = 0.0 0.0 1.0 +sampling.xy-domain.offsets = -63.45 0.0 63.45 -sampling.xz-domain.type = PlaneSampler -sampling.xz-domain.num_points = 256 128 -sampling.xz-domain.origin = 0.0 1280.0 0.0 -sampling.xz-domain.axis1 = 2550.0 0.0 0.0 -sampling.xz-domain.axis2 = 0.0 0.0 1270.0 +sampling.xz-domain.type = PlaneSampler +sampling.xz-domain.num_points = 256 128 +sampling.xz-domain.origin = 0.0 1280.0 0.0 +sampling.xz-domain.axis1 = 2550.0 0.0 0.0 +sampling.xz-domain.axis2 = 0.0 0.0 1270.0 diff --git a/docs/sphinx/walkthrough/turbines_inp.txt b/docs/sphinx/walkthrough/turbines_inp.txt index 360819ad67..df53f06b46 100644 --- a/docs/sphinx/walkthrough/turbines_inp.txt +++ b/docs/sphinx/walkthrough/turbines_inp.txt @@ -15,10 +15,10 @@ ABL.bndry_output_start_time = 7200.0 ABL.bndry_var_names = velocity temperature tke incflo.physics = ABL Actuator -io.restart_file = ../spinup/chk14400 +io.restart_file = ../spinup/chk14400 io.outputs = "actuator_src_term" turbulence.model = OneEqKsgsM84 # For neutral ABL, use "Smagorinsky" -TKE.source_terms = KsgsM84Src +TKE.source_terms = KsgsM84Src incflo.gravity = 0. 0. -9.81 # Gravitational force (3D) incflo.density = 1.225 # Reference density; make sure this agrees with OpenFAST values transport.viscosity = 1.0e-5 # Dynamic viscosity [N-s/m^2] @@ -31,17 +31,17 @@ transport.turbulent_prandtl = 0.3333 geometry.prob_lo = 0. 0. 0. # Lo corner coordinates geometry.prob_hi = 2560. 2560. 1280. # Hi corner coordinates amr.n_cell = 128 128 64 # Grid cells at coarsest AMRlevel -amr.max_level = 2 # Max AMR level in hierarchy +amr.max_level = 2 # Max AMR level in hierarchy geometry.is_periodic = 0 0 0 # Periodicity x y z (0/1) -xlo.type = mass_inflow -xlo.density = 1.225 -xlo.temperature = 290.0 +xlo.type = mass_inflow +xlo.density = 1.225 +xlo.temperature = 290.0 xlo.tke = 0.0 -xhi.type = pressure_outflow +xhi.type = pressure_outflow -ylo.type = slip_wall -yhi.type = slip_wall +ylo.type = slip_wall +yhi.type = slip_wall zlo.type = wall_model zhi.type = slip_wall @@ -94,21 +94,21 @@ sampling.output_frequency = 100 sampling.fields = velocity temperature #---- sample defs ---- -sampling.labels = xy-domain xz-domain - -sampling.xy-domain.type = PlaneSampler -sampling.xy-domain.num_points = 256 256 -sampling.xy-domain.origin = 0.0 0.0 91.0 -sampling.xy-domain.axis1 = 2550.0 0.0 0.0 -sampling.xy-domain.axis2 = 0.0 2550.0 0.0 -sampling.xy-domain.offset_vector = 0.0 0.0 1.0 -sampling.xy-domain.offsets = -63.45 0.0 63.45 - -sampling.xz-domain.type = PlaneSampler -sampling.xz-domain.num_points = 256 128 -sampling.xz-domain.origin = 0.0 1280.0 0.0 -sampling.xz-domain.axis1 = 2550.0 0.0 0.0 -sampling.xz-domain.axis2 = 0.0 0.0 1270.0 +sampling.labels = xy-domain xz-domain + +sampling.xy-domain.type = PlaneSampler +sampling.xy-domain.num_points = 256 256 +sampling.xy-domain.origin = 0.0 0.0 91.0 +sampling.xy-domain.axis1 = 2550.0 0.0 0.0 +sampling.xy-domain.axis2 = 0.0 2550.0 0.0 +sampling.xy-domain.offset_vector = 0.0 0.0 1.0 +sampling.xy-domain.offsets = -63.45 0.0 63.45 + +sampling.xz-domain.type = PlaneSampler +sampling.xz-domain.num_points = 256 128 +sampling.xz-domain.origin = 0.0 1280.0 0.0 +sampling.xz-domain.axis1 = 2550.0 0.0 0.0 +sampling.xz-domain.axis2 = 0.0 0.0 1270.0 #¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# # AVERAGING # @@ -131,56 +131,56 @@ averaging.stress.averaging_type = ReynoldsStress tagging.labels = T0_level_0_zone T1_level_0_zone T2_level_0_zone T0_level_1_zone T1_level_1_zone T2_level_1_zone # 1st refinement level -tagging.T0_level_0_zone.type = GeometryRefinement -tagging.T0_level_0_zone.shapes = T0_level_0_zone -tagging.T0_level_0_zone.level = 0 -tagging.T0_level_0_zone.T0_level_0_zone.type = box +tagging.T0_level_0_zone.type = GeometryRefinement +tagging.T0_level_0_zone.shapes = T0_level_0_zone +tagging.T0_level_0_zone.level = 0 +tagging.T0_level_0_zone.T0_level_0_zone.type = box tagging.T0_level_0_zone.T0_level_0_zone.origin = 520.0 1040.0 0.0 # -1D, -2D tagging.T0_level_0_zone.T0_level_0_zone.xaxis = 360.0 0.0 0.0 tagging.T0_level_0_zone.T0_level_0_zone.yaxis = 0.0 480.0 0.0 tagging.T0_level_0_zone.T0_level_0_zone.zaxis = 0.0 0.0 360.0 -tagging.T1_level_0_zone.type = GeometryRefinement -tagging.T1_level_0_zone.shapes = T1_level_0_zone -tagging.T1_level_0_zone.level = 0 -tagging.T1_level_0_zone.T1_level_0_zone.type = box +tagging.T1_level_0_zone.type = GeometryRefinement +tagging.T1_level_0_zone.shapes = T1_level_0_zone +tagging.T1_level_0_zone.level = 0 +tagging.T1_level_0_zone.T1_level_0_zone.type = box tagging.T1_level_0_zone.T1_level_0_zone.origin = 1160.0 1040.0 0.0 # -1D, -2D tagging.T1_level_0_zone.T1_level_0_zone.xaxis = 360.0 0.0 0.0 tagging.T1_level_0_zone.T1_level_0_zone.yaxis = 0.0 480.0 0.0 tagging.T1_level_0_zone.T1_level_0_zone.zaxis = 0.0 0.0 360.0 -tagging.T2_level_0_zone.type = GeometryRefinement -tagging.T2_level_0_zone.shapes = T2_level_0_zone -tagging.T2_level_0_zone.level = 0 -tagging.T2_level_0_zone.T2_level_0_zone.type = box +tagging.T2_level_0_zone.type = GeometryRefinement +tagging.T2_level_0_zone.shapes = T2_level_0_zone +tagging.T2_level_0_zone.level = 0 +tagging.T2_level_0_zone.T2_level_0_zone.type = box tagging.T2_level_0_zone.T2_level_0_zone.origin = 1800.0 1040.0 0.0 # -1D, -2D tagging.T2_level_0_zone.T2_level_0_zone.xaxis = 360.0 0.0 0.0 tagging.T2_level_0_zone.T2_level_0_zone.yaxis = 0.0 480.0 0.0 tagging.T2_level_0_zone.T2_level_0_zone.zaxis = 0.0 0.0 360.0 # 2nd refinement level -tagging.T0_level_1_zone.type = GeometryRefinement -tagging.T0_level_1_zone.shapes = T0_level_1_zone -tagging.T0_level_1_zone.level = 1 -tagging.T0_level_1_zone.T0_level_1_zone.type = box +tagging.T0_level_1_zone.type = GeometryRefinement +tagging.T0_level_1_zone.shapes = T0_level_1_zone +tagging.T0_level_1_zone.level = 1 +tagging.T0_level_1_zone.T0_level_1_zone.type = box tagging.T0_level_1_zone.T0_level_1_zone.origin = 580.0 1100.0 20.0 # -0.5D, -1.5D tagging.T0_level_1_zone.T0_level_1_zone.xaxis = 180.0 0.0 0.0 tagging.T0_level_1_zone.T0_level_1_zone.yaxis = 0.0 360.0 0.0 tagging.T0_level_1_zone.T0_level_1_zone.zaxis = 0.0 0.0 180.0 -tagging.T1_level_1_zone.type = GeometryRefinement -tagging.T1_level_1_zone.shapes = T1_level_1_zone -tagging.T1_level_1_zone.level = 1 -tagging.T1_level_1_zone.T1_level_1_zone.type = box +tagging.T1_level_1_zone.type = GeometryRefinement +tagging.T1_level_1_zone.shapes = T1_level_1_zone +tagging.T1_level_1_zone.level = 1 +tagging.T1_level_1_zone.T1_level_1_zone.type = box tagging.T1_level_1_zone.T1_level_1_zone.origin = 1220.0 1100.0 20.0 # -0.5D, -1.5D tagging.T1_level_1_zone.T1_level_1_zone.xaxis = 180.0 0.0 0.0 tagging.T1_level_1_zone.T1_level_1_zone.yaxis = 0.0 360.0 0.0 tagging.T1_level_1_zone.T1_level_1_zone.zaxis = 0.0 0.0 180.0 -tagging.T2_level_1_zone.type = GeometryRefinement -tagging.T2_level_1_zone.shapes = T2_level_1_zone -tagging.T2_level_1_zone.level = 1 -tagging.T2_level_1_zone.T2_level_1_zone.type = box +tagging.T2_level_1_zone.type = GeometryRefinement +tagging.T2_level_1_zone.shapes = T2_level_1_zone +tagging.T2_level_1_zone.level = 1 +tagging.T2_level_1_zone.T2_level_1_zone.type = box tagging.T2_level_1_zone.T2_level_1_zone.origin = 1860.0 1100.0 20.0 # -0.5D, -1.5D tagging.T2_level_1_zone.T2_level_1_zone.xaxis = 180.0 0.0 0.0 tagging.T2_level_1_zone.T2_level_1_zone.yaxis = 0.0 360.0 0.0 @@ -193,50 +193,50 @@ Actuator.labels = T0 T1 T2 Actuator.TurbineFastDisk.density = 1.225 # AirDens in OpenFAST models -Actuator.T0.type = TurbineFastDisk +Actuator.T0.type = TurbineFastDisk Actuator.T0.openfast_input_file = T0_OpenFAST/NREL-2p8-127.fst -Actuator.T0.base_position = 640.0 1280.0 0.0 -Actuator.T0.rotor_diameter = 126.9 -Actuator.T0.hub_height = 86.5 -Actuator.T0.num_points_blade = 64 -Actuator.T0.num_points_tower = 12 -Actuator.T0.epsilon = 5.0 5.0 5.0 -Actuator.T0.epsilon_tower = 5.0 5.0 5.0 -Actuator.T0.openfast_start_time = 0.0 -Actuator.T0.openfast_stop_time = 99999.0 -Actuator.T0.nacelle_drag_coeff = 0.0 -Actuator.T0.nacelle_area = 0.0 -Actuator.T0.yaw = 0.0 -Actuator.T0.output_frequency = 10 - -Actuator.T1.type = TurbineFastDisk +Actuator.T0.base_position = 640.0 1280.0 0.0 +Actuator.T0.rotor_diameter = 126.9 +Actuator.T0.hub_height = 86.5 +Actuator.T0.num_points_blade = 64 +Actuator.T0.num_points_tower = 12 +Actuator.T0.epsilon = 5.0 5.0 5.0 +Actuator.T0.epsilon_tower = 5.0 5.0 5.0 +Actuator.T0.openfast_start_time = 0.0 +Actuator.T0.openfast_stop_time = 99999.0 +Actuator.T0.nacelle_drag_coeff = 0.0 +Actuator.T0.nacelle_area = 0.0 +Actuator.T0.yaw = 0.0 +Actuator.T0.output_frequency = 10 + +Actuator.T1.type = TurbineFastDisk Actuator.T1.openfast_input_file = T1_OpenFAST/NREL-2p8-127.fst -Actuator.T1.base_position = 1280.0 1280.0 0.0 -Actuator.T1.rotor_diameter = 126.9 -Actuator.T1.hub_height = 86.5 -Actuator.T1.num_points_blade = 64 -Actuator.T1.num_points_tower = 12 -Actuator.T1.epsilon = 5.0 5.0 5.0 -Actuator.T1.epsilon_tower = 5.0 5.0 5.0 -Actuator.T1.openfast_start_time = 0.0 -Actuator.T1.openfast_stop_time = 99999.0 -Actuator.T1.nacelle_drag_coeff = 0.0 -Actuator.T1.nacelle_area = 0.0 -Actuator.T1.yaw = 0.0 -Actuator.T1.output_frequency = 10 - -Actuator.T2.type = TurbineFastDisk +Actuator.T1.base_position = 1280.0 1280.0 0.0 +Actuator.T1.rotor_diameter = 126.9 +Actuator.T1.hub_height = 86.5 +Actuator.T1.num_points_blade = 64 +Actuator.T1.num_points_tower = 12 +Actuator.T1.epsilon = 5.0 5.0 5.0 +Actuator.T1.epsilon_tower = 5.0 5.0 5.0 +Actuator.T1.openfast_start_time = 0.0 +Actuator.T1.openfast_stop_time = 99999.0 +Actuator.T1.nacelle_drag_coeff = 0.0 +Actuator.T1.nacelle_area = 0.0 +Actuator.T1.yaw = 0.0 +Actuator.T1.output_frequency = 10 + +Actuator.T2.type = TurbineFastDisk Actuator.T2.openfast_input_file = T2_OpenFAST/NREL-2p8-127.fst -Actuator.T2.base_position = 1920.0 1280.0 0.0 -Actuator.T2.rotor_diameter = 126.9 -Actuator.T2.hub_height = 86.5 -Actuator.T2.num_points_blade = 64 -Actuator.T2.num_points_tower = 12 -Actuator.T2.epsilon = 5.0 5.0 5.0 -Actuator.T2.epsilon_tower = 5.0 5.0 5.0 -Actuator.T2.openfast_start_time = 0.0 -Actuator.T2.openfast_stop_time = 99999.0 -Actuator.T2.nacelle_drag_coeff = 0.0 -Actuator.T2.nacelle_area = 0.0 -Actuator.T2.yaw = 0.0 -Actuator.T2.output_frequency = 10 +Actuator.T2.base_position = 1920.0 1280.0 0.0 +Actuator.T2.rotor_diameter = 126.9 +Actuator.T2.hub_height = 86.5 +Actuator.T2.num_points_blade = 64 +Actuator.T2.num_points_tower = 12 +Actuator.T2.epsilon = 5.0 5.0 5.0 +Actuator.T2.epsilon_tower = 5.0 5.0 5.0 +Actuator.T2.openfast_start_time = 0.0 +Actuator.T2.openfast_stop_time = 99999.0 +Actuator.T2.nacelle_drag_coeff = 0.0 +Actuator.T2.nacelle_area = 0.0 +Actuator.T2.yaw = 0.0 +Actuator.T2.output_frequency = 10