# BC for the stream function in the computing domain

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

 July 26, 2000, 13:21 BC for the stream function in the computing domain #1 Pierre Forges Guest   Posts: n/a July 26, 2000 Dear colleagues, ====> NOTATIONS (x,y) Cartesian coordinates. (r,s) Curvilinear coordinates. P Stream Function P is defined as follow: u = dP/dy = dP/dr dr/dy + dP/ds ds/dy [1] v =-dP/dx =-dP/dr dr/dx - dP/ds ds/dx [2] ====> INLET BOUNDARY CONDITION For x=0 and y in [0,ymax] u=1 v=0 Now if I go to the computing domain, I have: 1=dP/dr dr/dy + dP/ds ds/dy 0=dP/dr dr/dx + dP/ds ds/dx Two equations with two unknowns: dP/dr and dP/ds Solving for dP/ds and integrating from s=0 to s>0 yields: P(r=0,s)= P(0,0) + Int_{0}^{y(s)} [dr/dx] /[dr/dx ds/dy - dr/dy ds/dx] dy I have used my 2 BC, namely u=1 and v=0, to get a Dirichlet BC for the stream function P. I think it's the way to do it. Let's see now another boundary condition. ====> OUTLET BOUNDARY CONDITION For x=xmax and y in [0,ymax] du/dx=0 dv/dx=0 In the computing domain I have: du/dr dr/dx + du/ds ds/dx=0 dv/dr dr/dx + dv/ds ds/dx=0 Using eqs [1] and [2] yields: [d^2P/dr^2 dr/dy + d^2P/drds ds/dy] dr/dx + [d^2P/drds dr/dy + d^2P/ds^2 ds/dy] ds/dx=0 [3] [d^2P/dr^2 dr/dx + d^2P/drds ds/dx] dr/dx + [d^2P/drds dr/dx + d^2P/ds^2 ds/dx] ds/dx=0 [4] From these two equations, how can I get my boundary condition for the stream function P? I cannot isolate dP/ds from eqs [3] and [4]! If it was possible I would have been able to integrate from s=0 to s>0 as I did before. Same remark applies when I have the following BC: du/dy=0 and v=0. It's a BC for a symmetry line. Thank you so much in advance for replying to my question. My very best regards, Dr. Pierre Forges UAE University Mech. Eng. Dept. pforges@uaeu.ac.ae

 July 26, 2000, 14:47 Re: BC for the stream function in the computing do #2 John C. Chien Guest   Posts: n/a (1). The only reason of using the coordinate transformation is that after the transformation, the boundary condition specification will be easier. (2). For example, initially the wall is a curve, that is y=y(x). After the coordinate transformation, this becomes Eta=Eta1=a constant, using the body aligned coordinates. (3). In other words, the original wall boundary condition is Psi=a constant along y=y(x), now it becomes Psi=a constant along Eta=Eta1=a constant. (3). The same is true for the inlet condition. If u=1 along yinlet=yinlet(x), then after transformation, it becomes u=1 along Xi=0. And if you can set v_eta=0, then the flow will be uniform at Xi=0 in Eta direction. (4). So, the principle of using the coordinate transformation is to make the specification of the boundary conditions easier in the transformed coordinates. (5). In terms of the simple English, the inlet boundary condition is: at the inlet plane, velocity is uniform and is along the streamwise direction. (velocity has only Xi component) This inlet plane could be positioned originally at an angle to the x-y coordinate system and was difficult to specify the inlet boundary condition. After the transformation, it becomes easier. (6). In your case, you can keep the transformation the same as the original coordinates in the inlet and the exit region!!!! And use general transformation elsewhere. In other words, there is no transformation in the inlet and the exit regions. That is in those regions, Xi=x, Eta=y. Then those coordinate transformation factors will have value of either 0 or 1. In this way, u=1 at x=0, is the same as u=1 at Xi=0 or =Xi_zero (can be non-zero value). (7). the idea is: don't make the simple problem (u=1 at x=0) complicated through the coordinate transformation. as what you have just illustrated in your message. You are making the problem more difficult, not easier to solve.

 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 gschaider OpenFOAM Installation 120 December 2, 2009 11:23 ComFlu Main CFD Forum 1 October 25, 2009 14:25 ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50 Fred Main CFD Forum 10 January 26, 1999 13:58 Willard Lee Main CFD Forum 1 January 20, 1999 12:03

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