Class Mountains
The Mountains class is defined by the following equations:
(1) w_i,aa+(1/eps)*w_i,33=0 where 0<eps<<1 and a=1,2 and w_i is the extension of the surface velocity where w_i(x_3=1)=v_i.
(2) Integration of topography PDE using Taylor-Galerkin upwinding to stabilize the advection terms
H^(t+dt)=H^t+dt*w_3+w_hat*dt*[(div(w_hat*H^t)+w_3)+(dt/2)+H^t],
where w_hat=w*[1,1,0], dt<0.5*d/max(w_i), d is a characteristic element size; H(x_3=1)=lambda (?)
|
|
__init__(self,
domain,
v,
eps=0.01,
z=1)
Sets up the level set method. |
|
|
|
|
getDomain(self)
Returns the domain. |
|
|
|
|
getDt(self)
Returns the time step value. |
|
|
|
|
getH(self)
Returns the mesh size. |
|
|
|
|
| setTolerance(self,
tolerance=0.001) |
|
|
|
|
update(self,
u=None,
H_t=None,
dt=None,
verbose=False)
Sets a new W and updates the H function. |
|
|
__init__(self,
domain,
v,
eps=0.01,
z=1)
(Constructor)
|
|
Sets up the level set method.
- Parameters:
domain - the domain where the mountains is usedeps - the parameter for (1)z - the height of the box
|
update(self,
u=None,
H_t=None,
dt=None,
verbose=False)
|
|
Sets a new W and updates the H function.
- Parameters:
|