kiva的入口条件设定

zhzw

对于kiva，好就好在它提供了开发的源代码，只是程序太过庞大，你无法全部读通，而有时程序本身又没有提供你所需求的功能，需自行开发，难。
本人现在想处理多个不同入口条件的问题，不知如何处理，切向速度如何给定？请教各位大侠。

计算流体力学

你是否有pc机上调通的程序，能否给发一个过来，谢谢！！！， 我也正在调是该程序， 已读懂了绝大部分，但有一些算法不明白，或你能否给提供下面的文章
C.W.Hirt, A.A.Amsden, and J.L.cook. An Arbitrary Lagrangian-Eulerian commputing method for all flow speeds. J Comput Phys, 14.227(1974)

merryfish

你有源程序吗？

zhzw

[这个贴子最后由zhzw在 2003/05/19 09:58pm 编辑]

你是否有pc机上调通的程序，能否给发一个过来，谢谢！！！， 我也正在调是该程序， 已读懂了绝大部分，但有一些算法不明白，或你能否给提供下面的文章
C.W.Hirt, A.A.Amsden, and J.L.cook. An Arbitrary Lagrangian-Eulerian commputing method for all flow speeds. J Comput Phys, 14.227(1974)

这篇文章我手头就有，但是复印的，不知如何提供？文字文字

devydevyhu

您好
   我现在想做采用颗粒轨道模型的三维非平衡化学反应
   请问kiva能实现修改成这样吗？
   请问您可以将kiva发给我吗？
   或者我们可以交换也可以是有偿的
多谢

kiva

kiva改成多入口条件比较容易实现，需要改动的部分有rinput.f   setup.f  bccin.f 及 bc.f 很容易改动，至于切向速度的问题，在kiva网格生成时，系统定义了Ghost Cell（this is my view）,可以利用它来实现

luhao

Applies to:
KIVA3-V (032497)
Description:
Enables multiple inlet conditions through the use of regions.
Each region has the same inlet condition. So in order to get
multiple inlets it is necessary to define new regions.
There are six variables that has to be modified. They are:
rgamin, spdin0, velin, roin0, pcc and pcceff.
The routines involved are: setup, bcpexd, bcresp, bcpgrad, bcccfl, ccflux, setupbc, piston
and, of course, comkiva.i.
The patch_comkiva_1.1-1.2 applies on comkiva.i patched by
patch_comkiva-1.1 obtained from Detailed Chemistry
HOWTO:
Download the gzipped tarfile: patches_1.0-1.2.tar.gz with all the necessary patches and files.
COMPOSITION
apply_all_1.0-1.2
apply_all_1.1.1.2
chem.dat
iprep.mult
itape5.mult_1.1-1.2
patch_bcccfl_1.0-1.2
patch_bcpexd_1.0-1.2
patch_bcpgrad_1.0-1.2
patch_bcresp_1.0-1.2
patch_ccflux_1.0-1.2
patch_piston_1.0-1.2
patch_setup_1.0-1.2
patch_setupbc_1.0-1.2
patch_comkiva_1.0-1.2
patch_comkiva_1.1-1.2
Remember to create a backup-copy before you start.
unzip and untar the 'tar.gz'-file in the source directory.
Apply the patches...
make kiva3v and keep your fingers crossed. Run kiva3v...
Comments:
The supplied itape5 file has three major differences compared to the original itape5.
They are related with the velin, nregin0 and pcc identifiers.
  velin 0.00000e+0 0.0000e+0 0.0000e+0 0.0000e+0
nregin0    1    2    3    4
  numpcc    1
     0.0   1.00e+6   1.00e+6   1.00e+6   1.00e+6
As you probably have guessed the velocities for the different inlets are set
by velin and the inflow/initial conditions for the different regions
are set by the nregin0 identifier.
Since the inlet temperature is implicitly set by the pressure at the inlet,
which is the first value of the pcc-array. This had to be modified too.
Therefore the statement in subroutine setup
        if(bndpres.eq.1.0) then
          read ( 5,9340) id(1),numpcc
          write(12,9340) ' numpcc ',numpcc
was changed to
        if(openbnd.eq.1.0) then
          read ( 5,9340) id(1),numpcc
          write(12,9340) ' numpcc ',numpcc
So the pcc arrays are read even if there are no 'pure' pressure boundaries.
Also... The supplied itape5 file has no chemistry in it.

zhzw

多谢各位！