D03PLF Example Program Results
Example 1
Method Parameters:
Number of mesh points used = 201
Relative tolerance used = 0.250E-03
Absolute tolerance used = 0.100E-04
Integration Results:
Global error is less than 100 times the local error tolerance.
Integration Statistics:
Number of time steps (nearest 50) = 150
Number of function evaluations (nearest 100) = 1400
Number of Jacobian evaluations (nearest 20) = 20
Number of iterations (nearest 100) = 400
Example 2
Problem Parameter and initial conditions:
gamma = 1.400
e(x<0.5,0) = 2.500 e(x>0.5,0) = 0.250
rho(x>0.5,0) = 1.000 rho(x>0.5,0) = 0.125
Method Parameters:
Number of mesh points used = 141
Relative tolerance used = 0.500E-03
Absolute tolerance used = 0.500E-02
Solution
t x density velocity pressure
0.100 0.0000 1.0000 0.0000 1.0000
0.100 1.0000 -0.0000 1.0000
0.200 1.0000 -0.0000 1.0000
0.300 1.0000 -0.0000 1.0000
0.400 0.8668 0.1665 0.8188
0.500 0.4299 0.9182 0.3071
0.600 0.2969 0.9274 0.3028
0.700 0.1250 0.0000 0.1000
0.800 0.1250 -0.0000 0.1000
0.900 0.1250 -0.0000 0.1000
1.000 0.1250 0.0000 0.1000
0.200 0.0000 1.0000 0.0000 1.0000
0.100 1.0000 -0.0000 1.0000
0.200 1.0000 -0.0000 1.0000
0.300 0.8718 0.1601 0.8253
0.400 0.6113 0.5543 0.5022
0.500 0.4245 0.9314 0.3014
0.600 0.4259 0.9277 0.3030
0.700 0.2772 0.9272 0.3031
0.800 0.2657 0.9276 0.3032
0.900 0.1250 -0.0000 0.1000
1.000 0.1250 0.0000 0.1000
Number of time steps (nearest 50) = 150
Number of function evaluations (nearest 50) = 400
Number of Jacobian evaluations (nearest 1) = 1
Number of iterations (nearest 1) = 2