The uses of a weighting factor along with a time step in a single-step trapezoidal method to solve a first-order parabolic system have been systematically studied. The weighting factors are used in two main types: constants and variables. The most commonly used constant weighting factors can be defined by the ratio of the Fibonacci sequence. Among them, the optimal weighting factor is 0.618, resulting in a balance between the overall accuracy and efficiency. With the finite element formulation, the space and time dimensions can be discretized separately. For the time discretization only, there exists a zero-error dimensionless time step if a weighting factor is within the range of 0.5-1.0. By taking advantage of the zero-error condition, the weighting factor can be correlated with a time step. The influence of spatial dimensions is lumped into a nonzero eigenvalue of the system. Through validity tests of two benchmark linear problems, the variable weighting factor for a single-step trapezoidal method is shown to be accurate, efficient, and stable. The relevant features have been captured.
[1]
Suhas V. Patankar,et al.
A NEW FINITE-DIFFERENCE SCHEME FOR PARABOLIC DIFFERENTIAL EQUATIONS
,
1978
.
[2]
W. L. Wood.
Practical Time-Stepping Schemes
,
1990
.
[3]
Optimal Weighting Factor for Single-Step Trapezoidal Method
,
2002
.
[4]
Some New Finite-Difference Schemes For Parabolic Differential Equations
,
1982
.
[5]
OPTIMAL EXPONENTIAL DIFFERENCE SCHEME FOR SOLVING PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS
,
1988
.
[6]
Kumar K. Tamma,et al.
THE TIME DIMENSION AND A UNIFIED MATHEMATICAL FRAMEWORK FOR FIRST-ORDER PARABOLIC SYSTEMS
,
2002
.
[7]
C. W. Gear.
NUMERICAL INTEGRATION OF STIFF ORDINARY DIFFERENTIAL EQUATIONS. Report No. 221.
,
1967
.
[8]
J. Z. Zhu,et al.
The finite element method
,
1977
.
[9]
T. Belytschko,et al.
Computational Methods for Transient Analysis
,
1985
.
[10]
J. Crank,et al.
A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type
,
1947
.