Changeset 62 for palm/trunk/DOC/app/chapter_2.0.html
- Timestamp:
- Mar 13, 2007 2:52:40 AM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/app/chapter_2.0.html
r54 r62 1 1 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> 2 <html> 3 <head> 4 <meta content="text/html; charset=windows-1252" 5 http-equiv="CONTENT-TYPE"> 6 <title>PALM chapter 2.0</title> 7 <meta content="StarOffice 7 (Win32)" name="GENERATOR"> 8 <meta content="Marcus Oliver Letzel" name="AUTHOR"> 9 <meta content="20040719;14534028" name="CREATED"> 10 <meta content="20041117;10385730" name="CHANGED"> 11 <meta content="parallel LES model" name="KEYWORDS"> 12 <style> 13 <!-- 14 @page { size: 21cm 29.7cm } 15 --> 16 </style> 17 </head> 18 <body lang="en-US" dir="ltr"> 19 <h2 style="line-height: 100%;"><font size="4">2.0 Basic techniques of 2 <html><head> 3 <meta content="text/html; charset=windows-1252" http-equiv="CONTENT-TYPE"><title>PALM chapter 2.0</title> <meta content="StarOffice 7 (Win32)" name="GENERATOR"> <meta content="Marcus Oliver Letzel" name="AUTHOR"> <meta content="20040719;14534028" name="CREATED"> <meta content="20041117;10385730" name="CHANGED"> <meta content="parallel LES model" name="KEYWORDS"> <style> 4 <!-- 5 @page { size: 21cm 29.7cm } 6 --> 7 </style></head> 8 9 <body style="direction: ltr;" lang="en-US"><h2 style="line-height: 100%;"><font size="4">2.0 10 Basic techniques of 20 11 the LES model and its parallelization </font> 21 </h2> 22 <p style="line-height: 100%;">LES modelsgenerally permit the12 </h2><p style="line-height: 100%;">LES models 13 generally permit the 23 14 simulation of turbulent flows, whereby those eddies, that carry the 24 15 main energy are resolved by the numerical grid. Only the … … 28 19 simulated directly (they are explicitly resolved) and their effects are 29 20 represented by the advection terms. </p> 30 <p style="font-style: normal; line-height: 100%;">PALM is based on the 21 <p style="font-style: normal; line-height: 100%;">PALM is 22 based on the 31 23 non-hydrostatic incompressible Boussinesq equations. It contains a 32 24 water cycle with cloud formation and takes into account infrared 33 25 radiative cooling in cloudy conditions. The model has six prognostic 34 quantities in total u,v,w, liquid water potential temperature 35 <font face="Thorndale, serif">Θ</font><sub>l </sub>(BETTS, 26 quantities in total – u,v,w, liquid water potential 27 temperature 28 <font face="Thorndale, serif">Θ</font><sub>l 29 </sub>(BETTS, 36 30 1973), total water content q and subgrid-scale turbulent kinetic energy 37 31 e. The … … 45 39 1993). The water cycle is closed by using a simplified version of 46 40 KESSLERs scheme (KESSLER, 1965; 1969) to parameterize precipitation 47 processes (M ÜLLER and CHLOND, 1996). Incompressibility is41 processes (MÜLLER and CHLOND, 1996). Incompressibility is 48 42 applied by means of a Poisson equation for pressure, which is solved 49 43 with a direct method (SCHUMANN and SWEET, 1988). The Poisson equation … … 58 52 horizontal directions. At the lower surface, either temperature/ 59 53 humidity or their respective fluxes can be prescribed. </p> 60 <p style="font-style: normal; line-height: 100%;">The advection terms 54 <p style="font-style: normal; line-height: 100%;">The 55 advection terms 61 56 are treated by the scheme proposed by PIACSEK and WILLIAMS (1970), 62 57 which conserves the integral of linear and quadratic quantities up to … … 66 61 with the third-order Runge-Kutta scheme. A second-order Runge-Kutta 67 62 scheme, a leapfrog scheme and an Euler scheme are also implemented.</p> 68 <p style="line-height: 100%;">By default, the time step is computed 63 <p style="line-height: 100%;">By default, the time step is 64 computed 69 65 with respect to the different criteria (CFL, diffusion) and adapted 70 66 automatically. In case of a non-zero geostrophic 71 67 wind the coordinate system can be moved along with the mean wind in 72 68 order to maximize the time step (Galilei-Transformation). </p> 73 <p style="font-style: normal; line-height: 100%;">In principle a model 69 <p style="font-style: normal; line-height: 100%;">In 70 principle a model 74 71 run is carried out in the following way: After reading the control 75 72 parameters given by the user, all prognostic variables are … … 90 87 corrected with the help of the pressure solver. Following this, all 91 88 diagnostic turbulence quantities including possible 92 Prandtl-layer quantities are computed. At the end of a time89 Prandtl-layer–quantities are computed. At the end of a time 93 90 step the data output requested by the user is made 94 91 (e.g. statistic of analyses for control purposes or profiles and/or 95 92 graphics data). If the given end-time was reached, binary data maybe 96 93 be saved for restart. </p> 97 <p style="font-style: normal; line-height: 100%;">The model is based 94 <p style="font-style: normal; line-height: 100%;">The 95 model is based 98 96 on the originally non-parallel LES model which has been operated at the 99 97 institute since 1989 … … 109 107 Users can choose between a two- and a one-dimensional domain 110 108 decomposition. A 1D-decomposition is preferred on machines with a 111 slow network interconnection. In case of a 1D-decomposition, the 109 slow network interconnection. In case of a 1D-decomposition, 110 the 112 111 grid points along x direction are 113 112 distributed among the individual processors, but in y- and z-direction 114 113 all respective grid points belong to the same PE. </p> 115 <p style="line-height: 100%;">The calculation of central differences or 114 <p style="line-height: 100%;">The calculation of central 115 differences or 116 116 non-local arithmetic operations (e.g. global 117 117 sums, FFT) demands communication and an appropriate data exchange … … 126 126 y-direction, the data which lie distributed on the individual central 127 127 processing elements, have to be collected and/or relocated before. 128 This happens by means of the routine <tt>MPI_ALLTOALLV</tt>. Certain 128 This happens by means of the routine <tt>MPI_ALLTOALLV</tt>. 129 Certain 129 130 global operations like e.g. the search for absolute maxima or minima 130 131 within the 3D-arrays likewise require the employment of special MPI 131 132 routines (<tt>MPI_ALLREDUCE</tt>). </p> 132 <p style="line-height: 100%;">Further details of the internal model 133 <p style="line-height: 100%;">Further details of the 134 internal model 133 135 structure are described in the <a href="../tec/index.html">technical/numerical 134 136 documentation</a>. <br> 135 137 </p> 136 <hr><font color="#000080"><font color="#000080"><br> 137 <a href="chapter_1.0.html"><font color="#000080"><img width="32" 138 height="32" border="2" align="bottom" name="Grafik1" src="left.gif"></font></a><a 139 href="index.html"><font color="#000080"><img width="32" height="32" 140 border="2" align="bottom" name="Grafik2" src="up.gif"></font></a><a 141 href="chapter_3.0.html"><font color="#000080"><img width="32" 142 height="32" border="2" align="bottom" name="Grafik3" src="right.gif"></font></a><br> 143 </font></font><br> 144 <p style="line-height: 100%;"><span style="font-style: italic;">Last 145 change: </span>14/04/05 (SR)<font color="#000080"><font color="#000080"><br> 146 </font></font></p> 147 </body> 148 </html> 138 <hr><font color="#000080"><font color="#000080"><br><a href="chapter_1.0.html"><font color="#000080"><img name="Grafik1" src="left.gif" align="bottom" border="2" height="32" width="32"></font></a><a href="index.html"><font color="#000080"><img name="Grafik2" src="up.gif" align="bottom" border="2" height="32" width="32"></font></a><a href="chapter_3.0.html"><font color="#000080"><img name="Grafik3" src="right.gif" align="bottom" border="2" height="32" width="32"></font></a><br> 139 </font></font><br><p style="line-height: 100%;"><span style="font-style: italic;">Last 140 change: </span>$Id$<font color="#000080"><font color="#000080"><br> 141 </font></font></p></body></html>
Note: See TracChangeset
for help on using the changeset viewer.