/[escript]/trunk/finley/test/python/convection.py
ViewVC logotype

Diff of /trunk/finley/test/python/convection.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2444 by gross, Wed May 27 08:45:55 2009 UTC revision 2445 by gross, Fri May 29 03:23:25 2009 UTC
# Line 52  C_P=1                           # heat c Line 52  C_P=1                           # heat c
52  K=1.                            # thermal conductivity  K=1.                            # thermal conductivity
53  CHI=0.                          # Taylor Quinny coefficient  CHI=0.                          # Taylor Quinny coefficient
54  MUE=None                        # elastic shear modulus  MUE=None                        # elastic shear modulus
55  TAU_Y=3*10**(2.5)               # Drucker-Prager cohesion factor  TAU_Y=5*10**(2.5)               # Drucker-Prager cohesion factor
56  BETA=0                          # Drucker-Prager friction factor  BETA=0                          # Drucker-Prager friction factor
57  TAU_0=2*10**(2.5)               # transition stress  TAU_0=2*10**(2.5)               # transition stress
58  N=3                             # power for power law  N=3                             # power for power law
# Line 65  ETA_N0=1.                       # viscos Line 65  ETA_N0=1.                       # viscos
65    
66  T_TOL=1.e-4                     # tolerance temperature transport  T_TOL=1.e-4                     # tolerance temperature transport
67  TOL2=0.1                        # tolerance for flow update at a timestep. (large value will do only one correction step)  TOL2=0.1                        # tolerance for flow update at a timestep. (large value will do only one correction step)
68  FLOW_TOL=1.e-5                  # tolerance for inconcompressible flow solver  FLOW_TOL=1.e-4                  # tolerance for inconcompressible flow solver
69  FLOW_SUB_TOL=1.e-10             # sub-tolerance for inconcompressible flow solver  FLOW_SUB_TOL=1.e-8             # sub-tolerance for inconcompressible flow solver
70  NUSSELT_FN="nusselt.csv"  NUSSELT_FN="nusselt.csv"
71  VERBOSE=True  VERBOSE=True
72  DT_VIS=T_END/500                # time distane between two visulaization files  DT_VIS=T_END/500                # time distane between two visulaization files
# Line 250  else: Line 250  else:
250    p=Ra*(x2[DIM-1]-0.5*x2[DIM-1]**2-0.5)    p=Ra*(x2[DIM-1]-0.5*x2[DIM-1]**2-0.5)
251    nusselt_file=FileWriter(NUSSELT_FN,append=False)    nusselt_file=FileWriter(NUSSELT_FN,append=False)
252    
253    p_last=p
254  x=dom.getX()  x=dom.getX()
255  #  #
256  #   set up heat problem:  #   set up heat problem:

Legend:
Removed from v.2444  
changed lines
  Added in v.2445

  ViewVC Help
Powered by ViewVC 1.1.26