1 | <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> |
---|
2 | <html><head> |
---|
3 | |
---|
4 | <meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"> |
---|
5 | |
---|
6 | |
---|
7 | |
---|
8 | <title>PALM chapter 3.9</title><meta name="GENERATOR" content="StarOffice 7 (Win32)"> |
---|
9 | |
---|
10 | <meta name="AUTHOR" content="Siegfried Raasch"> |
---|
11 | |
---|
12 | <meta name="CREATED" content="20040809;13460943"> |
---|
13 | |
---|
14 | <meta name="CHANGED" content="20041112;15085727"> |
---|
15 | |
---|
16 | <meta name="KEYWORDS" content="parallel LES model"> |
---|
17 | |
---|
18 | <style> |
---|
19 | <!-- |
---|
20 | @page { size: 21cm 29.7cm } |
---|
21 | --> |
---|
22 | </style></head> |
---|
23 | <body style="direction: ltr;" lang="en-US"> |
---|
24 | |
---|
25 | <h3 style="line-height: 100%;">3.9 Setup for turbulent inflow</h3> |
---|
26 | |
---|
27 | <p style="line-height: 100%;">Starting from version 3.5 |
---|
28 | PALM includes a so-called turbulence recycling method which allows a |
---|
29 | turbulent inflow with non-cyclic horizontal boundary conditions. The |
---|
30 | method follows the one described by Lund et al. (1998, J. Comp. Phys., <span style="font-weight: bold;">140</span>, 233-258), modified by Kataoka and Mizuno (2002, Wind and Structures, <span style="font-weight: bold;">5</span>, 379-392). The method is switched on by setting the initial parameter <a href="chapter_4.1.html#turbulent_inflow">turbulent_inflow</a> = <span style="font-style: italic;">.TRUE.</span>.</p><p style="line-height: 100%;">The turbulent signal A'(y,z) to be imposed at the left inflow boundary is taken from the same simulation at a fixed distance x<sub>r</sub> from the inflow (given by parameter <a href="chapter_4.1.html#recycling_width">recycling_width</a>): A'(y,z) = A(x<sub>r</sub>,y,z) - <span style="font-weight: bold;">A(z)</span>, where <span style="font-weight: bold;">A(z)</span> |
---|
31 | is the horizontal average between the inflow boundary and the recycling |
---|
32 | plane. The turbulent quantity A'(y,z) is then added to a mean inflow |
---|
33 | profile a(z) which is the horizontal (and temporal) average as |
---|
34 | calculated from a precursor run. This is done after each timestep.</p><p style="line-height: 100%; font-weight: bold;">So far (09/16/08), a turbulent signal is imposed to all prognostic variables except humidity (passive scalar) and salinity!</p><p style="line-height: 100%;">In |
---|
35 | order to avoid a continous growth of the turbulent layer with time, the |
---|
36 | height up to which turbulence is recycled is limited by parameter <a href="chapter_4.1.html#inflow_damping_height">inflow_damping_height</a>. By default, this height is calculated and taken from the precursor run. </p>The following tries to give a short recipe about how to create a proper setup for using the turbulence recycling method:<br><br><ul><li>A |
---|
37 | precursor run with cyclic horizontal boundary conditions has to be |
---|
38 | carried out in order to create a (quasi-) stationary state. Since the |
---|
39 | horizontal average from this precursor run is used as the mean inflow |
---|
40 | profile for the main run, <span style="font-weight: bold;">the |
---|
41 | wall-normal velocity component must point |
---|
42 | into the domain at every grid point and its magnitude should be large |
---|
43 | enough in order to guarantee an inflow even if a turbulence signal is |
---|
44 | added</span>.<br></li><li>The |
---|
45 | main run requires from the precursor run the mean profiles to |
---|
46 | be used at the inflow. For this, the horizontally and temporally |
---|
47 | averaged mean profiles as provided with the standard PALM output are |
---|
48 | used. The user has to set parameters <a href="chapter_4.2.html#dt_data_output_pr">dt_data_output_pr</a>, <a href="chapter_4.2.html#averaging_interval">averaging_interval</a>, |
---|
49 | etc. for the precursor run appropriately, so that an output is done <span style="font-weight: bold;">at |
---|
50 | the end</span> of the precursor run. The profile information is then contained |
---|
51 | in the restart (binary) file created at the end of the precursor run |
---|
52 | and can be used by the main run. <span style="font-weight: bold;">It is very important that the mean profiles at the end of the precursor run are in a stationary or quasi-stationary state</span>, because otherwise it may not be justified to use them as constant profiles at the inflow. <span style="font-weight: bold;">Also, turbulence at the end of the precursor run should be fully developed. </span>Otherwise, the main run would need an additional spinup-time at the beginning to get the turbulence to its final stage.<br></li><li>In order to read the binary data from the precursor run, parameter <a href="chapter_4.1.html#initializing_actions">initializing_actions</a> = <span style="font-style: italic;">'read_data_for_recycling'</span> has to be set in the main run. .... |
---|
53 | set bc_lr = 'dirichlet/radiation' ... |
---|
54 | conserve_volume_flow = .TRUE. ? ....</li><li>The main run allows to use a larger domain along x and y than in the precursor run ...</li><li>Recycling width should be equal to the domain size of the precursor run (or larger?)</li></ul><pre style="line-height: 100%;"><font style="font-size: 10pt;" size="2">#</font><br><font style="font-size: 10pt;" size="2">RUN_CONTROL<span style="font-weight: bold;">_O</span></font><font style="font-size: 10pt;" size="2"> out:loc:tr d3<span style="font-weight: bold;">o</span># ~/palm/current_version/JOBS/$fname/MONITORING <span style="font-weight: bold;">_</span></font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">o</span>_</font><font style="font-size: 10pt;" size="2">rc</font><br><font style="font-size: 10pt;" size="2">HEADER</font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">_O</span> out:loc:tr d3<span style="font-weight: bold;">o</span># ~/palm/current_version/JOBS/$fname/MONITORING <span style="font-weight: bold;">_</span></font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">o</span>_</font><font style="font-size: 10pt;" size="2">header</font><br><font style="font-size: 10pt;" size="2">PLOT1D_PAR<span style="font-weight: bold;">_O</span></font><font style="font-size: 10pt;" size="2"> out:loc:tr pr<span style="font-weight: bold;">o</span># ~/palm/current_version/JOBS/$fname/OUTPUT <span style="font-weight: bold;">_</span></font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">o</span>_</font><font style="font-size: 10pt;" size="2">pr_par</font><br><font style="font-size: 10pt;" size="2">PLOT1D_DATA</font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">_O</span> out:loc:tr pr<span style="font-weight: bold;">o</span># ~/palm/current_version/JOBS/$fname/OUTPUT <span style="font-weight: bold;">_</span></font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">o</span>_</font><font style="font-size: 10pt;" size="2">pr_in</font><br></pre><br><hr> |
---|
55 | <p style="line-height: 100%;"><br> |
---|
56 | |
---|
57 | <font color="#000080"><font color="#000080"><a href="chapter_3.8.html"><font color="#000080"><img style="border: 2px solid ; width: 32px; height: 32px;" alt="" src="left.gif" name="Grafik1"></font></a><a href="index.html"><font color="#000080"><img src="up.gif" name="Grafik2" align="bottom" border="2" height="32" width="32"></font></a><a href="chapter_4.0.html"><font color="#000080"><img style="border: 2px solid ; width: 32px; height: 32px;" alt="" src="right.gif" name="Grafik3"></font></a></font></font></p> |
---|
58 | |
---|
59 | <p style="line-height: 100%;"><i>Last |
---|
60 | change: </i> $Id: chapter_3.9.html 220 2008-12-18 07:00:36Z maronga $</p> |
---|
61 | |
---|
62 | </body></html> |
---|