1 0 GfsAxi GfsBox GfsGEdge {} {
Time { end = 0 }
PhysicalParams { L = 50 }
AdvectionParams { scheme = none }
ApproxProjectionParams { tolerance = 1e-10 }
Refine 4
Refine (LEVEL + 1./50.*(x*x + y*y)*(4. - LEVEL))
Global {
}
Solid (ellipse (0., 0., A0, A0))
Init {} {
U = U0
Phi = {
double r = sqrt (cx*cx + cy*cy);
return U0*A0*A0*A0*cx/(2.*r*r*r);
}
}
OutputErrorNorm { start = end } { awk '{ print LEVEL " " $7 " " $9}' } { v = U } {
s = (dx("Phi") + 1.)
}
OutputSimulation { start = end } sim-LEVEL.gfs
}
GfsBox {
left = Boundary { BcDirichlet U U0 }
right = Boundary { BcDirichlet U U0 }
}