Report from Simulation, class Manager

by Kåre Dump

Report file

/auto/tor/NO/NO/Work/Diffpack-3.5.00/NO/doc/Book/src/app/Pipeflow/Verify/SIMULATION-report.html

generated

Mon Feb 5 15:03:37 2001.

Created by user

unknown user

running the program

/auto/tor/NO/NO/Work/Diffpack-3.5.00/NO/doc/Book/src/app/Pipeflow/../app
(last modified , bytes)

on the i386-linux host

ampc-218.si.sintef.no.

PURPOSE:

Contents

• 1 Momentum and energy equation administration - simulation #1
  
  • 1.1 Menu items that have been varied in the multiple loop
  • 1.2 Parameters in constitutive relations
  • 1.3 Grid and iteration parameters
  
  
• 2 Momentum equation - simulation #1
  
  • 2.1 Problem dependent parameters
  • 2.2 Summary of the equation solver
  • 2.3 Summary of the nonlinear equation solver
    • 2.3.1 Nonlinear solution strategy 1 (gives values from last switch)
    • 2.3.2 Nonlinear solution strategy 2 (gives values from last switch)
    
    
  • 2.4 Visualization of the velocity field
  
  
• 3 Energy equation - simulation #1
  
  • 3.1 Other parameters
  • 3.2 Summary of the equation solver
  • 3.3 Summary of the nonlinear equation solver
    • 3.3.1 Nonlinear solution strategy 1 (gives values from last switch)
    • 3.3.2 Nonlinear solution strategy 2 (gives values from last switch)
    
    
  • 3.4 Visualization of the temperature field
  
  
• 4 Momentum and energy equation administration - simulation #2
  
  • 4.1 Menu items that have been varied in the multiple loop
  • 4.2 Parameters in constitutive relations
  • 4.3 Grid and iteration parameters
  
  
• 5 Momentum equation - simulation #2
  
  • 5.1 Problem dependent parameters
  • 5.2 Summary of the equation solver
  • 5.3 Summary of the nonlinear equation solver
    • 5.3.1 Nonlinear solution strategy 1 (gives values from last switch)
    • 5.3.2 Nonlinear solution strategy 2 (gives values from last switch)
    
    
  • 5.4 Visualization of the velocity field
  
  
• 6 Energy equation - simulation #2
  
  • 6.1 Other parameters
  • 6.2 Summary of the equation solver
  • 6.3 Summary of the nonlinear equation solver
    • 6.3.1 Nonlinear solution strategy 1 (gives values from last switch)
    • 6.3.2 Nonlinear solution strategy 2 (gives values from last switch)
    
    
  • 6.4 Visualization of the temperature field
  
  
• 7 Momentum and energy equation administration - simulation #3
  
  • 7.1 Menu items that have been varied in the multiple loop
  • 7.2 Parameters in constitutive relations
  • 7.3 Grid and iteration parameters
  
  
• 8 Momentum equation - simulation #3
  
  • 8.1 Problem dependent parameters
  • 8.2 Summary of the equation solver
  • 8.3 Summary of the nonlinear equation solver
    • 8.3.1 Nonlinear solution strategy 1 (gives values from last switch)
    • 8.3.2 Nonlinear solution strategy 2 (gives values from last switch)
    
    
  • 8.4 Visualization of the velocity field
  
  
• 9 Energy equation - simulation #3
  
  • 9.1 Other parameters
  • 9.2 Summary of the equation solver
  • 9.3 Summary of the nonlinear equation solver
    • 9.3.1 Nonlinear solution strategy 1 (gives values from last switch)
    • 9.3.2 Nonlinear solution strategy 2 (gives values from last switch)
    
    
  • 9.4 Visualization of the temperature field
  
  
• 10 Momentum and energy equation administration - simulation #4
  
  • 10.1 Menu items that have been varied in the multiple loop
  • 10.2 Parameters in constitutive relations
  • 10.3 Grid and iteration parameters
  
  
• 11 Momentum equation - simulation #4
  
  • 11.1 Problem dependent parameters
  • 11.2 Summary of the equation solver
  • 11.3 Summary of the nonlinear equation solver
    • 11.3.1 Nonlinear solution strategy 1 (gives values from last switch)
    • 11.3.2 Nonlinear solution strategy 2 (gives values from last switch)
    
    
  • 11.4 Visualization of the velocity field
  
  
• 12 Energy equation - simulation #4
  
  • 12.1 Other parameters
  • 12.2 Summary of the equation solver
  • 12.3 Summary of the nonlinear equation solver
    • 12.3.1 Nonlinear solution strategy 1 (gives values from last switch)
    • 12.3.2 Nonlinear solution strategy 2 (gives values from last switch)
    
    
  • 12.4 Visualization of the temperature field
  
  
• 13 User annotations
  

*

1 Momentum and energy equation administration - simulation #1

1.1 Menu items that have been varied in the multiple loop

• power law exponent : 1
• alpha : 0

1.2 Parameters in constitutive relations

• Reference viscosity (mu_0) : 1.000000
• Reference temperature (T_0) : 0.000000
• Coefficient for temperature in viscosity (alpha) : 0.000000
• Power law exponent (n in viscosity expression) : 1.000000
• Initial guess used : Boundary_layer

1.3 Grid and iteration parameters

• input file containing grid : P=PreproBox|d=2 [0,1]x[0,1]|d=2 elm_tp=ElmTensorProd2 div=[4,4], grading=[1,1]
• casename : SIMULATION_m01
• plot directory : SIMULATION-report_figures/sec1
• no of nodes : 25
• no of elements : 4
• element type : ElmTensorProd2
• No of Gauss-Seidel iterations between momentum and energy equation: 2
• Number of continuation iterations : 0
• Accumulated number of Gauss-Seidel internal iterations : 6 ( 3 3 )
• Max outer Gauss-Seidel iterations : 10
• Tolerance in outer Gauss-Seidel iteration: 0.000010
• Outer Gauss-Seidel loop converged!
• Difference measure in Gauss-Seidel iteration : 0.000000

2 Momentum equation - simulation #1

2.1 Problem dependent parameters

• Pressure gradient : 1.000000

2.2 Summary of the equation solver

2.3 Summary of the nonlinear equation solver

• Strategy of solution : 1
• Number of linearizing strategies : 2
• Number of switches : 1
• Maximum allowed nonlinear difference : 0.000100

2.3.1 Nonlinear solution strategy 1 (gives values from last switch)

• nonlinear equation solver : NewtonRaphson
• eps : 3.7622e-17 (max: 1e-06)
• no of iterations : 1 (max: 10)
• type of stopping criterion : 1

2.3.2 Nonlinear solution strategy 2 (gives values from last switch)

• nonlinear equation solver : SuccessiveSubst
• eps : 1e+20 (max: 1e-06)
• no of iterations : 0 (max: 10)
• type of stopping criterion : 1

2.4 Visualization of the velocity field

• Plot of the "velocity" field:
  • simres2mtv -f SIMULATION_m01 -s -n 'velocity' -b -o '%comment="2D n=25 ElmTensorProd2 G " contstyle=2 nsteps=40' -G 0
  • plotmtvps SIMULATION-report_figures/sec2/velocity.ps -plotall -colorps SIMULATION-report_figures/sec2/velocity.mtv
  • Further experimentation: plotmtv [options] SIMULATION-report_figures/sec2/velocity.mtv

PostScript figure
The velocity field in run no. 1

Interactive plotmtv display of the figure above.

• Curve plot made by the following Unix command:
  • curveplot gnuplot -f .SIMULATION_m01.curve.map -r 'Cross' 's1' '.' -ps SIMULATION-report_figures/sec2/s1.ps

PostScript figure
Plot of the field along line 1, run no. 1

L_2 norm : 1.63254e-01 *

3 Energy equation - simulation #1

• casename : SIMULATION_m01
• plot directory : SIMULATION-report_figures/sec3
• no of nodes : 25
• no of elements : 4
• element type : ElmTensorProd2

3.1 Other parameters

3.2 Summary of the equation solver

3.3 Summary of the nonlinear equation solver

• Strategy of solution : 1
• Number of linearizing strategies : 2
• Number of switches : 1
• Maximum allowed nonlinear difference : 0.000100

3.3.1 Nonlinear solution strategy 1 (gives values from last switch)

• nonlinear equation solver : NewtonRaphson
• eps : 9.58909e-17 (max: 1e-06)
• no of iterations : 1 (max: 10)
• type of stopping criterion : 1

3.3.2 Nonlinear solution strategy 2 (gives values from last switch)

• nonlinear equation solver : SuccessiveSubst
• eps : 1e+20 (max: 1e-06)
• no of iterations : 0 (max: 10)
• type of stopping criterion : 1

3.4 Visualization of the temperature field

• Plot of the "temperature" field:
  • simres2mtv -f SIMULATION_m01 -s -n 'temperature' -b -o '%comment="2D n=25 ElmTensorProd2 G " contstyle=2 nsteps=40' -G 0
  • plotmtvps SIMULATION-report_figures/sec3/temperature.ps -plotall -colorps SIMULATION-report_figures/sec3/temperature.mtv
  • Further experimentation: plotmtv [options] SIMULATION-report_figures/sec3/temperature.mtv

PostScript figure
The temperature field in run no. 1

Interactive plotmtv display of the figure above.

*

4 Momentum and energy equation administration - simulation #2

4.1 Menu items that have been varied in the multiple loop

• power law exponent : 0.5
• alpha : 0

4.2 Parameters in constitutive relations

• Reference viscosity (mu_0) : 1.000000
• Reference temperature (T_0) : 0.000000
• Coefficient for temperature in viscosity (alpha) : 0.000000
• Power law exponent (n in viscosity expression) : 0.500000
• Initial guess used : Boundary_layer

4.3 Grid and iteration parameters

• input file containing grid : P=PreproBox|d=2 [0,1]x[0,1]|d=2 elm_tp=ElmTensorProd2 div=[4,4], grading=[1,1]
• casename : SIMULATION_m02
• plot directory : SIMULATION-report_figures/sec4
• no of nodes : 25
• no of elements : 4
• element type : ElmTensorProd2
• No of Gauss-Seidel iterations between momentum and energy equation: 2
• Number of continuation iterations : 0
• Accumulated number of Gauss-Seidel internal iterations : 25 ( 22 3 )
• Max outer Gauss-Seidel iterations : 10
• Tolerance in outer Gauss-Seidel iteration: 0.000010
• Outer Gauss-Seidel loop converged!
• Difference measure in Gauss-Seidel iteration : 0.000000

5 Momentum equation - simulation #2

5.1 Problem dependent parameters

• Pressure gradient : 1.000000

5.2 Summary of the equation solver

5.3 Summary of the nonlinear equation solver

• Strategy of solution : 1
• Number of linearizing strategies : 2
• Number of switches : 1
• Maximum allowed nonlinear difference : 0.000100

5.3.1 Nonlinear solution strategy 1 (gives values from last switch)

• nonlinear equation solver : NewtonRaphson
• eps : 5.05932e-07 (max: 1e-06)
• no of iterations : 1 (max: 10)
• type of stopping criterion : 1

5.3.2 Nonlinear solution strategy 2 (gives values from last switch)

• nonlinear equation solver : SuccessiveSubst
• eps : 5.07302e-07 (max: 1e-06)
• no of iterations : 10 (max: 10)
• type of stopping criterion : 1

5.4 Visualization of the velocity field

• Plot of the "velocity" field:
  • simres2mtv -f SIMULATION_m02 -s -n 'velocity' -b -o '%comment="2D n=25 ElmTensorProd2 G " contstyle=2 nsteps=40' -G 0
  • plotmtvps SIMULATION-report_figures/sec5/velocity.ps -plotall -colorps SIMULATION-report_figures/sec5/velocity.mtv
  • Further experimentation: plotmtv [options] SIMULATION-report_figures/sec5/velocity.mtv

PostScript figure
The velocity field in run no. 2

Interactive plotmtv display of the figure above.

• Curve plot made by the following Unix command:
  • curveplot gnuplot -f .SIMULATION_m02.curve.map -r 'Cross' 's1' '.' -ps SIMULATION-report_figures/sec5/s1.ps

PostScript figure
Plot of the field along line 1, run no. 2

L_2 norm : 3.34052e-02 *

6 Energy equation - simulation #2

• casename : SIMULATION_m02
• plot directory : SIMULATION-report_figures/sec6
• no of nodes : 25
• no of elements : 4
• element type : ElmTensorProd2

6.1 Other parameters

6.2 Summary of the equation solver

6.3 Summary of the nonlinear equation solver

• Strategy of solution : 1
• Number of linearizing strategies : 2
• Number of switches : 1
• Maximum allowed nonlinear difference : 0.000100

6.3.1 Nonlinear solution strategy 1 (gives values from last switch)

• nonlinear equation solver : NewtonRaphson
• eps : 5.0969e-07 (max: 1e-06)
• no of iterations : 1 (max: 10)
• type of stopping criterion : 1

6.3.2 Nonlinear solution strategy 2 (gives values from last switch)

• nonlinear equation solver : SuccessiveSubst
• eps : 1e+20 (max: 1e-06)
• no of iterations : 0 (max: 10)
• type of stopping criterion : 1

6.4 Visualization of the temperature field

• Plot of the "temperature" field:
  • simres2mtv -f SIMULATION_m02 -s -n 'temperature' -b -o '%comment="2D n=25 ElmTensorProd2 G " contstyle=2 nsteps=40' -G 0
  • plotmtvps SIMULATION-report_figures/sec6/temperature.ps -plotall -colorps SIMULATION-report_figures/sec6/temperature.mtv
  • Further experimentation: plotmtv [options] SIMULATION-report_figures/sec6/temperature.mtv

PostScript figure
The temperature field in run no. 2

Interactive plotmtv display of the figure above.

*

7 Momentum and energy equation administration - simulation #3

7.1 Menu items that have been varied in the multiple loop

• power law exponent : 1
• alpha : 1

7.2 Parameters in constitutive relations

• Reference viscosity (mu_0) : 1.000000
• Reference temperature (T_0) : 0.000000
• Coefficient for temperature in viscosity (alpha) : 1.000000
• Power law exponent (n in viscosity expression) : 1.000000
• Initial guess used : Boundary_layer

7.3 Grid and iteration parameters

• input file containing grid : P=PreproBox|d=2 [0,1]x[0,1]|d=2 elm_tp=ElmTensorProd2 div=[4,4], grading=[1,1]
• casename : SIMULATION_m03
• plot directory : SIMULATION-report_figures/sec7
• no of nodes : 25
• no of elements : 4
• element type : ElmTensorProd2
• No of Gauss-Seidel iterations between momentum and energy equation: 4
• Number of continuation iterations : 0
• Accumulated number of Gauss-Seidel internal iterations : 18 ( 8 10 )
• Max outer Gauss-Seidel iterations : 10
• Tolerance in outer Gauss-Seidel iteration: 0.000010
• Outer Gauss-Seidel loop converged!
• Difference measure in Gauss-Seidel iteration : 0.000000

8 Momentum equation - simulation #3

8.1 Problem dependent parameters

• Pressure gradient : 1.000000

8.2 Summary of the equation solver

8.3 Summary of the nonlinear equation solver

• Strategy of solution : 1
• Number of linearizing strategies : 2
• Number of switches : 1
• Maximum allowed nonlinear difference : 0.000100

8.3.1 Nonlinear solution strategy 1 (gives values from last switch)

• nonlinear equation solver : NewtonRaphson
• eps : 6.79459e-17 (max: 1e-06)
• no of iterations : 2 (max: 10)
• type of stopping criterion : 1

8.3.2 Nonlinear solution strategy 2 (gives values from last switch)

• nonlinear equation solver : SuccessiveSubst
• eps : 1e+20 (max: 1e-06)
• no of iterations : 0 (max: 10)
• type of stopping criterion : 1

8.4 Visualization of the velocity field

• Plot of the "velocity" field:
  • simres2mtv -f SIMULATION_m03 -s -n 'velocity' -b -o '%comment="2D n=25 ElmTensorProd2 G " contstyle=2 nsteps=40' -G 0
  • plotmtvps SIMULATION-report_figures/sec8/velocity.ps -plotall -colorps SIMULATION-report_figures/sec8/velocity.mtv
  • Further experimentation: plotmtv [options] SIMULATION-report_figures/sec8/velocity.mtv

PostScript figure
The velocity field in run no. 3

Interactive plotmtv display of the figure above.

• Curve plot made by the following Unix command:
  • curveplot gnuplot -f .SIMULATION_m03.curve.map -r 'Cross' 's1' '.' -ps SIMULATION-report_figures/sec8/s1.ps

PostScript figure
Plot of the field along line 1, run no. 3

L_2 norm : 1.63379e-01 *

9 Energy equation - simulation #3

• casename : SIMULATION_m03
• plot directory : SIMULATION-report_figures/sec9
• no of nodes : 25
• no of elements : 4
• element type : ElmTensorProd2

9.1 Other parameters

9.2 Summary of the equation solver

9.3 Summary of the nonlinear equation solver

• Strategy of solution : 1
• Number of linearizing strategies : 2
• Number of switches : 1
• Maximum allowed nonlinear difference : 0.000100

9.3.1 Nonlinear solution strategy 1 (gives values from last switch)

• nonlinear equation solver : NewtonRaphson
• eps : 9.53561e-17 (max: 1e-06)
• no of iterations : 2 (max: 10)
• type of stopping criterion : 1

9.3.2 Nonlinear solution strategy 2 (gives values from last switch)

• nonlinear equation solver : SuccessiveSubst
• eps : 1e+20 (max: 1e-06)
• no of iterations : 0 (max: 10)
• type of stopping criterion : 1

9.4 Visualization of the temperature field

• Plot of the "temperature" field:
  • simres2mtv -f SIMULATION_m03 -s -n 'temperature' -b -o '%comment="2D n=25 ElmTensorProd2 G " contstyle=2 nsteps=40' -G 0
  • plotmtvps SIMULATION-report_figures/sec9/temperature.ps -plotall -colorps SIMULATION-report_figures/sec9/temperature.mtv
  • Further experimentation: plotmtv [options] SIMULATION-report_figures/sec9/temperature.mtv

PostScript figure
The temperature field in run no. 3

Interactive plotmtv display of the figure above.

*

10 Momentum and energy equation administration - simulation #4

10.1 Menu items that have been varied in the multiple loop

• power law exponent : 0.5
• alpha : 1

10.2 Parameters in constitutive relations

• Reference viscosity (mu_0) : 1.000000
• Reference temperature (T_0) : 0.000000
• Coefficient for temperature in viscosity (alpha) : 1.000000
• Power law exponent (n in viscosity expression) : 0.500000
• Initial guess used : Boundary_layer

10.3 Grid and iteration parameters

• input file containing grid : P=PreproBox|d=2 [0,1]x[0,1]|d=2 elm_tp=ElmTensorProd2 div=[4,4], grading=[1,1]
• casename : SIMULATION_m04
• plot directory : SIMULATION-report_figures/sec10
• no of nodes : 25
• no of elements : 4
• element type : ElmTensorProd2
• No of Gauss-Seidel iterations between momentum and energy equation: 4
• Number of continuation iterations : 0
• Accumulated number of Gauss-Seidel internal iterations : 41 ( 28 13 )
• Max outer Gauss-Seidel iterations : 10
• Tolerance in outer Gauss-Seidel iteration: 0.000010
• Outer Gauss-Seidel loop converged!
• Difference measure in Gauss-Seidel iteration : 0.000000

11 Momentum equation - simulation #4

11.1 Problem dependent parameters

• Pressure gradient : 1.000000

11.2 Summary of the equation solver

11.3 Summary of the nonlinear equation solver

• Strategy of solution : 1
• Number of linearizing strategies : 2
• Number of switches : 1
• Maximum allowed nonlinear difference : 0.000100

11.3.1 Nonlinear solution strategy 1 (gives values from last switch)

• nonlinear equation solver : NewtonRaphson
• eps : 5.48522e-07 (max: 1e-06)
• no of iterations : 1 (max: 10)
• type of stopping criterion : 1

11.3.2 Nonlinear solution strategy 2 (gives values from last switch)

• nonlinear equation solver : SuccessiveSubst
• eps : 5.40124e-07 (max: 1e-06)
• no of iterations : 9 (max: 10)
• type of stopping criterion : 1

11.4 Visualization of the velocity field

• Plot of the "velocity" field:
  • simres2mtv -f SIMULATION_m04 -s -n 'velocity' -b -o '%comment="2D n=25 ElmTensorProd2 G " contstyle=2 nsteps=40' -G 0
  • plotmtvps SIMULATION-report_figures/sec11/velocity.ps -plotall -colorps SIMULATION-report_figures/sec11/velocity.mtv
  • Further experimentation: plotmtv [options] SIMULATION-report_figures/sec11/velocity.mtv

PostScript figure
The velocity field in run no. 4

Interactive plotmtv display of the figure above.

• Curve plot made by the following Unix command:
  • curveplot gnuplot -f .SIMULATION_m04.curve.map -r 'Cross' 's1' '.' -ps SIMULATION-report_figures/sec11/s1.ps

PostScript figure
Plot of the field along line 1, run no. 4

L_2 norm : 3.34157e-02 *

12 Energy equation - simulation #4

• casename : SIMULATION_m04
• plot directory : SIMULATION-report_figures/sec12
• no of nodes : 25
• no of elements : 4
• element type : ElmTensorProd2

12.1 Other parameters

12.2 Summary of the equation solver

12.3 Summary of the nonlinear equation solver

• Strategy of solution : 1
• Number of linearizing strategies : 2
• Number of switches : 1
• Maximum allowed nonlinear difference : 0.000100

12.3.1 Nonlinear solution strategy 1 (gives values from last switch)

• nonlinear equation solver : NewtonRaphson
• eps : 8.27916e-07 (max: 1e-06)
• no of iterations : 1 (max: 10)
• type of stopping criterion : 1

12.3.2 Nonlinear solution strategy 2 (gives values from last switch)

• nonlinear equation solver : SuccessiveSubst
• eps : 1e+20 (max: 1e-06)
• no of iterations : 0 (max: 10)
• type of stopping criterion : 1

12.4 Visualization of the temperature field

• Plot of the "temperature" field:
  • simres2mtv -f SIMULATION_m04 -s -n 'temperature' -b -o '%comment="2D n=25 ElmTensorProd2 G " contstyle=2 nsteps=40' -G 0
  • plotmtvps SIMULATION-report_figures/sec12/temperature.ps -plotall -colorps SIMULATION-report_figures/sec12/temperature.mtv
  • Further experimentation: plotmtv [options] SIMULATION-report_figures/sec12/temperature.mtv

PostScript figure
The temperature field in run no. 4

Interactive plotmtv display of the figure above.

13 User annotations

1. Open the HTML file in an editor: emacs SIMULATION-report.html&
2. Search for "User annotations"
3. Replace this text by your own (conclusions, comments)
4. Add links to movies, PostScript-files etc. with additional information (below are a couple of examples)

   Plot of ... (NOTE: The link won't work unless you put a sensible filename in here...)

   Movie ... (NOTE: The link won't work unless you put a sensible filename in here...)