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

Homework Help: Velocity in microchannel with temporal temperature variation

  1. Oct 30, 2017 #1
    a microchannel of length 2L and width h in the thermal cycling region. the temperature profile upload_2017-10-30_14-21-20.png .....(1)
    the cyclic temperature profile leads to a time dependent density upload_2017-10-30_14-22-42.png ....(2)
    using the mass conservation equation i.e. upload_2017-10-30_14-23-38.png ...(3)
    and momentum balance equation i.e. upload_2017-10-30_14-25-21.png ....(4)
    we have to find the exact solution of u, v & p

    The attempt at a solution
    equ(1) → upload_2017-10-30_14-29-10.png
    equ(2) → upload_2017-10-30_14-29-35.png
    length scales are normalised as upload_2017-10-30_14-30-20.png
    so equ(2) becomes upload_2017-10-30_14-30-55.png
    no slip & no penetration boundary condition at the walls (y=0,h) ; and constant pressure at channel entrance(x=-L) and exit(x=L). for small reynolds number eq(3) becomes upload_2017-10-30_14-37-9.png
    this equation normalised as upload_2017-10-30_14-45-34.png
    as upload_2017-10-30_14-49-1.png
    and from upload_2017-10-30_14-49-22.png
    we can get upload_2017-10-30_14-51-27.png
    & upload_2017-10-30_14-53-56.png is it the right process?

  2. jcsd
  3. Oct 31, 2017 #2
    What is the exact statement of the problem?
  4. Nov 10, 2017 #3
  5. Nov 10, 2017 #4
    I'm a bit confused over the problem statement. It says, "There are no spatial temperature gradients present along the length of the micro channel," and yet the temperature equation they provide indicates otherwise. I'm having trouble articulating what is happening physically in this problem. Are you supposed to assume that the length scale of the periodic temperature variations (lambda) is small compared to the distance 2L?

    Please articulate what you interpret is happening in the problem. Is it that the density variations are causing fluid flow pumping in the x direction?
  6. Nov 15, 2017 #5
    sorry this is the real physical problem IMG_20171115_141125.jpg IMG_20171115_141217.jpg IMG_20171115_145801.jpg IMG_20171115_145838.jpg IMG_20171026_151233.jpg
  7. Nov 15, 2017 #6
    The two fluid flow problems you have presented here are extremely different, and I would approach them entirely differently. I am assuming that you really want to solve the problem involving the moving temperature wave, indicated in the first post. I would assume that the spatial wavelength ##\lambda## is small compared to the length of the tube 2L, so that, for all intents and purposes, the tube could be assumed infinitely long. I would then change the frame of reference of the observer such that he is moving along with the temperature wave at the speed c. In that case, he would observe a density variation that is a function only of x. This then would reduce the problem to a steady state flow problem.

    I'm prepared to reveal more hints on how to implement the approach I described above involving a moving observer if what I said is not sufficiently clear.
    Last edited: Nov 16, 2017
  8. Nov 16, 2017 #7
    Sir, for the second problem is my process right? Or I have to solve in another way. More preciously if this process is wrong, how I can get the exact solution of velocities?
  9. Nov 16, 2017 #8
    As I said, in my judgment, because of the combined x and t variation, you are not going to get it solved the way you are doing it. But, the method I was explaining will work.

    From the frame of reference of an observer moving to the left at velocity c (i.e., moving coordinate system), you have $$T=T_0+\Delta T \cos{2\pi \frac{X}{\lambda}}$$ where X now replaces ##x+ct##. So now the density is a function only of X. As reckoned from this moving coordinate system, the continuity equation becomes: $$\frac{\partial (\rho U)}{\partial X}+\frac{\partial (\rho V)}{\partial y}=0$$where U and V are the velocity components reckoned from the frame of reference of the moving observer.

    Boundary conditions on the flow are now: ##U=+c## and V= 0 at y=h and y=0.

    Last edited: Nov 18, 2017
  10. Nov 19, 2017 #9
    It is too bad that @big dream has not returned. I find this a very interesting problem. Is there anyone else out there who is interested in pursuing the solution to this problem?
  11. Nov 20, 2017 #10
    I was sick. so now should I put stream function? or I have to normalise continuity equation first.
  12. Nov 20, 2017 #11
    Before doing anything else, I need you to prove to yourself that, if we make the following transformation of variables,

    t = T
    y = Y
    x = X - ct
    v = V(X,Y)
    u = U(X,Y) - c
    p = P(X)

    our equations reduce to:
    $$\frac{\partial (\rho U)}{\partial X}+\frac{\partial (\rho V)}{\partial Y}=0\tag{1}$$
    $$\frac{\partial P}{\partial X}=\mu \frac{\partial ^2U}{\partial Y^2}\tag{2}$$
    subject to the boundary conditions ##U=c## and ##V=0## at ##Y=0,h##. These are equivalent to the equations that an observer who is traveling at a constant speed c in the negative x direction would write. This observer would conclude that, as reckoned from his frame of reference, the system is at steady state, with all parameters functions of X and Y only.
    Last edited: Nov 20, 2017
  13. Nov 20, 2017 #12
    Yes, sir I got it. Now ?
  14. Nov 20, 2017 #13
    OK. Starting with Eqn. 2, in post #11, solve for U in terms of dP/dX using the boundary conditions at y = 0, h on U. What do you obtain?
  15. Nov 21, 2017 #14
    boundary condition: U= c at Y= 0,h
  16. Nov 21, 2017 #15
    Close, but this doesn't satisfy the boundary condition at Y=h. Please try once more.
  17. Nov 21, 2017 #16
  18. Nov 21, 2017 #17
    Much better. Now please take the integrated average of this between 0 and h to get the average axial velocity ##\bar{U}(X)##
  19. Nov 21, 2017 #18
  20. Nov 21, 2017 #19
    OK. Now divide by h to get Ubar.

    After that, please integrate the continuity equation (Eqn. 1 in post #11) between Y=0 and Y = h. What do you get?

    Unfortunately, I have to leave for a few hours now, but I'll be back later.
  21. Nov 21, 2017 #20
    Should I put the value of U and then evaluate?
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted