Help parabolic inlet using groovyBC
1 Attachment(s)
Hi ,
I've been experimenting with different boundary conditions and am trying to set up velocity profile using groovyBC. At the inlet i require a power law condition given by U=2.5*(z/zr)^0.299... Code:
inlet " ^^ " I am also a little confused on how to specify my "y" coordinate varying from 40 to 100 in groovy BC, I think i've done it incorrectly. Any help would be much appreciated. Thanks in advance, Regards , Abhinav 
I'm also now trying to use groovy BC to model a parabolic inlet velocity profile in a pipe flow.
Any help please? 
pts() returns a vector and you are using it as a scalar. Try something like this instead:
Code:
inlet 
Hi Aidan,
Thanks for the prompt reply. I tried implementing "dir=vector (0,1,0)" as suggested, but then when i run i get the following... Quote:
Regards, Abhinav 
1 Attachment(s)
Hi,
Thanks for the reply. I am trying to put a parabolic velocity profile at a pipe inlet. u(r) = umean * (1 (r/R)^2) My domain is a half pipe (see attached Picture). Pipe is oriented along the xdirection. Code:
inlet It's supposed to be 0.178 m/s at the inlet centre. I think I haven't got the transformation from the cylindrical to cartesian coordinate correctly. Can someone please help on the transformation. Thank you very much. Stefano 
groovy BC
Hi foamers,
I've been trying out this new boundary condition describing inflow vertical velocity profile given as: u(z)=u*/k*log((zdo)/zo) i've used groovyBC at my inlet wall as follows... Code:
inlet Quote:
regards, Abhinav 
I could get to run my case without any error.The problem was a small mistake at my end in defining my points.(sorry)
Thank you, Abhinav 
1 Attachment(s)
Hello Stefano and Abhinav,
I am trying to set parabolic inlet velocity using the groovyBC on a radial mesh which I generated using the guidelines from this website http://openfoamwiki.net/index.php/Si...nical_diffuser A snapshot of my mesh is attached . My flow is in the Z direction and it is a radial mesh. My doubt is that to apply the parabolic inlet do I have to perform transformation. Secondly I dont get the syntax of the groovyBC as used by Stefano. Some advice would be really helpful bye regards K.Suresh kumar 
1 Attachment(s)
Hi Kumar,
I just used the groovyBC syntax as given in the groovyBC Wiki (http://openfoamwiki.net/index.php/Contrib_groovyBC). I made a bit of mistake in my earlier post. The setup for the 0/U file supposed to be: Code:
inlet Un is the centreline velocity at the inlet c is the position vector of the centre of the inlet patch. The parabolic profile is just simply given by: U(r) = Un * (1(r/R)^2), where R is the radius of the inlet patch. Just put the above expression onto your 0/U for the inlet patch, set rpm to zero, change the value of c and Un. It's already in the cylindrical coordinate. It should work straight away  no additional transformation required. Be warned, epsilon will blow up in less than 10 iterations (if you use kepsilon and steady RANS, and starting the solution from zero). You need to start the run with laminar RAS model, then switch to kepsilon after the initial instabilities near the wall dissapear. Good luck. Please let me know if you need further info. Best Regards, Stefano 
Hi Stefano,
Thanks for the clear explanation of the use of the syntax. I will try the groovyBC in the format that you have specified. Actually I am still having some issues with the mesh . Once I am done with that, I will straight away use the groovyBC in the way you have specified. I will get back to you, If I have any errors. Thanks again. bye regards K.Suresh kumar 
1 Attachment(s)
Hi Stefano,
I had to change the mesh, because I want to compare the results with another solver FS3D. This time I am using a cartesian mesh and the solver interFoam. So I first had to define the volume fraction in a circular way from the center. Now I want to put a parabolic velocity profile on it. I started with your formulation. internalField uniform (0 0 0); boundaryField { inletWall { type groovyBC; variables "rpm=0.0;Un=20;c=vector(2.1e3,2.1e3,0);n=sum(normal())/mag(sum(normal()));p=pos()c;r=mag(p)+1.0e10;R=max(r);xt=vector(n.y,n.x,0);xT=xt/mag(xt);yt=vector(n.x*n.z,n.y*n.z,n.x*n.x+n.y*n.y);yT=yt/mag(yt);"; valueExpression "(alpha1+internalField(alpha1))*Un*normal()*(1pow(r/R,2)) + (rpm*pi/30)*((p & yT)*xT  (p & xT)*yT)"; value uniform (0 0 20); timelines (); } atmosphere { type fixedValue; value uniform (0 0 0); } } But I am not getting the parabolic velocity profile.The definition of c is the location of the centre of the round inlet. And I want to have a max velocity of 20 m/s. Could you give me any suggestions on how to set a parabolic velocity profile for this kind of problem regards K.Suresh kumar 
Quote:
I don't understand why you put "(alpha1+internalField(alpha1))" in your value expression. This is a U BC isn't it? I don't think you need that alpha1. But again, I have not run interFoam before. Another suggestion, it may be simpler to put the centre of the inlet face at 0,0,0. It makes things simpler because you don't have to worry about rounding. The other thing, is the patch "atmosphere" your outlet? You cannot have both fixedValue for velocity on both inlet and outlet. If you fixed the velocity at inlet, then you need to use zeroGradient for velocity at the outlet boundary. You can then use zeroGradient pressure at Inlet, and fixedValue pressure at outlet. Kind Regards, Stefano 
Hi Stefano,
Thanks for the suggestions. Actually I solved the problem for setting the parabolic inlet. You are right about the velocity at the outlet. I will change it. Your explanation was very useful. regards K.Suresh kumar 
information in parabolic inlet
Hello foamers
I 'm using open foam (IcoFoam)to solver my problem.I want use the parabolic velocity inlet for a cylinder.I understand stephane use a same my problem. Stephane ,I want your help.I have attache the geometry and the condition inlet. Thank 
information in parabolic inlet
1 Attachment(s)
Hello Stefano
I send you my geomtry .you can find it on the attachment.I work with IcoFoam( laminar staionnary) and i want a parabolic velocity inlet. Thank you . 
GroovyBC doubt
Hello,
I am doing simulation on a rotor blade. My flow is a pressure based flow. I have my Velocity at inlet which is given by direction cosines dir_cos1 = 4.04 dir_cos2 = EXPR{57.2958*(acos(0.07045*(z/sqrt(z*z+y*y))))} dir_cos3 = EXPR{57.2958*(acos(0.07045*(y/sqrt(z*z+y*y))))} Can someone tell me how to give this using groovyBC. Thanks a lot !! 
Quote:

All times are GMT 4. The time now is 17:28. 