Changes between Version 314 and Version 315 of doc/app/initialization_parameters


Ignore:
Timestamp:
Mar 20, 2017 10:27:18 AM (8 years ago)
Author:
hellstea
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • doc/app/initialization_parameters

    v314 v315  
    13571357Allowed values are '' 'dirichlet' '' (p(k=nz+1)= 0.0) or '' 'neumann' '' (p(k=nz+1)=p(k=nz)).\\\\
    13581358Simultaneous 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 \\\\
    13601359In 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.
     1360Note 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.
    13621361}}}
    13631362|----------------
     
    13851384}}}
    13861385{{{#!td style="vertical-align:top"
    1387 'initial_gradient'
     1386'initial_gradient'\\\\
     1387('nested' for nest domains)
    13881388}}}
    13891389{{{#!td
     
    14211421}}}
    14221422{{{#!td style="vertical-align:top"
    1423 'neumann'
     1423'neumann'\\\\
     1424('nested' for nest domains)
    14241425}}}
    14251426{{{#!td
     
    14561457}}}
    14571458{{{#!td style="vertical-align:top"
    1458 'initial_gradient'
     1459'initial_gradient'\\\\
     1460('nested' for nest domains)
    14591461}}}
    14601462{{{#!td
     
    14671469(up to k=nz the prognostic equation for the scalar concentration is solved).
    14681470
    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.
     1471When 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.\\\\
     1472In 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.\\\\
    14701473}}}
    14711474|----------------