next up previous contents
Next: Numerical examples: linear dynamic Up: Time Discretized Operators: Type Previous: Explicit representations of time

On the design of an L-stable second-order accurate explicit time discretized operator representation

In this section, the design of a stabilized second-order accurate explicit representation of the form tex2html_wrap_inline33868 (in a nonlinear sense) of the time discretized operator pertaining to the Type 2 classification which inherits unconditional stable features is presented, in which it is not necessary to solve coupled system of equations at every time step except the starting time step (as preliminary efforts; we hope to disseminate more recent advances in the near future).

In fact, if we select tex2html_wrap_inline34308 , tex2html_wrap_inline34310 and m=1, tex2html_wrap_inline34314 in Algorithm 4, a new explicit time discretized operator representation is given by the following:

  Alr10300

Since the above representation of two-field d-v form may become somewhat computationally intensive, an alternate d-form is described next.

d-form:

  eqnarray10358

  eqnarray10376

  eqnarray10406

Prop10422

Rem10430

Spectral properties

The curves of spectral radius of the L-stable second-order accurate stabilized explicit time discretized operator representation with different artificial damping ratios are shown in Figure 43 (a - b). It is evident that: 1) The spectral radii are less than unity for dynamic problems with various physical damping ratio cases; 2) The spectral radius tends to zero as tex2html_wrap_inline27631 tends to infinity, i.e. it is L-stable; 3) The spectral radius is reduced as the artificial damping ratio is increased; and 4) The spectral radius closely mimics the theoretical solution very well for lower frequency.

   figure10437
Figure 43: Spectral radius of second-order accurate explicit representation

Amplitude and phase errors:

These are obtained as follows:

eqnarray10444

The comparison of dissipation and dispersion characteristics with the central difference method are given in Figures 44(a) and 44(b).

   figure10456
Figure 44: Dissipation and dispersion for second-order accurate explicit representation

Extension of the explicit time discretized representation for nonlinear dynamic systems

A general form of the second-order nonlinear dynamic system is considered as follows:

equation10467

Alr10476

Accuracy

Following Hoff and Taylor [42], the central difference method is second-order accurate for nonlinear dynamic problems, and hence the suggested method is second-order accurate for nonlinear dynamic problems.


next up previous contents
Next: Numerical examples: linear dynamic Up: Time Discretized Operators: Type Previous: Explicit representations of time

Xiangmin Zhou
Wed Sep 6 20:47:30 CDT 2000

The views and opinions expressed in this page are strictly those of the page author.
The contents of this page have not been reviewed or approved by the University of Minnesota.