-
-
Notifications
You must be signed in to change notification settings - Fork 36
/
Copy pathCONTROLS.html
6 lines (5 loc) · 23.8 KB
/
CONTROLS.html
1
2
3
4
5
6
<html><head><link rel="stylesheet" type="text/css" href="style.css"/></head><body> <H2> <BR> *CONTROLS </H2> <P> Keyword type: step <P> This option is used to change the iteration control parameters. It should only be used by those users who know what they are doing and are expert in the field. A detailed description of the convergence criteria is given in Section 6.10. There are two, mutually exclusive parameter: PARAMETERS and RESET. The RESET parameter resets the control parameters to their defaults. The parameter PARAMETERS is used to change the defaults. It can take the value TIME INCREMENTATION, FIELD, LINE SEARCH, NETWORK, CFD or CONTACT. If the TIME INCREMENTATION value is selected, the number of iterations before certain actions are taken (e.g. the number of divergent iterations before the increment is reattempted) can be changed and effect of these actions (e.g. the increment size is divided by two). The FIELD parameter can be used to change the convergence criteria themselves. <P> LINE SEARCH can be used to change the line search parameters (only for face-to-face penalty contact). The line search parameter scales the correction to the solution calculated by the Newton-Raphson algorithm such that the residual force is orthogonal to the correction. This requires the solution of a nonlinear equation, and consequently an iterative procedure. In CalculiX this procedure is approximated by a linear connection between: <UL> <LI>the scalar product of the residual force from the last iteration with the solution correction in the present iteration (corresponds to a line search parameter of zero) and </LI> <LI>the scalar product of the residual force in the present iteration with the solution correction in the present iteration (corresponds to a line search parameter of one). </LI> </UL> For details of the line seach algorithm the reader is referred to [103]. <P> With the NETWORK parameter the convergence criteria for network iterations can be changed. The parameters <SPAN CLASS="MATH"><B><IMG WIDTH="23" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2190.png" ALT="$ c_{1t}$"></B></SPAN>, <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2191.png" ALT="$ c_{1f}$"></B></SPAN> and <SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2192.png" ALT="$ c_{1p}$"></B></SPAN> express the fraction of the mean energy balance, mass balance and element balance the energy balance residual, the mass balance residual and the element balance residual is not allowed to exceed, respectively. The parameters <SPAN CLASS="MATH"><B><IMG WIDTH="23" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2201.png" ALT="$ c_{2t}$"></B></SPAN>, <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2202.png" ALT="$ c_{2f}$"></B></SPAN>, <SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2203.png" ALT="$ c_{2p}$"></B></SPAN> and <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2204.png" ALT="$ c_{2a}$"></B></SPAN> is the fraction of the change in temperature, mass flow, pressure and geometry since the beginning of the increment the temperature, mass flow, pressure and geometry change in the actual network iteration is not allowed to exceed, respectively. The same applies to the parameters <SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2195.png" ALT="$ a_{1t}$"></B></SPAN>, <SPAN CLASS="MATH"><B><IMG WIDTH="27" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2196.png" ALT="$ a_{1f}$"></B></SPAN>, <SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2197.png" ALT="$ a_{1p}$"></B></SPAN>, <SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2207.png" ALT="$ a_{2t}$"></B></SPAN>, <SPAN CLASS="MATH"><B><IMG WIDTH="27" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2208.png" ALT="$ a_{2f}$"></B></SPAN>, <SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2209.png" ALT="$ a_{2p}$"></B></SPAN> and <SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2212.png" ALT="$ a_{2a}$"></B></SPAN>, except that they are absolute values and not fractions, e.g. the mean enery balance residual should not exceed <SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2195.png" ALT="$ a_{1t}$"></B></SPAN> etc. Therefore they have appropriate units. <P> With the CFD parameter the maximum number of iterations in certain fluid loops can be influenced. A fluid calculation within CalculiX is triggered at the start of a new mechanical increment. This increment is subdivided into fluid increments based on the physical fluid properties. For each fluid increment iterations are performed. Usually, iterations are performed until convergence of the fluid increment or until the maximum allowed number of iterations is reached. This is the first parameter <SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="15" ALIGN="BOTTOM" BORDER="0" SRC="img2257.png" ALT="$ iitt$"></B></SPAN> (“transient”). In fluid calculations the unknowns in the equation systems are the quantities (velocity..) at the element centers. The values at the face centers and the gradients are calculated based on these element center quantities. In case the mesh is not orthogonal, iterations have to be performed. The number of these iterations is expressed by <SPAN CLASS="MATH"><B><IMG WIDTH="29" HEIGHT="31" ALIGN="MIDDLE" BORDER="0" SRC="img2258.png" ALT="$ iitg$"></B></SPAN> (“geometry”) and <SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="31" ALIGN="MIDDLE" BORDER="0" SRC="img2259.png" ALT="$ iitp$"></B></SPAN> (taking non-orthogonality into account in the pressure correction equation, “pressure”). This is the second and third parameter. For a perfectly rectangular grid these values can be set to zero. Finally, the parameter <SPAN CLASS="MATH"><B><IMG WIDTH="22" HEIGHT="31" ALIGN="MIDDLE" BORDER="0" SRC="img2260.png" ALT="$ jit$"></B></SPAN> specifies how many coupled pressure-temperature iterations have to be performed. For incompressible flow the default value of 1 should not be changed. For inviscid compressible flow this value may have to be increased up to 4, whereas for viscid compressible flow this value has rarely to be changed. <P> Finally, the CONTACT parameter is used to change defaults in the face-to-face penalty contact convergence algorithm (cf. Section 6.10.2). This relates to <UL> <LI>the maximum relative difference in number of contact elements to allow for convergence (delcon). The corresponding absolute difference, which may not be exceeded is defined as the number of contact elements in the previous iteration times delcon. </LI> <LI>the fraction of contact elements which is removed in an aleatoric way before repeting an increment in case of a local mimimum in the solution (alea) </LI> <LI>the integer factor by which the normal spring stiffness (in case of linear pressure-overclosure) and stick slope are reduced in case of divergence or too slow convergence (kscalemax) </LI> <LI>the maximum number of iterions per increment (itf2f). </LI> </UL> <P> First line: <UL> <LI>*CONTROLS </LI> <LI>Enter the PARAMETERS parameter and its value, or the RESET parameter. </LI> </UL> <P> There are no subsequent lines if the parameter RESET is selected. <P> Following lines if PARAMETERS=TIME INCREMENTATION is selected: <P> Second line: <UL> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="18" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2177.png" ALT="$ I_0$"></B></SPAN> iteration after which a check is made whether the residuals increase in two consecutive iterations (default: 4). If so, the increment is reattempted with <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2261.png" ALT="$ D_f$"></B></SPAN> times its size. </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2178.png" ALT="$ I_R$"></B></SPAN> iteration after which a logarithmic convergence check is performed in each iteration (default: 8). If more than <SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2174.png" ALT="$ I_C$"></B></SPAN> iterations are needed, the increment is reattempted with <SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2262.png" ALT="$ D_C$"></B></SPAN> its size. </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2263.png" ALT="$ I_P$"></B></SPAN> iteration after which the residual tolerance <!-- MATH $R^\alpha_p$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2264.png" ALT="$ R^\alpha_p$"></B></SPAN> is used instead of <!-- MATH $R^\alpha_n$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2265.png" ALT="$ R^\alpha_n$"></B></SPAN> (default: 9). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2174.png" ALT="$ I_C$"></B></SPAN> maximum number of iterations allowed (default: 16). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="20" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2266.png" ALT="$ I_L$"></B></SPAN> number of iterations after which the size of the subsequent increment will be reduced (default: 10). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2267.png" ALT="$ I_G$"></B></SPAN> maximum number of iterations allowed in two consecutive increments for the size of the next increment to be increased (default: 4). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="20" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2268.png" ALT="$ I_S$"></B></SPAN> Currently not used. </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2269.png" ALT="$ I_A$"></B></SPAN> Maximum number of cutbacks per increment (default: 5). A cutback is a reattempted increment. </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="20" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2270.png" ALT="$ I_J$"></B></SPAN> Currently not used. </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2271.png" ALT="$ I_T$"></B></SPAN> Currently not used. </LI> </UL> <P> Third line: <UL> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2261.png" ALT="$ D_f$"></B></SPAN> Cutback factor if the solution seems to diverge(default: 0.25). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2262.png" ALT="$ D_C$"></B></SPAN> Cutback factor if the logarithmic extrapolation predicts too many iterations (default: 0.5). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2272.png" ALT="$ D_B$"></B></SPAN> Cutback factor for the next increment if more than <SPAN CLASS="MATH"><B><IMG WIDTH="20" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2266.png" ALT="$ I_L$"></B></SPAN> iterations were needed in the current increment (default: 0.75). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="27" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2273.png" ALT="$ D_A$"></B></SPAN> Cutback factor if the temperature change in two subsequent increments exceeds DELTMX (default: 0.85). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2274.png" ALT="$ D_S$"></B></SPAN> Currently not used. </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="29" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2275.png" ALT="$ D_H$"></B></SPAN> Currently not used. </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2276.png" ALT="$ D_D$"></B></SPAN> Factor by which the next increment will be increased if less than <SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2267.png" ALT="$ I_G$"></B></SPAN> iterations are needed in two consecutive increments (default: 1.5). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="30" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2277.png" ALT="$ W_G$"></B></SPAN> Currently not used. </LI> </UL> <P> Following line if PARAMETERS=FIELD is selected: <P> Second line: <UL> <LI><!-- MATH $R^{\alpha}_n$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2278.png" ALT="$ R^{\alpha}_n$"></B></SPAN> Convergence criterion for the ratio of the largest residual to the average force (default: 0.005). The average force is defined as the average over all increments in the present step of the instantaneous force. The instantaneous force in an increment is defined as the mean of the absolute value of the nodal force components within all elements. </LI> <LI><!-- MATH $C^{\alpha}_n$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2279.png" ALT="$ C^{\alpha}_n$"></B></SPAN> Convergence criterion for the ratio of the largest solution correction to the largest incremental solution value (default: 0.01). </LI> <LI><!-- MATH $q^{\alpha}_0$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="20" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2280.png" ALT="$ q^{\alpha}_0$"></B></SPAN> Initial value at the start of a new step of the time average force (default: the time average force from the previous steps or 0.01 for the first step). </LI> <LI><!-- MATH $q^{\alpha}_u$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2281.png" ALT="$ q^{\alpha}_u$"></B></SPAN> user-defined average force. If defined, the calculation of the average force is replaced by this value. </LI> <LI><!-- MATH $R^{\alpha}_p$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2282.png" ALT="$ R^{\alpha}_p$"></B></SPAN> Alternative residual convergence criterion to be used after <SPAN CLASS="MATH"><B><IMG WIDTH="21" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2263.png" ALT="$ I_P$"></B></SPAN> iterations instead of <!-- MATH $R^\alpha_n$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2265.png" ALT="$ R^\alpha_n$"></B></SPAN> (default: 0.02). </LI> <LI><!-- MATH $\epsilon^{\alpha}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="19" HEIGHT="15" ALIGN="BOTTOM" BORDER="0" SRC="img2283.png" ALT="$ \epsilon^{\alpha}$"></B></SPAN> Criterion for zero flux relative to <!-- MATH $q^{\alpha}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="20" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2284.png" ALT="$ q^{\alpha}$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="37" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2285.png" ALT="$ 10^{-5}$"></B></SPAN>). </LI> <LI><!-- MATH $C^{\alpha}_{\epsilon}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2286.png" ALT="$ C^{\alpha}_{\epsilon}$"></B></SPAN> Convergence criterion for the ratio of the largest solution correction to the largest incremental solution value in case of zero flux (default: <SPAN CLASS="MATH"><B><IMG WIDTH="37" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img925.png" ALT="$ 10^{-3}$"></B></SPAN>). </LI> <LI><!-- MATH $R^{\alpha}_l$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="30" ALIGN="MIDDLE" BORDER="0" SRC="img2287.png" ALT="$ R^{\alpha}_l$"></B></SPAN> Convergence criterion for the ratio of the largest residual to the average force for convergence in a single iteration (default: <SPAN CLASS="MATH"><B><IMG WIDTH="37" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2288.png" ALT="$ 10^{-8}$"></B></SPAN>). </LI> </UL> <P> Following line if PARAMETERS=LINE SEARCH is selected: <P> Second line: <UL> <LI>not used. </LI> <LI><!-- MATH $s^{ls}_{max}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="37" HEIGHT="38" ALIGN="MIDDLE" BORDER="0" SRC="img2289.png" ALT="$ s^{ls}_{max}$"></B></SPAN> Maximum value of the line search parameter (default: 1.01). </LI> <LI><!-- MATH $s^{ls}_{min}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="36" HEIGHT="38" ALIGN="MIDDLE" BORDER="0" SRC="img2290.png" ALT="$ s^{ls}_{min}$"></B></SPAN> Minimum value of the line search parameter (default: 0.25). </LI> <LI>not used. </LI> <LI>not used. </LI> </UL> <P> Following line if PARAMETERS=NETWORK is selected: <P> Second line: <UL> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="23" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2190.png" ALT="$ c_{1t}$"></B></SPAN> (default: <!-- MATH $5 \cdot 10^{-7}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="56" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2291.png" ALT="$ 5 \cdot 10^{-7}$"></B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2191.png" ALT="$ c_{1f}$"></B></SPAN> (default: <!-- MATH $5 \cdot 10^{-7}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="56" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2291.png" ALT="$ 5 \cdot 10^{-7}$"></B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2192.png" ALT="$ c_{1p}$"></B></SPAN> (default: <!-- MATH $5 \cdot 10^{-7}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="56" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2291.png" ALT="$ 5 \cdot 10^{-7}$"></B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="23" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2201.png" ALT="$ c_{2t}$"></B></SPAN> (default: <!-- MATH $5 \cdot 10^{-7}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="56" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2291.png" ALT="$ 5 \cdot 10^{-7}$"></B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2202.png" ALT="$ c_{2f}$"></B></SPAN> (default: <!-- MATH $5 \cdot 10^{-7}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="56" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2291.png" ALT="$ 5 \cdot 10^{-7}$"></B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2203.png" ALT="$ c_{2p}$"></B></SPAN> (default: <!-- MATH $5 \cdot 10^{-7}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="56" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2291.png" ALT="$ 5 \cdot 10^{-7}$"></B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="25" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2204.png" ALT="$ c_{2a}$"></B></SPAN> (default: <!-- MATH $5 \cdot 10^{-7}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="56" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2291.png" ALT="$ 5 \cdot 10^{-7}$"></B></SPAN>). </LI> </UL> <P> Third line: <UL> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2195.png" ALT="$ a_{1t}$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="33" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2198.png" ALT="$ 10^{20}$"></B></SPAN>[M][L]<SPAN CLASS="MATH"><B><IMG WIDTH="11" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img6.png" ALT="$ ^2$"></B></SPAN>/[t]<SPAN CLASS="MATH"><B><IMG WIDTH="11" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img7.png" ALT="$ ^3$"></B></SPAN>; unit in SI: Watt). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="27" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2196.png" ALT="$ a_{1f}$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="33" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2198.png" ALT="$ 10^{20}$"></B></SPAN>[M]/[t]; unit in SI: kg/s). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2197.png" ALT="$ a_{1p}$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="33" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2198.png" ALT="$ 10^{20}$"></B></SPAN>[-]; dimensionless). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2207.png" ALT="$ a_{2t}$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="33" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2198.png" ALT="$ 10^{20}$"></B></SPAN>[T]; unit in SI: K). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="27" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2208.png" ALT="$ a_{2f}$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="33" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2198.png" ALT="$ 10^{20}$"></B></SPAN>[M]/[t]; unit in SI: kg/s). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2209.png" ALT="$ a_{2p}$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="33" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2198.png" ALT="$ 10^{20}$"></B></SPAN>[M]/([t]<SPAN CLASS="MATH"><B><IMG WIDTH="11" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img6.png" ALT="$ ^2$"></B></SPAN>[L]); unit in SI: Pa). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2212.png" ALT="$ a_{2a}$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="33" HEIGHT="18" ALIGN="BOTTOM" BORDER="0" SRC="img2198.png" ALT="$ 10^{20}$"></B></SPAN>[L]; unit in SI: m). </LI> </UL> <P> Here, [M], [L], [T] and [t] are the units for mass, length, temperature and time. <P> Following line if PARAMETERS=CFD is selected: <P> Second line: <UL> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="26" HEIGHT="15" ALIGN="BOTTOM" BORDER="0" SRC="img2257.png" ALT="$ iitt$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="20" HEIGHT="14" ALIGN="BOTTOM" BORDER="0" SRC="img2292.png" ALT="$ 20$"></B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="29" HEIGHT="31" ALIGN="MIDDLE" BORDER="0" SRC="img2258.png" ALT="$ iitg$"></B></SPAN> (default: <SPAN CLASS="MATH"><B>0</B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="31" ALIGN="MIDDLE" BORDER="0" SRC="img2259.png" ALT="$ iitp$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="12" HEIGHT="13" ALIGN="BOTTOM" BORDER="0" SRC="img1508.png" ALT="$ 1$"></B></SPAN>). </LI> <LI><SPAN CLASS="MATH"><B><IMG WIDTH="22" HEIGHT="31" ALIGN="MIDDLE" BORDER="0" SRC="img2260.png" ALT="$ jit$"></B></SPAN> (default: <SPAN CLASS="MATH"><B><IMG WIDTH="12" HEIGHT="13" ALIGN="BOTTOM" BORDER="0" SRC="img1508.png" ALT="$ 1$"></B></SPAN>). </LI> </UL> <P> Following line if PARAMETERS=CONTACT is selected: <P> Second line: <UL> <LI>delcon (<SPAN CLASS="MATH"><B><IMG WIDTH="29" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2293.png" ALT="$ \ge 0$"></B></SPAN>; default: <SPAN CLASS="MATH"><B><IMG WIDTH="40" HEIGHT="14" ALIGN="BOTTOM" BORDER="0" SRC="img2294.png" ALT="$ 0.001$"></B></SPAN>). </LI> <LI>alea (<SPAN CLASS="MATH"><B><IMG WIDTH="29" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2295.png" ALT="$ 0 \le$"></B></SPAN> alea <SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2296.png" ALT="$ \le 1$"></B></SPAN>; default: <SPAN CLASS="MATH"><B><IMG WIDTH="24" HEIGHT="14" ALIGN="BOTTOM" BORDER="0" SRC="img2297.png" ALT="$ 0.1$"></B></SPAN>). </LI> <LI>kscalemax (<SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2298.png" ALT="$ \ge 1$"></B></SPAN>, integer; default: <SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="14" ALIGN="BOTTOM" BORDER="0" SRC="img2251.png" ALT="$ 100$"></B></SPAN>). </LI> <LI>itf2f (<SPAN CLASS="MATH"><B><IMG WIDTH="28" HEIGHT="29" ALIGN="MIDDLE" BORDER="0" SRC="img2298.png" ALT="$ \ge 1$"></B></SPAN>, integer; default: <SPAN CLASS="MATH"><B><IMG WIDTH="20" HEIGHT="14" ALIGN="BOTTOM" BORDER="0" SRC="img2299.png" ALT="$ 60$"></B></SPAN>). </LI> </UL> <P> <PRE>
Example:
*CONTROLS,PARAMETERS=FIELD
1.e30,1.e30,0.01,,0.02,1.e-5,1.e-3,1.e-8
</PRE> <P> leads to convergence in just one iteration since nearly any residuals are accepted for convergence (<!-- MATH $R^{\alpha}_n=10^{30}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="76" HEIGHT="33" ALIGN="MIDDLE" BORDER="0" SRC="img2300.png" ALT="$ R^{\alpha}_n=10^{30}$"></B></SPAN> and <!-- MATH $C^{\alpha}_n=10^{30}$ --> <SPAN CLASS="MATH"><B><IMG WIDTH="76" HEIGHT="33" ALIGN="MIDDLE" BORDER="0" SRC="img2301.png" ALT="$ C^{\alpha}_n=10^{30}$"></B></SPAN>. <P> <P><P> <BR> Example files: beammrco. <P> </body></html>