Changeset 2100 for palm

Timestamp:

Jan 5, 2017 4:40:16 PM (8 years ago)

Author:

suehring

Message:

Prevent extremely large SGS-velocities by limiting TKE in stochastic terms.

File:

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

palm/trunk/SOURCE/lpm_advec.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Prevent extremely large SGS-velocities in regions where TKE is zero, e.g. 
+! at the begin of simulations and/or in non-turbulent regions. 
 ! 
 ! Former revisions:
@@ -1578 +1579 @@
     REAL(wp) ::  v_sgs   !< subgrid-scale velocity component 
 
+!-- At first, limit TKE to a small non-zero number, in order to prevent 
+!-- the occurrence of extremely large SGS-velocities in case TKE is zero, 
+!-- (could occur at the simulation begin).
+    e_n = MAX( e_n, 1E-20_wp )
 !
 !-- Please note, terms 1 and 2 (drift and memory term, respectively) are