source: palm/trunk/TESTS/cases/rans_tkee/MONITORING/rans_tkee_rc @ 4102

Last change on this file since 4102 was 4102, checked in by suehring, 5 years ago

Bugfix, set Neumann boundary conditions for the subgrid TKE at vertical walls instead of implicit Dirichlet conditions that always act as a sink term for the subgrid TKE. Therefore, add new data structure for vertical surfaces and revise the setting of the boundary grid point index space. Moreover, slightly revise setting of boundary conditions at upward- and downward facing surfaces. Finally, setting of boundary conditions for subgrid TKE and dissipation (in RANS mode) is now modularized. Update test case results.

File size: 11.9 KB
RevLine 
[3541]1
2
3
41D run control output:
5------------------------------
6
7ITER.   HH:MM:SS    DT      UMAX   VMAX    U*   ALPHA   ENERG.
8-------------------------------------------------------------
9      0 00:00:00    0.01    3.00   0.10  0.100   45.0     4.45
10   3325 12:00:04   12.30    3.36   0.87  0.127   26.2     4.48
11
12 ******************************    --------------------------------------------
[4102]13 * PALM 6.0  Rev: 3915M       *    atmosphere - run with 1D - prerun
[3541]14 ******************************    --------------------------------------------
15
[4101]16 Date:               2019-07-17    Run:       rans_tkee__gfortran_default__4   
[4102]17 Time:                 17:23:23    Run-No.:   00
[4101]18 Run on host:        gfortran_d
[3541]19 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
20 ------------------------------------------------------------------------------
21
22 Numerical Schemes:
23 -----------------
24
25 --> Use the TKE-e turbulence closure (RANS mode).
26 --> Use the boussinesq approximation for the model equations.
27 --> Solve perturbation pressure via multigrid method (w-cycle)
28     number of grid levels:                    4
29     Gauss-Seidel red/black iterations:        2
30     fixed number of multigrid cycles:           2
31     perturbation pressure is calculated at every Runge-Kutta step
32 --> Momentum advection via Wicker-Skamarock-Scheme 5th order
33 --> Scalar advection via Wicker-Skamarock-Scheme 5th order
34 --> Loop optimization method: cache
35 --> Time differencing scheme: runge-kutta-3       
36 ------------------------------------------------------------------------------
37
38
39 Run time and time step information:
40 ----------------------------------
41
42 Timestep:             variable     maximum value: 20.000 s    CFL-factor: 0.90
43 Start time:              0.000 s
44 End time:               50.000 s
45
46
47 Computational grid and domain size:
48 ----------------------------------
49
50 Grid length:      dx =      10.000 m    dy =      10.000 m
51                dz(1) =      10.000 m
52                dz(2) =     999.000 m
53
54 Domain size:       x =    320.000 m     y =    320.000 m  z(u) =    965.000 m
55
56
57
58 Number of gridpoints (x,y,z):  (0:  31, 0:  31, 0:  97)
59 Subdomain size (x,y,z):        (    16,     16,     98)
60
61
62    Characteristic levels of the geo. wind component ug:
63
64       Height:         0.0     0.0  m
65       ug:            3.00    3.00  m/s
66       Gradient:    ------    0.00  1/100s
67       Gridpoint:        0       0
68
69    Characteristic levels of the geo. wind component vg:
70
71       Height:         0.0     0.0  m
72       vg:            0.00    0.00  m/s
73       Gradient:    ------    0.00  1/100s
74       Gridpoint:        0       0
75
76
77 Topography information:
78 ----------------------
79
80 Topography: read_from_file                         
81 Topography grid definition convention:
82 cell center (scalar grid points)
83
84
85
86 Boundary conditions:
87 -------------------
88
89                     p                    uv                                  pt
90
91 B. bound.: p(0)     = p(1)   | uv(0)     = -uv(1)                | pt(0)     = pt(1)           
92 T. bound.: p(nzt+1) = 0      | uv(nzt+1) = ug(nzt+1), vg(nzt+1)  | pt(nzt+1) = pt(nzt) + dpt/dz_
93
94                     e
95
96 B. bound.: e(0)     = e(1)                                                                     
97 T. bound.: e(nzt+1) = e(nzt) = e(nzt-1)                                                         
98
99 Bottom surface fluxes are used in diffusion terms at k=1
100       Predefined constant heatflux:    0.001000 K m/s
101
102
103    Constant flux layer between bottom surface and first computational u,v-level:
104
105       z_mo =   5.00 m   z0 = 0.1000 m   z0h = 0.10000 m   kappa = 0.40
106       Rif value range:     -20.00 <= rif <= 20.00
107
108
109 Lateral boundaries:
110       left/right:  dirichlet/radiation
111       north/south: cyclic             
112
113       use_cmax: T
114       pt damping layer width =     0.00 m, pt damping factor = 0.0000
115
116
117 Initial profiles:
118 ----------------
119
120    Characteristic levels of the initial temperature profile:
121
122       Height:           0.0     0.0  m
123       Temperature:   300.00  300.00  K
124       Gradient:      ------    0.00  K/100m
125       Gridpoint:          0       0
126
127
128 List output:
129 -----------
130
131    1D-Profiles:
132       Output every             9999999.90 s
133
134
135 Data output:
136 -----------
137
138
139    1D-Profiles:
140       Output format: netCDF 64bit offset                     
141
142
143       Profile: u, w"u", w*u*, wu, u*2, v, w"v", w*v*, wv, v*2, w, w*2,     
144
145              : e, e*, km, l,                                               
146       Output every                10.00 s
147       Time averaged over           0.00 s
148       Averaging input every        0.00 s
149
150    Time series:
151       Output format: netCDF 64bit offset                     
152
153       Output every                 1.00 s
154
155 ------------------------------------------------------------------------------
156
157
158 Physical quantities:
159 -------------------
160
161    Geograph. latitude  :   latitude  = 52.5 degr
162    Geograph. longitude :   longitude = 13.3 degr
163    Angular velocity    :   omega  = 0.729E-04 rad/s
164    Coriolis parameter  :   f      =  0.000116 1/s
165                            f*     =  0.000089 1/s
166
167    Day of the year at model start :   day_init      =     172
168    UTC time at model start        :   time_utc_init =     0.0 s
169
170    Gravity             :   g      =  9.8 m/s**2
171
172    Reference state used in buoyancy terms: initial_profile
173
174
175 Cloud physics quantities / methods:
176 ----------------------------------
177
178
179
180 LES / Turbulence quantities:
181 ---------------------------
182
183    initial TKE is prescribed as  0.10E+00 (m/s)**2
184    TKE is not allowed to fall below  0.10E-02 (m/s)**2
185    Mixing length is limited to 1.80 * z
186    RANS-mode constants: c_0 =   0.55000
187                         c_1 =   1.44000
188                         c_2 =   1.92000
189                         c_3 =   1.44000
190                         c_4 =   0.00000
191                         sigma_e    =   1.00000
192                         sigma_diss =   1.30000
193
194
195 Actions during the simulation:
196 -----------------------------
197
198
199
200 1D-Model parameters:
201 -------------------
202
203    Simulation time:                    43200.0 s
204    Run-controll output every:          43200.0 s
205    Vertical profile output every:      86400.0 s
206    Mixing length calculation:         blackadar           
207    Dissipation calculation:           prognostic         
208
209
210
[4101]211 *** dynamic module disabled
[3541]212
213 ------------------------------------------------------------------------------
214
215
216
217Run-control output:
218------------------
219
220RUN  ITER. HH:MM:SS.SS    DT(E)     UMAX     VMAX     WMAX     U*    W*      THETA*     Z_I     ENERG.   DISTENERG    DIVOLD     DIVNEW     UMAX(KJI)    VMAX(KJI)    WMAX(KJI)   ADVECX   ADVECY   MGCYC
221---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
[4102]222  0      0 00:00:00.00   0.0200A   3.3592   0.8728   0.0000  0.127  0.00   0.000E+00      0.  0.442E+01  0.442E-02  0.000E+00  0.000E+00   20  29  29    6  29  29    0  29  29    0.000    0.000       0
223  0      1 00:00:00.02   0.0400A   3.3802   1.7430  -1.2281  0.153  0.00  -0.390E-02     10.  0.444E+01  0.460E-02  0.155E-03  0.701E-04   19   7  19    2  14  13    2  26  29    0.000    0.000       2
224  0      2 00:00:00.06   0.0800A   3.3832   1.7620  -1.3655  0.151  0.20  -0.421E-02     10.  0.444E+01  0.493E-02  0.874E-04  0.359E-04   19   6  21    2  14  13    2  26  29    0.000    0.000       2
225  0      3 00:00:00.14   0.1600A   3.3844   1.7515   1.3704  0.150  0.00  -0.428E-02      0.  0.444E+01  0.496E-02  0.977E-04  0.323E-04   19   5  22    2  14  13    2  13   1    0.000    0.000       2
226  0      4 00:00:00.30   0.3200A   3.3845   1.7361   1.3743  0.150  0.00  -0.429E-02      0.  0.444E+01  0.456E-02  0.521E-04  0.687E-05   19   5  24    2  14  13    2  13   1    0.000    0.000       2
227  0      5 00:00:00.62   0.6400A   3.3838   1.7102   1.3741  0.149  0.00  -0.432E-02      0.  0.444E+01  0.437E-02  0.385E-04  0.278E-04   19   4  26    2  14  13    2  13   1    0.000    0.000       2
228  0      6 00:00:01.26   1.2800A   3.3833   1.6639   1.3660  0.148  0.00  -0.437E-02      0.  0.444E+01  0.430E-02  0.704E-04  0.332E-04   19   1   0    2  14  13    2  13   1    0.000    0.000       2
229  0      7 00:00:02.54   2.5600A   3.3827   1.5852   1.3342  0.147  0.39  -0.449E-02     30.  0.444E+01  0.420E-02  0.618E-04  0.165E-04   19   1   0    2  14  13    2  13   1    0.000    0.000       2
230  0      8 00:00:05.09   2.6600A   3.3826   1.4669   1.2489  0.145  0.00  -0.484E-02     20.  0.444E+01  0.421E-02  0.144E-03  0.499E-04   20   3  31    2  14  13    2  13   1    0.000    0.000       2
231  0      9 00:00:07.75   2.6200D   3.3825   1.4248   1.2091  0.142  0.00  -0.540E-02      0.  0.444E+01  0.399E-02  0.137E-03  0.424E-04   20   3   0    2  17   0    2  13   1    0.000    0.000       2
232  0     10 00:00:10.37   2.6600A   3.3822   1.4058   1.2176  0.140  0.00  -0.536E-02      0.  0.444E+01  0.380E-02  0.971E-04  0.324E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
233  0     11 00:00:13.03   2.6600A   3.3818   1.3920   1.2409  0.138  0.00  -0.554E-02      0.  0.444E+01  0.372E-02  0.912E-04  0.361E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
234  0     12 00:00:15.69   2.6600A   3.3818   1.3841   1.2506  0.137  0.00  -0.573E-02      0.  0.444E+01  0.366E-02  0.867E-04  0.300E-04   19   6  24    2  17   0    1  12   2    0.000    0.000       2
235  0     13 00:00:18.35   2.6600A   3.3820   1.3790   1.2520  0.135  0.00  -0.581E-02      0.  0.444E+01  0.362E-02  0.864E-04  0.282E-04   19   6  23    2  17   0    1  12   2    0.000    0.000       2
236  0     14 00:00:21.01   2.6600A   3.3822   1.3757   1.2484  0.134  0.00  -0.617E-02      0.  0.444E+01  0.358E-02  0.862E-04  0.277E-04   19   6  23    2  17   0    1  12   2    0.000    0.000       2
237  0     15 00:00:23.67   2.6600A   3.3824   1.3733   1.2423  0.133  0.00  -0.644E-02      0.  0.444E+01  0.353E-02  0.800E-04  0.249E-04   19   6  24    2  17   0    1  12   2    0.000    0.000       2
238  0     16 00:00:26.33   2.6600A   3.3825   1.3716   1.2356  0.132  0.00  -0.633E-02      0.  0.444E+01  0.348E-02  0.755E-04  0.250E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
239  0     17 00:00:29.00   2.6600A   3.3827   1.3703   1.2296  0.131  0.00  -0.648E-02      0.  0.444E+01  0.345E-02  0.734E-04  0.242E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
240  0     18 00:00:31.66   2.6600A   3.3829   1.3692   1.2245  0.130  0.00  -0.650E-02      0.  0.444E+01  0.342E-02  0.719E-04  0.232E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
241  0     19 00:00:34.32   2.6600A   3.3830   1.3683   1.2206  0.129  0.00  -0.661E-02      0.  0.444E+01  0.339E-02  0.709E-04  0.228E-04   19   7  24    2  17   0    1  12   2    0.000    0.000       2
242  0     20 00:00:36.98   2.6600A   3.3832   1.3675   1.2176  0.128  0.00  -0.648E-02      0.  0.444E+01  0.336E-02  0.692E-04  0.224E-04   19   7  24    2  17   0    1  12   2    0.000    0.000       2
243  0     21 00:00:39.64   2.6600A   3.3834   1.3668   1.2153  0.128  0.00  -0.649E-02      0.  0.444E+01  0.334E-02  0.674E-04  0.221E-04   19   8  24    2  17   0    1  12   2    0.000    0.000       2
244  0     22 00:00:42.29   2.6600A   3.3836   1.3662   1.2136  0.127  0.00  -0.658E-02      0.  0.444E+01  0.331E-02  0.661E-04  0.218E-04   19   8  24    2  17   0    1  12   2    0.000    0.000       2
245  0     23 00:00:44.95   2.6600A   3.3837   1.3657   1.2122  0.127  0.00  -0.664E-02      0.  0.444E+01  0.329E-02  0.651E-04  0.215E-04   19   8  24    2  17   0    1  12   2    0.000    0.000       2
246  0     24 00:00:47.61   2.6600A   3.3839   1.3652   1.2111  0.126  0.00  -0.661E-02      0.  0.444E+01  0.327E-02  0.642E-04  0.212E-04   19   8  24    2  17   0    1  12   2    0.000    0.000       2
247  0     25 00:00:50.27   2.6600A   3.3841   1.3647   1.2101  0.126  0.00  -0.654E-02      0.  0.444E+01  0.325E-02  0.634E-04  0.210E-04   19   8  25    2  17   0    1  12   2    0.000    0.000       2
Note: See TracBrowser for help on using the repository browser.