# --- Simulation time control parameters --- time.stop_time = 18800 # Max (simulated) time to evolve [s] time.max_step = -1 time.fixed_dt = 0.1 # Fixed timestep size (in seconds). If negative, then time.cfl is used time.cfl = 0.95 time.checkpoint_interval = 2500 time.checkpoint_start = 36000 incflo.physics = ABL Actuator # List of physics models to include in simulation. incflo.verbose = 0 io.check_file = /lustre/scratch/asim/latest_run/Unstable/chk io.restart_file = /lustre/scratch/asim/latest_run/Unstable/chk36000 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 = 1.642e-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. # Lo corner coordinates geometry.prob_hi = 2560. 640 1280. # Hi corner coordinates amr.n_cell = 256 64 128 # Grid cells at coarsest AMRlevel amr.max_level = 2 # Max AMR level in hierarchy geometry.is_periodic = 0 0 0 xlo.type = mass_inflow xlo.density = 1.225 xlo.temperature = 0.0 xlo.tke = 0.0 xhi.type = pressure_outflow ylo.type = mass_inflow ylo.density = 1.225 ylo.temperature = 0.0 ylo.tke = 0.0 yhi.type = pressure_outflow zlo.type = wall_model zlo.temperature_type = wall_model zlo.tke_type = zero_gradient zhi.type = slip_wall zhi.temperature_type = fixed_gradient zhi.temperature = 0.003 # --- ABL parameters --- ICNS.source_terms = BoussinesqBuoyancy CoriolisForcing ActuatorForcing BodyForce ABLMeanBoussinesq ABL.bndry_io_mode = 1 ABL.bndry_file = /lustre/scratch/asim/latest_run/Unstable/bndry_file.nc ABL.bndry_planes = xlo ylo ABL.bndry_var_names = velocity temperature tke ABL.bndry_output_format = native incflo.velocity = 7.0 0.0 0.0 ABLForcing.abl_forcing_height = 90.0 ABL.kappa = 0.4 ABL.normal_direction = 2 ABL.surface_roughness_z0 = 0.5 ABL.reference_temperature = 298 ABL.surface_temp_rate = 0.0 ABL.surface_temp_flux = 0.08 # 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 ABL.wall_shear_stress_type = local 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 BoussinesqBuoyancy.reference_temperature = 298 BodyForce.magnitude = 0.0002772541892390438 0.0010086798505818944 0.0 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 = 298 295.5 305.803 304.5 303.5 time.plot_interval = 4500 io.plot_file = /lustre/scratch/asim/latest_run/Unstable/plt io.KE_int = -1 #---- 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 = 700 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 = 750 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 #---- actuator defs ---- Actuator.labels = T0 Actuator.T0.type = TurbineFastDisk Actuator.T0.openfast_input_file = Turbine/NREL-2p8-127.fst Actuator.T0.base_position = 800 320.0 0.0 Actuator.T0.rotor_diameter = 127.0 Actuator.T0.hub_height = 90.0 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 = 10000.0 Actuator.T0.fllc = false Actuator.T0.nacelle_drag_coeff = 1.0 Actuator.T0.nacelle_area = 8.0 Actuator.T0.yaw = 270.0 Actuator.T0.output_frequency = 1 Actuator.T0.num_blades = 3 Actuator.T0.use_tip_correction = true Actuator.T0.use_root_correction = true Actuator.TurbineFastDisk.density = 1.225 #---- postprocessing defs ---- io.post_processing_directory = /lustre/scratch/asim/latest_run/Unstable/ incflo.post_processing = sampling # --- Sampling parameters --- sampling.output_frequency = 25 sampling.fields = velocity temperature tke #---- 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 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 #---- extra params ---- #== END AMR-WIND INPUT ==