# --- Simulation time control parameters --- time.stop_time = 28900 # Max (simulated) time to evolve [s] time.max_step = -1 time.fixed_dt = -1 # Fixed timestep size (in seconds). If negative, then time.cfl is used time.cfl = 0.95 incflo.physics = ABL #Actuator # List of physics models to include in simulation. incflo.verbose = 0 time.plot_interval = 4000 # Steps between plot files time.checkpoint_interval = 3000 # Steps between checkpoint files io.plot_file = "/lustre/scratch/asim/Wind_Energy_Cases/Stable/No_Turbine/plt" io.check_file = "/lustre/scratch/asim/Wind_Energy_Cases/Stable/No_Turbine/chk" ABL.bndry_file = /lustre/scratch/asim/Wind_Energy_Cases/Stable/bndry_file.native ABL.bndry_io_mode = 1 # 0 = write, 1 = read ABL.bndry_planes = ylo xlo ABL.bndry_output_start_time = 28000 ABL.bndry_var_names = velocity temperature tke io.restart_file = /lustre/scratch/asim/Wind_Energy_Cases/Stable/chk56000 ######## need to change incflo.use_godunov = true incflo.godunov_type = weno_z turbulence.model = OneEqKsgsM84 TKE.source_terms = KsgsM84Src nodal_proj.mg_rtol = 1e-06 nodal_proj.mg_atol = 1e-12 mac_proj.mg_rtol = 1e-06 mac_proj.mg_atol = 1e-12 diffusion.mg_rtol = 1e-06 diffusion.mg_atol = 1e-12 temperature_diffusion.mg_rtol = 1e-10 temperature_diffusion.mg_atol = 1e-13 incflo.gravity = 0.0 0.0 -9.81 # Gravitational acceleration vector (x,y,z) [m/s^2] incflo.density = 1.225 # Fluid density [kg/m^3] transport.viscosity = 1e-05 # Fluid dynamic viscosity [kg/m-s] transport.laminar_prandtl = 0.7 # Laminar prandtl number transport.turbulent_prandtl = 0.3333 # Turbulent prandtl number # --- Geometry and Mesh --- geometry.prob_lo = 0.0 0.0 0.0 geometry.prob_hi = 2560.0 640.0 1280.0 amr.n_cell = 256 64 128 # Number of cells in x, y, and z directions amr.max_level = 2 geometry.is_periodic = 0 0 0 xlo.type = mass_inflow xlo.density = 1.225 xlo.temperature = 295.0 xlo.tke = 0.0 xhi.type = pressure_outflow ylo.type = mass_inflow ylo.density = 1.225 ylo.temperature = 295.0 ylo.tke = 0.0 yhi.type = pressure_outflow zlo.type = wall_model zhi.type = slip_wall zhi.temperature_type = fixed_gradient zhi.temperature = 0.003 #¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# # PHYSICS # #.......................................# ICNS.source_terms = BoussinesqBuoyancy CoriolisForcing BodyForce ABLMeanBoussinesq #ActuatorForcing ##--------- Additions by calc_inflow_stats.py ---------# ABL.wall_shear_stress_type = "local" ABL.inflow_outflow_mode = true ABL.wf_velocity = 3.8732822942681806 0.27257386624677665 ABL.wf_vmag = 3.895320565681261 ABL.wf_theta = 298.4848218926262 BodyForce.magnitude = 0.00030728511300339714 0.0007657185887876884 0.0 BoussinesqBuoyancy.read_temperature_profile = true BoussinesqBuoyancy.tprofile_filename = avg_theta.dat # --- ABL parameters --- incflo.velocity = 7.0 0.0 0.0 ABLForcing.abl_forcing_height = 90 CoriolisForcing.latitude = 36.57236455833334 CoriolisForcing.rotational_time_period = 86164.0900027328 CoriolisForcing.north_vector = 0.0 1.0 0.0 CoriolisForcing.east_vector = 1.0 0.0 0.0 ABL.normal_direction = 2 BoussinesqBuoyancy.reference_temperature = 295 ABL.temperature_heights = 0.0 90 600.0 700.0 1700.0 # Make sure the top height here goes above the domain height ABL.temperature_values = 295 299 301 302 303 ABL.perturb_temperature = true ABL.cutoff_height = 50 ABL.perturb_velocity = true ABL.perturb_ref_height = 50 ABL.Uperiods = 4.0 ABL.Vperiods = 4.0 ABL.deltaU = 1.0 ABL.deltaV = 1.0 ABL.kappa = .40 ABL.surface_roughness_z0 = 0.03 ABL.reference_temperature = 295 ABL.surface_temp_rate = -0.04 ABL.surface_temp_flux = 0.0 # Surface temperature flux [K-m/s] ABL.mo_beta_m = 16.0 # Monin-Obukhov Beta m parameter ABL.mo_gamma_m = 5.0 # Monin-Obukhov Gamma m parameter ABL.mo_gamma_h = 5.0 # Monin-Obukhov Gamma h parameter ABL.random_gauss_mean = 0.0 ABL.random_gauss_var = 1.0 #---- tagging defs ---- tagging.labels = T0_level_0_zone T0_level_1_zone 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 = 500.0 150.0 0.0 tagging.T0_level_0_zone.T0_level_0_zone.xaxis = 600.0 0.0 0.0 tagging.T0_level_0_zone.T0_level_0_zone.yaxis = 0.0 350.0 0.0 tagging.T0_level_0_zone.T0_level_0_zone.zaxis = 0.0 0.0 300.0 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 = 700.0 210.0 20.0 tagging.T0_level_1_zone.T0_level_1_zone.xaxis = 200.0 0.0 0.0 tagging.T0_level_1_zone.T0_level_1_zone.yaxis = 0.0 220.0 0.0 tagging.T0_level_1_zone.T0_level_1_zone.zaxis = 0.0 0.0 200.0 #---- postprocessing defs ---- incflo.post_processing = sampling averaging sampling.type = Sampling sampling.output_frequency = 2000 sampling.fields = velocity temperature tke averaging.type = TimeAveraging averaging.averaging_window = 60.0 averaging.averaging_start_time = 30000.0 #---- averaging defs ---- averaging.labels = means stress averaging.means.averaging_type = ReAveraging averaging.means.fields = velocity averaging.stress.averaging_type = ReynoldsStress averaging.stress.fields = velocity #---- sample defs ---- sampling.labels = xy xz yz sampling.xy.type = PlaneSampler sampling.xy.num_points = 256 256 sampling.xy.origin = 0.0 0.0 0.0 sampling.xy.axis1 = 2560.0 0.0 0.0 sampling.xy.axis2 = 0.0 640.0 0.0 sampling.xy.normal = 0.0 0.0 1.0 sampling.xy.offsets = 25 26.5 80 90 145 153 160 170 180 190 sampling.xz.type = PlaneSampler sampling.xz.num_points = 256 256 sampling.xz.origin = 0.0 0.0 0.0 sampling.xz.axis1 = 2550.0 0.0 0.0 sampling.xz.axis2 = 0.0 0.0 1270.0 sampling.xz.normal = 0.0 1.0 0.0 sampling.xz.offsets = 257 265 320 375 383 sampling.yz.type = PlaneSampler sampling.yz.num_points = 256 256 sampling.yz.origin = 0.0 0.0 0.0 sampling.yz.axis1 = 0.0 640.0 0.0 sampling.yz.axis2 = 0.0 0.0 1280.0 sampling.yz.normal = 1.0 0.0 0.0 sampling.yz.offsets = 800 910 927 1020 1054 1130 1181 1240 1308 1350 1435 1460 1562 1570 1680 1689 1790 1816 1900 1943 2010 2070 2120 2197 2230 2324 2451 2500 #¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨# # TURBINES # #.......................................# #Actuator.labels = T1 # #Actuator.T1.type = TurbineFastDisk #Actuator.T1.openfast_input_file = OpenFAST3p3_NREL2p8-127.bugfix/NREL-2p8-127.fst #Actuator.T1.base_position = 800 320 0.0 #Actuator.T1.rotor_diameter = 127 #Actuator.T1.hub_height = 89 #Actuator.T1.num_points_blade = 50 #Actuator.T1.num_points_tower = 12 #Actuator.T1.epsilon = 10.0 10.0 10.0 #Actuator.T1.epsilon_tower = 5.0 5.0 5.0 #Actuator.T1.openfast_start_time = 0.0 #Actuator.T1.openfast_stop_time = 10000 #Actuator.T1.nacelle_drag_coeff = 0.0 #Actuator.T1.nacelle_area = 0.0 #Actuator.T1.yaw = 270.0 #Actuator.T1.output_frequency = 1000 #Actuator.TurbineFastDisk.density = 1.225 #####