multiphaseInterFoam: timestep error by simulating a co-extrusion nozzle
1 Attachment(s)
Hey FOAMers,
I am a new OpenFOAM 2.1.1. user and have a problem which I try to solve since 10 days. Here my case: I want to simulate a coextrusion nozzle and a Jet-Breakup at the nozzle outlet with the multiphaseInterFoam solver but the timesteps decrease and the Simulation stops after 30 to 80 s (clocktime). Here a cut of the logfile: Code:
Courant Number mean: 7.84266e-07 max: 0.107767 A Image of teh geometry is attached: Inlet1 is Oil and Inlet2 Water, the Internalfield is filled with air. Tried several kinds of Boundary conditions for U (velocity), p_rgh (totalpressure, fixedValue and zeroGradient for inlets and outlet) and alphas (fixedvalue and zeroGradient for inlets) but none of them worked, the current 0 folder looks as follows: U: Code:
dimensions [0 1 -1 0 0 0 0]; Code:
dimensions [1 -1 -2 0 0 0 0]; Code:
dimensions [0 0 0 0 0 0 0]; alphawater: Code:
dimensions [0 0 0 0 0 0 0]; alphaoil: Code:
dimensions [0 0 0 0 0 0 0]; How mentioned I use the multiphaseInterFoam solver of the dambreak tutorial, without any changes. In the fvSolution file i played with the cAlpha ( 0 to 2) without any result. fvSolution: Code:
solvers fvSchemes: Code:
FoamFile Controldict: Code:
application multiphaseInterFoam; Maybe someone can help me to find the problem. Thank you in advance, Leon |
problem solved. Need to put a small field of the phases directly at the inlets and everything works fine.
|
hi, did you mean need a setfieldsDict to initialize the internalfield? like:
... regions ( boxToCell { box ( 0 0 -1 ) ( 0.1461 0.292 1 ); fieldValues ( volScalarFieldValue alphawater 1 volScalarFieldValue alphaoil 0 ... |
Hi, Quatschinsky:
From tutorial of the solver like dambreak4phase, I found the file "alphaair" write like this: boundaryField { leftWall { type alphaContactAngle; thetaProperties ( ( water air ) 90 0 0 0 ( oil air ) 90 0 0 0 ( mercury air ) 90 0 0 0 ( water oil ) 90 0 0 0 ( water mercury ) 90 0 0 0 ( oil mercury ) 90 0 0 0 ); value uniform 0; } .... but yours did not like this, can u tell me whether there are some default setting in it? |
Hey salame
Yeah, exactly i mean setting Fields with the setFieldsdict. If you just set the Value of the phase at the Inlet equal 1 in the alpha files it was not working but with running the setFieldsdict it is. Just the first 5 to 10 cells at the inlet need to be field with the phase also. You just need alphacontactangle if the angle between the wall and the phases is not equal 90°. If it is equal 90° you can also take zeroGradient. Regards, Leon |
Quote:
by the way, i think a extend question. If the two fluids are injected one by one, how set the setfieldsDict? is it ok that i reset the setfieldsDict file after the first fluid finish injection? then how can i determine the revelant "volScalarFieldValue"? |
buiding problem
Hi all,
I tried to make multiphaseInterFoam solver in my user directory but got the error as follows : SOURCE=mymultiphaseInterFoam.C ; g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-100 -I../myinterFoam -ImultiphaseMixture/lnInclude -I/opt/openfoam220/src/transportModels -I/opt/openfoam220/src/transportModels/incompressible/lnInclude -I/opt/openfoam220/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam220/src/turbulenceModels/incompressible/turbulenceModel -I/opt/openfoam220/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam220/src/OpenFOAM/lnInclude -I/opt/openfoam220/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linuxGccDPOpt/mymultiphaseInterFoam.o mymultiphaseInterFoam.C:36:33: fatal error: mymultiphaseMixture.H: No such file or directory compilation terminated. make: *** [Make/linuxGccDPOpt/mymultiphaseInterFoam.o] Error 1 wmake libso in the multiphaseMixture inside the directory worked fine. Can anyone please let me know where the error is ? (Hopefully it kind of belongs to the same thread) Thanks |
Quote:
I met the problem quite same with you. I can not understand " put a small field of the phases directly at the inlets". Can you give me a more detailed explanation. Thank you in advance!:) |
All times are GMT -4. The time now is 18:26. |