]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fixed docs of step-34.
authorheltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 30 Jan 2013 16:57:46 +0000 (16:57 +0000)
committerheltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 30 Jan 2013 16:57:46 +0000 (16:57 +0000)
git-svn-id: https://svn.dealii.org/trunk@28197 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-34/doc/intro.dox
deal.II/examples/step-34/doc/results.dox
deal.II/examples/step-34/parameters.prm

index 71cfc948e94bf52511d1e7de092c79626c3b4d8a..420350d0b679264c1fedc0b08286cea2e6e4698f 100644 (file)
@@ -1,19 +1,11 @@
 <br>
 
-<i>This program was contributed by Luca Heltai.
-</i>
+<i>This program was contributed by Luca Heltai (thanks to Michael
+Gratton for pointing out what the exact solution should have been in
+the threedimensional case).  </i>
 
 <a name="Intro"></a>
 
-<h1> Warning! </h1>
-
-The three dimensional solution generated by this example program is
-wrong. There is a factor of $1/2$ that appears, and we haven't figured
-out where it comes from yet. For the moment this issue is "solved" by
-multiplying the wind function in the parameter file by two.
-
-If you think you have spotted the mistake, please let us know.
-
 <h1>Introduction</h1>
 
 <h3> Irrotational flow </h3>
@@ -666,3 +658,16 @@ files in the current directory and an object of type HyperBallBoundary
 will then be attached to the triangulation to allow mesh refinement
 that respects the continuous geometry behind the discrete initial
 mesh.
+
+For a sphere of radius $a$ translating at a velocity of $U$ in the $x$ direction, the potential reads
+\[ 
+\phi = -\frac{1}{2}U \left(\frac{a}{r}\right)3 r \cos\theta
+
+\] see, e.g. J.N. Newman, \emph{Marine Hydrodynamics}, 1977,
+pp. 127. For unit speed and radius, and restricting $(x,y,z)$ to lie
+on the surface of the sphere, 
+\[ \phi = -x/2.\] In the test problem,
+the flow is $(1,1,1)$, so the appropriate exact solution on the
+surface of the sphere is the superposition of the above solution with
+the analogous solution along the $y$ and $z$ axes, or \[ \phi =
+\frac{1}{2}(x + y + z) \]
\ No newline at end of file
index d534ed9800ae77b75ef7cf38a2475f1ced4270be..e4c92e08bb1df660cebf8603a00aef5ae088b373 100644 (file)
@@ -32,7 +32,7 @@ subsection Exact solution 3d
 
   # Separate vector valued expressions by ';' as ',' is used internally by the
   # function parser.
-  set Function expression = x+y+z   # default: 0
+  set Function expression = .5*(x+y+z)   # default: 0
 
   # The name of the variables as they will be used in the function, separated
   # by ','.
@@ -76,7 +76,7 @@ subsection Wind function 3d
 
   # Separate vector valued expressions by ';' as ',' is used internally by the
   # function parser.
-  set Function expression = 2; 2; 2 # default: 0; 0; 0
+  set Function expression = 1; 1; 1 # default: 0; 0; 0
 
   # The name of the variables as they will be used in the function, separated
   # by ','.
index 5e12caf422f9eb48ae23d0b519150e8248730cc3..60edd7c35ed36e24eaf4f06c906e43fef272d0b9 100644 (file)
@@ -27,7 +27,7 @@ subsection Exact solution 3d
 
   # Separate vector valued expressions by ';' as ',' is used internally by the
   # function parser.
-  set Function expression = x+y+z   # default: 0
+  set Function expression = .5*(x+y+z)
 
   # The name of the variables as they will be used in the function, separated
   # by ','.
@@ -71,7 +71,7 @@ subsection Wind function 3d
 
   # Separate vector valued expressions by ';' as ',' is used internally by the
   # function parser.
-  set Function expression = 2; 2; 2 # default: 0; 0; 0
+  set Function expression = 1; 1; 1
 
   # The name of the variables as they will be used in the function, separated
   # by ','.

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.