|
发表于 2002-9-2 15:11:49
|
显示全部楼层
求教!
由于word软件的原因我没有办法看到你的公式,提供下面阻力源项的添加方法,希望能给你一些帮助。
Q1. Question: - How do I set resistance terms, such as the Ergun Equation?
Answer:
Introduction
The ERGUN equation represents the pressure drop associated with flow through a packed bed, or more generally, a porous medium.
In PHOENICS, pressure drops are applied as sinks in the momentum equation. The units of the momentum equation are those of force, so the pressure drop must also be turned into a force.
The ERGUN equation can be written as:
dp/dx = a*U1 + b*U12
where a=µ/K and b is an empirical constant. Here, µ is the dynamic viscosity of the fluid flowing through the porous media and K is an empirical constant called the permeability.
The pressure drop can be turned into a force, as follows:
Fx = Area * dp
= Area * dp/dx * dx
Now Area = dy * dz;
so dx * Area = Volume
Hence:
Fx = Volume * (a*U1 + b*U12)
Implementation
The velocities used in the equation are superficial velocities, so
This can be introduced as two patches:
PATCH(ERGUN1, VOLUME, ........
COVAL(ERGUN1, U1, a, 0.0)
PATCH(ERGUN2, VOLUME, ........
COVAL(ERGUN2, U1, -b, 0.0)
Note that the -ve COefficient in patch ERGUN2 is activating a quadratic source - as stated in POLIS.
In the VR-Editor, the above can be achieved by giving an object the type User-Defined, and setting two patches with the Attributes as described above.
|
|