2dgl
所属分类:数值算法/人工智能
开发工具:Fortran
文件大小:217KB
下载次数:4
上传日期:2016-01-26 10:46:09
上 传 者:
JoeChiang
说明: 计算流体力学方法 二维Fortran程序源代码
(BOOK and Code Computational methods for fluid dynamics 2DGL)
文件列表:
2dgl (0, 2015-09-08)
2dgl\mg (0, 2015-09-25)
2dgl\mg\bound.inc (616, 1997-10-06)
2dgl\mg\bound.ing (516, 1997-10-06)
2dgl\mg\caffa.f (75762, 1998-08-31)
2dgl\mg\charac.inc (132, 1997-10-06)
2dgl\mg\circl.ing (168, 1997-10-07)
2dgl\mg\coef.inc (114, 1997-10-06)
2dgl\mg\examples (0, 2015-09-08)
2dgl\mg\examples\buocyl (0, 2015-09-25)
2dgl\mg\examples\buocyl\buocyl.cin (662, 1998-09-07)
2dgl\mg\examples\buocyl\buocyl.gin (1544, 1998-09-07)
2dgl\mg\examples\buocyl\buocyl.pin (308, 1998-09-07)
2dgl\mg\examples\cav45l (0, 2015-09-25)
2dgl\mg\examples\cav45l\cav45l.cin (773, 1997-12-23)
2dgl\mg\examples\cav45l\cav45l.gin (896, 1998-08-30)
2dgl\mg\examples\cav45l\cav45l.pin (316, 1997-10-13)
2dgl\mg\examples\cav90l (0, 2015-09-25)
2dgl\mg\examples\cav90l\cav90l.cin (758, 1999-03-30)
2dgl\mg\examples\cav90l\cav90l.gin (822, 1998-08-30)
2dgl\mg\examples\cav90l\cav90l.pin (318, 1999-03-30)
2dgl\mg\examples\chanel (0, 2015-09-25)
2dgl\mg\examples\chanel\chanel.cin (10044, 1998-09-07)
2dgl\mg\examples\chanel\chanel.gin (1146, 1997-12-24)
2dgl\mg\examples\chanel\chanel.pin (307, 1997-11-13)
2dgl\mg\examples\concil (0, 2015-09-25)
2dgl\mg\examples\concil\concil.cin (819, 1998-09-10)
2dgl\mg\examples\concil\concil.gin (684, 1998-08-30)
2dgl\mg\examples\concil\concil.pin (434, 1998-09-07)
2dgl\mg\examples\concil\concil.pin1 (506, 1998-09-07)
2dgl\mg\examples\damdif (0, 2015-09-25)
2dgl\mg\examples\damdif\damdif.cin (3102, 1998-09-07)
2dgl\mg\examples\damdif\damdif.gin (1104, 1998-09-07)
2dgl\mg\examples\damdif\damdif.pin (309, 1998-09-07)
2dgl\mg\examples\hotcol (0, 2015-09-25)
2dgl\mg\examples\hotcol\hotcol.cin (812, 1998-08-30)
2dgl\mg\examples\hotcol\hotcol.gin (820, 1998-08-30)
2dgl\mg\examples\hotcol\hotcol.pin (318, 1997-10-12)
2dgl\mg\float.inc (75, 1998-08-30)
2dgl\mg\geo.inc (109, 1997-10-07)
... ...
CAFFA-MG: Multigrid Version of the CAFFA-code
=============================================
The CAFFAC.F version of the code (see directory "sg") was extended
by implementing the multigrid acceleration for the outer iterations
(multigrid acceleration for inner iterations does not bring much;
it would only make sense for the pressure-correction equation and
would reduce the computing time by at most a factor of 2; see directory
"solvers" for a multigrid version of SIP, which could be directly
implemented in CAFFA-codes instead of single-grid SIP solver).
The method is described in Chap. 11, Sect. 11.1; a detailed description
can also be found in a paper:
Lilek, Z., Muzaferija, S., Peric, M.: Efficiency and accuracy
aspects of a full-multigrid SIMPLE algorithm for three-dimensional
flows, Numer. Heat Transfer, Part B, Vol. 31, pp. 23--42 (1997).
The program includes a lot of comments which describe all the
extensions; users familiar with the single-grid version and after
reading Sect. 11.1 will have no problems in following the steps.
In the directory "examples" there are several test cases set up for
the multigrid computations: 'cav45l' (lid-driven cavity with inclined
walls at 45 deg.), 'cav90l' (lid-driven square cavity), 'hotcol'
(square cavity with one wall hot and the other cold), 'chanel'
(complex channel with a strongly non-orthogonal grid), 'concil'
(buoyancy-driven flow between concentric cylinders, inner hot and
outer cold; only one half computed), etc. The parameters
used in these test calculations have not been optimized (number
of iterations on each grid level, under-relaxation factors etc.).
The speed-up compared with the single-grid version is remarkable.
Note that for unsteady problems one usually needs only few outer
iterations per time step, so multigrid acceleration does not lead
to a significant reduction of computing time. Steady, elliptic
problems may be speeded-up by as much as 100 times (at sixth level)!
I guess that an additional modification is needed in the code for
O- and C-grids (such test problems have not been tried so far); the
new version will replace this one after the correction is applied.
I appreciate the help provided by Martin Schmid (who found and
corrected most of the bugs I made in programming the code extensions)
and Samir Muzaferija (who helped find the most difficult ones).
M. Peric, Hamburg, 19***
peric@schiffbau.uni-hamburg.de
近期下载者:
相关文件:
收藏者: