Application of Program Generation Technology in Solving Heat and Flow Problems

Based on a new DIY concept for software development, an automatic program-generating technology attached on a software system called as Finite Element Program Generator （FEPG） provides a platform of developing programs, through which a scientific researcher can submit his special physico-mathematical problem to the system in a more direct and convenient way for solution. For solving flow and heat problems by using finite element method, the stabilization technologies and fraction-step methods are adopted to overcome the numerical difficul- ties caused mainly due to the dominated convection. A couple of benchmark problems are given in this paper as examples to illustrate the usage and the superiority of the automatic program generation technique, including the flow in a lid-driven cavity, the starting flow in a circular pipe, the natural convection in a square cavity, and the flow past a circular cylinder, etc. They are also shown as the verification of the algorithms.

New Splitting Methods for Convection-Dominated Diffusion Problems and Navier-Stokes Equations

We present a new splitting method for time-dependent convention-dominated diffusion problems.The original convention diffusion system is split into two sub-systems:a pure convection system and a diffusion system.At each time step,a convection problem and a diffusion problem are solved successively.A few important features of the scheme lie in the facts that the convection subproblem is solved explicitly and multistep techniques can be used to essentially enlarge the stability region so that the resulting scheme behaves like an unconditionally stable scheme;while the diffusion subproblem is always self-adjoint and coercive so that they can be solved efficiently using many existing optimal preconditioned iterative solvers.The scheme can be extended for solving the Navier-Stokes equations,where the nonlinearity is resolved by a linear explicit multistep scheme at the convection step,while only a generalized Stokes problem is needed to solve at the diffusion step and the major stiffness matrix stays invariant in the time marching process.Numerical simulations are presented to demonstrate the stability,convergence and performance of the single-step and multistep variants of the new scheme.