hiuluom |
June 7, 2017 06:37 |
Hi all,
I solved the above problems by adding phic at inlet and outlet BC for k and epsilon
k
Code:
inlet_pipe
{
type turbulentIntensityKineticEnergyInlet;
value $initialField;
phi phic;
alpha alphac;
intensity 0.03;
}
Code:
outlet
{
type inletOutlet;
phi phic;
alpha alphac;
inletValue $internalField;
value $internalField;
}
epsilon
Code:
inlet_pipe
{
type turbulentMixingLengthDissipationRateInlet;
value $initialField;
phi phic;
alpha alphac;
mixingLength 0.042;
}
Code:
outlet
{
type inletOutlet;
phi phic;
alpha alphac;
inletValue $internalField;
value $internalField;
}
But I get a new error with that:
Code:
Courant Number mean: 5.18008e-007 max: 0.0244965
Time = 2e-005
Evolving kinematicCloud
Solving 3-D cloud kinematicCloud
Cloud: kinematicCloud
Current number of parcels = 0
Current mass in system = 0
Linear momentum = (0 0 0)
|Linear momentum| = 0
Linear kinetic energy = 0
Rotational kinetic energy = 0
model1:
number of parcels added = 0
mass introduced = 0
Parcels absorbed into film = 0
New film detached parcels = 0
Parcel fate (number, mass) : patch inlet_pipe
- escape = 0, 0
- stick = 0, 0
Parcel fate (number, mass) : patch outside_wall_pipe
- escape = 0, 0
- stick = 0, 0
Parcel fate (number, mass) : patch inside_wall_pipe_1
- escape = 0, 0
- stick = 0, 0
Parcel fate (number, mass) : patch symmetry_pipe
- escape = 0, 0
- stick = 0, 0
Parcel fate (number, mass) : patch symmetry_2
- escape = 0, 0
- stick = 0, 0
Parcel fate (number, mass) : patch outside_wall
- escape = 0, 0
- stick = 0, 0
Parcel fate (number, mass) : patch outlet_ambient
- escape = 0, 0
- stick = 0, 0
Parcel fate (number, mass) : patch inside_wall_pipe_2
- escape = 0, 0
- stick = 0, 0
DILUPBiCG: Solving for Ucx, Initial residual = 1, Final residual = 1.25761e-009, No Iterations 2
DILUPBiCG: Solving for Ucy, Initial residual = 1, Final residual = 9.96601e-010, No Iterations 2
DILUPBiCG: Solving for Ucz, Initial residual = 1, Final residual = 1.65935e-006, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.00719945, No Iterations 18
GAMG: Solving for p, Initial residual = 0.157115, Final residual = 0.00130879, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0174056, Final residual = 0.000154434, No Iterations 3
time step continuity errors : sum local = 1.30207e-007, global = -1.40036e-008, cumulative = -1.40036e-008
GAMG: Solving for p, Initial residual = 0.00302915, Final residual = 1.90974e-005, No Iterations 5
GAMG: Solving for p, Initial residual = 0.00119684, Final residual = 6.68533e-006, No Iterations 3
GAMG: Solving for p, Initial residual = 0.000211614, Final residual = 7.64506e-007, No Iterations 5
time step continuity errors : sum local = 6.32305e-010, global = -7.29831e-011, cumulative = -1.40765e-008
DILUPBiCG: Solving for epsilon, Initial residual = 0.000435878, Final residual = 1.6674e-008, No Iterations 1
--> FOAM FATAL ERROR:
request for volVectorField U from objectRegistry region0 failed
available objects of type volVectorField are
3
(
Uc_0
Uc
cloudVolSUSu
)
From function objectRegistry::lookupObject<Type>(const word&) const
From a message OF, I do not know "request for volVectorField U" where BCs are or solution?
The velocity Uc as the same Goldschmidt OF tutorial. I hope foamers give me a recommendation to solve problem.
Thanks
|