Ignore:
Timestamp:
Jun 12, 2019 11:52:39 AM (2 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/time_integration.f90

    r4017 r4022  
    2525! -----------------
    2626! $Id$
     27! Call synthetic turbulence generator at last RK3 substep right after boundary
     28! conditions are updated in offline nesting in order to assure that
     29! perturbations are always imposed
     30!
     31! 4017 2019-06-06 12:16:46Z schwenkel
    2732! Mass (volume) flux correction included to ensure global mass conservation for child domains.
    2833!
     
    12401245             CALL nesting_offl_bc
    12411246!
    1242 !--       Impose a turbulent inflow using synthetic generated turbulence,
    1243 !--       only once per time step.
    1244           IF ( use_syn_turb_gen  .AND.  time_stg_call >= dt_stg_call  .AND.                        &
    1245              intermediate_timestep_count == intermediate_timestep_count_max )  THEN
     1247!--       Impose a turbulent inflow using synthetic generated turbulence.
     1248          IF ( use_syn_turb_gen  .AND.                                                             &
     1249               intermediate_timestep_count == intermediate_timestep_count_max )  THEN
    12461250             CALL cpu_log( log_point(57), 'synthetic_turbulence_gen', 'start' )
    12471251             CALL stg_main
Note: See TracChangeset for help on using the changeset viewer.