# How to get BC for Psi in the computing domain?

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

 July 26, 2000, 12:56 How to get BC for Psi 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 00 yields: P(r=0,s)= 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 0

 July 26, 2000, 20:10 Re: How to get BC for Psi in the computing domain? #2 John C. Chien Guest   Posts: n/a (1). The coordinate transformation is not done UNTIL the coordinate transformation factors are determined or specified. You are the one who must specify these coordinate transformation factors. (2). So, you can arbitrarily create a new coordinate transformation (just like a structured mesh) and have all non-zero transformation factors. (3). But, who would do thing like that. (4). If your original inlet condition is nice and clean in x-y coordinate system, then keep the transformed coordinates exactly the same as the x-y coordinates by specifying the proper transformation factors in that region. In simple sentense, no stretching and rotation in the inlet region when creating the new coordinates. In this way, the new is identical to the old. But you are still using the new one, except that the transformation factors in those regions are such that nothing has changed. (5). This is the way to create the new coordinates. If you do not follow this rule, you will end up with a general distorted mesh over the inlet, and will have to use interpolation method to specify the inlet conditions in the new coordinates. So, you have to generate the mesh (or the cew coordinates) to simply the work by specifying the transformation factors.(or how the mesh should be stretched or rotated). (6). In other words, in your message, you have not even finished your coordinate transformation, because the transformation factors are not specified yet. You are the one who must supply the values of ds/dx(at i,j), ds/dy(at i,j), dr/dx(at i,j), and dr/dy(at i,j). (if s=x at i=1, and r=y at i=1, then ds/dx=1 at i=1, dr/dy=1 at i=1. ds/dy=dr/dx=0. Or you can let s=const*x, r=const2*y, the ds/dy and dr/dx are still zeros.)

 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 HMR CFX 5 October 10, 2016 05:57 TWaung CFX 4 May 1, 2012 03:14 Kiat110616 CFX 4 April 3, 2011 22:43 mike CFX 1 March 19, 2008 08:22 Pierre Forges Main CFD Forum 1 July 26, 2000 14:47

All times are GMT -4. The time now is 20:52.