Changes between Version 22 and Version 23 of doc/tec/noncyclic

Timestamp:

Feb 22, 2021 4:48:11 PM (4 years ago)

Author:

wagner

Comment:

--

doc/tec/noncyclic

@@ -1 +1 @@
 = Non-cyclic lateral boundary conditions =
 
-Figure 1 shows the grid structure for non-cyclic boundary conditions at the left/right boundary '''LB/RB''' ([../../app/inipar/#bc_lr bc_lr]) and figure 2 for non-cyclic boundary conditions at the north/south boundary '''NB/SB''' ([../../app/inipar/#bc_ns bc_ns]).
+Figure 1 shows the grid structure for non-cyclic boundary conditions at the left/right boundary '''LB/RB''' ([../../app/initialization_parameters/#bc_lr bc_lr]) and figure 2 for non-cyclic boundary conditions at the north/south boundary '''NB/SB''' ([../../app/initialization_parameters/#bc_ns bc_ns]).
 The indices (i,j,k) represent the directions (x,y,z).
 The model domain extends from -1:nx+1 in the x-direction, from -1:ny+1 in the y-direction and from 0:nzt+1 in the z-direction.
@@ -20 +20 @@
 In case of a Dirichlet condition, the values at i = 0 (j = 0) are taken from i = -1 (j = -1).
 In case of a radiation boundary condition, the solution of the Sommerfeld equation overwrites the prognostic values at i = 0 (j = 0). \\
-For non-cyclic lateral boundary conditions, the parameter [../../app/inipar/#psolver psolver] has to be set to '' 'multigrid' '' because the default FFT-solver can only be applied for cyclic boundary conditions. \\\\
+For non-cyclic lateral boundary conditions, the parameter [../../app/initialization_parameters/#psolver psolver] has to be set to '' 'multigrid' '' because the default FFT-solver can only be applied for cyclic boundary conditions. \\\\
 
 === Inflow boundary ===
@@ -30 +30 @@
  }}}
 t denotes the time, Δt the time step and s,,init,, the initialization profile of the scalar quantities which is constant in time.
-The quantities at the inflow are set by the initial vertical profiles (see [../../app/inipar/#initializing_actions initializing_actions]).
+The quantities at the inflow are set by the initial vertical profiles (see [../../app/initialization_parameters/#initializing_actions initializing_actions]).
 A Neumann condition is used for the subgrid-scale turbulent kinetic energy e (here e.g. for a left-right flow):
 {{{
@@ -51 +51 @@
 \end{cases} . \quad (4)$
  }}} 
-d,,f,, is a damping factor to control the damping intensity, and d,,w,, is the width of the relaxation region extending from the inflow. Quantities d,,f,, and d,,w,, can be set with parameters [../../app/inipar/#pt_damping_factor pt_damping_factor] and [../../app/inipar/#pt_damping_width pt_damping_width], respectively.
+d,,f,, is a damping factor to control the damping intensity, and d,,w,, is the width of the relaxation region extending from the inflow. Quantities d,,f,, and d,,w,, can be set with parameters [../../app/initialization_parameters/#pt_damping_factor pt_damping_factor] and [../../app/initialization_parameters/#pt_damping_width pt_damping_width], respectively.
 Both parameters have to be set by the user and must be adjusted case-by-case, because both parameters depend on the numerical and physical conditions, so that application of universal default values is not possible. 
 So far, we have experience with gravity waves in case of cold air outbreaks, which grow in amplitude up to quite extreme values, if no damping is applied.
-In the respective simulations, we used typical values for [../../app/inipar/#pt_damping_factor pt_damping_factor] of 0.05 and for [../../app/inipar/#pt_damping_width pt_damping_width] of 25 km in order to prevent the gravity waves from growing.
+In the respective simulations, we used typical values for [../../app/initialization_parameters/#pt_damping_factor pt_damping_factor] of 0.05 and for [../../app/initialization_parameters/#pt_damping_width pt_damping_width] of 25 km in order to prevent the gravity waves from growing.
 
 === Outflow boundary ===
@@ -61 +61 @@
 For the scalar quantities, Neumann boundary conditions are used at the outflow boundary which is the simplest way.
 For the velocity components, a Neumann condition would require to be considered in the solution of the Poisson equation for perturbation pressure, which has not been realized so far, because it requires some technical effort. 
-Instead, PALM offers a radiation boundary conditions for the velocity components, which are not in conflict with the pressure solver (see [../../app/inipar/#bc_lr bc_lr] and [../../app/inipar/#bc_ns bc_ns]).
+Instead, PALM offers a radiation boundary conditions for the velocity components, which are not in conflict with the pressure solver (see [../../app/initialization_parameters/#bc_lr bc_lr] and [../../app/initialization_parameters/#bc_ns bc_ns]).
 For the radiation condition, the Sommerfeld radiation equation is solved at the outflow
 {{{
@@ -146 +146 @@
 
 PALM offers the possibility of a mass flux correction at the outflow (e.g. Tian, 2004).
-If parameter [../../app/inipar/#conserve_volume_flow conserve_volume_flow] is set true, the mass flux at the inflow and outflow is calculated by:
+If parameter [../../app/initialization_parameters/#conserve_volume_flow conserve_volume_flow] is set true, the mass flux at the inflow and outflow is calculated by:
 {{{
 #!Latex
 $\dot{m} = \sum_{k=1}^{nz-1} \Delta z(k) \sum_{l=0}^{nx_i} \psi(l,k) \Delta x_i \; . \quad (12)$
 }}}
-where Δx,,i,, and ψ is equal to Δy (Δx) and u (v) in case of [../../app/inipar/#bc_lr bc_lr] ([../../app/inipar/#bc_ns bc_ns]). 
+where Δx,,i,, and ψ is equal to Δy (Δx) and u (v) in case of [../../app/initialization_parameters/#bc_lr bc_lr] ([../../app/initialization_parameters/#bc_ns bc_ns]). 
 The correction factor for the outflow velocity, which is necessary due to different mass fluxes at inflow and outflow, can be calculated by
 {{{