|
[Sponsors] |
Using FTCS to solve a reduced form of Navier-Stokes eqn |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 30, 2015, 14:57 |
Using FTCS to solve a reduced form of Navier-Stokes eqn
|
#1 |
New Member
Justin
Join Date: Jan 2015
Posts: 28
Rep Power: 11 |
I'm relatively new to CFD and took on research in CFD. I'm a mechanical engineering major...so I'll do my best in describing the given problem.
The problem regards two parallel plates extended to infinity at a distance of h apart. The fluid between the plates has a given density and kinematic viscosity. The upper plate is stationary and the lower plate is suddenly set in motion with a constant velocity of 40 m/s. The spacing h is 4 cm. A constant streamwise pressure gradient of dp/dx is imposed within the domain at the instant motion starts. A spacial size of 0.001 m is specified. The reduced form of the Navier-Stokes equation is Using finite differences, I get the following I'm confused about the term. I'm supposed to solve for velocity using the FTCS explicit scheme with 1) , 2) Pa/m, and 3) Pa/m. So how to I account for the pressure gradient? Any help is SUPER appreciated. And no, this is not my research itself. This is something my professor gave me as a project in preparation for our research. I don't need an outright answer...just some guidance. Thanks! Last edited by DA6righthand; March 30, 2015 at 15:05. Reason: Accidentally posted before finishing post |
|
March 30, 2015, 15:46 |
|
#2 |
Member
Alex
Join Date: Jan 2014
Posts: 54
Rep Power: 12 |
Take a look at pressure-velocity coupling. 12 Steps to Navier-Stokes gives a good introduction and introduces this method, but when I looked at these lectures I felt they were seriously lacking depth.
Or...if you have a temperature boundary condition too and assume IC flow, you can use the equation of state P = pRT and solve the energy equation with momentum. This would make more sense physically, since you will lose energy to heat, but will introduce another level of complexity to your problem. You also need at least one boundary condition for your pressure. Dale Andersons book is really thorough and probably gives insight into your problem, but it isn't for beginners. |
|
March 30, 2015, 15:49 |
|
#3 |
Member
Join Date: Jul 2013
Posts: 55
Rep Power: 12 |
||
March 30, 2015, 16:21 |
|
#4 |
New Member
Justin
Join Date: Jan 2015
Posts: 28
Rep Power: 11 |
The form of the NS equation I gave in the original post is what the book gave me. Actually, I'm using the Tannehil/Anderson book to learn from. The given problem is from an older edition of the book (I have the most recent edition).
As I said, I'm relatively new to CFD. I struggled through chapter 2 of the Anderson book (I taught myself basic PDE methods like separation of variables via my differential equations book) but my professor felt I had a good enough grasp on the material to move on to chapter 3. I'm currently in chapter 3 of Anderson which I find much easier than chapter 2 because of a numerical analysis course I previously took. My first reaction to the problem in my first post was to let be constant but I'm completely unsure of that. But if I use the central difference, I'll end up with pressure and velocity terms mixed together which leads me to confusion when trying to use von Neumann analysis to determine the stability criteria. Also, I pretty much typed the problem from the book verbatim so everything about the problem is given above. Hopefully I didn't confuse you guys... |
|
March 30, 2015, 16:25 |
|
#5 |
Member
Alex
Join Date: Jan 2014
Posts: 54
Rep Power: 12 |
Quote:
|
|
March 30, 2015, 16:54 |
|
#6 |
Senior Member
Join Date: Jul 2009
Posts: 354
Rep Power: 18 |
If you are imposing the pressure gradient then simply treat it as a source term. The problem description certainly makes it sound like that is how you should proceed.
|
|
April 1, 2015, 18:33 |
|
#7 |
New Member
Justin
Join Date: Jan 2015
Posts: 28
Rep Power: 11 |
So here's where I'm at. From what has been discussed thus far and from the problem statement, we can say pressure varies linearly in the x-direction and hence is a constant. Using the FTCS scheme yields
And then solving for the time is Now I need to determine the stability criteria via von Neumann analysis. Because we only have an equation that varies in time and with height , my first reaction is to let . I know that after I substitute the exponentials into the discretized PDE, the amplification factor and that . I'm at a loss of what to do with the term though after trying to find . I end up with complex exponentials which depend on y and t in the equation for which doesn't seem correct. |
|
April 2, 2015, 03:06 |
|
#8 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71 |
the last term is just a number...
a suggestion, the factor dt*mu/dy^2 = alpha, therefore the Von Neumann analysis should be performed checking the maximum of |G| for a range of values for alpha, dt and wavenumber. |
|
April 2, 2015, 04:14 |
|
#9 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71 |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Filtered navier stokes equation..LES:: | Palani Velladurai | Main CFD Forum | 7 | September 6, 2013 02:51 |
Newbie:Viscoelasticity and Navier stokes equation | Rajil Saraswat | Main CFD Forum | 2 | June 9, 2003 07:21 |
1-d navier stokes equation | Robert | Main CFD Forum | 0 | December 20, 2002 02:14 |
help: I am trying to solve Navier Stokes compressible and viscid flow | Jose Choy | Main CFD Forum | 2 | May 18, 2000 05:45 |
any better way to solve 3D stokes flow? | Yangang Bao | Main CFD Forum | 0 | September 29, 1999 10:26 |