Changeset 4122 for palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

Timestamp:

Jul 26, 2019 1:11:56 PM (5 years ago)

Author:

schwenkel

Message:

Implement reset method as bottom boundary condition

File:

• palm/trunk/SOURCE/lagrangian_particle_model_mod.f90 (modified) (8 diffs)

palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implement reset method as bottom boundary condition
+! 
+! 4121 2019-07-26 10:01:22Z schwenkel
 ! Implementation of an simple method for interpolating the velocities to
 ! particle position
@@ -173 +176 @@
 ! Description:
 ! ------------
+!> The embedded LPM allows for studying transport and dispersion processes within
+!> turbulent flows. This model including passive particles that do not show any
+!> feedback on the turbulent flow. Further also particles with inertia and
+!> cloud droplets ca be simulated explicitly.
 !>
+!> @todo test lcm
+!>       implement simple interpolation method for subgrid scale velocites
+!> @note <Enter notes on the module>
+!> @bug  <Enter bug on the module>
 !------------------------------------------------------------------------------!
  MODULE lagrangian_particle_model_mod
@@ -253 +264 @@
     USE NETCDF
 #endif
-
-
-     USE arrays_3d,                                                             &
-        ONLY:
-
-    USE indices,                                                               &
-        ONLY:  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
-
-    USE kinds
-
-    USE pegrid
 
     IMPLICIT NONE
@@ -1074 +1074 @@
        CASE ( 'reflect' )
           ibc_par_b = 2
+
+       CASE ( 'reset' )
+          ibc_par_b = 3
 
        CASE DEFAULT
@@ -4371 +4374 @@
     INTEGER(iwp) ::  i2             !< grid index (x) of current particle position
     INTEGER(iwp) ::  i3             !< grid index (x) of intermediate particle position
+    INTEGER(iwp) ::  index_reset    !< index reset height
     INTEGER(iwp) ::  jr             !< dummy for sorting algorithmus
     INTEGER(iwp) ::  j1             !< grid index (y) of old particle position
@@ -4379 +4383 @@
     INTEGER(iwp) ::  k3             !< grid index (z) of intermediate particle position
     INTEGER(iwp) ::  n              !< particle number
+    INTEGER(iwp) ::  particles_top  !< maximum reset height
     INTEGER(iwp) ::  t_index        !< running index for intermediate particle timesteps in reflection algorithmus
     INTEGER(iwp) ::  t_index_number !< number of intermediate particle timesteps in reflection algorithmus
@@ -4417 +4422 @@
     REAL(wp) ::  prt_y          !< current particle position (y)
     REAL(wp) ::  prt_z          !< current particle position (z)
+    REAL(wp) ::  ran_val        !< location of wall in z
+    REAL(wp) ::  reset_top      !< location of wall in z
     REAL(wp) ::  t_old          !< previous reflection time
     REAL(wp) ::  tmp_t          !< dummy for sorting algorithmus
@@ -4487 +4494 @@
                      particles(n)%rvar3 < 0.0_wp )  THEN
                    particles(n)%rvar3 = -particles(n)%rvar3
+                ENDIF
+             ELSEIF ( ibc_par_b == 3 )  THEN
+!
+!--             Find reset height. @note this works only in non-strechted cases
+                particles_top = INT( pst(1) / dz(1) )
+                index_reset = MINLOC( prt_count(nzb+1:particles_top,j,i), DIM = 1 )
+                reset_top = zu(index_reset)
+                iran_part = iran_part + myid
+                ran_val = random_function( iran_part )
+                particles(n)%z       = reset_top *  ( 1.0  + ( ran_val / 10.0_wp) )
+                particles(n)%speed_z = 0.0_wp
+                IF ( curvature_solution_effects ) THEN
+                   particles(n)%radius = particles(n)%aux1
+                ELSE
+                   particles(n)%radius = 1.0E-8
                 ENDIF
              ENDIF