Monday, 17 June 2013

Further Swirl Investigations

I have spent the last few weeks gaining a deeper understanding of the processes governing swirl atomisation. I want to report on that now, together with giving a bit more information on the gas centred swirl coaxial concept.

I have been using the following principal references in my study of swirl atomisation:-


  • Theory and Practice of Swirl Atomizers, Yuri I. Khavkin, Taylor & Francis, 2004
  • Atomization and Sprays, Arthur H. Lefebvre, CRC Press, 1988
  • An Appraisal of Swirl Atomizer Inviscid Flow Analysis Part 1 & 2, Dr. John J. Chinn, Journal of Atomization and Sprays, 1993

My studies of swirl theory led me to believe that if the swirl chamber and outlet orifice were of the same diameter then the coefficient of discharge would be unity for all practical purposes.

Here is a photograph of the swirl inducer used in the trials:-




As previously mentioned, this swirler was formed by generating a two start, 3mm pitch metric trapezoidal thread on a 12mm diameter brass bar. The finished turned diameter of the swirler was 11.55mm. Due to the particular geometry of the profile tool used to cut this thread, it was not possible to depth it fully with a two start helix. Hence various precision measuring instruments and techniques were used to determine the key passage dimensions. The total equivalent flow area of the passages was calculated, using the equivalent diameter as defined in the last post. The total equivalent diameter of the swirl inducer above was 2.426 x 10^-3 metres.

Where EA = Equivalent Area, this gave:-


EA = 4.62 x 10^-6 metre square (0.007 inch square)

The swirler was made a transition fit in a short cylindrical body. The extremity of the inducer protruded approximately 1mm from the base of the body. In this way any swirl effects on the flow were eliminated. The measured flow would be the baseline level through the swirl inducer passages.

The assembled unit was tested using mains water at a pressure of 500kPa (72psi). The first test was a simple check of pressure drop. This was measured as 70kPa (10 psi) using an in line gauge.

The theoretical flow rate of the swirl inducer was then calculated using the standard pipe flow relation and the equivalent area:-

m = A (2p deltaP)^0.5  (1)

Where:-

m = mass flow
A= flow area
p = fluid density
deltaP = Pressure drop

Cd was not included as at this stage an ideal flow rate was being calculated. Substituting the known values in (1):-

m = 4.62 x 10^-6 (2 x 1000 x 70 x 10^3)^0.5

m = 4.62 x 10^-6 x 11.832 x 10^3

m = 0.0547 kg/sec (0.120 lb/sec)

The actual flow rate was then measured by timed discharge into an 8 litre bucket. The bucket was filled in an average time of 143.4 seconds, giving a flow rate of:-

8/143.4 = 0.0557 kg/sec (0.122 lb/sec)

This is a reasonable agreement given the limits of accuracy of both the flow area and flow rate measurement methods.

Reading Khavkin, I found vindication for my choice of equivalent diameter to calculate the swirler flow rate. Khavkin gives the following relation for swirl atomiser equivalent diameter, De:-


De = [(4{m/(2p deltaP)^0.5})/pi)]^0.5  (2)

Substituting the relevant values into (2) shows an equivalent diameter De of 2.433 x 10^-3 metres. This compares well with the calculation shown above.

Next the swirl inducer was fitted into a body having a 15mm long x 11.5mm diameter swirl chamber. The swirl chamber was parallel all the way to the outlet orifice, creating a so called "open" swirl atomiser.

This set up was trialled and the by now familiar thin swirling cone was produced:-






The break up length gives some indication of the swirl angular velocity. Due to the 15mm length of the swirl chamber, friction has reduced the velocity to the extent that the liquid viscosity and surface tension is able to resist the centrifugal disturbing force and produce a longer cone.

The discharge time for 8 litres with this set up was 146 seconds. This gives a flow rate of 0.0547 kg/sec. This compares very well with the calculated value for the swirler and the baseline figure without interference from the swirl chamber. If anything, a Cd value of around 0.97 to 0.98 is being displayed. This is so close to 1 that it is tempting to write this off to measurement and timing errors. That said, for all practical purposes the Cd can be considered to be unity.

The next trial was to investigate the length of the swirl chamber on the break up length and stability of the swirling cone. It should be emphasised here that a swirling external cone is not important in terms of the operation of a gas centred swirl coaxial injector. Rather the presence and condition of the cone is a metric of the quality of the internal swirling liquid film, which is an important parameter for a gas centred swirl device. The swirl chamber dimensions were reduced in 2mm increments by moving the swirl inducer towards the exit. It was found that the absolute minimum swirl chamber length for a stable cone was about 4mm. At 2mm swirl length the cone showed signs of instability (tearing). This could have been a symptom of shorter break up due to the reduced friction, but I believe that this explanation would produce an intact shorter cone, rather than the tear displayed.

The shorter cone produced by the 4mm swirl length:-



And the cone produced by the 2mm swirl length:-



The shorter cone length can be seen here, as can the tear in the fabric of the cone.

The flow rate for all of the swirl chamber length trials were of the order of 0.055 kg/sec.

It has been mentioned that the quality of the internal swirling liquid film is an important performance parameter for the gas centred swirl injector. There are various empirically derived film thickness relations in Lefebvre. Performance of swirl atomisers is highly dependant on specific geometry. None of the methods in Lefebvre related to an open swirl atomiser. Not surprisingly they did not give very sensible results when tried. 

According to Dr. John Chinn of UMIST, the air core within a swirl atomiser adjusts itself to permit the maximum flow rate. This is known as the principle of maximum flow, and is analogous to the flow of water over a weir.

This means that the exit diameter of the swirl atomiser is reduced by the diameter of the air core and effectively becomes an annulus. It was shown earlier that a swirl inducer equivalent diameter of 2.426 x 10^-3 metres was required to flow 0.055 kg/sec. As the same flow rate is exiting the atomiser, then the area of this annulus must be of the order of 4.62 x 10^-6 metres square. The air core in a standard swirl atomiser tapers outwards towards the exit diameter. In an open swirl atomiser the air core is parallel.

The photograph below shows the body with the swirl inducer pressed into position. It has been highlighted to show the air core annulus concept:-



The swirling liquid film thickness may be estimated as follows. If the flow area is A1 and the air core area A2, then the air core area A2 may be found by subtracting A1 from the total exit area Ae:-

Given that De is 11.5 x 10^-3 m,


Ae = 0.1038 x 10^-3 metres square

A1 = 4.62 x 10^-6 metres square

Hence A2 = Ae - A1 

A2 = 0.09918 x 10^-3 metres square

This means the air core diameter Dac is 11.23 x 10^-3 m. Subtracting this from the exit diameter De and dividing by two gives the thickness t of the liquid film:-


t = (D - Dac)/2

t = (11.5 x10^-3 - 11.23 x 10^-3)/2

t = 0.135 x 10^-3 m

t ~135 microns

Obviously I have no practical way of measuring this dimension apart from the assumption that particle size reflects film thickness to some degree. I had speculated in the previous post that the film thickness was of the order of 180 microns based on particle size. What does seem clear is that if the exit diameter is increased then film thickness should decrease. That would have a positive effect on atomisation and would also shorten the cone due to the thinner sheet having less resistance to inertial forces. That said, the friction from the larger swirl chamber wall surface area would tend to retard the angular velocity, somewhat negating this effect.

I trialled this by opening the last 10mm of the swirl chamber to 13mm. The swirler was installed to make the total swirl chamber length 15mm, as in the first trial which gave a cone length of 52mm. Here is a photograph of the 13mm De unit:-



It can be seen that the cone length has approximately halved. In addition, tearing can be seen on the cone surface suggesting it has been thinned. The discharge rate was still of the order of 0.055kg/sec.

The approximate swirling liquid film thickness for the 13mm De value is 115 micron, so approximately half of that when the De was 11.5mm.

For an effective gas centred swirl coaxial injector, a stable, thin swirling liquid film is essential. When this thin film is sheared by the central gas flow, very small droplets result. The possibility of swirling the central gas flow has occurred to me. I tried running compressed air at 690 kPa (100psi) through the swirl test unit. This produced a strong vacuum at the centre of the swirler, presumably where the air core was being generated. I suspect that in a production device this vacuum would disrupt the swirling film and lead to coalescence and droplet size increase. It may also lead to recirculation and attendant flow instabilities.

To conclude, the investigations so far carried out on swirl atomisation have shown the following:-


  • The Cd of the open swirl atomiser can be considered unity for all practical purposes
  • The equivalent diameter of the swirl inducer can be used to predict output flow
  • With the current swirl inducer, a minimum swirl length of 4mm is required for a stable film
  • Increasing swirl length increases stability of the film but increases cone length, decreasing swirl length decreases stability and decreases cone length
  • For the current swirl inducer it appears possible to approximate the liquid film thickness
  • Increasing the exit diameter of the swirler appears to decrease the liquid film thickness   
  • Halving the film thickness appears to halve the cone length (tentative)
These trials have provided some excellent information, all of which is highly pertinent to the design of a gas centred swirl coaxial injector. The design of this device can now commence. Keep watching this space.















6 comments:

  1. Hi,

    I just wanted you to know that there are folks like me who are fascinated by your posts, and really enjoy reading about your journey.

    I drop by about once a month to see if you've posted something new. And while your mathematics for atomization may be a bit beyond my reach, I look at them anyway, hoping that some of your knowledge will rub off...

    Please keep posting. I can't wait to see you build an engine and test it.

    Best,

    Bennett Dawson
    Vermont, USA

    ReplyDelete
  2. Dear Bennett,

    Firstly my apologies for the delay in replying to your comment. Thank you very much for your interest in my project and your encouraging words. My "knowledge" of atomisation is based on a solid engineering background, careful observation and common sense. So far it seems to be working. The project is taking far longer than I ever anticipated. I suppose I could have simply forged ahead and built a thrust chamber, only to stand and watch it melt/explode before my eyes. Rather I chose to build from the ground up; the result may turn out to be the same but I feel I will have learned a little more.

    Thanks again, British Reaction Research.

    ReplyDelete
  3. Note that Professor Yuriy Khavkin died May 17, 2013 at

    ReplyDelete
    Replies
    1. Yuriy was working on his books in New Hope, Minnesota at his daughter's house. I was assisting in some of the English. As an Electrical Engineer I could understand some of what Yuiy would explain to me in the process

      Delete
  4. We found a water fountain spraying with the configuration closed as you show and I took the picture for his book for Yuriy. We were working on the design of a new Turbine concept light weight - Yuriy died March 2013 in Brooklyn Center, Minnesota near family at age 92.- Yuriy knew Swirl Aomizers and now I know more than I did. Roy Yeater

    ReplyDelete
  5. Dear Roy,

    Thank you for your comments and for informing me of the death of Professor Khavkin. It is gratifying to know that he spent his last days engaged in work he loved and close to his family. His intellectual legacy is a sound one, and I send my heartfelt condolences to his family.

    Regards,

    British Reaction Research.

    ReplyDelete