55 |
DIM=2 |
DIM=2 |
56 |
H=1. |
H=1. |
57 |
L=4*H |
L=4*H |
58 |
THETA=0.5 # time stepping THETA=0.5 cranck nicolson |
USE_BACKWARD_EULER=False |
59 |
TOL=1.e-4 |
TOL=1.e-4 |
60 |
PERTURBATION=0.1 |
PERTURBATION=0.1 |
61 |
DT=1.e-4 |
DT=1.e-4 |
180 |
# |
# |
181 |
# set up heat problem: |
# set up heat problem: |
182 |
# |
# |
183 |
heat=TemperatureCartesian(dom,theta=THETA) |
heat=TemperatureCartesian(dom,useBackwardEuler=USE_BACKWARD_EULER) |
184 |
heat.setTolerance(TOL*extratol) |
heat.setTolerance(TOL*extratol) |
185 |
|
|
186 |
fixed_T_at=whereZero(x[DIM-1])+whereZero(H-x[DIM-1]) |
fixed_T_at=whereZero(x[DIM-1])+whereZero(H-x[DIM-1]) |