Simple Box Model, Inlet/Outlet Flow

Hello, this is sort of a follow-up to a previous post.

I made a 10x10x10 box with 1-inch thick walls. There is one inlet and one outlet opposite each other.

I am trying to do a simple fluid flow in FLUENT through the inlet to the outlet. My wall temperature at 320K and the air inside initially at 300K.

 

No external boundary conditions. Wall-Fluid Boundary Coupled. Standard initialization at 300K, Box patched at 320K.

Inlet set to 0.5 m/s velocity magnitude. 

Using k-epsilon viscous model.

 

Everything appears fine in my set-up, but when I run the calculation my energy residual goes through the roof, setting insane temperature readings. Any idea what is causing the non-convergence?

Thank you

 

«1

Comments

  • DrAmineDrAmine GermanyForum Coordinator
    edited January 22
    Can you share some pictures of BCS? Please check the material properties. Do not patch the solid with any temperature
  • DrAmineDrAmine GermanyForum Coordinator
    edited January 22
    Moving to Fluids category
  • DrAmineDrAmine GermanyForum Coordinator
    edited January 22
    Moving to Fluids category
  • edited January 22
    I cannot get any screenshots at the moment, but i have the wall-box-fluid boundary coupled, and the inlet set to 0.5 m/s. That is the only change I've made.

    And why not patch the box to a temperature? That is the only way I have been advised to set initial conditions. Is there another way to do that?

    Thanks
  • DrAmineDrAmine GermanyForum Coordinator
    edited January 22

    Just for debugging I was mentioning that. If you are doing steady run than initialization is not really relevant: it is only a guess. Please add screenshots as soon as you can and check the material properties.

  • edited January 22

    I left all boundary conditions defaulted except for the inlet velocity. The materials are Air for the fluid and aluminum for the box.

    Here are pictures:

    A short summary of what I am doing:

    I simple heat transfer via temperature difference between room (320K) and incoming air (300K, air in room also initially 300K).

    There is one inlet and one outlet, left at 0 gauge pressure.

     

    Thanks

  • DrAmineDrAmine GermanyForum Coordinator
    edited January 22

    Can you share other information like residual plot and solution methods: I know the case is simple and I would probably recreate it within 10 minutes but we want to know what is going wrong here. Can you monitor the temperature inside the solid shell?

  • DrAmineDrAmine GermanyForum Coordinator
    edited January 22

    I cannot reproduce the issue on simple pipe with shell similar to your case: also check your mesh!

  • edited January 22

    When I initialize with the patch and check the countours, I get 300K inside the box and 320K in the shell just as I wanted, but it still wont converge.

    Here is my residuals. This time I did hybrid initialization with the settings previously screen-shot:

    Inside Fluid Contour:

    Box Interior Contour:

     

    Thanks

     

  • RobRob UKForum Coordinator
    edited January 22

    Please create a plane through the domain and plot temperature & flow velocity on that. 

    As an aside, never plot on default:interior or the fluid:interior type labels. You're OK here as it's a small model, when you display on this surface you're displaying on every facet in the domain/model. With 10-1000M cells you'll generally cause the graphics card to become very upset & go on strike. 

  • DrAmineDrAmine GermanyForum Coordinator
    edited January 22

    Can you show the grid and where one can see the BCs. It looks like there are some singularities in the corners. Moreover as it is steady-state and there is no heat flux defined at the solid earlier or later the solid will have the same temperature as the fluid. So think better in realistic boundary conditions.

  • edited January 22

    I get the same results, but here is what that looks like:

  • DrAmineDrAmine GermanyForum Coordinator
    edited January 22

    The same result as?

    Switch on coupled solver with pseudo-transient and run the case for longer time. Double check your mesh and use isothermal boundary condition at the outer wall of the solid first. There is definitely something wrong in your case which you need to figure out.

  • edited January 22
    Same result meaning I get the same temperature values. Sorry for confusion.

    Outside boundary should be isothermal because I didn't define it.

    As far as meshing, I didn't do any fine meshing details or editing, but would that be the cause of this non-convergence?

    Thanks
  • RobRob UKForum Coordinator
    edited January 22

    Can you plot velocity as well?  With negative temperature (in Kelvin) there's something very wrong with the result. 

  • edited January 22

    Here is a picture of the cross section with the velocity streamlines:

  • edited January 22

    Just ran the model with the Couple Solver and Pseudo Transient turned on. Here are the residuals and post fluent plane contours:

    Same issues as before, but temperature seems more even although still unrealistic.

     

    Thanks

  • edited January 22

    I am also getting the message "temperature limited to 5.000e3 in xxx cells" and "temperature limited to 1.000e0 in xxx cells" repeatedly with iterations as I run the calculation. Not sure if this gives any more information considering the energy equations are failing to converge. 

  • RobRob UKForum Coordinator
    edited January 23

    Plot vectors and mesh on the same plane. Residuals suggest the solution is not converged, so it's now a case of figuring out why. 

  • edited January 23
    Yes, why is the solution not converging. Any idea?
Sign In or Register to comment.