# the source term

 Register Blogs Members List Search Today's Posts Mark Forums Read

 April 26, 2002, 08:59 the source term #1 neil Guest   Posts: n/a In the manual of CFX4, the default source term of the transport equation is multiplied by volume fraction. So When I add other terms as source term, should I divide the volume fraction? If that is right, how about the momentum equation? Does anyone know it?

 April 26, 2002, 10:01 Re: the source term #2 cfd guy Guest   Posts: n/a Hi neil, Pay attention on units you're using. Suppose that you want to modify your V Velocity Equation. If you use a source term like this: -m*V/K with m and K your own parameters. The resulting unit must be Force/Volume -> [N/m3]. Take a look at the transport equations and you will check the units. Continuity: Mass flow rate / volume -> [Kg/(s*m3)] About linearization, I suggest you to take a look at Patankar's book (Numerical Heat Transfer and Fluid Flow, 1980). Good luck, cfd guy

 May 9, 2002, 01:55 Re: the source term #3 zhu Guest   Posts: n/a Hi, cfd guy the question puzzles me, when I want to add user scalar transport equation , the source term doesn't deal with volume fraction , i don't know whether the source term is divided by volume fraction. would you clear me out? thanks in advance. zhu

 May 9, 2002, 07:49 Re: the source term #4 cfd guy Guest   Posts: n/a Hi zhu, I don't know your specific scalar transport equation but I can say that depending the equation you have, you can't divide it by volume fraction or any parameter. In my previous post I said that care must be taken with units. Volume fraction only matters in a two-phase flow (at least), so you'll have two scalar transport equation, one for each phase. Perhaps your question is about linearization. Did you take a look at specific books? Regards, cfd guy

 May 10, 2002, 02:52 Re: the source term #5 zhu Guest   Posts: n/a hi, cfd guy, thanks i think my question isn't about linearization. i want to know whether the source term of two-phase flow is multiplied by volume fraction in program. for example, user scalar transport equation of A phase is as follow. a(r*density*scal)/a(t) + v.[r*density*U*scal-dif*v(scal)]=r*(B-P)+C*(scal of B phase-scal of A phase) r,scal, dif,B,P is volume fraction ,user scalar,diffusivity of A phase, respectively. B-gravity, t--time. v--laplacian functor. C--interphase exchange coefficient In mathematical models chapter of the cfx use guid ,the generic scalar advection-diffusion equation takes the form is given. In the equation the source term S is multiplied by volume fraction.when the source term S is added in subroutine usrsrc,whether the source terms S is automatically multiplied by volume fraction during run the program? so,in subroutine usrsrc,after linearization, whether i specify the source term S is equal to [r*(B-P)+C*(scal of B phase-scal of A phase)], or [r*(B-P)+C*(scal of A phase-scal of B phase)]/r? i don't know. looking forward to your post. zhu

 May 10, 2002, 13:41 Re: the source term #6 cfd guy Guest   Posts: n/a Hi zhu, I think I finally understood your question, sorry. It's a good question. The equation you mentioned is Eq. 3-328 from the online manual (PDF File, CFX-4.3 and 4.4 versions) and 12.10.1 of the Solver manual.1 from the printed manual (CFX-4.2 Version). If you take a look at the first term on right hand side (RHS) of the 3-328 equation you'll see r(alpha)*S(alpha), where r and S are volume fraction and source term of phase alpha, repectively. I believe, I'm not sure, that the code multiplies r(alpha) with the source terms (SP and SU), ie. if you have a source like this: phi(alpha)*r(alpha)/K(alpha), where phi and K are the scalar variable and a specific constant of alpha phase, you can put this source term in two ways: SP = -1/K(alpha) ----> (Remember that r(alpha) is multiplied later) or SU = phi(alpha)/K(alpha) ----> (You should take care using this implementation. Remember the linearization). Well, I believe that the code does like I mentioned above. I may be wrong. Commercial codes are "black boxes" and, unfortunately, you can't have full control of them. Your question also applies with the Volume. Once I had the similar doubt, however my problem was with volume of cells, not the volume fraction, when I was applying some source terms into an equation. If you're still in doubt you should contact some CFX office. I'm also interested in their answer. Good luck! cfd guy

 May 16, 2002, 21:54 Re: the source term #7 zhu Guest   Posts: n/a hi, cfd guy, thanks, I consult CFX office, they say that volume fraction isn't considered in sourse terms, in the program code, the source term s isn't multiplied by volume fraction.that is ,if a source like this: phi(alpha)*r(alpha)/K(alpha), where phi and K are the scalar variable and a specific constant of alpha phase, SU = phi(alpha))*r(alpha)/K(alpha) . p.s. i have a question about volume of cells. whether the source terms is multiplied the volume of cell. such as , convection-diffusion=s, i want add one term phi(alpha)*r(alpha)/K(alpha) in the s, whether the term multiplys volume of cell? thanks. zhu

 May 17, 2002, 09:09 Re: the source term #8 cfd guy Guest   Posts: n/a Hi zhu, I'm glad you had an answer, after all this question seems to be easy to obtain an answer, but as I said before the cfd codes are black boxes. All I can say about volume of cells is that you must satisfy the units of your transport equation. We'll keep the same example used before, but we're gonna add this source into the V velocity equation. Source Term: phi(alpha)*r(alpha)/K(alpha) Note that "phi" has dimension of L*T^-1 (velocity, m/s in SI units). r is dimensionless (volume fraction). And SUPPOSE that K has dimension of T*M^-1 (s/Kg). Thus, the source term will have dimension of Kg.m/s^2 [Newton]. This dimension must be guaranteed for the momentum equations. However SUPPOSE again that K has dimension of T*L^3*M^-1 (s.m3/Kg). In this case your source term must be multiplied by the volume of cell. SU = phi(alpha)*r(alpha)*VOL(CELL)/K(alpha) or SP = -r(alpha)*VOL(CELL)/K(alpha) Well, I believe that's it. Kind regards, cfd guy

 June 9, 2002, 01:36 Re: the source term #9 neil Guest   Posts: n/a Hi,cfd guy and Zhu. I am glad to see your discussions about the problem. However, I still have some ideas about the first question. If the volume fraction should not be divided by user, is it suitable for the momentum equations or general scalar equations or both? Thanks neil

 June 9, 2002, 01:39 Re: the source term #10 neil Guest   Posts: n/a I mean if the volume fraction should not be divided by user in both the momentum equation and user scalar equations.

 July 15, 2002, 17:10 plotting radial velocity #11 hamza Guest   Posts: n/a hi anyone i am using CFX5.5 to simulate an coaxial jet mixer. i want to plot radial velocity profile at specific axial point. i observe that i can plot radail profile at the inlet, outlet, and the defult but i can not plot at specific point for example if i have pipe x=2m long and diameter of .05m and i want plot radial profile at point x=0.5m. if you don't mind explain how to do that in cfx5.5 in detail. thank you.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post zwdi FLUENT 14 June 27, 2017 15:40 GGerber OpenFOAM Installation 54 April 24, 2015 16:02 hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24 jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51 Greg Perkins FLUENT 0 October 13, 2000 23:03

All times are GMT -4. The time now is 12:36.