Changeset 4022 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

Jun 12, 2019 11:52:39 AM (5 years ago)

Author:

suehring

Message:

Synthetic turbulence generator: Revise bias correction of imposed perturbations (correction via volume flow can create instabilities in case the mean volume flow is close to zero); Introduce lower limits in calculation of coefficient matrix, else the calculation may become numerically unstable; Impose perturbations every timestep, even though no new set of perturbations is generated in case dt_stg_call /= dt_3d; Implement a gradual decrease of Reynolds stress and length scales above ABL height (within 1 length scale above ABL depth to 1/10) rather than an abrupt decline; Bugfix in non-nested case: use ABL height for parametrized turbulence; Offline nesting: Rename subroutine for ABL height calculation and make it public; Change top boundary condition for pressure back to Neumann

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (2 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Call synthetic turbulence generator at last RK3 substep right after boundary
+! conditions are updated in offline nesting in order to assure that 
+! perturbations are always imposed 
+! 
+! 4017 2019-06-06 12:16:46Z schwenkel
 ! Mass (volume) flux correction included to ensure global mass conservation for child domains.
 ! 
@@ -1240 +1245 @@
              CALL nesting_offl_bc
 !
-!--       Impose a turbulent inflow using synthetic generated turbulence, 
-!--       only once per time step.
-          IF ( use_syn_turb_gen  .AND.  time_stg_call >= dt_stg_call  .AND.                        &
-             intermediate_timestep_count == intermediate_timestep_count_max )  THEN
+!--       Impose a turbulent inflow using synthetic generated turbulence.
+          IF ( use_syn_turb_gen  .AND.                                                             &
+               intermediate_timestep_count == intermediate_timestep_count_max )  THEN
              CALL cpu_log( log_point(57), 'synthetic_turbulence_gen', 'start' )
              CALL stg_main