[ff3d-users] problem with hydrostatic sample

[Robert Li]

Dear Stephane Del Pino,
  I found a hydrostatic sample on the internet that
you wrote on Mar.,21,2003. The title is Domain
Decomposition for the Couplex Exercise.

  I have some problem with it. 

1. When I run it, the error is as following
Generating Surface Meshes
./solver\Domain.hpp:138: Trying to access Object in
case of R^3 domain

2.I don't clear about the concepts of computational
domain and the fictitious domain. I think
computational domain is domain D, fictitious domain is
mesh, right?

3.We usually only set boundary conditions on border of
computational domain D. Here you also set boundary
condition on Mesh border, except for computational
domain<1,0,0>. Is this necessary or right?

Thanks!

Robert


>>>>>>>>
/*camera {location<-10000,-10000,0> look_at
<1900,250,150>}
light_source { <2, 4, -3> color rgb <1,1,0>}
background{ color rgb <1,1,1>}
*/
box { <18440,244,2000>, < 21680,264,4000>
pigment { color rgb <1,0,0> }}


>>>>>>>>>>>
//The mesh is uniform with 240  60  30 points.
double N = 30;
vector n = (8*N, 2*N, N);
vector a = (0,0,0);
//vector b = (25000,695*L,6000);
vector b = (25000,695,6000);
scene S=pov("dechet2.pov");
mesh Mesh = structured(n,a,b);
domain D= domain(S);
function dogger = one(y<=199.99);
function clay =one(y <= (295.01+
55.01*x/25000))*one(y>199.99);
function
limestone=one(y<=594.99)*one(y>(295.01+55.01*x/25000));
function marl = one(y>594.99);
function K = (3.1536e-5 * marl + 6.3072 * limestone
                        + 3.1536e-6*clay + 25.2288*dogger
)*(1-one(<1,0,0>));

function penal = 1-one(<1,0,0>);

solve(H) in D by Mesh
        
cg(maxiter=900,epsilon=1E-10),krylov(precond=ichol),memory(matrix=none)

{
 pde(H) 
        -div(K*grad(H))=0;
        dnu(H)=0 on <1,0,0>;
        H =289* dogger + 310* limestone
                + clay* (289 + (y-200)*21.0/150)
                + marl * (310 + (y - 595)*0.3) on Mesh xmax;
        H =286*dogger + 200*limestone
                + clay * (286 - (y-200)*86/95)
                + marl * (200 - (y - 595)*0.2) on Mesh xmin;
        H =180 + 160 * x/25000 on Mesh ymax;
};

//plot(medit,"H",H*penal,Mesh);
save(medit,"H",H*penal,Mesh,dos);
save(medit,"H",Mesh,dos);

__________________________________
Do you Yahoo!?
Yahoo! SiteBuilder - Free web site building tool. Try it!
http://webhosting.yahoo.com/ps/sb/