Commit 9dc33b41 authored by Alexander Kissinger's avatar Alexander Kissinger
Browse files

New section for implicit spatial discretization methods shortly explaining the...

New section for implicit spatial discretization methods shortly explaining the box and cc method. Reviewed by Lena.

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@10237 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent e371876c
This diff is collapsed.
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!-- Created with Inkscape (http://www.inkscape.org/) -->
<svg
xmlns:dc="http://purl.org/dc/elements/1.1/"
xmlns:cc="http://creativecommons.org/ns#"
xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"
xmlns:svg="http://www.w3.org/2000/svg"
xmlns="http://www.w3.org/2000/svg"
xmlns:sodipodi="http://sodipodi.sourceforge.net/DTD/sodipodi-0.dtd"
xmlns:inkscape="http://www.inkscape.org/namespaces/inkscape"
width="1052.3622"
height="744.09448"
id="svg2"
sodipodi:version="0.32"
inkscape:version="0.48.3.1 r9886"
version="1.0"
sodipodi:docname="cc_disc.eps"
inkscape:output_extension="org.inkscape.output.svg.inkscape">
<defs
id="defs4">
<marker
inkscape:stockid="Arrow1Sstart"
orient="auto"
refY="0.0"
refX="0.0"
id="Arrow1Sstart"
style="overflow:visible">
<path
id="path3538"
d="M 0.0,0.0 L 5.0,-5.0 L -12.5,0.0 L 5.0,5.0 L 0.0,0.0 z "
style="fill-rule:evenodd;stroke:#000000;stroke-width:1.0pt;marker-start:none"
transform="scale(0.2) translate(6,0)" />
</marker>
<marker
inkscape:stockid="Arrow1Mstart"
orient="auto"
refY="0.0"
refX="0.0"
id="Arrow1Mstart"
style="overflow:visible">
<path
id="path3532"
d="M 0.0,0.0 L 5.0,-5.0 L -12.5,0.0 L 5.0,5.0 L 0.0,0.0 z "
style="fill-rule:evenodd;stroke:#000000;stroke-width:1.0pt;marker-start:none"
transform="scale(0.4) translate(10,0)" />
</marker>
<marker
inkscape:stockid="Arrow1Mend"
orient="auto"
refY="0.0"
refX="0.0"
id="Arrow1Mend"
style="overflow:visible;">
<path
id="path3727"
d="M 0.0,0.0 L 5.0,-5.0 L -12.5,0.0 L 5.0,5.0 L 0.0,0.0 z "
style="fill-rule:evenodd;stroke:#000000;stroke-width:1.0pt;marker-start:none;"
transform="scale(0.4) rotate(180) translate(10,0)" />
</marker>
<inkscape:perspective
sodipodi:type="inkscape:persp3d"
inkscape:vp_x="0 : 526.18109 : 1"
inkscape:vp_y="0 : 1000 : 0"
inkscape:vp_z="744.09448 : 526.18109 : 1"
inkscape:persp3d-origin="372.04724 : 350.78739 : 1"
id="perspective10" />
</defs>
<sodipodi:namedview
id="base"
pagecolor="#ffffff"
bordercolor="#666666"
borderopacity="1.0"
inkscape:pageopacity="0.0"
inkscape:pageshadow="2"
inkscape:zoom="2.96"
inkscape:cx="236.49949"
inkscape:cy="354.0982"
inkscape:document-units="px"
inkscape:current-layer="layer1"
showgrid="true"
inkscape:window-width="1280"
inkscape:window-height="963"
inkscape:window-x="0"
inkscape:window-y="17"
inkscape:window-maximized="0"
showguides="true"
inkscape:guide-bbox="true">
<inkscape:grid
type="xygrid"
id="grid2385"
visible="true"
enabled="true" />
<sodipodi:guide
orientation="0,1"
position="380,520"
id="guide4263" />
<sodipodi:guide
orientation="1,0"
position="360,400"
id="guide4265" />
<sodipodi:guide
orientation="0,1"
position="260,320"
id="guide4267" />
<sodipodi:guide
orientation="1,0"
position="160,420"
id="guide4269" />
</sodipodi:namedview>
<metadata
id="metadata7">
<rdf:RDF>
<cc:Work
rdf:about="">
<dc:format>image/svg+xml</dc:format>
<dc:type
rdf:resource="http://purl.org/dc/dcmitype/StillImage" />
<dc:title></dc:title>
</cc:Work>
<cc:Work
rdf:about="">
<dc:format>image/svg+xml</dc:format>
<dc:type
rdf:resource="http://purl.org/dc/dcmitype/StillImage" />
</cc:Work>
</rdf:RDF>
</metadata>
<g
inkscape:label="Ebene 1"
inkscape:groupmode="layer"
id="layer1">
<rect
style="fill:#00a3ff;fill-opacity:1;fill-rule:nonzero;stroke:none"
id="rect5186"
width="195.92877"
height="198.20135"
x="260"
y="122.09448"
ry="2.4040923" />
<g
id="g3686">
<rect
ry="1.2165556"
y="124.64518"
x="59.951752"
height="399.42368"
width="399.94101"
id="rect2383"
style="opacity:1;fill:none;fill-opacity:1;fill-rule:nonzero;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1" />
<path
id="path2389"
d="M 260,123.33807 L 258.8744,525.88098"
style="fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:3;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1;stroke-miterlimit:4;stroke-dasharray:none" />
<path
id="path2391"
d="M 60,324.09448 L 460,324.09448"
style="fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:3;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1;stroke-miterlimit:4;stroke-dasharray:none" />
</g>
<g
id="g5054" />
<flowRoot
xml:space="preserve"
id="flowRoot3695"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#000000;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
transform="translate(-487.86472,28.082659)"><flowRegion
id="flowRegion3697"><rect
id="rect3699"
width="242.61003"
height="29.817987"
x="440.4931"
y="25.751904" /></flowRegion><flowPara
id="flowPara3701"
style="font-size:24px">FV mesh</flowPara></flowRoot> <path
style="fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:3;stroke-linecap:butt;stroke-linejoin:miter;marker-end:url(#Arrow1Mend);stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1"
d="M 78.265175,87.068982 L 103.04528,117.84893"
id="path3711" />
<g
id="g5069" />
<flowRoot
xml:space="preserve"
id="flowRoot5158"
style="font-size:24px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#000000;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
transform="translate(28.86722,-96.7741)"><flowRegion
id="flowRegion5160"><rect
id="rect5162"
width="20.330421"
height="32.528717"
x="295.46921"
y="284.62631"
style="font-size:24px" /></flowRegion><flowPara
id="flowPara5164">i</flowPara></flowRoot> <path
sodipodi:type="arc"
style="fill:#000000;fill-opacity:1;fill-rule:nonzero;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;stroke-dashoffset:0"
id="path5188"
sodipodi:cx="261.58511"
sodipodi:cy="329.3533"
sodipodi:rx="13.553634"
sodipodi:ry="13.553634"
d="m 275.13875,329.3533 a 13.553634,13.553634 0 1 1 -27.10727,0 13.553634,13.553634 0 1 1 27.10727,0 z"
transform="translate(98.75746,-105.25882)" />
<g
id="g5398"
transform="translate(-415.05597,-263.03423)">
<flowRoot
transform="translate(364.22784,7.914998)"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#000000;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
id="flowRoot5362"
xml:space="preserve"><flowRegion
id="flowRegion5364"><rect
style="fill:#000000;fill-opacity:1"
y="493.35226"
x="509.61661"
height="54.214523"
width="66.412804"
id="rect5366" /></flowRegion><flowPara
id="flowPara5368">x</flowPara></flowRoot> <flowRoot
transform="translate(299.66602,81.480568)"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#000000;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
id="flowRoot5370"
xml:space="preserve"><flowRegion
id="flowRegion5372"><rect
style="fill:#000000;fill-opacity:1"
y="444.55917"
x="608.55817"
height="28.462643"
width="29.818018"
id="rect5374" /></flowRegion><flowPara
style="font-size:20px;fill:#000000;fill-opacity:1"
id="flowPara5376">ij</flowPara></flowRoot> </g>
<path
style="fill:none;stroke:#000000;stroke-width:3.95269799;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;marker-start:url(#Arrow1Sstart);marker-end:none"
d="m 359.00741,368.24356 -0.18282,-35.05902 0,0 0,0 0,0"
id="path2552"
inkscape:connector-curvature="0" />
<g
id="g4730"
transform="translate(-151.29942,-244.81931)">
<flowRoot
transform="translate(-9.52871,93.700078)"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#000000;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
id="flowRoot4694"
xml:space="preserve"><flowRegion
id="flowRegion4696"><rect
style="fill:#000000;fill-opacity:1"
y="493.35226"
x="509.61661"
height="54.214523"
width="66.412804"
id="rect4698" /></flowRegion><flowPara
style="font-weight:bold"
id="flowPara4700">n</flowPara></flowRoot> <flowRoot
transform="translate(-74.09053,167.26565)"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#000000;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
id="flowRoot4702"
xml:space="preserve"><flowRegion
id="flowRegion4704"><rect
style="fill:#000000;fill-opacity:1"
y="444.55917"
x="608.55817"
height="28.462643"
width="29.818018"
id="rect4706" /></flowRegion><flowPara
style="font-size:20px;fill:#000000;fill-opacity:1"
id="flowPara4708">ij</flowPara></flowRoot> <flowRoot
transform="translate(-76.92541,135.39594)"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#000000;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
id="flowRoot4710"
xml:space="preserve"><flowRegion
id="flowRegion4712"><rect
style="fill:#000000;fill-opacity:1"
y="444.55917"
x="608.55817"
height="28.462643"
width="29.818018"
id="rect4714" /></flowRegion><flowPara
style="font-size:20px;fill:#000000;fill-opacity:1"
id="flowPara4716" /></flowRoot> <flowRoot
xml:space="preserve"
id="flowRoot3215"
style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:times;-inkscape-font-specification:times"
transform="translate(-208.70058,302.81931)"><flowRegion
id="flowRegion3217"><rect
id="rect3219"
width="331.08109"
height="304.05405"
x="547.97296"
y="46.12151" /></flowRegion><flowPara
id="flowPara3221"></flowPara></flowRoot> </g>
<rect
style="fill:none;stroke:#0000ff;stroke-width:8;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;stroke-dashoffset:0"
id="rect5124"
width="199.23842"
height="197.88306"
x="260"
y="125.76817"
ry="1.2165556" />
<flowRoot
xml:space="preserve"
id="flowRoot5192-9"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#0000ff;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
transform="translate(-59.350824,-359.68247)"><flowRegion
id="flowRegion5194-6"><rect
id="rect5196-0"
width="83.980377"
height="54.214508"
x="509.61661"
y="493.35226"
style="fill:#0000ff;fill-opacity:1" /></flowRegion><flowPara
id="flowPara5198-1">cv</flowPara></flowRoot> <flowRoot
xml:space="preserve"
id="flowRoot5200-6"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#0000ff;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
transform="matrix(1,0,0,0.91106788,-103.91264,-246.15951)"><flowRegion
id="flowRegion5202-4"><rect
id="rect5204-4"
width="29.818018"
height="28.462643"
x="608.55817"
y="444.55917"
style="fill:#0000ff;fill-opacity:1" /></flowRegion><flowPara
id="flowPara5206-4"
style="font-size:20px;fill:#0000ff;fill-opacity:1">i</flowPara></flowRoot> <rect
style="fill:none;stroke:#0000ff;stroke-width:8;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;stroke-dashoffset:0"
id="rect5124-2"
width="199.23842"
height="197.88306"
x="260.7616"
y="324.2114"
ry="1.2165556" />
<flowRoot
xml:space="preserve"
id="flowRoot5158-2"
style="font-size:24px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#000000;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
transform="translate(30.52465,109.67227)"><flowRegion
id="flowRegion5160-0"><rect
id="rect5162-5"
width="20.330421"
height="32.528717"
x="295.46921"
y="284.62631"
style="font-size:24px" /></flowRegion><flowPara
id="flowPara5164-7">j</flowPara></flowRoot> <path
sodipodi:type="arc"
style="fill:#000000;fill-opacity:1;fill-rule:nonzero;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;stroke-dashoffset:0"
id="path5188-2"
sodipodi:cx="261.58511"
sodipodi:cy="329.3533"
sodipodi:rx="13.553634"
sodipodi:ry="13.553634"
d="m 275.13875,329.3533 a 13.553634,13.553634 0 1 1 -27.10727,0 13.553634,13.553634 0 1 1 27.10727,0 z"
transform="translate(98.86125,94.74118)" />
<text
xml:space="preserve"
style="font-size:10px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;text-align:start;line-height:125%;letter-spacing:0px;word-spacing:0px;writing-mode:lr-tb;text-anchor:start;fill:#000000;fill-opacity:1;stroke:none;font-family:times;-inkscape-font-specification:times"
x="775"
y="430.58096"
id="text3289"
sodipodi:linespacing="125%"><tspan
sodipodi:role="line"
id="tspan3291"
x="775"
y="430.58096" /></text>
<flowRoot
xml:space="preserve"
id="flowRoot5192-7"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#0000ff;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
transform="translate(-59.609375,-33.11923)"><flowRegion
id="flowRegion5194-7"><rect
id="rect5196-7"
width="83.980377"
height="54.214508"
x="509.61661"
y="493.35226"
style="fill:#0000ff;fill-opacity:1" /></flowRegion><flowPara
id="flowPara5198-8">cv</flowPara></flowRoot> <flowRoot
xml:space="preserve"
id="flowRoot5200-8"
style="font-size:40px;font-style:normal;font-weight:normal;text-align:center;text-anchor:middle;fill:#0000ff;fill-opacity:1;stroke:none;font-family:Bitstream Vera Sans"
transform="matrix(1,0,0,0.91106788,-104.17119,80.403726)"><flowRegion
id="flowRegion5202-5"><rect
id="rect5204-3"
width="29.818018"
height="28.462643"
x="608.55817"
y="444.55917"
style="fill:#0000ff;fill-opacity:1" /></flowRegion><flowPara
id="flowPara5206-8"
style="font-size:20px;fill:#0000ff;fill-opacity:1">j</flowPara></flowRoot> <flowRoot
xml:space="preserve"
id="flowRoot3332"
style="fill:black;stroke:none;stroke-opacity:1;stroke-width:1px;stroke-linejoin:miter;stroke-linecap:butt;fill-opacity:1;font-family:times;font-style:normal;font-weight:normal;font-size:10px;line-height:125%;letter-spacing:0px;word-spacing:0px;-inkscape-font-specification:times;font-stretch:normal;font-variant:normal;text-anchor:start;text-align:start;writing-mode:lr"><flowRegion
id="flowRegion3334"><rect
id="rect3336"
width="117.56757"
height="82.432434"
x="843.24323"
y="375.17557" /></flowRegion><flowPara
id="flowPara3338"></flowPara></flowRoot> <path
style="fill:none;stroke:#000000;stroke-width:2.76887465;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;marker-start:none;marker-end:url(#Arrow1Mend)"
d="m 474,266.14954 c -11.90811,-0.54837 -28.83138,3.83247 -42.3991,6.56559 -13.42199,2.70377 -29.8628,5.09037 -41.20783,11.37009 -21.99582,12.17512 -21.72525,14.81164 -26.69354,19.1751 -4.55198,3.99783 -3.51158,9.51744 -4.69953,14.20754"
id="path7526"
sodipodi:nodetypes="csssc"
inkscape:connector-curvature="0" />
<path
sodipodi:type="arc"
style="fill:#000000;fill-opacity:1;fill-rule:nonzero;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;stroke-dashoffset:0"
id="path5188-2-6"
sodipodi:cx="261.58511"
sodipodi:cy="329.3533"
sodipodi:rx="13.553634"
sodipodi:ry="13.553634"
d="m 275.13875,329.3533 a 13.553634,13.553634 0 1 1 -27.10727,0 13.553634,13.553634 0 1 1 27.10727,0 z"
transform="translate(-101.13875,94.741184)" />
<path
sodipodi:type="arc"
style="fill:#000000;fill-opacity:1;fill-rule:nonzero;stroke:#000000;stroke-width:3;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none;stroke-dashoffset:0"
id="path5188-4"
sodipodi:cx="261.58511"
sodipodi:cy="329.3533"
sodipodi:rx="13.553634"
sodipodi:ry="13.553634"
d="m 275.13875,329.3533 a 13.553634,13.553634 0 1 1 -27.10727,0 13.553634,13.553634 0 1 1 27.10727,0 z"
transform="translate(-101.24254,-105.25882)" />
<path
sodipodi:type="arc"
style="fill:#000000;fill-opacity:1;fill-rule:nonzero;stroke:none"
id="path7528"
sodipodi:cx="666.89191"
sodipodi:cy="460.98639"
sodipodi:rx="8.7837839"
sodipodi:ry="10.135136"
d="m 675.67569,460.98639 a 8.7837839,10.135136 0 1 1 -17.56757,0 8.7837839,10.135136 0 1 1 17.56757,0 z"
transform="translate(-307.40753,-136.44093)" />
</g>
</svg>
......@@ -137,7 +137,7 @@ unknowns in the system description. The
available primary variables are, e.\ g., saturations, mole/mass
fractions, temperature, pressures, etc.
\input{box}
\input{spatialdiscretization}
\section{Available models}
......
\section{Box method - A short introduction}\label{box}
\section{Implicit Spatial Discretization Schemes}\label{spatialdiscretization}
For the spatial discretization the so called BOX-method is used which unites the advantages of the finite-volume (FV) and finite-element (FE) methods.
For the implicit models there are two spatial discretization schemes (BOX and Cell Centered Finite Volume Method) available which are shortly introduced
in this section.
\subsection{Box method - A short introduction}\label{box}
The so called BOX-method is unites the advantages of the finite-volume (FV) and finite-element (FE) methods.
First, the model domain $G$ is discretized with a FE mesh consisting of nodes i and corresponding elements $E_k$. Then, a secondary FV mesh is constructed by connecting the midpoints and barycenters of the elements surrounding node i creating a box $B_i$ around node i (see Figure \ref{pc:box}a).
......@@ -118,3 +123,23 @@ The consideration of the time discretization and inserting $W_j = 1$ finally lea
\label{eq:discfin}
V_i \frac{\hat u_i^{n+1} - \hat u_i^{n}}{\Delta t} + \int_{\partial B_i} F(\tilde u^{n+1}) \cdot \mathbf n \; d{\varGamma}_{B_i} - V_i \: q^{n+1} \: = 0
\end{equation}
\subsection{Cell Centered Finite Volume Method - A short introduction}\label{cc}
\begin{figure} [h]
\centering
\includegraphics[width=0.4\linewidth,keepaspectratio]{EPS/cc_disc}
\caption{\label{pc:cc} Discretization of the Cell Centered Finite Volume Method}
\end{figure}
The cell centered finite volume method uses the elements of the grid as control volumes.
For each control volume all discrete values are determined at the element/control volume center (see Figure~\ref{pc:cc}).
The mass or energy fluxes are evaluated at the integration points ($x_{ij}$), which are located at the midpoints of the control
volume faces. This is a two point flux approximation since the flux between the element/control volume centers $i$ and $j$ is calculated
only with information from these two points. In contrast the BOX-method uses a multi-point flux approximation where all nodes of the
element influence the flux between two specific nodes. \\
Neumann boundary conditions are applied at the boundary control volume faces and Dirichlet boundary conditions at the boundary control volumes. \\
The cell centered finite volume method is robust and mass conservative but should only be applied for structured grids
(the control volume face normal vector ($n_{ij}$) should be parallel to the direction of the gradient between the two element/control
volume centers).
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment