# contribution of convective source in total source?

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 July 27, 2008, 20:31 contribution of convective source in total source? #1 Antony Guest   Posts: n/a Hi all, I just read a Phoenics document concerning source. I can't grasp the understanding of convective source of a conserved variable (PHI) as a result of setting COVAL(Patch_Name,PHI,C,Vphi) : S_phi = T*C*(V_phi - Phi_p) + max(0,Sp)*(V_phi - Phi_p) where Sp is mass source imposed to the same patch. Suppose I already set a positive mass source Sp in a patch X. Now I want to set the total source of a conserved variable (Phi) in the patch X as S_phi = 20*(10 - Phi_p). I thought this was simply done by setting COVAL(X,Phi,20,10). But having read a Phoenics document about convective source, I am confused. Given the mass source Sp, does it mean I need to set COVAL (X,Phi,20-Sp/T,10) so that in total I will have : S_phi = T*(20- Sp/T)*(10 - phi_p) + Sp*(10 - phi_p) = T*20*(10 - phi_P) ? My understanding prior to this was that, after defining a mass source S, we also need to set the value (let's name it VS_Phi) of Phi transported by the source S, thus giving the total source : S_phi = T*20(10 - phi_p) + S*VS_Phi Now I don't understand the importance of setting this VS_phi if the total source is defined by T*C*(V_phi - Phi_p) + max(0,Sp)(V_phi - Phi_p) I also tried to work out this source with the Phoenics Mathematical Basis, but still cannot catch how this convective source T*C*(V_phi - Phi_p) + max(0,Sp)*(V_phi - Phi_p) come about in the equation. Can somebody release my deep confusion ? Regards, Antony

 July 29, 2008, 08:01 Re: contribution of convective source in total sou #2 Mike Guest   Posts: n/a The source of Phi associated with a mass source patch is given by: S_phi = [ C_phi*T + max(0,Sp)]*(V_phi - Phi_p) where Sp = C_p*T*(V_p - P_p) and T is a geometric multiplier, such as area. As an example, the PATCH settings would be: PATCH(X,WEST,1,1,1,NY,1,NZ,1,LSTEP) COVAL(X,P1,FIXFLU,mf_in) COVAL(X,U1,ONLYMS,u1_in) COVAL(X,C1,ONLYMS,c1_in) so that T=Aw (west area), C_phi=ONLYMS=0.0, and Sp = mf_in*Aw with mf_in>0.0 Thus,with V_c1=c1_in and V_u1=u1_in, we have S_C1 = mf_in*Aw*(c1_in - c1_p) and S_u1 = mf_in*Aw*(u1_in - u1_p) The coeffcient C_phi would be greater than zero if one wished to account for the diffusive inflow of material.

 July 29, 2008, 23:26 Re: contribution of convective source in total sou #3 Antony Guest   Posts: n/a Thank's a lot, Mike. A nice example. Your first sentence makes me doubt my understanding. With the source defined by : S_phi = [ C_phi*T + max(0,Sp)]*(V_phi - Phi_p) Which one is true in the following ? 1. S_phi is the total source of Phi in a cell designated by the patch 2. S_phi is a source of Phi associated with a mass source patch only (not the total source of Phi) ? If I have the following setting (with mf_in > 0) : PATCH(X,WEST,1,1,1,NY,1,NZ,1,LSTEP) COVAL(X,P1,FIXFLU,mf_in) COVAL(X,C1,ONLYMS,c1_in) PATCH(Y,CELL,1,1,1,NY,1,NZ,1,LSTEP) COVAL(Y,C1,Cov,c2_in) which one is true in the following ? 1. S_C1 = mf_in*Aw*(c1_in - c1_p) + [Cov + mf_in*Aw]*(c2_in - c1_p) 2. S_C1 = mf_in*Aw*(c1_in - c1_p) + Cov*(c2_in - c1_p) Regards, Antony

 July 29, 2008, 23:35 Re: contribution of convective source in total sou #4 Antony Guest   Posts: n/a Sorry, this is the correct post. I forgot to preview the previous message. ----------------------------------------------------------- Thank's a lot, Mike. A nice example. Your first sentence makes me doubt my understanding. With the source defined by : S_phi = [ C_phi*T + max(0,Sp)]*(V_phi - Phi_p) Which one is true in the following ? 1. S_phi is the total source of Phi in a cell designated by the patch 2. S_phi is a source of Phi associated with a mass source patch only (not the total source of Phi) ? Likewise, if I have the following setting (with mf_in > 0) : PATCH(X,WEST,1,1,1,NY,1,NZ,1,LSTEP) COVAL(X,P1,FIXFLU,mf_in) COVAL(X,C1,ONLYMS,c1_in) PATCH(Y,CELL,1,1,1,NY,1,NZ,1,LSTEP) COVAL(Y,C1,Cov,c2_in) which one is true in the following ? 1. S_C1 = mf_in*Aw*(c1_in - c1_p) + [Cov + mf_in*Aw]*(c2_in - c1_p) 2. S_C1 = mf_in*Aw*(c1_in - c1_p) + Cov*(c2_in - c1_p) Regards, Antony

 July 30, 2008, 05:15 Re: contribution of convective source in total sou #5 mike Guest   Posts: n/a S_phi is a source of Phi associated with a mass source patch only. If you have the following settings (with mf_in > 0) : PATCH(X,WEST,1,1,1,NY,1,NZ,1,LSTEP) COVAL(X,P1,FIXFLU,mf_in);COVAL(X,C1,ONLYMS,c1_in) and PATCH(Y,CELL,1,1,1,NY,1,NZ,1,LSTEP) COVAL(Y,C1,Cov,c2_in) Then PATCH X creates a "mass-flow" source for C1 givcen by: S_C1_X = mf_in*Aw*(c1_in - c1_p) and PATCH Y creates a C1 source given by: S_C1_Y = Cov*(c2_in - c1_p) Consequently, the total source of C1 is given by (2) at the bottom of your posting, ie: S_C1_TOT = S_C1_X + S_C1_Y = mf_in*Aw*(c1_in - c1_p) + Cov*(c2_in - c1_p) S_phi defines total source for phi

 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 NH CFX 1 December 16, 2008 18:24 Li CFX 2 July 28, 2008 02:33 Kumar CFX 0 May 12, 2008 09:33 Ramcharan CFX 0 February 13, 2008 08:27 Mehdi Main CFD Forum 1 February 18, 2004 00:40

All times are GMT -4. The time now is 16:46.

 Contact Us - CFD Online - Top