The CHAM LOGO

Click here for general information about Simulation Scenarios.

Brief description of scenario:
SOLVERS

by Brian Spalding and Alexey Ginevsky, January 2013

Contents

  1. Purpose
  2. Menu items (SP)
  3. Menu items (USP)
  4. Some results (SP)
  5. Some results (USP)

1. Purpose

The Solvers SimScene contains set of tests for all the linear-equation solvers which have been introduced into Earth.

The tests concern steady heat conduction through a medium with temperature-independent conductivity and no convection. Such tasks can, in principle be colved within a single sweep. Very often however more sweeps are required as a consequence of round-off errors

. Each task can be solved using various linear equation solvers at the same conditions. Tasks use as structured as unstructured grid.

After running a task user is able to see parameters of solver as dependencies maximum residual, maximum corrections, sum of residuals, sum of corrections and internal parameters of solver vs iteration number.

User can see the dependence temperature near the cold and hot objects and in the middle of domain from sweep to sweep. Besides user is able to analyze convergence due to nett-sources of objects and sum of them, which mean the balance of heat fluxes for these tasks. Last method is more correct for physical purposes.

The full list of solvers is:

The most of solvers can utilise pre-conditioners. Solvers from package SPARSEKIT can use following pre-conditioners each:

Conjugate Residual Solver (KIVA-II) can be used with all pre-conditioners from SPARSEKIT, with solvers PHSOL and MADI as pre-conditioners and with point-by-point and Jackobi solvers as pre-conditioners.

The main purpose of creation SimScene SOLVERS is preparing the tool for testing of linear equation solvers. Well know fact is that rate of convergence depends on quality of linear equation solver. From one hand non-linear equations of hydrodynamics tasks do not require exact solution of linear system on each sweep, because linear eqautions has not exact coefficient, but from other hand, better performance of linear equation solver leads to fast convergence even for non-linear tasks.

Test case

The first task for testing of solvers is case "2D conductivity with body in xz plane". User is free for choice to use some body into domian or remain domain empty.

This is conductivity task with constant conductivity, where cold temperature is set on the west side and hot temperature is set on east side. This task has not exact solution if some body is placed in the centre of domain. Convergence is not always good, if conductivity of body is strongly different from conductivity of domain material (air). Therefore, it is good test for convergence of linear equation solvers.

As a example the view of domain with cylinder as body is demonstrated on the following figure:

Implementations of this case for structured grid and ustructured grid have some differencies. Therefore, descriptions of menu items are different. See Menu items (SP) and Menu items (USP) accordingly.