虚心求教，一个简单的管内两相流的问题,内祥）

fhcloud

我做的是这样的一个简单的两相流问题：水平放置的圆管（直径25mm，厚1mm）内的两相是氮气和微米级的小水滴（或叫做水雾），氮气和水雾的体积百分比为：0.9：0.1
氮气的参数：初速设为1m/s（可以改变）,密度1.25kg/m3,
水雾的参数：初速为0，密度是24.4kg/m3，水雾直径是5微米（实际水雾直径不止是5um，还有更小的和更大些的，我不知道如何模拟这种多直径的水雾，也想请教大家）
圆管可以认为是没摩擦的光管，长1.5m，这两相从光管的一头流进另一头流出（参数暂时就这么多，不够或有多请各位指导）。我现在想模拟一下管中的两相速度分布和管中的不同位置的水的浓度分布情况。
以下是卧参照“我整理出的多相流计算资料”中的第一篇用phoenics 3.5版做的模拟设定，结果很失败，希望各位能够指出错误。
1.首先，我用cylindrical－polar坐标画了一个图形：包括3个object，分别是圆管的壁面（blockage），一个进口（inlet），一个出口（outlet），也有3个方向：z方向为管长方向（也就是两相在管中的主流方向），y是径向，x是旋转角的方向（可能表达不对）。
2.Within VR Editor.
Click on 'Main menu', 'Output', and toggle Monitor display from ASCII to GRAPHICS. Set the Monitor update frequency to 1. Click 'Top menu' and 'OK'.
Activate solution of the second phase.
Click on 'Main menu'.
Click on 'Models'.
Click on 'One_phase', and from the list of multi-phase models select 'IPSA Full'. Click 'OK'.

3.Set interphase linkages
Click on 'Settings' next to IPSA_FULL.
Click on Interphase Drag law 'User', and from the list of Drag laws select 'Standard drag curve'.
Click 'OK'.
Set the particle diameter to 5.0E-6 m, and the minimum slip velocity to 1.0E-3 m/s.
Click 'Previous Panel'.
4.Set properties of the second phase.
Click on 'Properties'.
Click on 'Set Phase 2 properties'.
Set the density of phase 2 to 24.4 kg/m3.
5.Set initial values
Click on 'Initialisation'
Set initial value of W1to 13 ,W2  to 0
Set initial value of R1 (phase 1 volume fraction) to 0.9, and that of R2 (phase 2 volume fraction) to 0.1. You will have to use the < > buttons next to Variable to scroll the list of variables.
6.Set relaxation parameters
Click on 'Numerics'.
Set the Gloabal convergence criterion to 0.01%.
Click on 'Relaxation control', and set the DTFALS values for U1, U2，W1 ,W2 , V1andW2 all to 0.001.
Click 'Previous panel', 'Top Menu', 'OK'.

7.Set the inlet values for the two phases
Click twice on the inlet object (called INLET). Click on 'Attributes'.
Set the inlet density to 1.125. (This is actually the phase 1 density, 1.25, multiplied by the phase 1 volume fraction, 0.9).
Click on Setting values for 'Phase 1' to switch to Phase 2. Set the inlet density for phase 2 to 2.44 (this is the phase 2 density, 24.4, multiplied by the phase 2 volume fraction, 0.1)
Set the inlet velocity for phase 2 to 0m/s in the Z direction.
Click on 'OK', 'OK'.
8.Allow the second phase to leave the domain at the outlet.
Click twice on the outlet object (called OUTLET). Click on 'Attributes'.
Set the Ph 1 Coefficient to 1.0.
Set the Ph 2 Coefficient to 100.0 (this is the phase 1 coefficient multiplied by the phase 2 density).
Click on 'OK', 'OK'.

9.Running the Solver.
In the PHOENICS-VR environment, click on 'Run', 'Solver', then 'Local solver (Earth)' and 'yes'
when prompted to save data to files.

encle

结果很失败，??????
出现什么结果？？你预计会是什么结果？？
圆管的壁面（blockage）？？为什么不用plate，而用blockage

fhcloud

谢谢
我去试试
试完了把图片上传过来，你帮我看看

phoenics

你可以参考 library case p220 等
使用代数划移模型或GENTRA试试

fhcloud

上传不了东西阿