]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Finishing touches
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 10 Sep 2006 04:09:13 +0000 (04:09 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 10 Sep 2006 04:09:13 +0000 (04:09 +0000)
git-svn-id: https://svn.dealii.org/trunk@13884 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-24/doc/intro.dox
deal.II/examples/step-24/doc/results.dox

index 96dd500d918d1d2685b3417289fbec6ee89f7299..441c30ac719b9491ba48883cf3005a0c3cad4863 100644 (file)
@@ -186,3 +186,13 @@ G_{1}+(k\, \theta)G_{2}-(c_0k)^2\theta (\theta F^{n}+(1-\theta)F^{n-1}) \\
 Mv^n & = & -(c_0^2\,k\, \theta\, A+c_0B)\bar{p}^{n}+ G_2 - 
 c_0^2k(\theta F^{n}+(1-\theta)F^{n-1}) 
 @f}
+
+
+<h3>What the program does</h3>
+
+Compared to @ref step_23 "step-23", this programs adds the treatment of a
+simple absorbing boundary conditions. In addition, it deals with data obtained
+from actual experimental measurements. To this end, we need to evaluate the
+solution at points at which the experiment also evaluates a real pressure
+field. We will see how to do that using the VectorTools::point_value function
+further down below.
index e997942ef4cd02326c5149b830b8baecb3cb29c8..45fe566502688b9af76b39a561c9f970332d22ef 100644 (file)
@@ -103,6 +103,14 @@ figures:
 </tr>
 </table>
 
+Note that in the experimental data, the first signal (i.e. the left-most dark
+line) results from absorption at the tissue boundary, and therefore reaches
+the detectors first and before any of the signals from the interior. This
+signal is also faintly visible at the end of the traces, around 30 $\my s$,
+which indicates that the signal travelled through the entire tissue to reach
+detectors at the other side, after all the signals originating from the
+interior have reached them.
+
 As before, the numerical result better matches experimental ones by applying a
 bandwidth filter that matches the actual behavior of detectors (left) and by
 choosing a finer mesh (right):
@@ -118,7 +126,19 @@ choosing a finer mesh (right):
 </tr>
 </table>
 
-
+One of the important differences between the left and the right figure is that
+the curves look much less "angular" at the right. The angularity comes from
+the fact that while waves in the continuous equation travel equally fast in
+all directions, this isn't the case after discretization: there, waves that
+travel diagonal to cells move at slightly different speeds to those that move
+parallel to mesh lines. This anisotropy leads to wave fronts that aren't
+perfectly circular (and would produce sinusoidal signals in the stacked
+plots), but are bulged out in certain directions. To make things worse, the
+circular mesh we use (see for example @ref step_6 "step-6" for a view of the
+coarse mesh) is not isotropic either. The net result is that the signal fronts
+are not sinusoidal unless the mesh is sufficiently fine. The right image is a
+lot better in this respect, though artificts in the form of trailing spurious
+waves can still be seen.
 
 
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.