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

 April 18, 2002, 22:17 Re: about unsteady uds term #21 Greg Perkins Guest   Posts: n/a I think you'll have to make some modifications, since Fluent solves for the mass fraction of each species - not concentration directly. It thus makes, sense I think, to reformulate your equation in the form solved by Fluent. This is: d/dt(rho*Yi) + d/dxj(rho*Yi*uj) = diff + Si now your equation has a similar form - ie transient term, convection term, diffusion term and source term. To do the transformation, note that the concentration ci = rho*Yi/Wi where Yi = mass fraction, rho = density and Wi = species molecular weight. so in your case you'll get: d/dt(rho*Yi/Wi) + d/dxj(rho*Yi/Wi*uj) = diff + Si' where I've ignored doing the transform for the diffusion term and your source term for the absorption is Si' = F*dq/dt = F*k(q*-q)/Wi and I think these q's would be a fucntion of the concentration and hence mass fraction Yi. So lets just say for the moment that: Si' = f(Yi)/Wi If you multiply through by Wi, which is a constant you'll have the form given by Fluent. So all you need to so is to make sure you write a source term udf for the Si' term which will have the units kg/m3-s. You can do that fairly easily. As I said in my last post, this is quite easy once you've worked out what your term is by transforming it to something Fluent solves and understands. I don't think it would be advisable to solve this equation directly through a uds - and why bother anyway. Hope it helps Greg

 April 19, 2002, 00:31 Re: about unsteady uds term #22 Frank Guest   Posts: n/a Hi: Actually I have modified the equation, but I still cannot find any change. The following is my program, can you check for me: #include"udf.h" #define F 0.96 #define rho 1000 #define W 200 #define time_step 0.5 DEFINE_SOURCE(xmass_source1,cell,thread,ds,eqn1) {float k[2]={0.9,1.2}; float K[2]={5.9,5.1}; float b[2]={0.724,0.215}; real source1,phi1,phi_curr1,phi_old1,phi_balan1,time_di ff1; phi_balan1=K[0]*rho*C_YI(cell,thread,0)/W; (q*=KC) phi_old1=phi_1; phi_curr1=(1-time_step*k[0])*phi_old1+k[0]*time_step*phi_balan1; (dq/dt=k(q*-q) phi1=phi_curr1; time_diff1=(phi_curr1-phi_old1)/time_step; (dq/dt) source1=-F*time_diff1; return source1; } Where I use (q_curr-q_old)/time_step =k(q_curr-q_old)to represent dq/dt=k(q*-q) q*i=Kici. I do not know how to set q_old. At the beginning, q_old shoud be 0. Could you check this program for me and give me key suggestion?I feel so thankful for your any reply.

 April 20, 2002, 00:14 strange problem #23 Frank Guest   Posts: n/a Hi: I found why there is no change why I add source term. I print C_Yi(cell,thread,0), I found it is 0. It means that the species fraction cannot transfer to the program, so my program can not work. Could you tell what reason can cause this problem? Thank you very much

 May 2, 2002, 03:46 help me #24 Frank Guest   Posts: n/a Hi: Actually I have modified the equation, but I still cannot find any change. The following is my program, can you check for me: #include"udf.h" #define F 0.96 #define rho 1000 #define W 200 #define time_step 0.5 DEFINE_SOURCE(xmass_source1,cell,thread,ds,eqn1) {float k[2]={0.9,1.2}; float K[2]={5.9,5.1}; float b[2]={0.724,0.215}; real source1,phi1,phi_curr1,phi_old1,phi_balan1,time_di ff1; phi_balan1=K[0]*rho*C_YI(cell,thread,0)/W; (q*=KC) phi_old1=phi_1; phi_curr1=(1-time_step*k[0])*phi_old1+k[0]*time_step*phi_balan1; (dq/dt=k(q*-q) phi1=phi_curr1; time_diff1=(phi_curr1-phi_old1)/time_step; (dq/dt) source1=-F*time_diff1; return source1; } for component2, it is similar to the above program Where I use (q_curr-q_old)/time_step =k(q_curr-q_old)to represent dq/dt=k(q*-q) q*i=Kici. I do not know how to set q_old. At the beginning, q_old shoud be 0. Could you check this program for me and give me key suggestion?I feel so thankful for your any reply.

 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 wind Fluent UDF and Scheme Programming 5 June 21, 2013 05:39 Joseph CFX 14 April 20, 2010 15:45 kel85uk FLUENT 0 March 17, 2010 09:53 Ameya Durve FLUENT 0 February 24, 2009 14:52 Ryan, Lee FLUENT 4 October 18, 2004 10:20

All times are GMT -4. The time now is 11:03.