how to set up wall function in bubbleFoam

adouchihiitoko · October 17, 2010, 04:05

Hi
I try to calculate turbulent bubbly flow by bubbleFoam. I don't know how to set up wall function in k and epsilon file in 0 folder since "bubbleColumn" tutorial doesn't treat turbulent flow and wall function. I suppose that "zeroGradient" is unreasonable. Can anyone help me?

alberto · October 18, 2010, 15:51

Please, check the code. Wall-functions in bubbleFoam are coded directly in the solver and automatically applied when the turbulence model is active.

Best,

adouchihiitoko · October 24, 2010, 05:16

Thanks.
I found all files related turbulence e.g. wallfunctions.* are in \applications\solvers\multiohase\bubbleFoam.

Franz_J · November 4, 2010, 08:18

Hi Adouchihihiitoko
Did you manage to get bubbleFoam to solve with turbulence as the RASModel ? I get the following messages,

keyword laplacian(DepsilonEff,epsilon) is undefined in dictionary "/home/f/OpenFOAM/f-1.7.0/working2/sludge2D/test2/system/fvSchemes::laplacianSchemes"

file: /home/f/OpenFOAM/f-1.7.0/working2/test2/system/fvSchemes::laplacianSchemes from line 42 to line 45.

From function dictionary::lookupEntry(const word&, bool, bool) const
in file db/dictionary/dictionary.C at line 395.

l_r_mcglashan · November 4, 2010, 08:57

Just add the entry into system/fvSchemes, and be careful using the turbulence model, it's very basic! Remember it's not using the turbulence libraries that can be used for single phase flows, so the only entry it uses in the RAS dictionary is turbulence on/off.

If you look in kEpsilon.H:

Code:

    // Dissipation equation
    fvScalarMatrix epsEqn
    (
        fvm::ddt(beta, epsilon)
      + fvm::div(phib, epsilon)
      - fvm::laplacian
        (
            alphaEps*nuEffb, epsilon,
            "laplacian(DepsilonEff,epsilon)"
        )
      ==
         C1*beta*G*epsilon/k
       - fvm::Sp(C2*beta*epsilon/k, epsilon)
    );

You can see that for the laplacian discretisation the scheme is looked up in the fvSchemes dictionary.

Franz_J · November 5, 2010, 06:59

thanks Laurence
I did that and it's working. Now I'm running some test cases and it crashes after abut 50 time steps, Ur Courant number tends to diverge, but I can fiddle with it,
thanks again

October 17, 2010, 04:05	how to set up wall function in bubbleFoam	#1
adouchihiitoko New Member trister Join Date: May 2009 Posts: 3 Rep Power: 17	Hi I try to calculate turbulent bubbly flow by bubbleFoam. I don't know how to set up wall function in k and epsilon file in 0 folder since "bubbleColumn" tutorial doesn't treat turbulent flow and wall function. I suppose that "zeroGradient" is unreasonable. Can anyone help me?

October 18, 2010, 15:51		#2
alberto Senior Member Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,912 Rep Power: 36	Please, check the code. Wall-functions in bubbleFoam are coded directly in the solver and automatically applied when the turbulence model is active. Best, __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

November 4, 2010, 08:57		#5
l_r_mcglashan Senior Member Laurence R. McGlashan Join Date: Mar 2009 Posts: 370 Rep Power: 23	Just add the entry into system/fvSchemes, and be careful using the turbulence model, it's very basic! Remember it's not using the turbulence libraries that can be used for single phase flows, so the only entry it uses in the RAS dictionary is turbulence on/off. If you look in kEpsilon.H: Code: // Dissipation equation fvScalarMatrix epsEqn ( fvm::ddt(beta, epsilon) + fvm::div(phib, epsilon) - fvm::laplacian ( alphaEpsnuEffb, epsilon, "laplacian(DepsilonEff,epsilon)" ) == C1betaGepsilon/k - fvm::Sp(C2betaepsilon/k, epsilon) ); You can see that for the laplacian discretisation the scheme is looked up in the fvSchemes dictionary. __________________ Laurence R. McGlashan :: Website

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Pneumatic simulation - moving wall as a function of a pressure difference	jbmackay	OpenFOAM	0	September 22, 2010 16:51
Wall function formulation in CFX and Fluent	gravis	ANSYS	0	May 4, 2010 12:03
Compilation errors in ThirdPartymallochoard	feng_w	OpenFOAM Installation	1	January 25, 2009 07:59
Wall function problem in Fluent	mefpz	FLUENT	1	October 10, 2007 14:43
Env variable not set	gruber2	OpenFOAM Installation	5	December 30, 2005 05:27

October 24, 2010, 05:16		#3
adouchihiitoko New Member trister Join Date: May 2009 Posts: 3 Rep Power: 17	Thanks. I found all files related turbulence e.g. wallfunctions.* are in \applications\solvers\multiohase\bubbleFoam.

November 4, 2010, 08:18		#4
Franz_J New Member Franz Jacobsen Join Date: Oct 2009 Location: Brisbane Posts: 9 Rep Power: 17	Hi Adouchihihiitoko Did you manage to get bubbleFoam to solve with turbulence as the RASModel ? I get the following messages, keyword laplacian(DepsilonEff,epsilon) is undefined in dictionary "/home/f/OpenFOAM/f-1.7.0/working2/sludge2D/test2/system/fvSchemes::laplacianSchemes" file: /home/f/OpenFOAM/f-1.7.0/working2/test2/system/fvSchemes::laplacianSchemes from line 42 to line 45. From function dictionary::lookupEntry(const word&, bool, bool) const in file db/dictionary/dictionary.C at line 395.

November 5, 2010, 06:59		#6
Franz_J New Member Franz Jacobsen Join Date: Oct 2009 Location: Brisbane Posts: 9 Rep Power: 17	thanks Laurence I did that and it's working. Now I'm running some test cases and it crashes after abut 50 time steps, Ur Courant number tends to diverge, but I can fiddle with it, thanks again