Changes between Version 314 and Version 315 of doc/app/initialization_parameters
- Timestamp:
- Mar 20, 2017 10:27:18 AM (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/app/initialization_parameters
v314 v315 1357 1357 Allowed values are '' 'dirichlet' '' (p(k=nz+1)= 0.0) or '' 'neumann' '' (p(k=nz+1)=p(k=nz)).\\\\ 1358 1358 Simultaneous use of Neumann boundary conditions both at the top and bottom boundary ([#bc_p_b bc_p_b]) yields no consistent solution for the perturbation pressure in case that the multigrid method is used for solving the Poisson equation (see [#psolver psolver]), and should be avoided. Since at the bottom boundary the Neumann condition is a good choice (see [#bc_p_b bc_p_b]), in that case, a Dirichlet condition should be set at the top boundary.\\\\ 1359 \\\\1360 1359 In case of nested run, the default is not 'dirichlet', but 'neumann' instead.\\\\ 1361 Note that therefore the Poisson equation for perturbation pressure has all-Neumann boundary condition in nest-domains. This means that the solution is only determined up to a constant. Or in other words, the average value of the pressure is arbitrary. This implies that the convergence rate is usually somewhat lower than when Dirichlet condition is used at least on one of the boundaries. Users should pay attention to this and adjust the steering parameters of the multigrid solver carefully to achieve sufficient convergence without sacrificing too much computing time for the perturbation-pressure solution.1360 Note that therefore the Poisson equation for perturbation pressure has all-Neumann boundary condition in nest-domains. This means that the solution is only determined up to a constant. Or in other words, the average value of the pressure is arbitrary. This implies that the convergence rate is usually somewhat lower than when Dirichlet condition is used at least on one of the boundaries. Users using the nesting system should pay attention to this and adjust the steering parameters of the multigrid solver carefully to achieve sufficient convergence without sacrificing too much computing time for the perturbation-pressure solution. 1362 1361 }}} 1363 1362 |---------------- … … 1385 1384 }}} 1386 1385 {{{#!td style="vertical-align:top" 1387 'initial_gradient' 1386 'initial_gradient'\\\\ 1387 ('nested' for nest domains) 1388 1388 }}} 1389 1389 {{{#!td … … 1421 1421 }}} 1422 1422 {{{#!td style="vertical-align:top" 1423 'neumann' 1423 'neumann'\\\\ 1424 ('nested' for nest domains) 1424 1425 }}} 1425 1426 {{{#!td … … 1456 1457 }}} 1457 1458 {{{#!td style="vertical-align:top" 1458 'initial_gradient' 1459 'initial_gradient'\\\\ 1460 ('nested' for nest domains) 1459 1461 }}} 1460 1462 {{{#!td … … 1467 1469 (up to k=nz the prognostic equation for the scalar concentration is solved). 1468 1470 1469 When a constant scalar flux is used at the top boundary ([#top_scalarflux top_scalarflux]), bc_s_t = 'neumann' must be used, because otherwise the resolved scale may contribute to the top flux so that a constant value cannot be guaranteed. 1471 When a constant scalar flux is used at the top boundary ([#top_scalarflux top_scalarflux]), bc_s_t = 'neumann' must be used, because otherwise the resolved scale may contribute to the top flux so that a constant value cannot be guaranteed.\\\\ 1472 In case of nested run the default value of '''bc_s_t''' in the nest domains is not '' 'initial_gradient' '' but '' 'nested' '' instead. For the root domain of a nested run the default is '' 'initial_gradient' '' as usually.\\\\ 1470 1473 }}} 1471 1474 |----------------