Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Applicable domain of Stokes Flow

  1. Feb 17, 2010 #1
    I'm working on a calculation of flow through a rectangular duct and I'm assuming I'm in the Stokes' flow regime (Re<<1), but I also want to experiment on this system and I was wondering if anyone knows until what Re-number Stokes' flow is still a good approximation (and how good, i.e. in percentages error): Re<0.1, Re<0.01 ?!

    Does someone maybe have/know a scientific paper or book on this?

    Thanks in advance,

  2. jcsd
  3. Feb 17, 2010 #2
    That is one of the problems of fluid mechanics, knowing when you are outside of the equations known domain.

    Are you able to test your flow field in an experiment to back up the calculations? I am assuming that you are not, since then you probably wouldn't have asked the question.

    EDIT: In my copy of "Transport Phenomena" 2nd Edition by Bird, Stewart, and Lightfoot that this book "Microhydrodynamics: Principles and Selected Applications" that there is a thorough discussion of "creeping flow" problems.

    The book is fairly inexpensive. See link below.



    Last edited by a moderator: Apr 24, 2017
  4. Feb 17, 2010 #3

    Andy Resnick

    User Avatar
    Science Advisor
    Education Advisor

    For flow around a sphere, Stokes' law begins to fail around Re ~0.5 I don't have a decent (inexpensive) reference for this,
  5. Feb 18, 2010 #4
    I want to solve boundary layer flow and heat transfer using Finite difference method.
    Can anyone help me.
  6. Feb 18, 2010 #5
    Thanks for the answers Andy (gives me some clue about order of magnitude) and CFD (potentially interesting book!)

    This morning I found a journal article completely devoted to this topic in which criterion is proposed which matches well with a large number of experiments. Turns out that (working with porous structures like I do) Repore<0.8 for an accuracy of 99% or better, which is similar to the 0.5 Andy gave for flow around spheres (to which some porous media are very similar).

    For those who are interested:
    Comiti et al. / Chem. Eng. Sci. 55 (2000) 3057-3061
  7. Feb 18, 2010 #6
    Why is that in this thread? This thread has nothing to do with finite difference methods.

    Please create a new thread with your question.

    Also, provide enough information so that someone can actually help you. Such as, the specifics of your problem, any graphs/drawings of the system. Your question is so vague no one could possibly help you out.

  8. Feb 22, 2010 #7
    i want to solve convection boundary layer flow past a vertical plate problem using Finite difference method, that is, i have to solve momentum (NavierStokes) and energy equations. Most of the boundary layer problem is solved by similarity transform method. That is they convert PDE to ODE then they solve ODE by Runge Kutta Shooting method. I would like to solve the PDE directly.
    Momentum Eqn
    U_t+U U_x + V U_y = U_yy + Gr T
    Energy Eqn
    T_t+U T_x + V T_y = (1/Pr) T_yy
    Boundary Condns
    U=0; T=0; at X=0
    U -> 0; T -> 0; as Y -> infinity

    When i am solving the problem i could not get the correct velocity and temperature profiles as given in many papers. So i need help in this regard, particularly how to handle the infinite boundary condition in the computer code.
    Thanks in advance
  9. Feb 22, 2010 #8
    Two things:
    1. You're still posting in the wrong place, this topic is on the Applicable domain of Stokes Flow. This can severely limit the amount of response you get because some people just don't want to answer anymore if they see a question in a topic where it doesn't belong.

    2. That aside, I have two questions for you:
    a) By the looks of your equations you want to solve the boundary layer problem in 2D, how accurate do you want to have your solution? If accuracy is not of vital importance: a very good approximation (less than 1% error) is obtained when you do the similarity transform that you already mentioned and then solve the resulting ODE via an integral method (i.e. analytically but with an approximation).

    b) Are you coding the solver for the equations by hand? Because a standard CFD package such as Comsol or Fluent should be able to solve your set of equations easily and in that case the boundary conditions shouldn't be an issue at all.

    Final thing: maybe there is a moderator who can shift this thing to a new topic?!
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook