Simple Box Model, Inlet/Outlet Flow

  • 52 Views
  • Last Post 4 weeks ago
ilardimichael26 posted this 5 weeks ago

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

 

Order By: Standard | Newest | Votes
abenhadj posted this 5 weeks ago

Can you share some pictures of BCS? Please check the material properties. Do not patch the solid with any temperature

Best regards,

Amine

abenhadj posted this 5 weeks ago

Moving to Fluids category

Best regards,

Amine

abenhadj posted this 5 weeks ago

Moving to Fluids category

Best regards,

Amine

ilardimichael26 posted this 5 weeks ago

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

abenhadj posted this 5 weeks ago

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.

Best regards,

Amine

ilardimichael26 posted this 5 weeks ago

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

abenhadj posted this 5 weeks ago

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?

Best regards,

Amine

abenhadj posted this 5 weeks ago

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

Best regards,

Amine

ilardimichael26 posted this 5 weeks ago

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

 

rwoolhou posted this 5 weeks ago

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. 

abenhadj posted this 5 weeks ago

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.

Best regards,

Amine

ilardimichael26 posted this 5 weeks ago

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

abenhadj posted this 5 weeks ago

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.

Best regards,

Amine

ilardimichael26 posted this 5 weeks ago

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

rwoolhou posted this 5 weeks ago

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

ilardimichael26 posted this 5 weeks ago

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

ilardimichael26 posted this 5 weeks ago

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

ilardimichael26 posted this 5 weeks ago

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. 

rwoolhou posted this 5 weeks ago

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. 

ilardimichael26 posted this 4 weeks ago

Yes, why is the solution not converging. Any idea?

abenhadj posted this 4 weeks ago

Lets do a test: remove the shell solid zone at first, assign isothermal bc at the walls and use coupled solver and run: does it result in better behavior?

Best regards,

Amine

ilardimichael26 posted this 4 weeks ago

Okay, I suppressed the shell and re-meshed.

In Fluent I made sure the boundary from the fluid to outside was isothermal, and set my inlet to 0.5 m/s.

Here are screenshots:

Thanks

abenhadj posted this 4 weeks ago

Looks better but I still think your grid is not good. Know at the outer wall define shell conduction (in the wall panel) in Fluent and assign thickness and run again.

Best regards,

Amine

Close