# What should be the boundary condition for this case ?

Hello, I am trying to simulate this paper and I am stuck at this. I have tried piecewise polynomial function for specific heat and Ideal gas for density for air and for inlet condition I have used mass flow inlet and specified mass flow rate and temperature and for outlet pressure outlet. Still I am getting around 80K error. What should be the boundary conditions ? I have also tried applying gauge pressure with mass flow rate but solution diverges ? I have checked mesh quality.

Tagged:

## Answers

28MemberI think you have used correct boundary conditions. You need to find the error somewhere else.

Is your model 2D? The Swiss-roll recuperator is tubes with circular cross-sections, but if your model is a 2D and the pictures are showing your domain, note that the solver takes your model as a rectangular channel with theoretically infinite depth (perpendicular to the screen) so that gradients in the z-direction are zero.

8,855Forum CoordinatorIt looks reasonable. What's the near wall & free stream mesh like?

6MemberIn the paper they have specified that it is a 2D geometry so I have created a 2D geometry. The mesh looks like this. Also when I only provide mass flow rate the solution will converge but when I provide both mass flow rate and pressure Solution starts diverging.

8,855Forum CoordinatorMesh looks OK.

You shouldn't need to add a pressure to the inlets as it's not supersonic. Why a mass and not a velocity boundary? Is heat & mass conserved?

6MemberBecause in the paper they have given mass flow rate as Inlet Condition. Yes, heat and mass conserved.

8,855Forum CoordinatorOK, we tend to use velocity bc's over mass flow for incompressible flow. If mass & heat is conserved how well do cp etc compare with those used in the paper? Ie if you do a hand calc how far off on the properties would you be to account for 80K.

6MemberIn paper they haven't specified Cp value but have only told us that specific heat is temperature dependent. So I am trying to match outlet temperatures.

My simulation Results :

Output Cold : 812 K

Outlet Hot : 615 K

Reported Value in paper :

Output Cold : 759.3 K

Outlet Hot : 673.4 K

Here is the link to the paper : https://www.sciencedirect.com/science/article/abs/pii/S1359431108002895

8,855Forum CoordinatorSo, the solver is behaving and your problem now is to try and figure the material properties they've used. If it's temperature dependent you may struggle as you're trying to curve fit rather than find a single value.

6MemberI would like to know what could be the reasons behind this deviations and how could I replicate their results.

8,855Forum CoordinatorThe deviation is between a case in which all the fluxes balance, and a case where we don't know the material properties and (probably) all of the fluxes balance?

6MemberSo how can I overcome this problem ?