2

I'm designing a cooling system to cool nitrogen gas stream at 300°C with a mass flow rate of 15.6 g/s moving through an undulating channel in a cylindrical aluminum ring block with external fins. How much mass of aluminum or what special feature is needed to cool 15.6 g of nitrogen moving through and exiting to 25°C EVERY SECOND. It's to be noted that the block is immersed in a large pool of water, and that should be factored in as a cooling advantage. I'll love a method of calculation that allows flexibility with variables like mass flow rate, temperature, and channel diameter.

TechDroid
  • 213
  • 1
  • 8
  • One thing you'll need is the heat transfer coefficient(s) between the fluid and the block, which is a function of the channel's geometry, flow rate, flow regime, and temperature. You might find some rough estimates (order of magnitude) from tables in books or would need to do a simulation to get a more accurate value. Personally, I'd simulate the whole thing if you have access to the software. Otherwise, get ready for a bunch of differential equations. – user1318499 Feb 15 '19 at 06:08
  • Yeah, you made it more scary than I was imagining it. Just to clarify, doesn't the being immersed in so much water factor into the cooling effectiveness? – TechDroid Feb 15 '19 at 06:15
  • I believe he personally favor simulating the proposed scenario in real life or computationally than going through the complex math involved. Take it easy man. – TechDroid Feb 15 '19 at 17:29
  • @TechDroid Simulations of heat transfer coefficients do not give more accurate values than "rough estimates" from tables or books. This notion is seriously misplaced. Personally, I would sit down and understand the foundations to the problem before I would simulate anything (let alone "the whole thing"). Otherwise, get ready for a bunch of empirical guesses. – Jeffrey J Weimer Feb 15 '19 at 18:12
  • I have added an additional note to my answer to address the first part of your posting. By curtesy, you are to respond (by vote) whether what is now posted gives you what you need as an answer to your question. Alternatively, add further clarification to your question so that it can be address properly. – Jeffrey J Weimer Feb 18 '19 at 23:45
  • Please update your question to include the information that you have since given in the comments. – Jeffrey J Weimer Feb 19 '19 at 20:11
  • What happens when the body of water has heated up? – Solar Mike Feb 20 '19 at 11:58
  • It shouldn't, at least not by a considerable amount. It more like dropping a white hot Iron block into a lake and expect that to cause any noticable temperature different. However a subject of concern might be localize temperature increase around the heat exchanger, and I believe convection and other water dynamics happening should keep the immediate surrounding at the desired temperature. If you think not, I'll love to hear your thoughts. – TechDroid Feb 20 '19 at 13:38
  • I would prefer you to define the question appropriately, with values for the amount of water etc. – Solar Mike Feb 20 '19 at 14:07
  • I should probably close this discussion and open a fresh, more specific and detailed post. See you over there. – TechDroid Feb 20 '19 at 14:28
  • I am voting to close. The scope of the question continues to expand in comments, meaning its bounds are now becoming unclear. No disrespect is intended to the need for the OP to find a workable answer for the problem. All due respect is given to the demands to post a well-defined problem statement, to edit the problem statement as needed to keep the follow-up answers within bounds, and to put an end to the original question at a point where new questions that arise are clearly demanding their own thread. – Jeffrey J Weimer Feb 23 '19 at 20:52

2 Answers2

2

Cooling a Gas in a Channel (Tube)

Assume we have a cylindrical channel that has gas flowing in it. The channel is in a material that itself is cylindrical. We have a single pass tube heat exchanger or alternatively the equivalent of a 1-D extended system but with gas flow. The outside of the tube is exposed to air or water at a temperature $T_\infty$. The energy balance at a given point $z$ for the flow of the gas through this heat exchanger is basically as below.

$$ \dot{m}\tilde{C}_p \frac{dT}{dz} = - U P (T - T_\infty) $$

In this, $P$ is the perimeter area of the inside of the tube, $P = 2\pi r_{i}$. The factor $U$ is the overall heat transfer coefficient (based from the inside of the tube). It is a combination of the convection of the gas in the tube, the conduction through the tube, and the convection outside the tube.

Recast this in a dimensionless form with $\Theta = (T - T_\infty)/(T_h - T_\infty)$ and $Z = z/L$ where $T_h$ is the entering hot temperature and $L$ is the tube length. Allow $\beta = UPL/\dot{m}\tilde{C}_p$.

$$ d\ln(\Theta) = - \beta\ dZ $$

With the boundary condition $\Theta(Z=0) = 1$, the answer is $\ln(\Theta) = -\beta Z$. This gives the temperature profile of the gas in the tube. See this post for a comparable question and answer.

For the system proposed, $T_h = 300$ $^o$C. The value of $T_\infty$ should be BELOW 25 $^o$C. Otherwise, cooling could eventually require an infinite area heat exchanger (in essence an infinite length tube). Assume $T_\infty = 0$ $^o$C, and the goal is $\Theta = (25 - 0)/(300 - 0) = 0.083$. Use a 1 mm diameter tube with values of $\dot{m} =$ 16 g/s, $\tilde{C}_p = 1$ J/g for nitrogen, and $U = h_a = 50$ W/m$^2$ $^o$C. This latter assumption says that gas flow in the channel is the limiting value. The result is

$\beta = UPL/\dot{m}\tilde{C}_p = (50)\pi(0.001)L/((16)(1)) = 0.010 L$

$\ln(0.083) = -0.010 x \Rightarrow x = 250$

You need a 1 mm diameter tube that is 250 m long with walls that are cooled at 0 $^o$C.

Every 10x increase in heat transfer coefficient cuts the required length by 10x. Every 10x decrease in flow rate cuts the required length by 10x. Decreasing the external temperature $T_\infty$ gives a modest decrease in required length.

Flowing Gas Thru a Tube

The velocity of the gas through the tube is $v = \dot{V}/A$. Use $A = \pi D$ where $D$ is the diameter. The volumetric flow of an ideal gas is related to its mass flow as below with $T$ as temperature, $M$ as molar mass, and $p$ as pressure.

$$ \dot{V} = \dot{m} \frac{RT}{Mp} $$

Consider 1 bar at 300 $^o$C with nitrogen to obtain

$$ \dot{V} = \dot{m} \frac{(8.314)(573)}{(28)(101325)} = 0.00170 \dot{m}$$

Use the value 16 g/s to obtain $\dot{V} = 0.027$ m$^3$/s. In a 1 mm diameter tube, the required gas velocity becomes 8.6 m/s.

$$ v = 0.027 / ((\pi)(0.001)) = 8.6$$

At this point, we are left to find the required pressure drop to sustain 8.6 m/s through a 250 m long tube with a 1 mm internal diameter. This involves a review of the friction factor for fluids in pipe flow. The gas contracts as it cools. This demands further investigation too, because the basic calculations with incompressible fluids will fail with such cases.

Other Options

The analysis used a 1 mm diameter tube. Increasing the tube diameter will decrease the required length and decrease the required velocity. The issue at some point will be that a gas does not sustain a uniform temperature profile. It will cool along the wall of the tube but less at the center.

The analysis used only one tube. An alternative is to flow the same mass (volume) of gas through a multi-channel system. With $N$ tubes, the contact area increases by $N$. This means, the required tube length decreases by $N$. This alternative says, rather having the gas flow through one longer channel at higher speed, you have the gas flow through N channels at 1/N speed because area increases by N. The contact time calculated as area per volumetric flow is the same. The appreciation of why we flow gases through mesh grids to heat or cool them rapidly strikes home. As for the pressure drop, this issue should be addressed in a separate discussion.

A gas can be cooled by adiabatic expansion. An effective design may just be to allow the hot gas at high pressure to expand rapidly to lower pressure.

Finally, analysis at this problem shows an example of how a heat exchanger design equation is used for fluid flow in a pipe surrounded by a constant temperature fluid. This may be an easier approach to take from a practical standpoint.

Cooling a Gas Flowing in a Channel through a Block

The first part of the analysis above gives you the LENGTH OF A TUBE that is required to cool a gas to a given end temperature when the tube walls are held at a constant (lower) temperature.

The second part of the analysis above gives you the GAS VELOCITY that is required to pass the gas through the tube at a specific mass flow rate (g/s).

Assume that you will make your channel with a diameter $D$. Assume the block is cooler than your desired end temperature. Determine the length of channel that you need using the first part (the equation for $\Theta$ versus $Z$). Determine the gas velocity you need to support the desired mass flow rate using the second step. Finally, follow up with pressure drop calculations for gas flow in pipes. I suspect you will find out that you likely will not be able to get the velocity that you need in a single channel because the pressure drop that you will need across such a long channel is so high as to be physically unreachable. Also, any curves in the channel only act to increase the required pressure drop. So, a "curvy" channel will be far longer but require far higher pressure to force the gas through.

In the end, you may have to realize that a better option to cool the amount of gas you want as rapidly as you want is to run the gas through a multi-channel (porous) solid that is held at a lower temperature than you need.

Heating a Mass of Solid in a Given Time (Optional)

Admittedly, the first part of the question is likely not addressed in the above. Specifically, what time is needed for a mass of a solid to absorb a specific amount of heat (from whatever source).

The approach to answer this first question depends on making one or another of many different assumptions. When the gas at constant temperature is in full contact with the solid, the problem is an unsteady state convection + conduction problem. It may require analysis comparable to the use of Heisler charts. When the gas is in full contact with the solid but also changing its temperature, the Heisler chart analysis is replaced by an incremental or numerical integration approach (solve the Heisler chart at t, increment by dt and solve again). Finally, when the gas is flowing through the solid and the solid is not maintained at a constant temperature, one goes to the need for a Heisler chart or its foundations plus time integration plus position analysis.

Jeffrey J Weimer
  • 2,915
  • 6
  • 13
  • So you recommend adding more tube as well as meshes in the tubes? – TechDroid Feb 15 '19 at 16:46
  • The U in the analysis is meant for what material? – TechDroid Feb 15 '19 at 16:55
  • How did the Cp = 0.7 J/g be, is it an assumption or an actual value? My calc says 241.76 J/g using Q/m = c∆T. – TechDroid Feb 15 '19 at 17:10
  • "Every 10x increase in heat transfer coefficient cuts the required length by 10x. Every 10x decrease in flow rate cuts the required length by 10x. Decreasing the external temperature T∞ gives a modest decrease in required length." From the above statement, is the flow rate the mass flow rate or volumetric? Also in the whole analysis, the mass or thickness of aluminum surrounding the tube isn't factored, isn't that necessary? Just so you know, the tube/channel is bore through the aluminum block itself. – TechDroid Feb 15 '19 at 17:22
  • So I've thoroughly looked at the analysis and needless to say, it's of great help. All in all, the target of 15.6~16 g of nitrogen getting cooled from 300°C to 25°C is achieved in 1 SECOND. – TechDroid Feb 15 '19 at 17:46
  • Ok Jeffrey, you've come a long way and helped a lot, but since I'm just a regular dude looking for answers, things don't sink in as you might expect, so I'm really sorry if I'm asking supposedly obvious questions from your perspective. So to wrap it all up, your statement: "You need a 1 mm diameter tube that is 250 m long with walls that are cooled at 0 oC [to cool a nitrogen gas stream from 300°C to 25°C]", my question now is: the system I'm envisioning has a mass flow rate through the channel to be 15.6 g/s, and I need to cool all that in the same 1 sec, how does your analysis achieve that? – TechDroid Feb 19 '19 at 19:40
  • Just to reiterate Jeffrey, I'm trying to design an aluminum cooling device by boring tubular channel through a solid aluminum block and having fins attached to the outside in contact with the water. 15 g of nitrogen at 25°C must exit the channel every second, the nitrogen going in is 300°C. Now what I'm trying to figure out is how much mass or what geometry of the aluminum should I have to make that happen IN 1 SECOND. Remember the cooling aluminum block is immersed in a large body of water. So I'll love a method whereby I can change the mass flow rate or exit temp to tweak the design anytime. – TechDroid Feb 19 '19 at 19:53
  • I sincerely apologise for the misconception, I was just trying to make the question less personal and more beneficial to all. – TechDroid Feb 19 '19 at 21:03
  • So I used 20°C instead of 0 since it's more practical in the context of my design, then 5 mm diameter channel, and 46.8 g/s MFR. I would up with L~239.13m; V~0.079 m³/s; v~5.03 m/s. I took your advice and tried redesigning the system to run the gas through multiple channels N with length L each. While that solves the pressure drop problem, it raised another concern: won't that cut back on the time the gas spends in contact with the cooling system because the gas will have to be moving astonishingly fast at some point. We're talking about 4616 psi tapering driving pressure,... – TechDroid Feb 20 '19 at 06:53
  • ... it doesn't hit the rock bottom 0 psi, but come quite close relative to the initial driving pressure. The nitrogen is coming from a tank at 4616 psi initial pressure that tapers off really quickly, like 10 seconds release time. – TechDroid Feb 20 '19 at 06:58
  • I have added a note to my answer. Please credit my answer so far and perhaps close your question about cooling gas. Open a new one on pressure drop. Perhaps even open a new one on the actual design system itself. I am done using comments to continue to address an ever-expanding serious of questions that break off from the starting question. – Jeffrey J Weimer Feb 20 '19 at 13:52
  • It's been a great pleasure working with you and thanks so much for the answer. I think it requires me having at least 15 reputation to upvote your answer which I don't have. I'll do as you suggested. Thanks so much for the answer. – TechDroid Feb 20 '19 at 14:12
  • 1
    I'm glad this gives you a start on solving your problem! It gives me a basis for an interesting set of questions to pose in an undergraduate chemical engineering heat transfer course. :-) – Jeffrey J Weimer Feb 20 '19 at 14:32
  • @ Jeffrey, So I did post a new question focused more on the design, but was deemed a replicate, so I guess we'll have to continue here. Over there though, a user made a notable proposal of using an adiabatic cooling system instead, what's your take on this and how can it be applied to this system from your perspective? – TechDroid Feb 20 '19 at 18:47
  • So what do you say we explore that front. I've been doing some research and there are two ways about it: the Joule-Thomson cooling and the ideal gas adiabatic cooling. Both have their portion of influence in this scenario. As an expert what do you think? – TechDroid Feb 23 '19 at 15:43
  • I believe as the gas vents out of the tank it's stored, there's some real and ideal part to it based on the pressure it venting at a particular moment. – TechDroid Feb 23 '19 at 15:45
-1

Lets say you have Q per second calories of nitrogen heat to loose to the Aluminum for every 1 cm length of your pipe.

Then we have $$Q_{1\ cm\ length} = \frac{K*A(T_{hot}-T_{cold} )}{d(R^2_{aluminium}/R^2_{pipe})}\ $$

-K= thermal conductivity of the aluminum = $o.5 (cal/sec)/(cm^2 C/cm)$

-t = time, second

-R = radius

-T = temperature Celsius

-A = area of contact per one cm length of pipe.

-d = thickness of aluminum around pipe, approximate 0.8* one side, if the aluminum section is square.

So If we have the diameter of the pipe as D, we have A=pi*D and we can calculate d the thickness of aluminum.

If we have all the dimensions by plugging in the equation we can calculate the Q heat loss by nitrogen per second per each centimeter length of the pipe and if we multiply that by t we get total heat lss after some time.

Edit

After some comments I modify my answer as follows to include the length of the pipe and the time.

$$Q_{nitrogen} = \frac{K*A*L(T_{hot}-T_{cold} )}{d(R^2_{aluminium}/R^2_{pipe})}t\ $$ -L=length of pipe

This equation gives the Q loss by nitrogen after certain time t passing through the block. So if you need to loose more heat you have to make the pipe length L longer, or change other variables.

kamran
  • 21,714
  • 2
  • 18
  • 38
  • This is unbelievably helpful! This is great. Just to clarify, how does R²aluminum differs from R²pipe. And should the thickness of aluminum be measured from the center of the pipe to the edge of the block or otherwise? Does the equation allow to change the heat conductivity to something like J/s? – TechDroid Feb 15 '19 at 08:14
  • @TechDroid, R is from center, yes. R aluminium and pipe are substitude for a plane conductor thickness. of course you can change the units, but you be carful with dimensional integrity through out. disclaimer, this is a numerical method approximation. it gives a good insight into how the heat conduction works though. there is radiation and also the system is not linear. – kamran Feb 15 '19 at 08:24
  • The change in temperature in the equation must be for the desired temperature change in the nitrogen flow right? – TechDroid Feb 15 '19 at 08:25
  • So the aluminum being immersed in so much water means? – TechDroid Feb 15 '19 at 08:27
  • @TechDroid, this eq can have any one of varibles as either constant or hunted for. – kamran Feb 15 '19 at 08:28
  • We don't care. Water or air, we just assume T exterior of aluminum is that in eq. I got to go. – kamran Feb 15 '19 at 08:30
  • Sorry but I still don't get it, R pipe is the radius of pipe through the aluminum section, R aluminum is? – TechDroid Feb 15 '19 at 08:31
  • 0.8 of the side of square block cross section of aluminium. if the block side is a, R alum is 0.8a. – kamran Feb 15 '19 at 08:32
  • Just one more thing, (Thot - Tcold) is the change in temperature of the aluminum or nitrogen. I believe the pipe having a circular cross section makes it has a radius (Rpipe), how does an aluminum block have a radius. – TechDroid Feb 15 '19 at 10:37
  • The approach is not right. The approach calculates the steady state heat flow (W) across a block of aluminum with a thickness $d$ where the temperature difference is $T_h - T_c$ across the block and the area of the contact is fudged to be "area per length", meaning area per length of the block (or area per height of the block). This approach does not calculate the answer to the unsteady-state problem to cool a mass of gas from $T_h$ to $T_c$ is a given time $t$. – Jeffrey J Weimer Feb 15 '19 at 13:45
  • So I've you got a way out, a better approach? – TechDroid Feb 15 '19 at 15:01
  • Thanks for the update Kamran. So I've you taken a look at the input from Jeffrey, it seems like yours in a more broken down form. I'll love if you could merge both approach into one, or at least relate both so they can be used in different context to achieve the same result. Also, you're yet to answer my queries above on the Ralum variable and Temperature change. – TechDroid Feb 15 '19 at 18:19
  • @TechDroid The approach here cannot be "merged". It is simply wrong. It is a steady state mash-up to an unsteady state problem. It analyzes 1-D conduction through a slab when the system must instead be analyzed as heat exchanger or 1-D extended surface with convection. – Jeffrey J Weimer Feb 15 '19 at 18:30