-.. --.~'v 14...... .. 3; . $31,. 1.15? In )‘JII 'J‘.fu' . “ “5‘ ' I," #1” '.‘\‘3 l I? "1 l' ‘n " I . 1 ' l- ' w ‘ - "- I? n J . "1}, ‘n‘ AI I“ 111'“ r.'.'.‘ :‘muI‘I. I .I 1"3?‘ IirIL‘ .14 '«1. 3" '1 . In. 4-. thin??? ‘,t'-'.r".'-‘1 11‘ I .‘u.' . I 1.111.: I. '.:..,‘III' I y". > I I ._~ . 1- ., ." 0"“: III I .' .1 I . 1 I I1“ 'I ’I .1 . u'1I .I ‘1 ‘ ‘I 11 ' V " 11-w- ”11‘7"" :r-f '—1 ’1', 1‘ . . Ip-~ ' I." 1'. 'C'L ‘I. :1. '."' — I: f. ' -II- ‘F:?;?.—o‘4. v‘g " . “I1? . .‘ "m i 4.1““... , .1. A \ ,,‘ h 4 . -: ‘7 i "1311' 'M'UV .1111"?! 1' 'I'. "3 1| ’ '..I{ "U .f‘n'lsé'wi . 1-1' :I . .I ‘ «r 1le" , 3; - ":1 "a ‘ ‘I_ . . ‘. um" I .1 1. a - .3i .1 'I3‘ “LII ..... .1. ’fn kuw 'mw'exh IIIII ”w {”1" I «MI-11.1%, H n -. ' 91': .H‘ I? 3.??1?‘ 4?; .11 ':;11111f.'gl11?W'?|11 1IIIEI?' II . .I I11ItI.0I ,' ‘I. I.' II"_I-°I I I . I_ iII 'I.II “I '1"'§;;"i:."?? ?‘ L. " W :13 20‘? "\t I‘llgtlfi " '0, 11?.3 2‘ "3 w? t 2? ~ 3a; .1 " a: ‘V‘ 111 ."“'" .1 : ‘I.I 'I .1 111111.11; .‘ IIIII I>TI~. “:1". III .‘ '. ‘|.. Iffl 1.1% ' 1 '1? u ?',! 3114.1»? 1? .:';%%h' mw 1 I , ? ‘I'n: Wr'fiby~ I111311141? H III-I1II ”mug I1113'11I1I1'1 1 ll NY'I" -.I'.'IILI& J131Im' 2'131I 15.! x iv «M l I n 1"”.- l - ' H3”: ' I. ‘ ' I: '3'. . . - . , 1‘13: 'l: 3339"” 1 ”'3 ' ' 1‘1""? '- . ‘l ’ .‘u‘ " . ' ‘H I I ¢ I l I‘ :II “I! r II 1 I ' 1, L I 1 0 41 v I 1| 1 ‘ I 1 y 1 1. 1111.1, If: .1? . I I I l , '- I" .< ‘I 1111'II .1 .I> .19. I " '_ ' I ' I111“: I1{II-I11 \‘J'I 1 '33? JUL?" 1. :I‘ ? 'I.?ll.1;'1...11:' III?‘ I«1.-II.I 111?“ .' .~ A ., .. “ “In“. 111 1 II .43 [41,11 “111.1" 5 ul 11? 1 :1 1.1:: 1-, I - I 5 I11? I,' 1 .1 11/ a I I III I “W! .m? [21”??? I ’ ' . R '. MI: I I??' “ll?“ 3???? .I" "?????| W?!” I? ?I ll? M??? n? ’l I??? x? t‘ II'L‘: THE5‘5 lllllll Hill llllll lllllll/ llllllllll 3 1293 301070 3258* :‘ 2.1:? ”New ' 3,. in“ sitar-e - . U22"? as; I x _ / This is to certify that the thesis entitled Drying of Alfalfa Under Controlled Temperatures and Vapor Pressure Conditions. presented by Natalie Jo Wentz-Carroll has been accepted towards fulfillment of the requirements for M . S . degree inAgricultural Engineering 5 /7/ ///"' Major professorLS Date W 0-7 639 MS U is an Affirmative Action/Equal Opportunity Institution MSU LIBRARIES v RETURNING MATERIALS: Place in book drop to remove this checkout from your record. FINES will be charged if book is returned after the date stamped below. DRYING OF ALFALFA UNDER CONTROLLED TEMPERATURE AND VAPOR PRESSURE CONDITIONS BY Natalie Jo Wentz-Carroll A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Department of Agricultural Engineering 1984 ABSTRACT DRYING OF ALFALFA WITH CONTROLLED TEMPERATURE AND VAPOR PRESSURE BY Natalie Wentz-Carroll A controlled environment drying chamber was built and used to study the drying of alfalfa at various temperatures and vapor pressures. Several samples were treated with potassium carbonate to determine the effect this chemical had on the drying rate of alfalfa. Alfalfa samples were cut and placed in the drying chamber where the weight change was monitored to determine' the rate of drying. Drying rates were compared by calculating an average drying constant, an average mass transfer coefficient and by exponential curve fitting. Results show that drying constant and the mass transfer coefficient both can be used to model the drying curves of alfalfa. Drying curves were analyzed in three intervals, with cutoff points of 2.0 and 1.0 % moisture (dry basis). Drying temperatures had a direct relationship to the drying rate. Potassium carbonate treatments speeded the drying of alfalfa, espically for moistures above 1.0 %. This work is Dedicated to My Husband, Dale Michael Carroll ii ACKNOWLEDGMENTS Throughout the course of this graduate program the author has been very appreciative of the exceptional assistance provided by her major professor, Dr. C. Alan Rotz. His help, concern and consideration were unremitting. The author also wishes to express her gratitude for the time and involvement of Dr. R. Brook and Dr. J.W. Thomas throughout the duration of her research. The faculty and staff of the Michigan State Agricultural Department were always an infinite source of support, suggestion and concern for which the author thanks them. The assistance of Dr. R. Byler and D. Sprott was always useful and appreciated. iii TABLE OF CONTENTS LIST OF TABLES .0.00.000000000000000000000...0.0.0.0....v1 LIST OF FIGURES I.0..0....O...OOOOOOOOOOOOOOOOOOOOOOOOOOVj-ii CHAPTER PAGE 1. INTRODUCTION O...O0.0.0.0....OOOOOOOOOOOOOOOOOOOOOOO 1 PrOdUCtion 0..OOOOOOOOIOOOOOOIOOOOOOOOOOOOOOOOOOOOO. DrYing OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO0.0.00.0.0... NM 2. OBJECTIVES 0..0......0....OOOOOOOOOOOOOOOOOOOOOOO0.. S 3. LITERATURE REVIEW Plant Structure .................................... 6 Stomata and Guard Cells ........................ 8 Cuticle ....................................... 10 Water Loss ........................................ 11 Dry Matter and Respiration Losses ................. 13 Increasing Drying Rates ........................... 14 Mechanical Treatments to speed drying ......... 15 Chemical Treatments to speed drying ........... 16 Radiation ......................................... 18 Drying Analysis ................................... 18 4. EQUIPMENT Chamber and Air Handling 0.0.0000000000000000000000 22 Microcomputer and Control ......................... 29 SO£tware 0.0.00IOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO. 36 iv 5. DATA COLLECTION 0......OOOOOOOOIOOOOOOOOO00.0.00... 39 6. ANALYSIS OOOOOOOOOOOO0.0...OOOOOOOOOOOOOOOOOOOOOOO. 44 DrYing constant OI.0.0.0.000...OOOOOOOOOOOOOOOOOOOO‘ . . 45 Mass Transfer Coefficient ......................... Exponential curve Fitting O O O O O O O O O O O O O O O O O O O O O O O O 0 ga 7. RESULTS AND DISCUSSION ............................ 51 crop Maturity 0....0.0.0....OOOOOOOOOOOOOOOOOOOOOOO 56 DrYing constant 0 O O O O O O O O O O I O O O I O O O O O O O O O O O O O O O I O O O 59 M355 TranSfer COeffiCient I O O O O O O O O O O O O O O O O I O O O O O O O 63 Exponential Curve Fitting ......................... 66 Chemical Treatment O...OOOOOOOOOOOOOOIOO0.0.00.0... 66 8. SUMRY 00....COO....0.0...OOOOOOOOOOOOOOOOOOOOOOOO BIBLIOGRAPHY OI...OO...OOOOOOOOOOOOOOOOOOOOOOOOOOOOO... 80 APPENDI XES A. Computer Programs ............................. 86a B. Hourly Listing of Sample Data ................. 90 C. Drying Constants and Mass Transfer Coefficients 100 D . Plot of the Sample Moisture Content over Time 109 LIST OF TABLES Table 1. Regression Coefficients for the Temperature sensor calibration OOI.00......OOOOOOOOOOOOOOOOOOO. 2. Regression Coefficients to Obtain the Sample Mass used in Trials 1-25 OOOOOOOOOOOIOOOOOOOOOOOOO. 3. Alfalfa Maturity and Time of Cutting for all Trials 0......O...00.0.00...OOOOOOOOOOOOOOO00...... 4. Temperature Settings, Drying Temperatures, Vapor Pressures 'and Concentration of Moisture in the Drying Air for all Trials ......... 5. Volume of Water Displaced by Alfalfa at Different MOiSture contents OOOOOOOOOOOOOOOOOOOOOOO 6. Drying Constant and Mass Transfer Coefficient for Alfalfa Trials at Two Vapor Pressures and Different Temperatures .............. 7. Drying Constant and Mass Transfer Coefficient for Alfalfa Trials at Two Temperatures and Different Vapor Pressures ........ 8. Mean Drying Constant for Three Moisture Content Intervals with Different Temperatures and Vapor Pressures, and the ' Results of an Analysis of Variance on this Data .0.COO...O0......0..OOOOOOOOOOOOOOOOOO00...... 9. Mean Mass Transfer Coefficient for Three Moisture Content Intervals with Different Temperatures and Vapor Pressures, and the Results of an Analysis of Variance on this vi 10. 11. 12. 13. Data 0.00.0000...0.0.0.0...OOOOOOOOOCCOCOOOCOOO0.0. Coefficients of the Two Exponential Curve Fit to the Data, in Order of Increasing Temperatures 00......OOOOOOOOOOOOOOOOOOOOOOO0.0.0.. Coefficients of the Two Exponential Curve Fit to the Data, in Order of Increasing vapor Pressures 0.00.0.0...OOOOOOOOOOOOOOOOOOOOOOO. Mean Drying Constant and Mass Transfer Coefficient for Three Moisture Content Intervals at Four Environmental Conditions With Treated Alfalfa 0.0000000000000000000000...... Comparison of the Drying Constant and the Mass Transfer Coefficient Values for Treated ' Alfalfa Dried Simultanously with Control Alfalfa OOOOOOOOOOOOOOOOOOCOOOOOOOOOOOOOOOOOC0.0... CALCKHD, Program to Calculate the Drying Constant and the Mass Transfer Coefficient ........ BMDPBR, Program to Calculate a Two Exponential Best Fit Curve using BMDP (BiOMedical Data PaCkage) 00......OOOOOIOOOOOOOO... Listing of the Moisture Contents on an Hourly Basis for all Trials ....................... Output form CALCKHD, Drying Constants and Mass Transfer Coefficients ....................... vii 68 69 71 863 89 90 100 LIST OF FIGURES Figure Page 1. Plan View of the Equipment, Drying Chamber, Aminco Aire Unit, and the Microcomputer ............... 23 2. Cutout View of the Drying Chamber Showing Transducers, Sensors, Baffle and Alfalfa ”Olders 0.0.0.0000...0.0.0.000...OOOOOOOOOOOO0.00...... 24 3. Detail View of the Alfalfa Holders and TranSducer Mount 0.0.0.0....OOOOOCOOOOOOOOOIO00......O. 26 4. Block Diagram of the Process Controller and Data COlleCter .0.0000000000000000000000000000000000.0. 31 5. Plot of the Change in Moisture Content of the Actural Data over Time for Different Temperatures at Two Vapor Pressures ................... 57 6. Plot of the Change in Moisture Content of the Actual Data over time for Different Vapor Pressures at Two Temperatures ................... 58 7. Example of the Two Exponential Curve Fit (amp) for Tria151and2 OOOOOIOOOOOOOO0.0.0.0000...O. 67 8. Linear Regression Curves for Drying Constant Based on the Drying Temperature for Treated and Control Samples at Three Moisture Content Intervals ..................................... 73 9. Plot of the Two Exponential Equations of Best Fit Generated Data on Potassium carbonate 0.00.0.0.0.00000000000000000000000000......O. 75 D. Plot of the Actual Change in Moisture over viii Time for all Trials 0000......OOOOOOOOOOOOOOOOOOOOOOOO 109 ix CHAPTER 1 INTRODUCTION Alfalfa has been cultivated throughout recorded history (Bolton, 1962). Columella (approximately 60 AD) praised it for its longevity, soil improving qualities, ability to be cut four times a year, and for its feed value for cattle and horses. Alfalfa is now grown worldwide in the temperate zone and has been adapted to many soil types, rainfall amounts and a range of elevations. The major benefits of alfalfa over grasses are an increase in protein, vitamin and mineral content as well as greater yield, drouth resistance and nitrogen fixing capability. Compared to other forage legumes in the North Central area of the United States, alfalfa provides a higher yield. Production In 1981 alfalfa hay (including alfalfa mixes) was grown on 10.7 million ha in the United States yielding 7.6 million tonnes of hay (Agricultural Statictics, 1982). This was 44% of the total area harvested for hay (24.4 million ha) and 58% of the total number of tonnes produced. The production value was given as 9.2 billion dollars. In Michigan, alfalfa production exceeds that of all other forages combined. In 1981, 0.5 million ha were planted in alfalfa (79% of total hay area) which yielded 3.0 million tonnes (85% of all hay production). Michigan produces about 5% of the alfalfa grown in the United States. Michigan averaged 6 tonnes of hay and haylage per hectare in 1981 generally produced with a 3 cut system. (Michigan Agricultural Statistics, 1982). Drying When alfalfa is cut the plant generally contains 75-80% water. For safe storage as baled hay, the moisture content must be reduced to below 20%, wet basis. Usually this is accomplished by leaving the cut crop to dry in the field with occasional turning. When the desired moisture content is achieved, the alfalfa is baled and taken off the field. The drying process usually takes 2-4 days if no rain occurs. Rain on a cut plant causes leaching and microbial growth. This reduces protein and nutrient value, and increases dry matter loss. Respiration, the oxidation of organic compounds causing a release of energy, occurs in the uncut plant for maintenance and development. This process continues in the cut plant until it has dried to about 40% moisture content causing further dry matter losses as the plant continues to oxidize organic compounds (Wolf and Carson, 1973). The total losses during hay making have been estimated as high as 30%. At this rate 4.3 million tonnes of alfalfa were cut in 1981 to yield 3 million tonnes produced. This substantial loss (1.6 million dollars) would be decreased by reducing haying losses. Many methods to reduce these losses are being studied. Early investigation to speed drying and thereby reduce the chance of loss due to rain damage involved mechanical conditioning such as crushing, crimping, or laceration. These methods increased dry matter losses because increased mechanical manipulation caused greater leaf loss (Jones and Palmer, 1932; Kepner et al., 1960). Some chemicals, such as fusicosin, organic phosphates and potassium carbonate showed increased drying rate of forages and grasses in the laboratory (Turner, 1969; Harris, 1978; Tullberg and Angus, 1972; Wieghart et al., 1980). Investigation is currently being performed to determine the adapability of carbonate solutions to field use (Rotz et al., 1984). CHAPTER 2 OBJECTIVES The overall objective of this research was to quantify the affects of temperature and vapor pressure on the drying of alfalfa, with the intention of obtaining an equation based on envionmental factors to determine the necessary drying time for baling. The specific goals were: 1) To develop equipment for monitoring alfalfa drying rates in a drying chamber which allowed control of air temperature and vapor pressure. 2) To develop a procedure for data collection and analysis in order to compare drying rates under various treatments . 3) To quantify the affect of potassium carbonate solution on the drying rate of alfalfa under controlled conditions. CHAPTER 3 LITERATURE REVIEW Plant Structure The alfalfa plant is a perennial legume capable of fixing nitrogen (with rizobia bacteria). It is a dicot with a compound leaf consisting of three leaflets attached to the petiole that is attached to a true stem. The plant has a tap root which makes it more drouth resistant than grass. The tap root, however, increases the plant's subjectivity to winter heaving and therefore less winter hardy than grasses. Regrowth is by crown buds which allow alfalfa to be cut up to 10 times each year, with 3-4 cuts generally used.. As the crop matures the stem tissue strength is increased and the leaf to stem proportion decreases (Harris and Tullberg, 1980). The alfalfa stem is mainly a support and a conduction organ. Xylem transports water and dissolved materials up to 6 the leaves and phloem carries dissolved food material downward. The xylem and phloem are separated by the cambium. The plant stem is coverd by an epidermis, as is the leaf. The cuticle (or integument) is a multilayered membrane which covers the epidermal cells. Cutin is a waterproof waxy layer which covers the entire plant above ground. It is responsible for the tough, durable nature of the cuticle. The cutin density is related to many environmental conditions during plant growth, including temperature, radiant energy, humidity, and soil moisture availibility (Hadley, 1980). Hamilton (1975) adds air movement, stress conditions, elevation, soil nutrients and weathering to the list of factors which affect the wax density on a growing leaf. The cuticle makes the leaf nearly impermeable to water and carbon dioxide except by way of openings in the leaf called stomata which are surrounded and controlled by turgor of the guard cells. Turgor is controlled by water take-up or water loss by the guard cells (Zelitch, 1967). The plant Cuticle protects cellular tissues, absorbs ultraviolet light, conserves water and nutrients, participates in respiration, protects against infection, parisitic and herbicide damage and influences plant appearance (Hamilton, 1975). Water vapor must leave the plant through the stomata, through breaks in the surface, or by diffusion through the cuticle, (Mears and Roberts, 1970). Moisture moves as a vapor from regions of higher partial vapor pressure to regions of lower partial vapor pressure (Barre, 1938). The rate of diffusion is, therefore, proportional to the vapor pressure gradient and inversely proportional to the resistances of moisture movement in the plant material. Moisture diffusion through the cuticle is much slower than water loss through the stomata. Therefore, when the live plant is under water stress the stomata close. Stomata and guard Cells As was previously mentioned the stomata control uptake of CO2 as well as water .loss. The extent of stomatal opening is due to many different metabolic activities in the guard cells (Walker and Zelitch, 1963). Whitney et al. (1969) showed that some stomata do not change apperature with changing stimuli and seem "locked" in place. In the growing plant light is the major factor causing guard cell turgor to increase which opens the stoma. Zelitch (1961) found that 4300 dekalux of light for a 90 minute exposure fully opened the stomata of tobacco leaf discs floating on water. Covering the leaf disc closed the stomata in 30 minutes. Zelitch floated the discs on water as he found an abundent supply of water essential to obtain ’7‘ fully open stomata. He further stated that the substance controlling the guard cells was generated in or near these cells since covering half a leaf disc yielded no open stomata in the covered portion and 70% open stomata from the half in the light, within 90 minutes. Hense the controlling mechanism was affected by the light/dark conditions at each stomata and not for the leaf as a whole. Turner (1969) suggested potassium ions moving into guard cells could account for the increase in turger which opens the stomata. This would also indicate that stoma work independently of each other. Under field conditions stomata have been shown to close rapidly for the first two hours after a plant is cut to conserve water (Jones and Palmér, 1932; Pederson and Buchele, 1960). Oxygen is necessary for normal stomatal opening (Walker and Zelitch, 1963). Temperature and water availability also affect stomatal opening (Whitney, et al., 1969). Carbon dioxide causes closure of stomata in light (Staflet, 1957). As a plant wilts guard cells loose turgor which cause the stomatal closure in an effort to save internal water. (Zelitch, 1967; Kramer, 1963). The internal water balance and degree of water stress depends on the relative rates of water absorption and water loss (Kramer, 1963). Different methods have been used to indicate the relationship between the stomatal opening and water loss. 10 Staflelt (1957) suggests that guard cell width is a better indicator of water use than stomatal width. Meyer and Anderson (1939) stated that the rate of water diffusion is proportional to the perimeter of the stomata and not the area of the opening. Others, however, continue to use width of stomatal opening to indicate changing transpiration rates (Zelitch, 1961 and 1967; Walker and Zeltich, 1962; Schonherr, 1976; Whitney, et al., 1969). Whitney found that for alfalfa dried at high temperatures the degree of stomatal opening (partially vs. fully open) did not significantly influence drying rates although there was a significant difference in the drying rates of alfalfa leaves when the stomata were open as compared to when they were closed. Cuticle The cuticle plays the major role in water conservation by epidermal cells (Schieferstein and Loomis, 1956). It is the main factor affecting the deposition, distribution and retention of chemicals applied to forage as aqueous solutions or suspensions (Holloway, 1969). Other functions have been discussed previously. Once the stomata close cuticular resistance to water loss causes the reduced transpiration rates (Schieferstein and Loomis, 1956). For this reason much research has focused on reducing cuticular 11 resistance. Schonherr and Bukovac (1973) found that water permeability depends more on the state and distribution of soluable cuticular lipids than on their chemical composition. Hartel (1947 and 1951, translated from German and discussed in the paper by Schonherr and Bukovac, (1973)) observed that the maximum cuticular transpiration occured when pretreated with buffer ions at pH = 7. Schonherr (1976) found that water permeability was dependent on pH and cations of the buffer solution. The efficiency of cuticular waxes in limiting cuticular transpiration is a function of both the quantity and form of wax on a leaf surface (Denna, 1970). The amount of wax is higher in plants grown on soils low in phosphorus, potassium or nitrogen. Waxes are also denser and more compact if the plants are exposed to excessive winds. Waxes on the crown of the cuticle are often removed by weathering (as when two leaves rub together). Brushing leaves increases transpiration in the cuticular phase after the stomata have closed (Hall and Jones, 1961; Hamilton, 1975). Hamilton found that exposure to hexane vapor increased transpiration. Water Loss The amount of cutinization is also increased by continued plant water stress (Kramer, 1963). As mentioned, 12 water stress causes stomatal closure in an effort to retain water. Water loss appears to be linearly related to vapor pressure differences with the major limitation being cuticular resistance (Leshem et al., 1972; Mears and Roberts, 1970). Hence if the living plant is stressed, the amount of wax on the leaves increased in order to decrease water loss. Movement of water in a living plant is due to changes in the water potential. Water potential is the chemical potential of water in the plant or conceptually an expression of the Gibbs free energy of water in the plant (Slayter and Taylor, 1960; Merva, 1975). Under normal growing conditions, water is taken up by plant roots from a region of higher potential to a region of lower potential. It changes state from a liquid to a vapor both within the atmosphere outside the plant as well as within the plant. The rate of water uptake is limited by available water and the amount of water being transpired. Water vapor exits the plant primarily through stomata and to a smaller degree through the cuticle. The cut plant looses water fairly rapidly until the stomata close. Then drying is slowed as water must be moved through the resistive cuticle layer, or through breaks in the plant (Mears and Roberts, 1970). After the initial drying period is primarily from the stems through the leaves in the cut 13 plant (Harris and Tullberg, 1980). 251 Matter and Respiration Losses After afalfa is cut it continues to respire until it reaches 35-40% moisture content (Greenhill, 1959; Wolf and Carson, 1973). Most dry matter, nutrient and protein losses occur in the field, and these losses increase if the crop is damaged by rain (Shepherd et al., 1954). These researchers found a 38.5% loss of leaf dry matter occured for hay which was not damaged by rain. This value was increased to 47.3% for hay wetted by two showers and 74.5% if wetted by three showers. Klinner (1975) reported up to 30% dry matter losses to be expected during hay making. At any temperature dry matter losses are inversely related to the saturation deficit which is related to the drying rate. In other words, as the drying rate increases dry matter losses decrease (Greenhill, 1959). Greenhill also showed that as temperature increases dry matter losses increase. He sights a 7% dry matter loss due to continued respiration. Greenhill suggests that continued photosynthesis may offset respiration effects. Respiration is a linear function when related to moisture content but it increases exponentially with temperature to 25'C (Wood and Parker, 1971). Respiration ceases at tissue temperatures 14 (greater than°55 C (for 15 minutes) or at a moisture content less than 40% (Wolf and Carson, 1973). Microwave treatment for 30 seconds can also stop respiration (Priepke and Bruhn, 1970). Increasing Drying Rates Various treatments have been used to speed the drying rate of grasses and forages. Byers and Routley (1966) showed that steaming increased water movement from alfalfa. Whitney et al. (1969) showed that drying rates increased in grass specimans when leaves were steamed, cut, split, or exposed to petroleum vapor for 60 seconds. Exposure to petroleum ether vapor was shown to be a more effective treatment to increase the drying rate than dry heat, steam, or steam and petroleum ether vapor. The treatment has a similar effect to longitudinal splitting of the stem as studied on ryegrass (Harris et al., 1974). The effect of petroleum ether increased as the temperature of the ether increased since the waxes on a leaf were more easily removed with the hotter vapor. This study showed that holding leaves over hot petroleum ether vapor for 10 seconds was as effective as dipping for 60 seconds in a cold liquid petroleum ether solution. Harris and Thaine (1975) showed that mechanical and thermal treatments increased the drying rate in leaves and 15 to a greater extent in stems. Application of a flame to standing and cut plants also reduced drying time (Person and Sorenson, 1970). Mechanical Treatments To Speed Drying Mechanical field treatments such as crushing, crimping and maceration increase drying time but also normally increase dry matter losses (Jones and Palmer, 1932; Kepner et al., 1960; Shepperson, 1974). Mechanically treated plants initially dry faster since more cells are exposed to the air. Once the exposed cells are dry the drying the rates are the same as for untreated plants (Byers and Routley, 1966). Dry matter losses can be reduced by using less severe mechanical treatments such as tedders or turners (Shepperson, 1974). Chop length (25-100mm) had no observable effect on drying rates (Menzies and O'Callaghan, 1971) so minimal chopping is preferable. A method of turning after some drying has occured is essential since air movement through the swath is restricted by a seal formed by the top leaves as they wilt (Barrington et al., 1970). The rate of carotene losses is increased by crushing but the final carotene content is not affected since the plant dries faster (Fairbanks and Thierstein, 1966). Other investigation concerning mechanical treatment is not discussed here as the intended research does not involve this aspect of forage drying. 16 Chemical Treatments gg §pggg Drying Chemical treatments have also been used to alter drying rates. Herbicides caused no significant increase in drying under ideal conditions but increased the drying rate in unfavorable weather conditions when drying to 40% moisture content. Drying to 25% MC was the same for treated and control hay and treatment caused a loss in the green color of the dried product (Kennedy et al., 1954). Shepherd (1959) found that the herbicide ethylene dipyridelium dibromide did not affect subsequent plant growth. Dilute solutions of fusicosin increased cuticular as well as stomatal transpiration (Turner, 1969). Dipping in sodium azide increased the initial drying rate only (Tullberg and Angus, 1972). Tri-n-butyl phosphate doubled the drying rate of grass and reduced dry matter losses (Harris and May-Brown, 1976). Briphos compounds (orgnaic phosphates 16M, COGD and COZOD) also increased the drying rate of grass and reduced dry matter losses (Harris, 1978). The utilization of potassium carbonate to speed drying has been used in the grape industry in making raisins. In fact, Columella (60 AD) suggested the use of ashes and oil to increase the drying rate of grapes. More recently research has substantiated the use of KZCO3 to speed the drying of grapes in making raisins (Dudman and Grncarevic, 1962; Grncarevic et al., 1968). That this treatment does 17 not remove sufficient wax to speed drying (Dudman and Grncarevic 1962) so investigations have centered on finding other possible reasons for the increased drying rates. Increased drying rates may be due to a change from a hydrophobic to a hydrophylic character of the grape surface whereby the transport of water may occur in a liquid phase across the cuticle in emulsion filled air spaces rather than the slow vapor diffusion process (Dudman and Grncarevic, 1962; Chambers and Possingham, 1963). The effect of the potassium carbonate was reversed by washing the plants. In fact the drying rates slow to that of untreated plants when washed, indicating that the waxes are not permenantly altered (Grncarevic et a1. 1968). These researchers also found an artificial wax covered membrane system showed a similar reaction to drying when dipped in an oil and potassium carbonate solution. Potassuim carbonate (KZCO3) has been shown to increase the drying rate of alfalfa leaves and stems (Tullberg and Angus, 1972; Wieghart et al., 1980), affecting the stems more than the leaves (Tullberg and Angus, 1978; Johnson, 1983). This may be of benefit to reduce overdrying of leaves, which is a major problem in the drying of alfalfa. To reach an average moisture content of 20% for baling the leaves are generally dryer than 20%. The leaves then become brittle and are easily shattered and lost. This is a 18 paticular problem as the alfalfa leaves contain more protien than the stems. Radiation Increasing thermal radiation or the exposure time to this radiation increases moisture removal rates from plants (Person and Sorenson, 1962). These authors also found that at equal radiation intensities drying rates are dependent upon the wave length of radiation. Incident radiation is affected by clouds, rain and perhaps wind (VanElderen et al., 1972). Solar radiation increases the drying rate because leaf temperature is raised which affects the vapor pressure deficit (Hill et al., 1977). Ajibola et a1. (1980) have studied the solar spectral absorptivities for chopped, macerated and dewatered alfalfa in the field. However, there has not been extensive research of the effect of light as it relates to drying and its interrelation with various methods used to speed drying. Drying Analysis An exponentially decaying curve is generally used to model the drying of many agricultural products. This curve is generally of the form: MR = exp(-kt) where MR = (M-Me)/(Mo-Me) t - time (hour) k - drying constant M - the moisture content 19 Me - equilibrium moisture content Mo - initial moisture content for the interval In 1969 Whitney, et al. used an exponential curve of the form: (M-Me)/(Mo-Me) = B*exp(-kt) where k s C*exp(ST) C and S were determined from experimental data. Menzies and O'Callaghan (1971) used a similar equation and further noted that forages dried at temperatures under 80 C% dried in two or three distinct stages, with differing coefficients for each stage. These researchers found the drying constants were related to the temperature but not to the humidity of the air. The first cutoff point for the drying period was defined by MC = 0.83 + 0.58 Mo where MC - moisutre content Mo - initial moisture This relationship fit the data with a correlation coefficient of 0.81. The other cutoff moisture contents were not given in this paper. The equation to determine the drying constant was: k = Cl*exp(C2*t) t - temperature, degrees C C1,C2 - 0.031 and 0.024 for the first interval - 0.0078 and 0.037 for the second interval - 0.0027 and 0.046 for the third interval 20 Kemp et a1. (1972) used an equation similar to that of Whitney (with B = exp(C2) and k = -C1*LE): MR = exp(-C1*LE*t)exp(C2) where: LE - latent evaporation Cl - 0.1 C2 - 0.6 C1,C2 are coefficients determined from test data. A regression on data generated from this equation, and the actual test data gave a correlation coefficient of -0.96 for early bloom alfalfa and -0.98 for full bloom alfalfa. Substitution of k = 0.1(LE) into this equation gave correlation coefficients of 0.9 (early bloom) and 0.87 (full bloom). Hill (1976) also used a similar equation of the form: MR = e(-kt)*exp(0.8) where: k - 0.0079(VPD) + 0.098 VPD - vapor pressure deficit Temperatures from 19-31 C, and relative humidities from 50-93% were observed during testing. Hill warns that major errors occured when there were significant radiation effects, presumably due to the crop surface temperature being substantially higher than the ambient temperatures. He also found the greatest accuracy when the vapor pressure deficit was near 30 millibars below saturation. Hayhoe and Jackson (1974) used a similar equation to predict the final moisture content after drying for n days: I: * -* Mn Mo exp ( a sumPE) 21 where M - moisture at the end of the nth M2 - initial moisture at cutting a - weighting factor sumPE - the summation of the potential energy for each day of drying (mm) day .Harris and Thaine (1975, Harris and May-Brown, 1976) used the technique of plotting the relative water content versus time, where the relative water content was described as: ch = (Wt-Wd)/(Ws-Wd)*100 To summerize, most previous modeling of the drying of alfalfa, and other hays, is based upon an exponentially decaying curve. There have been various environmental conditions used as the independent variable in the models including ambient temperature, latent evaporation, vapor pressure deficit and the potential energy. CHAPTER 4 EQUIPMENT Chamber and Air Handling In order to monitor the drying of alfalfa in a controlled environment equipment was needed that would maintain desired temperature and vapor pressure conditions. The equipment used consisted of a drying chamber, sample holders, an Aminco Aire unit, weight sensors, and‘ a microcomputer to control the environment and collect data. Figure 1 gives a plan view of the equipment. All control systems, including the microcomputer and the rewiring of the Aminco Aire Unit were built and maintained by Richard Byler and are discussed in Chapter 4 of his Dissertation (Byler, 1983). The drying chamber, shown in Figure 2, was constructed of plywood. The floor of the chamber was 1.6 meters long with a plenum at each end. Each plenum was 0.3 meters deep 22 23 nousaaououofiz muomcmm wuaumumaamu « Auuouam .n xn wawamuvv umusaaououoaz mnu can .DHCD mufi< oucfie< .Hmnemzo wcwmua .ucwaafisvm mcu wo 3mw> swam .H muawfim nonfimso w:H%uQ mumwaon maaamm llllL _ L r A f mHHoo wcwumwm :8: 32 855 numb nouns « an muuufiom mewMH< ecu uflwm 24 Auuoumm .n he wcaamupv u 3mH> unous wnfimco wcfizua wnu mo 30cm u .mumuswmcmuh we“ .muomcom =Cmgwm agalceua=.~ .N wuswam 25 which gave a one meter length in the area containing the samples. The cross section of both the sampling area and the plenums was 30 x 50 cm2. Each plenum stood on a base, also made of plywood, at a height convient for changing and viewing the samples. Since the sample chamber was not supported by the base it was reinforced with angle iron. Between each plenum and the sample chamber a metal honey-comb was placed to create a uniform airflow by creating a pressure drop. An air baffle was placed in the intake plenum to disperse the air coming into the plenum from the Aminco unit. The holes in the baffle were located between 30-75 degrees off center, on both sides. There were no holes that would allow airflow directly onto the samples. The air entered and exited the plenums through the plenum floor and was circulated by a fixed speed squirral-cage fan. Airflow was measured with a hot wire anemometer at various points around the sample holders. The velocity varied from 0.1 m/s to 0.7 m/s with an average airflow of 130 m /h. Henderson and Pabis (1962) found that airflow rate variations in this range had an "insignificant" effect on drying. The reason to maintain a constant airflow was two-fold; to decrease variability between trials, and to circulate the air through the Aminco unit to maintain proper drying conditions. The sample chamber top was constructed with an aluminum 26 frame which fit over the plywood floor and plenums. The frame was covered with propafilm-C to allow light on the sample, viewing, and to reduce heat buildup. The inlet plenum was insulated with rigid foam 0.25 m thick to reduce heat loss. The sample chamber was also lined with rigid foam insulation except for 0.4 m of the chamber to allow viewing of the samples. The plywood was covered with several coats of polyeurthene to reduce moisture adsorption by the wood. To change samples the entire film covered aluminum frame was slid up and off the drying chamber. The alfalfa stems were held vertically in the chamber on two sample holders. Each sample holder was made in three parts, an aluminum rod, a hard plastic cross piece and 4 binder clips as shown in Figure 3. The binder clips held the alfalfa stems and fit over the cross piece. Each clip Figure 3. Detail View of the Alfalfa Holders and Transducer Mount (by D. Sprott) 27 held 7-8 stems. The cross piece had a hole in the center offset from the line of binder clips which fit on the aluminum rod protuding through a hole in the chamber floor. The rod was held upright by the transducer outside the chamber. Each transducer was wired to the microcomputer so the sample weight could be recorded at any time. The sensors were mounted on a cement block to allow shock absorption The cement block rested on two pieces of iron which were in turn mounted 14 cm high on a cross piece on the chamber base. The block did not Sufficiently damp out vibrations and additional filtering was necessary. A strain gage amplifier was used to do this. To allow adjustable placement of the sensors two tracks were built in the cement block. The sensors were then mounted on a base which could he slid along the track (Figure 3). This arrangment allowed easy positioning of the sensors directly under the holes in the chamber floor. Additionally this allowed a method of changing the number of sample holders. An Aminco Aire Conditing unit, model J4-5460, was used to control the environment. The unit was designed to condition up to 8.5 cubic meters of air per minute in cabinets of less than 1.13 cubic meter volume. The sample chamber had a volume of 0.24 cubic meters and the fan forced about 2.17 cubic meters of air per minute through the 28 system. This unit drew in air (in this case from the exit plenum of the drying chamber) and cooled it to the water temperature by spraying the air with a fine mist of water. The water temperature was maintained at the desired level with either a refrigeration unit or an electric heating unit as needed. The air was then heated to the desired temperature in the next portion of the Aminco by an electric heater. The heating and refrigeration units were controlled by the microcomputer rather than the original controller. The air was then moved into the drying (sample) chamber by the fan. Two ducts were used to connect the Aminco to the sample chamber. They were each 1 m long with a 0.127 m diameter. The tube leading from the air conditioning unit to the sample chamber was insulated with 0.03 m of foil-faced fiberglass insulation. The duct exiting the sample chamber was not insulated since the Aminco cooled the air to the water temperature so retaining heat would serve no purpose. Once steady state conditions were reached in the equipment they were easily held there as the load produced by the drying alfalfa was not great. There was some heat loss to the environment from the sample chamber, the amount of which varied dependent upon ambient and drying temperatures. This was only a minor problem, however, and with contol feedback using the drying chamber temperature 29 the test conditions were well maintained. The greatest load on the system was due to the systematic cooling, saturating and reheating of the circulating air in the Aminco. Microccomputer and Control The microcomputer was used to maintain constant drying conditions and to record the sample mass, and the system temperatures. Extensive revisions were performed by Byler (1983) to allow microcomputer control of the Aminco Aire unit. The relationship of the water temperature, the dry bulb temperature and the resultant relative humdidty was given graphically in the manufacturers catalog. This chart was used to determine the necessary water bath temperature to give the desired relative humidity at any dry bulb temperature. A formula was derived from this data by Byler (1983): DB = 77.3673-16.9764ln(RH)-0.082136*WT*ln(RH) +73.6944/(RH)+1.377*wt where DB = the dry bulb temperature,‘Celsius RH = the relative humidity, percent WT = the water bath temperature,'Celsius The mean square error for this equation was reported as less than 0.2. The original bymetalic thermoregulators of the Aminco Aire unit were replaced with integrated circuit precision 30 temperature sensors (LM 135) for better control. These sensors were wired to the microcomputer to be used for control and data output. For trials 1 to 16 a temperature sensor placed near the exit of the Aminco unit controlled the air temperature. There was some difficulty in maintaining the desired chamber temperature and vapor pressure with this equipment. After trial 16 a sensor near the samples was used to help control the desired temperature. This modification enhanced control since the heat loss that occured before the drying air reached the samples did not need to be estimated. Six analog temperature sensors were used. The microcomputer converted the analog signal to a digital signal. Two of these sensors were placed in the Aminco Aire unit, one located in the water bath (#2) and one near the electric heater (#4). The other four sensors were in the drying chamber, #3 was near the inlet and the other three (#‘s 5,6,7) were near the drying samples, one between the samples and one behind each sample. The numbers refer to the computer control number. The microcomputer as has been mentioned, was used to collect data and to control the sample environment. Figure 4 diagrams the information flow in the data acquisition and process control equipment. A more detailed discussion of the microcomputer may be found in the PhD Thesis by Byler (1983). The data collection program recorded data on a Auwflhm .m may umuumHHou mama EEEE 31 _| Law 30......- 12.5.0 0055‘ H“ l :65 .33.:0 .o:x::9ua.u.2 5.3.390 fist can uwHHouucoo mmwuopm msu mo Emuwmwo xuoam nun—.3... .a—con 23:2.th anus—m 3:0 3.3:”. .25. O i=3n>li 2.3.30 :u‘.> .c wusmwm nun—.3:— 32..“ 2.2.3 32 cassette tape for storage as well as on a hardcopy output, if desired. The temperature information obtained was used to control the air conditioning unit. National Semiconductor precision temperature sensors LM335 (National Semiconductor, 1980) were used. The voltage output was linear and operated in a range of -10 “C to + 100 °C with a nonlinearity over that range of 0.3 “C (National Semiconducter, 1980, pp. 9-22). Actual nonlinearity in the ranges used was expected to be 0.15 “C or better. All the thermocouples were calibrated by placing them in the water of the Aminco, varying the temperature with the microcomputer and recording the thermometer reading and the readings from the thermocouples as given to the microcomputer. A laboratory mercury thermometer marked in 0.1 °C was used. Linear regression curves were fit to this data. The thermocouples used in the chamber to record the drying temperature were also calibrated for a range of air temperatures from 10-50 ‘C. The linearity of these calibrations was very good with a correlation coefficient, r, greater than 0.9999 for all sensors. Table 1 shows the temperature sensor coefficients from the equation A*(digital number) + B. The conversion factors for sensor 4 are not precise since this sensor was used to control the air heater in the Aminco the temperature was not recorded. 33 Table 1. Regression Coefficients for the Temperature Sensor Calibration (from Byler, 1983) * * Sensor A B Number ” 2 0.003932 -18.85 3 0.003994 -l9.00 4 0.322 -15.94 5 0.003921 -19.05 6 0.003921 -18.79 7 0.003986 -18.52 *In the equation: A*(computer number) + B Table 2. Regression Coefficients to Obtain the Sample Mass Used in Experiments 1-25 ‘ Trials Sensor A* 8* r Number 0 1-9 ' 0 0.2215776 103.57 0.9963 1 0.2159596 129.23 0.9998 10-17 0 0.221799 190.251 0.9995 1 0.2136949 220.570 0.9998 18-25 0 0.0162765 61.409 0.9999 1 0.0157122 114.210 0.9999 ¥ . f * In the equation: mass = A*(computer number) + B r is the regression coefficient 34 There were two weight transducers used for the testing. They were strain gage transducers, model number 4850 (two pound capacity) built by GSE Incorporated. The rated nonlinearity was 0.02% (0.2 grams) of full scale and 0.01% (0.1 grams) typically. The sensors were numbered 0 and 1 in the microcomputer. The two sensors were used in order that two samples could be dried simultaneously. The operating temperature range for these transducers was -17.78 to 93.33 %C with a temperature effect of 0.00044%/C on rated output. The transducers were located outside the chamber to reduce temperature and humidity effects. This was espically important for tests at high humidities where the moisture conditions may have affected the transducers. The transducers were designed to be insensitive to other than vertical mass loading. The rated error per inch off-center loading at half capacity was 0.004% (0.04 gram) of full scale. To calibrate these transducers known masses (0 to 200 grams) were loaded on the sample holder and the corresponding digital output noted (from the microcomputer). A linear regression of the form A*(digital number) + B was performed on the weight sensor calibrations. The coefficients A and B are shown in Table 2. The correlation coefficient, r, was greater than 0.996 in each case. The coefficients changed after trials 9 and 17 because of system 36 changes. For trials 1 to 17 calibrations were performed after each test and these coefficients were used in subsequent analysis. The r value for any one of these trials was greater than 0.9999. It was not necessary, however, to enter different regression coefficients for each trial. As Table 2 shows by grouping the trials in three sets the resultant correlation coefficient was good enough to allow average coefficients to be used for each set. After trial 16 an excitation source for the weight transducer was used so there was less noise in the data. Also after trial 16 the weights were calibrated by taking a final reading, with the samples in place, and a reading without samples (0 9) since the linearity had been well established. One further note- although the weight transducers measured force and calculations were made using units of mass all data was collected in the same place that the calibrations were run so gravitational affects did not influence data collection. Software Digital control software and data acquisition software for trials 17-33 were written by Byler (1983) and the programs may be found in the appendicies of that reference with a discussion found in Chapter 4. This program output an initial data set, a data set every 10 minutes for four 37 hours, and a data set every half hour thereafter. A data set included the weights of the sample, the waterbath and chamber temperatures and the date and time of the data set. More frequent readings were taken at the beginning of a test since alfalfa dries faster at a high moisture content. To reduce system noise each weight reading was actually the average value of six readings taken over a one minute interval. The program was stored on tape and had to be retrieved from the tape after a system malfunction. Repeated reading from the tape caused the tape to age and become distorted, therefore, the program had to be reentered occasionally. Because of these problems a data collection program was stored in the microcomputer memory. It allowed the operator to choose any time interval desired (minutes) for data collection. Ten minute intervals were then used in trials 17-33. To calibrate the transducers a known weight was placed on the sample holder and allowed to rest for 5 minutes after which a data set was recorded for that weight. Calibrations after each test. included 0, 50 and 100 gram weights. This method was used for trials 1-16. After trial 16 a quicker method of calibration was used as the transducer linearity had been well established. A data set was recorded just before the sample was removed from the chamber to be dried. The alfalfa was immediately 38 weighed and dried. The sample holders were replaced and another data set taken now without any alfalfa. This method gave a regression line (two points) and allowed minimum disturbance of the transducers. The sample was weighed and these two points were used to determine the regression line. The original calibrations (to check linearity as well as find the coefficients) used many weights, from 0 to 200 grams. CHAPTER 5 DATA COLLECTION Before the alfalfa samples were gathered for each test the apparatus was checked to confirm that it was functioning' at the desired temperature and vapor pressure. Alfalfa was cut within 3 miles of the laboratory on the Michigan State farms and transported to the laboratory in a plastic bag to minimize drying during transport. A visual estimate of the percent bloom was made when the sample was cut. The alfalfa was placed in the chamber immediately so no storage was necessary. The top 22 cm of the plant was used in order to decrease variability in the product. Fifty grams of alfalfa were clipped to each holder. A 50 gram sample generally required 20-30 stems. Once the samples were in the chamber the data collection program was initiated. The data collection program put the sample weight and 39 40 temperature data on tape and printed this information, if desired. Appendix B gives the hourly moisture content data sets for all trials. The program transfered the data to the tape after 8 data sets had been taken and listed the track and block where it was stored on the tape. When the trial was completed the cassette tape was taken to another microcomputer system to transfer the data to the university mainframe computer (Cyber 750). This was done in two steps. First the data was read from the tape and put on a floppy disk (using the program TDCONTRL, Byler, 1983). For trials 1-16 the following procedure was then used. The data was transfered to the mainframe using CDCFAST (Byler, 1983). Once in the mainframe the moisture content for each data set was calculated from the sample mass which was determined from the regression equations obtained after each test. Analysis was performed on the change in moisture content over time. After trial 16, system changes in both software and hardware were made which necessitated minor data analysis software changes. The new software calculated the moisture content of the sample before the data was sent to the mainframe (RPAK, Byler, 1983). Once these files were transfered to the mainframe the programs were simplified because the moisture content, on which all subsequent analysis was based, no longer needed to be calculated. 41 The drying temperatures from the 4 chamber temperature sensors were averaged for each test to check the equipment control. In trials 17-33 control was very good, in the eariler trials (before the system changes) the actual temperatures were often very different from those desired. The high temperatures were particularly difficult to maintain when the control sensor was not near the samples. In order to have an idea of the maturity of the sample the leaf-stem ratio and internodal lengths were determined from remaining alfalfa as soon as a test was initiated. For the leaf-stem ratio, a 25 9 sample was used (the top 22 cm of each stem only). The leaves were removed and weighed and the remaining stems were weighed. The ratio was then calculated. To measure internodal lengths, six 22 cm alfalfa stems were used. To get a representative value per trial the uppermost four internode lengths were averaged for the six stems. The leaf-stem ratio, internodal length and percent bloom and initial moisture contents for each trial are listed in Table 3. The trials generally lasted for 60 hours however this time length was occasionally shortened due to equipment malfunction. In order for the alfalfa attain a moisture content below 30 percent in the drying chamber the alfalfa was allowed to dry for 60 hours. Once the trial was 42 Table 3. Alfalfa Maturity and Time of Cutting for all Trials L/S INL BLM MC A,B Trial cut/year ratio (cm) (%) (%db) untreated ' 1 3/82 1.4 3.5 pre 4,33 4,25 2 3/82 1.6 2.2 pre 2,35 2,77 3 3/82 1.4 3.0 10 3,33 3,50 4 3/82 1.5 2.0 10 3,37 3,42 5 3/82 2.1 3.3 10 2,53 2,51 5 3/82 2.3 2.1 10 2.84 3.07 7 3/82 1.9 2-2 50 2.94 2.98 3 4/83 1.4 3.3 pre 3,34 3,23 9 4/82 2.1 3.1 pre 3.57 3.69 10 4/82 1.7 2-9 pre 3.42 3.00 11 4/82 1.9 3.0 pre 3.23 3.22 12 4/82 2.0 2.1 pre 3.38 3.23 13 4/82 1.8 2.5 pre 3.12 3.01 14 4/82 1.7 1.4 pre 2.85 2.96 15 4/82 2.2 1.6 pre 2.55 2.69 16 4/82 1.8 2.2 pre 2.73 2.70 17 1/83 1.2 2.3 pre 4.76 5.10 18 1/83 1.0 2.6 pre 7.25 6.24 19 1/83 1.0 3.2 pre 5.82 6.18 20 1/83 0.9 3.8 pre 6.02 5.98 21 1/83 1.1 3.8 pre 6.02 5.98 22 1/83 1.4 4.1 pre 2.75 3.64 23 1/83 1.3 na pre 1.27 1.19 24 1/83 NA na 5 3.17 2.99 25 1/83 1.6 3.3 pre 3.08 3.32 treated 25 2/83 1.7 3.32 pre 3,55 3,72 27 2/83 0.9 na pre 4.08 4.34 28 3/83 NA na pre 2.50 3.19 29 3/83 NA na pre 4,29 4,33 30 3/83 NA na pre 4,43 3,13 31 3/83 NA na pre 3,50 3,10 32 3/83 NA na pre 3,07 3,04 33 3/83 NA na pre 3,03 3,42 INL - average of top 4 internodal lengths (6 stems/trial) na - information not available BLM - bloom VP - vapor pressure L/S ratio - leaf/stem ratio MC - initial moisture content (dry basis) 43 completed the samples were removed from the chamber and weighed. They were placed in paper bags and oven dried at 103 ”C for at least 24 hours. Weighing the oven dried sample gave the final dry weights to be used in moisture content calculations. Trials 1-16 were performed in late summer and fall of 1982 while trials 17-33 were performed during the summer of 1983. In trials 26-33 both samples were sprayed with a 2.76% potassium carbonate solution. With this conditioning treatment, two trials were performed at each of 4 drying conditions. Samples in trials 1-25 were not treated, except in three cases where sample B was sprayed with a solution to speed drying. The potassium carbonate solution was used in trials 2 and 16 while a commercially available drying solution was used in trial 6. CHAPTER 6 ANALYSIS The calculated change in moisture content, with time, was analyzed to determine variations in drying rate for different environments. Three methods were used. The first method calculated a drying constant which was used as a measure of the drying rate. The second method of analysis required the calculation of the mass transfer coefficient. The third involved analysis of the coefficients of a best fit exponential curve for each trial. With the analysis of the drying constant and mass transfer coefficient each drying curve was divided into three sections dependent upon moisture content. A drying .constant and mass transfer coefficient were determined for moistures greater than 2 (dry basis), between 2 and 1, and for a moisture content less than 1. This was done to separate the drying process into a fast drying period (before stomatal closure), a 44 45 medium drying period (closed stomata), and a slow drying period (closed stomata, high resistance to water loss). It was beyond the scope of this study to determine at what moisture contents these periods might occur so the cutoff points are somewhat arbitrary. Drying Constant The first procedure (drying constant analysis) has been used and described by Priepke and Bruhn (1970). The drying curve was assumed to fit an exponential model where the drying constant was the slope of the natural log of the moisture ratio as a function of time. The drying constant, k, was determined as: k = -(1/t)*ln((M-Me)/(Mo-Me)) length of time interval(h) initial moisture content for the interval (db) M a final moisture content for the interval (db) Me a equilibrium moisture content (db) where: t Mo This equation is often written in exponential form: MR . e-kt where: MR - the moisture ratio = (M-Me)/(Mo-Me) A drying constant, k, was calculated for each data collection interval and an average value determined for the 3 ranges previously mentioned. The drying constants for each range were then compared using the SPSS statictical package (Nie et al., 1975) with temperature and vapor 46 pressure as the independent variables. In order to calculate the drying constant an equilibrium moisture content (Me) had to be assumed. Values for the Me of alfalfa dried under various environmental conditions were obtained from Bakker-Arkema (1962). Mass Transfer Coefficient The second method of analysis was based on the calculation of the mass transfer coefficient. This coefficient was calculated using the equation (Holman, 1976): m = HD(Cp- Ca) 7 where: m a diffuse mass flux H a mass transfer coefficient C = concentration of moisture in the plant C a concentration of moisture in the a air = initial moisture content for the M1 time interval M2 = final moisture content for the time interval t1 = initial time t2 = final time The laws of Dalton and Avogadro state that the molal composition of a mixture is proportional to the distribution of partial pressures. This gives the molal humidity, generally labeled f, as the mass of water vapor in mols per 47 one mol of air: f s Pv/Pa Pv - vapor pressure Pa a air pressure a Pt - Pv (total pressure - vapor pressure) So the concentration of moisture in the air was determined from the vapor pressure and the equation: Ca: Pv/(101.3 - Pv) where: Pv = vapor pressure Table 4 shows the Ca values for each trial. To determine the value of Cp the volume of alfalfa at various moisture contents was needed. Samples of 50 grams of alfalfa at moisture contents of 0.2 to 4.88 were immersed in water to determine their volumes (Table 5). A linear regression was used giving the equation: Vol = 17.67*MC + 27.52 where: MC = moisture content The r2 value was 0.83 for this regression. The concentration of moisture in the plant could then be determined using the equation: Cp= MC*DM/Vol where: MC = moisture content DM - dry mass of the alfalfa Vol = the volume of alfalfa at MC The values of the drying constant, k, the mass transfer coefficient, HD, and the correlation coefficient were 48 Table 4. Temperature Settings, Drying temperature, Vapor Pressure and Concentration of Moisture in the Drying Air for all Trials T c wa DB TRIAL (‘C) 9v a (kPa) kPa water (‘C) ('C) kPa air 1 26 1.08 0.011 5.0 35 2 26 1.08 0.011 5.0 35 3 30 1.91 0.019 15.0 35 4 25 1.87 0.019 15.0 27 5 26 1.88 0.019 15.0 31 6 29 1.92 0.019 15.0 40 7 23 1.26 0.013: 8.0 25 8 24 1.55 0.016 12.0 25 9 25 1.94 0.020 15.5 25 10 21 1.88 0.019 18.0 25 11 26 2.25 0.023 18.0 26 12 21 1.36 0.014 12.5 30 13 30 2.12 0.021 17.0 40 14 22 1.24 0.012 7.5 25 15 17 0.91 0.009 2.0 18 16 27 1.67 0.017 12.5 40 17 26 1.68 0.017 25.0 40 18 40 4.79 0.050 30.0 40 19 33 2.55 0.026 19.3 33 20 33 1.90 0.019 14.0 33 21 33 1.55 0.016 10.0 33 22 33 1.30 0.013 6.5 33 23 38 1.58 0.016 10.0 38 24 38 1.88 0.019 13.5 38 25 38 1.13 0.011 3.0 38 ggeatEd sampléfs 5.76 0.060 14.5 35 27 30 4.24 0.044 10.0 30 28 40 7.37 0.078 19.0 40 29 25 3.17 0.032 6.0 25 30 25 3.17 0.032 6.0 25 31 30 4.24 0.044 10.0 30 32 35 5.76 0.060 14.5 35 33 40 7.03 0.075 19.0 40 C0- concentration of moisture in the air P - vapor pressure, WB - wet bulb tempera DB - dry bulb temperature t - temperature, dry bulb kPa (kilo Pascals) ture 49 Table 5. Volume of Water Displaced by Alfalfa at Different Moisture Contents Sample Moisture Mass Content Volume (9) (db) (ml) # m 51.9 4.88 100 25.68 2.01 60 14.70 0.72 50 9.72 0.14 30 9.87 0.15 30 9.61 0.11 25 10.46 0.20 35 10.27 0.20 25 50 determined using the program CALCKHD (Appendix A). This program used the moisture content and time for each data set to calculate the drying rate (as described by k and HD) for each time interval. The output for this program is seen in Appendix C. Exponential Curve Fitting The third method of analysis was exponential curve fitting. The statistical software package BMDP, program P3R (nonlinear regression), (Jennrich, 1981) was used to fit an .exponential curve to the data. This program, BMDP3R, is given in Appendix A. A model utilizing two exponents was used: MC = P1*exp(P2*t) + P3*exp(P4*t) + P5 \ where: Px = coefficients (x=l,2,3,4) P5 = equilibrium moisture content (from Bakker-Arkema, 1962) t = time This equation was similar to the equation used to calculate the drying constant where: (-kt) (-kt) 140- M = (M0- Me)e 4’ Me = (P1 + P3)e + P5 This comparison assumes -k = P2 = P4 (the exponents). Models utilizing 3 and 4 exponents were also used. These models did not fit the data better than the 2 exponential model so no further analysis was done. CHAPTER 7 RESULTS AND DISCUSSION The drying rate of alfalfa was studied by analyzing the change in mass, with time, of cut alfalfa plants for samples dried under different temperatures and vapor pressures. All mass loss was assumed to be due to loss of water from the plant as the plant dried. As was mentioned in the literature review, Greenhill (1959) reported a 7% dry matter loss due to continued respiration during drying but suggested that continued photosynthesis may offset this loss. A graph of the change in moisture content, with time, is given for each trial in Appendix D. These curves show the exponential nature of the drying alfalfa plant. Table 6 shows values for the drying constants, k, and the mass transfer coefficients, HD, for trials with the same vapor pressure conditions (different temperatures). The values given are for transducer A from 3 trials at each of 51 Table 6. Drying Constants and Mass Transfer Coefficients for Alfalfa Trials at Two Vapor Pressures and Different Temperatures (Transducer A) Vapor Pressure = 1.25 kPa Vapor Pressure 2 1.88 kPa Trial 7 22 25 10 5 24 Temp 23 33 38 21 26 38 INL 2.2 4.1 3.3 2.9 3.3 na L/S 1.9 1.4 1.6 1.7 2.1 na Bloom 50 pre 10 pre 10 5 k1 0.089 0.137 0.201 0.041 0.053 0.190 k2 0.067 0.120 0.201 0.035 0.056 0.222 k3 0.033 0.049 0.044 na 0.041 0.072 HD1 0.515 1.020 1.203 0.339 0.271 1.241 HDZ 0.278 0.607 0.909 0.187 0.209 1.058 HD3 0.097 0.146 0.161 na 0.097 0.294 INL - internodal length (cm) L/S - leaf/stem ratio Bloom - percent na - information not available k - drying constant k1 - average k value for moisture contents greater than 2.0 k2 - average k value for 2.0 > moisture > 1.0 k3 - average k value for moisture < 1.0 mass transfer coefficient value for moisture > 2.0 value for 2.0 > moisture > 1.0 value for moisutre < 1.0 HD - HD1 - average H HD2 - average HD HD3 - average HD 53 two vapor pressures. Generally k and HD increase with increasing temperatures, indicating faster drying at higher temperatures. A linear regression, with temperature as the independent variable and drying constant as the dependent variable gave an R2 value greater than 0.82 for each of the three moisture ranges. Regression with HD as the dependent 2 value less than 0.65 for all three variable gave an R moisture ranges. Table 6 also shows that there was a decrease in the magnitude of k and HD at lower moisture contents. This would be expected due to the greater resistance of the plant to water loss at lower moisture contents. Table 7 shows the k and HD values for trials at different vapor pressures and two temperatures. In this table no pattern seems to exist between the drying constant nor the mass transfer coefficient with increasing vapor pressure. As before, however, there was a decrease in the magnitude of k and HD with decreasing moisture contents. These two results suggest that drying time was directly related to the ambient temperature while not directly affected by the vapor pressure in the range used in these trials. It is possible that other factors, including maturity, weathering and water stress during growing, may have masked the effect of the vapor pressure on drying. The change in moisture content, over time, was plotted 54 Table 7. Drying Constants and Mass Transfer Coefficients for Alfalfa Trials at Two Temperatures and Different Vapor Pressures (Transducer A) Temperature . 26 C Trial 1 l7 5 ll VP 1.08 1.68 1.88 2.25 INL 3.5 2.3 3.3 3.0 L/S 1.4 1.2 2.1 1.9 Bloom pre pre 10 pre k1 0.139 0.257 0.053 0.042 k2 0.052 0.164 0.056 0.050 k3 0.039 0.072 0.041 0.050 HD1 1.272 2.577 0.271 0.282 HD2 0.354 1.027 0.209 0.210 HD3 0.185 0.310 0.097 0.159 Temperature s 33 C Trial 22 21 20 19 VP 1.30 1.55 1.90 2.55 INL 4.1 3.8 3.8 3.2 L/S 1.4 1.1 0.9 1.0 Bloom pre pre pre pre k1 0.137 0.114 0.102 0.209 k2 0.120 0.078 0.051 0.026 k3 0.049 0.030 0.026 0.062 HDl 1.020 1.334 1.242 2.713 HDZ 0.607 0.544 0.358 0.961 HD3 0.146 0.155 0.137 0.372 INL - internodal length (cm) L/S - leaf/stem ratio Bloom - percent k - drying constant k1 - average k value for moisture contents greater than 2.0' k2 - average k value for 2.0 > moisture > 1.0 k3 - average k value for moisture < 1.0 HD - mass transfer coefficient HD1 - average H value for moisture > 2.0 HDZ - average HD value for 2.0 > moisture > 1.0 HD3 - average HD value for moisutre < 1.0 56 for the conditions of Table 6 in Figure 5. This figure shows drying occurs more rapidly as temperatures increased for a vapor pressure of 1.25 or 1.88 kPa. The change in moisture over time at four vapor pressures for temperatures of 26 and 33 °C are given in Figure 6 (refer to Table 7). The drying rate due to changes in vapor pressure were not consistant. At 26 “C the fastest drying was for 1.68 kPa. The alfalfa dried with vapor pressures of 1.30 kPa and 2.55 kPa dried at a similar rate which was faster than the alfalfa dried at 1.55 kPa or 1.90 kPa. The environment with the greatest vapor pressure (2.25 kPa) had the slowest drying at 26 “C but the fastest at 33 “C. This apparant interaction deserves further study. Crop Maturity Information collected on the leaf-stem ratio and percent bloom indicate that the drying rate may be affected by the physical properties of the plant when it was cut. The moisture content (of the mature alfalfa plant was a product of many variables including the soil fertility and acidity, available moisture during growing, wind and other environmental conditions. In general at high percent bloom the k values are lower than for trials using less mature alfalfa. For example, trials 7 and 14 were performed under similar conditions except that trial 14 was in pre-bloom and fi 57 VHPOR PRESSURE : 1.25 kPa 4-00 - l NTENT (DB) C3 . 00 uf.’ (KN :3 p. .‘2 23 °0 .00 TIME IN HOURS 8 .45“) g VAPOR PRESSURE = 1-88 1.1». "' 0 ~8 A 345: z- 0 m"1 m 38 “C p. 2: C3 00 m9 (KN-1 :3 p. 8 252’ ‘13-‘30 ~00 50-0 00 715.00 TIHE IN HOURS Figure 5. Plot of the Change in Moisture Content of the Actual Data Over Time for Different Temperatures at Two Vapor Pressures g 58 a” TEHPERHTURE = 26% C3 U - A 1.08 kPa i-g 0 1-68 kPa EV; m [-88 kPa 2E 22'1255 kPa C3 (now ufi- age: :3 p_. 2 8'4 - ~ r31.00 25.00 55.00 75.00 TIHE IN HOURS 8 505° TEHPERHTURE = 33°C H A 1.30 kPa P8 0 1055 kPa E; m 1 ~90 kPa E Q 2'55 kPa CD 00 In? tum“ :3 p. 8 as: " c11.00 25.00 55.00 75.00 TIME IN HOURS Figure 6. Plot of the Change in Moisture Content of the Actual Data Over Time for Different Vapor Pressures at Two Temperatures 59 trial 7 was at 50% bloom. The K values for trial 7 are lower than those for trial 14. This was not unexpected since the older plant has a lower moisture content when cut . and will not release water as easily as a plant at a higher moisture. This effect was also observed in trials on grasses, by Menzies and O'Callaghan (1971). Menzies and O'Callaghan stated that generally the maturity of the grass did not have much effect on the drying characteristics of the grass. Wieghart (1979) observed, however, that older alfalfa dried faster than the less mature alfalfa when cut. This author suggested the faster drying was due to a greater amount of weathering sustained by the older plant's cuticle, and a subsequent decrease in cuticle effectiveness. Higher leaf-stem ratios generally yield lower k values, which would be expected as this is another test of the level of maturity. The internodal lengths did not noticably affect the drying constant. Table 6 shows that there may be a correlation in these three parameters. As the percent bloom increased the leaf/stem ratio also increased and the internodal length decreased. Crop maturity was not a variable of major concern in this study. However, due to the possible affect of maturity on drying an attempt was made to keep the variabilities discussed here to a minimum. 6O Drying Constant The mean drying constant, k, for 4 temperatures, 4 vapor pressures, and 3 moisture content levels are given in Table 8. These values are separated into 2 groups, those for transducer A and those values for samples dried on transducer B. This was done to check that these two transducers gave the same results. There was a general pattern of increasing drying constant magnitude (as shown by the means) with increasing temperatures. As mentioned previously this relationship was confirmed by a linear regression. The values of k at temperatures from 31 to 35 ’8 are lower than the values in the previous group (ZS-30°C). This was not expected and may be due to trials being performed in different years. Trials in the 31-35 °C temperature range were run early in 1983 (a total of 4 trials) and the trials performed at 26-30 %C and 1.31-l.8 kPa were run in 1982 (a total of 8 trials). The differences in the material, as well as the variability in the '83 trials (less repetitions) may have Caused the lower k values at higher temperatures. The analysis of variance shows that the differences in the transducers (or sides of the chamber) are negligable (Table 8). This was particularly true for sample moisuures greater than 1.0 (a=0.82). For sample moisture less than 61 Table 8. Mean Drying Constants for 3 Moisture Content Intervals with Different Temperatures and Vapor Pressures, and the Results of an analysis of variance on this data Transducer A Transducer B Temperature kl k2 k3 k1 k2 k3 17-25 C .083 .067 .042 .085 .071 .044 26-30 C .164 .112 .051 .148 .104 .059 31-35 C .140 .094 .042 .146 .085 .047 36-40 C .209 .192 .056 .243 .224 .065 Vapor Pressure 0.7-1.3 kPa .130 .103 .040 .132 .108 .046 l.31-1.8 kPa .162 .133 .053 .160 .123 .062 1.81-2.0 kPa .130 .080 .042 .093 .081 .046 2.01-4.5 kPa .207 .133 .062 .221 .154 .061 overall mean .141 .108 .048 .142 .111 .053 Signifiance level from the analysis of variance k1 k2 k3 Temperature .021 .001 .031 Vapor Pressure .124 .106 .008 Transducer .954 .820 .140 Interactions ' kl k2 k3 temperature-vapor pressure .722 .270 .001 temperature-transducer ,92} 726. 720 vapor pressure-transducer .978 “911 '43 temperature-vapor pressure-transducer .992 '979 I823 R k - drying constant kl - average k value for moisture > 2.0 k2 - average k value for 2.0 > moisture > 1.0 R3 - average k value for moisture < 1.0 62 1.0 the possible difference in the transducers has a value of a-0.l4 (86% confidence that the means of the two samples are not equal). Of greater interest was the effect of temperature, vapor pressure and their interactions. For moisture contents greater than 2.0 the analysis shows a 98% confidence that the means are not equal for the effect of temperature and a 88% confidence that the vapor pressure effect was significantly different. For moisture contents from 1.0 to 2.0 the level of confidence increases to 99.9% for temperature and 99.4% for vapor pressure. The interaction was not significant. For moisture contents less than 1.0 temperature and vapor pressure effects had a significant interaction (99.9% confidence). This may be due to the increased resistance to further water loss by the plant at low moisture contents. When the plant is severely stressed the internal and external environments, particularly the difference in moisture, becomes a governing factor in further water loss. This analysis indicates that initial drying was less dependent on environmental conditions which would be expected when stomata are open. Then as stomata close and plant resistance to water loss increases the environmental conditions, singularly and combined, play a larger role in drying. Regression analysis was used to determine the closeness 63 of fit of the model based on k to the actual data. The predicted moisture content was calculated using the average drying constant for each trial. The predicted values were compared to the actual values. The coefficients of determination (r2) are given in Appendix C, with the k and HD values. The fit of the theoretical curves were very good. Of the 99 values (33 trials with 3 moisture intervals per trial) there were 6 cases where r2 was less than 0.9 and 3 cases in which r2 was less than 0.8. For the interval of moisture content greater than 2.0 the lowest correlation value was r2 = 0.95. The interval with moisture contents less than 1.0 gave the poorest fit possibly due to the equipment noise having a significant effect upon the data and causing inaccurate results. Mass Transfer Coefficient Table 9 gives the results of an analysis of variance on the mass transfer coefficient with the mean HD values for the 4 temperatures, 4 vapor pressures and 2 transducer levels. The mean HD increases with increasing temperature which was also shown in Table 7. As was the case with k an increasing vapor pressure does not show an increase in the mean mass transfer coefficient values. However HD does increase with increasing vapor pressure, except in the l.81-2.0 kPa range. Trials in all vapor pressure ranges 64 Table 9. Mean Mass Transfer Coefficient for Three Moisture Content Intervals with Different Temperatures and Vapor Pressures, and the Results of an analyais of variance on this data Transducer A Transducer B Temperature HD1 HDZ HD3 HD1 HD2 HD3 17-25 C 0.524 0.304 0.125 0.533 0.314 0.117 26-30 C 1.121 0.528 0.188 1.176 0.526 0.195 31-35 C 1.577 0.617 0.203 1.628 0.563 0.224 36-40 C 1.854 1.092 1.007 2.402 1.437 0.769 Va or Pressure 0.7-1.3 kPa 0.836 0.467 0.134 0.875 0.499 0.140 0.31-1.8 kPa 1.177 0.659 0.191 1.263 0.628 0.236 1.81-2.0 kPa 0.732 0.388 0.162 0.702 0.389 0.175 2.01-4.5 kPa 2.308 0.924 1.045 2.738 1.203 04123 Overall mean 1.120 0.561 0.308 1.247 0.621 0.284 Signifiance level from the analysis of variance HD1 HD2 HD3 Temperature .001 .001 .001 Vapor Pressure .001 .001 .001 Transducer .631 .555 .001 Interactions Hnl H02 HD3 Temperature-vapor pressure .001 .005 .001 Temperature-transducer .779 .292 .001 vapor pressure-transducer .821 .532 .001 temperature-vapor pressure-transducer .539 .290 .001 HD - mass transfer coefficient - average value for moisture > 2 - average value for 2.0 > moisture > 1.0 1 'HD3 - average value for moisture < 1.0 65 were performed during both 1982 and 1983 so it was unclear what might have caused this, other than natural sampling variability. Similar to the drying constant analysis the effect of the trandsucer was not significant until the low moisture content levels (MC<1.0) are reached. Unlike the analysis of k, however, temperature and vapor pressure have a significant (confidence 99.9%) affect on HD at all moisture contents. The interaction of temperature and vapor pressure was significant for all three moisture content groups also. This result reflects the key point of the mass transfer coefficient analysis. This analysis is based on the moisture in the plant and in the air and consequently the air vapor pressure interaction with drying temperature is important. In modeling the drying of alfalfa with the mass transfer coefficient (non-exponential in nature) the drying temperature and vapor pressure are important at all moisture contents. This may be because this model does not assume the natural exponential nature of drying alfalfa. The correlation coefficients for the mass transfer coefficinets are in Appendix C as previously mentioned. Of the 99 regression coefficients there were 14 values with 2 values of r < 0.9 and 7 values less than 0.8. The fit was better at the highest moisture content interval (greater 66 than 2.0) as with k. Again the poorest fit was seen at the lowest moisture contents, probably due to system noise causing inaccurate results. Exponential Curve Fitting The results of the exponential curve fitting were inconclusive. Figure 7 shows the curves generated from this analysis for trials 1 and 2. The curves fit the data quite well but the problem of analysis lies in the lack of a means of comparison of the equations. Table 10 lists the calculated coefficients in order of increasing temperatures and Table 11 lists these coefficients in order of increasing vapor pressure. There are no definite trends among the variables with the change in environmental conditions. In fact trials with the same temperature and vapor pressure, such as trials 1 and 2, had markedly different coefficients. This method of analysis modeled alfalfa drying very well but could not be used to compare trials. Because there was no noticable patterns in the coefficients this analysis was not performed on trials 10-33. Chemical Treatment The eight (8) trials performed with potassium carbonate C U. 3.75 I p Q A Y 2.‘ 0 U i f 4 -¢ ‘0‘- Lo - 3. 2.25 V Z 5 T 1." J a 07- 00 ‘0000.0000.0000’0000‘0000’0000'0 000000007: 0 0 0 0 0 O 0 0 0 0 . . 0 0 0 0 0 O 0 . 000000 . 0 0 000000 ‘000000000’000090000’0000’0000‘ :20: :70. {‘0' . A a g. - 9‘ a... fit... . ' C Time (hours) 1 - transducer A 0’0000‘0000‘0000'00009000090000°0 . O of 0 0' 0 0' 0 0' 0 0 0 o 0 ' 0 0 ' 0 O P 0 0 . 0 0 0 0 0 ‘3 0 0 ' 0 0 ' 0 0 . 0 o -0 0 0 .9 0 . 000' . ‘00000: : g , . ”0000000000000 . O O 0.000090000’0000’0000°0000°0000°0 15 0E 7' 0 " 9' 5 Time Lhours) Trial 2 — transducer A 0 - actual (observed) data point P - predicted data point * - observed and predicted data points coincide Figure 7. l and 2 67 0" AOL. -‘VI>0U ‘ 0'0 avg-040mm 8 Trial 10" I D u... 0..OOO...00......OOO0.0000......O O O .- O 0' . 0° 0 0‘ g o.“ . o o . 0 . . . 00 0 0 ' . . 0 0 C0 . 0 9 . 0 P . 0 ;. . . 00‘ o 00 o . 0000‘ . ‘00000= . 0 3000000000 5. . 0’000090000’0000°0000‘0000‘0000‘0 120‘ 370; :20. ‘ .. Q: a g ‘ .p. O ..-0 .-.U 3.. . Time, (hours) 1 - transducer B 0‘0000‘0000‘0000‘0000‘0000’0000’0 0 O O O: O of . 0' 0 O- O o; o 0' 0 0° 0 0' 0 0‘ 0 O. o 0' 0 . S . 0 O ‘ 0 0 O 0 0 " . C0: 0 . 0000 , , 3000000000000000000000 . O o o 0‘000000000‘0000900009000090000’0 IL 0: Tims.(hours) Trial 2 - transducernB Example of the Two Exponential Curve Fit (BMDP) for Trials 68 Table 10. Coefficients of the Two Exponential Curve, Fit to the Data, in Order of Increasing Temperature Trial t VP P1 P2 P3 P4 P5 j Transducer A 7 23 1.3 2.61 -0.07 0.06 -0.00 0.10 8 24 1.6 5.00 -0.04 -1. 82 -0.03 0.11 4 25 1.9 2.09 -0.10 0.89 -0.03 0.13 9 26 1.9 0.56 -l.83 2. 78 -0.04 0.12 2 26 1.0 2.42 -0.11 0. 20 -0.01 0.08 1 26 1.1 1.90 -0.26 1.87 -0.04 0.08 5 26 1.9 2.48 -0. 06 0.10 -0.00 0.12 6 29 1.9 0.48 -0. 84 2.11 -0.08 0.11 3 30 1.9 3.14 -0. 08 0.11 -0.00 0.09 ansducer 7 23 1.3 2.66 -0.07 0.12 -0.00 0.10 8 24 1.6 3.13 -0.05 0.04 -0.00 0. 09 4 25 1.9 0.44 -l.79 2.88 -0.06 0.13 9 25 1.9 0.58. -2.76 2. 95 -0.05 0.12 2 26 1.0 2.29 -0.82 0. 48 -0.06 0. 08 1 26 1.1 2.36 -0.17 1. 39 -0. 03 0. 08 5 26 1.9 2.45 -2.67 0.05 -0. 00 0.12 6 29 1.9 2.47 -0.29 0.42 -0. 06 0.09 3 30 1.9 3.18 -0.09 0.19 -0. 00 0.10 t - temperature, dry bulb Vp - vapor pressure, kPa Px - Coefficients in the equation: Plexp(P2*time) + P3exp(P4*time) + P5 P5 - the equilibrium moisture content, fixed, from Bakker-Arkema (1962) Table 11. 69 Coefficients of the Two Exponential Curve, Fit to the Data, in Order of Increasing Vapor Pressure Trial VP t P1 P2 P3 P4 P5 Transducer A 2 1.0 26 2.42 -0.11 0.20 -0.01 0.08 1 1.1 26 1.90 -0.26 1.87 -0.04 0.08 7 1.3 23 2.61 -0.07 0.06 -0.00 0.10 8 1.6 24 5.00 -0.04 -1.82 -0.03 0.11 4 1.9 25 2.09 -0.10 0.89 -0.03 0.13 5 1.9 26 2.48 -0.06 0.10 -0.00 0.12 9 1.9 25 0.56 -1.83 2.79 -0.04 0.12 6 1.9 29 0.48 -0.84 2.11 -0.08 0.11 3 1.9 30 3.14 -0.08 0.11 -0.00 0.09 Transducer B 2 1.0 25 2.29 -0.82 0.48 -0.06 0.08 l 1.1 26 2.36 -0.17 1.39 -0.03 0.08 7 1.3 23 2.66 -0.07 0.12 -0.00 0.10 8 1.6 24 3.13 -0.05 0.04 -0.00 0.09 4 1.9 24 0.44 -1.79 2.88 -0.06 0.13 5 1.9 26 2.45 -0.06 0.05 -0.00 0.12 9 1.9. 25 0.58 -2.67 2.95 -0.05 0.12 6 1.9 29 2.47 -0.29 0.42 -0.05 0.09 3 1.9 30‘ 3.18 -0.09 0.19 -0.00 0.10 t - temperature, dry bulb Vp - vapor pressure, kPa Px - Coefficients in the equation: Plexp(P2*time) + P3exp(P4*time) + P5 P5 - the equilibrium moisture content, fixed, from Bakker-Arkema (1962) 70 solution sprayed on both alfalfa samples were compared using the drying constant and mass transfer coefficient analysis described previously. The mean values of k and HD, for the three levels of moisture, are shown in Table 12. As with the untreated samples the magnitude of both k and HD decrease with decreasing moisture content of the plant or with increasing drying temperatures. The magnitudes of the average drying constant or mass transfer coefficient at each level of moisture content are, in every case, greater for the treated product than the untreated alfalfa, indicating faster drying with the potassium carbonate solution. The Chemically treated alfalfa was dried under conditions of higher vapor pressures than the untreated alfalfa so a direct comparison was not possible. Linear regression was used to analyze each moisture content interval for both the potassium carbonate treatment and the control. Drying temperature was the independent variable and the drying constant the dependent variable. The linear regression fit the treated samples best with r2 greater than 0.84 for moisture contents above 1.0. The coefficient of determination for the treated samples at moisture contents below 1.0 and for all the control sample moisture intervals was less than 0.41. In spite of this poor fit the regression lines indicate the effect of the potassium carbonate solution to speed drying. The slope of 71 Table 12. Mean Drying Constants and 'Mass Transfer Coefficients for Three Moisture Content Intervals at Four Environmental Conditions with Treated Alfalfa Condition A B C D mean k1 0.116 0.147 0.236 0.422 0.230 k2 0.108 0.118 0.241 0.399 0.216 k3 0.112 0.072 0.158 0.154 0.108 mean 0.112 0.112 0.212 0.325 HDl 0.922 1.234 1.699 2.923 1.694 H02 0.580 0.676 1.252 1.923 . 1.108 HD3 0.485 0.315 0.466 0.344 0.402 mean 0.662 0.742 1.139 1.730 Conditions A - temperature = 25 C, vapor pressure - 3.17 kPa B - temperature = 30 C, vapor pressure - 4.24 kPa C - temperature = 35 C, vapor pressure = 5.76 kPa D - temperature - 40 C, vapor pressure = 7.03-7.37 kPa k - drying constant kl - average k value for moisture > 2.0 k2 - average k value for 2.0 > moisture > 1.0 k3 - average k value for moisture < 1.0 - mass transfer coefficient H H31 - average value for moisture > 2 HD2 - average value for 2.0 > moisture > 1.0 HDB - average value for moisture < 1.0 72 the line for untreated. samples with moisture contents greater than 2.0 was 0.007. The slope of the line for the treated samples, with moisture content greater than 2.0 was 0.020. This is a 186% increase in the slope of this regression line. The second moisture content interval (1.0-2.0) also had an increased drying constant as indicated by the increased slope for the treated samples. The slope of the line for the treated samples was 0.020 and the slope of the line determined by the untreated samples was 0.006 (233% increase). The final moisture content interval (less than 1.0) was also affected by the chemical treatment, the slope of the regression line increasing from 0.0006 to 0.006 due to chemical treatment. A graph of these regression lines is given in Figure 8. The effect of the chemical treatment is greater at higher drying temperatures as indicated by the increased slope of these regression lines. At drying temperatures above 35 ’C the effect of potassium carbonate is significant. At 25 ’C the treated samples are observed to have lower drying constants than the untreated samples. This is not actually the case, it is due to the linear model. Exponential models were not used in this analysis as they did not increase the Coefficient of determination. There were two cases in which the sample on thansducer B was sprayed with a potassium carbonate solution (trials 2 73 In Control sample, I - 1 0 Control sample, I - 2 o ‘ Control sample, I - 3 f O Treated sample, I - 1 :5q .0 Treated sample, I - 2 f, Treated sample, I-l O T P-c:q I: c J p. 03 . 2c: ‘ 1. o9: _ Uo" '4 . (.3 z -- - ha m— / O CDCD‘A” .7 + _+_ + A O 9 I 925. 00 30 .00 35 .00 4b .00 45 .00 TEMPERHTURE C I - 1: moisture contents greater than 2.0 (dry basis) I - 2: moisture contents greater than 1.0 and less than 2.0 (dry basis) I - 3 moisture contents less than 1.0 (dry basis) Figure 8. Linear Regression Curves for Drying Constant Based on the Drying Temperature for Treated and Control Samples at Three Moisture Content Intervals 74 and 16, trial 6 used a similar, commercially available solution). These treatments significantly increased both the k and HD values as shown in Table 13. The. increased drying rate was higher at the higher moisture contents. The treated samples achieved lower moisture contents much faster than the untreated and continued to dry somewhat faster. This may been seen graphically in Appendix D. The curve fitting package, BMDP-P3R, was also utilized, but in a slightly different manner. Two trials, at each of 4 drying conditions, with two samples per trial were performed giving a total of 4 drying curves for each drying condition (temperature and vapor pressure set). The P3R package was used on these four curves to give a 2 exponential equation of best fit. The BMDP 2 exponential equations were as shown: At 25 C: 3.69exp(-0.12t) + 0.18exp(-4.06t) + 0.18 At 30 C: 2.40exp(-0.10t) + 0.98exp(-0.36t) + 0.30 At 35 C: 0.23exp(-0.00t) + 3.19exp(-0.22t) + 0.06 At 40 C: 2.75exp(-0.44t) + 0.25exp(-0.01t) - 0.07 A plot of these four equations is given in Figure 9 and 'indicates that as the temperature increases the initial slope of the drying curves are steeper. The final moisture content was lower at the higher temperatures, except for the sample dried at 25°C. Note that these trials were performed with a relative humidity of 35% rather than a constant vapor pressure. The drying rate increased as the ambient 75 Table 13. Comparison of the Drying Constant and the Mass Transfer Coefficient Values for Treated Alfalfa Dried Simultanously with Control Alfalfa Trial 2 Trial 16 A g A B k1 0.131 0.577 0.213 0.399 k2 0.093 0.625 0.180 0.530 k3 0.025 0.042 0.073 0.070 HD1 0.705 3.091 1.095 2.059 HDZ 0.383 2.560 0.727 2.101 H03 0.114 0.133 0.209 0.191 k - drying constant K1 - average k value for moisture contents greater than 2.0 K2 - average k value for 2.0 > moisture > 1.0 K3 - average k value for moisture < 1.0 H - mass transfer coefficient H 1 - average HD value for moisture > 2.0 HD2 - average HD value for 2.0 > moisture > 1.0 H 3 - average HD value for moisutre < 1.0 76 _ 4.00 CD¢31CDCII (JICOCWH73C7 go O'Ik(3|§ 0- JBCUH3350 3.00 2.00 MOISTURE CONTENT (081 1,00 .. —~— .00 15.00 20.00 35.00 40.00 50.00 TIME IN HOURS Figure 9. Plot of the Two Exponential Equations of Best Fit Generated from Potassium Carbonate Treatment Data cp.00 temperature increased. 77 CHAPTER 8 SUMMARY Three methods were used to analyze and compare the drying of alfalfa under controlled conditions. A drying chamber was built and' used with equipment developed by Byler (1983) to monitor the drying of alfalfa under controlled conditions. The drying chamber held two samples which could be simultaniously monitored by a microprocessor which was connected to the weight transducers which held the samples. A data collection procedure was established and analyis techniques were tried. Data was transfered to the mainframe computer where average drying constant and mass transfer coefficients were determined for each trial. Curve fitting was also used on some trials. 78 79 Both methods of analysis used, determination of a drying constant, k, and a mass transfer coefficient, HD, fit the data well, particularly at moisture contents greater than 1.0 (dry basis). The drying constant increased with an increase in ambient temperature. The drying temperature was the driving force in the drying constant analysis, and both the drying temperature and the vapor pressure were important factors in the mass transter coefficient analysis. Potassium carbonate increased the drying rate of cut alfalfa. This effect was greater at ambient temperatures above 35 ’8 where this solution would be particularly useful to speed drying. The increased drying rate was observed at all temperatures but might not be economically feasible for ambient temperatures below 30'C. B I BLI OGRAPHY BIBLIOGRAPHY Ajibola O., Koegel R., and Bruhn H. D.: Radiant energy and its relation to forage drying. Transactions of the ASAE 23(5):1297-1300, 1980. Bakker-Arkema F. W., Hall C. W. and Benne E. J.: Equilibrium Moisture Content of Alfalfa, MI Ag Exp Sta Quat Bul 44:492-496, 1962. Barre H.J.: Vapor Pressures in studying moisture transfer problems. Agricultural Engineering 19(6):247-249,1938. Barrington G. P. and Bruhn H. D.: Effect of mechanical forage-harvesting devices on field curing rates and relative harvesting losses. Transactions of the ASAE 13:874-878,1970. Bolton J. L.: Alfalfa; Botony, Cultivation and Utilization, London, England, Leonard Hill TBooksT Limited, Interscience Publishers, Inc., 1962. Byers G.L. and Routly D.G.: Alfalfa drying - overcoming natural barriers. Agricultural Engineering, Sept., 1966. Byler R.K.: Parameter estimation methodology in selected moisture desorption models. (Unpublished Dissertation, Michigan State University, 1983). Chambers T.C. and Possingham J.V.: Studies of the fine structure of the way layer of sultana grapes. Aust J Biol Sci 16:818-825, 1963. Columella L.J.M.: Husbandry 1g Twelve Books, translated by Pliney et al., London, printed for A. Miller, M.CC.XLV (1745). pp. 520-521 (Chapter XVI - of making dried raisins or raisins of the sun.) @ 60 A.D. 80 81 Dudman W.F. and Grncarevic M.: Determination of the surface waxy substances of grapes. J Sci Food Agric 13:221-224, 1962 vanElderen E., vanHoven S.P.J.H. and Kroeze G.H.: Short wave radiation in a wheat crop during the harvest period. J Agric Engng Res 17:94-98, 1972. Fairbanks 6.8. and Thierstein C.E.: Performance of hay conditioning machines. Transactions of the ASAE:182-184,1966. Greenhill W.L.: The respiration drift of harvested pasture plants during drying. J Sci Food Agric 10:495-501. 1959. Grncarevic M., Radler F. and Possingham J.V.: The dipping effect causing increased drying of grapes demonstrated with an artificial cuticle. Am J Enol 19:27-29, 1968. Hadley N.F.: Surface waxes and integumentary permeability. American Scientist 68:546-553, 1980. Hall D.M. and Jones R.L.: Physiological significance of surfage wax on leaves. Nature 191:95-96, 1961. Hamilton J.W.: Transpiration control by native clover epicuticular wax. Advancing Frontiers of Plant Sci 30:175-187, 1975. Hammer P.C. and Day C.L.: Thin layer hay drying. Transactions of the ASAE 1967. Harris C.E.: The effect of organic phosphates on the drying rate of grass leaves and dry matter losses during drying. J Agric Sci, Camb. 91:185-189, 1978. Harris C.E. and May-Brown R.: The effect of tri-n-butly phospate on the drying rate and respiration rate of grass leaves measured in the laboratory. J Agric Sci, Camb. 86:531-535, 1976. Harris C.E. and Thaine R.: The effect of thermal and mechinical treatment on the drying rates of leaves and flowering stem internodes of Italian rye grass (Lolium nullifloriun) measured in the laboratory. J Agric Sci, Camb. 85:325-329, 1975. 82 Harris C.E., Thaine R. and Marjatta Sarisaco H.I.: Effectiveness of some mechanical, thermal and chemical laboratory treatments on the drying rates of leaves and stem internodes of grass. J Agric Sci, Camb. 83:353-358, 1974. Harris C.E. and Tullberg J.N.: Pathways of water loss from legumes and grasses cut for conservation. Grass and Forage Sci 35:1-11, 1980. Hayhoe H. N. and Jackson L. P.: Weather effects on Hay Drying Rates, Can J Plant Sci, 54:479-484, 1974. Hearle W.L. and Paterson H.: Thhe design and operation of systems with special reference to drying in the bale. J Proc Inst Agric Engrs 17:39-50, 1961. Hill J.D., Ross I.J. and Barfield B.J.: Use of vapor pressure deficit to predict drying time for alfalfa hay. Transactions of the ASAE 20(2):372-374, 1977. Holloway P.J.: Chemistry 'of leaf waxes in relation to wetting. Holman, J.P.: Heat Transfer McGraw-Hill Book Co., New York, New York, 1976. J Sci Food Agric 20:124-128, 1969. Jennrich R.: Nonlinear REgression, P3R, Abstract. BMDP Statistical Software, University of California Press, Berekley, CA, 1981. Jones T.N. and Palmer L.O.: Field curing of hay as influenced by plant physiological reactions. Agricultural Engineering 199-200, 1932. Kaldy M.S., Hanna M.R. and Smoliak 8.: Influence of drying methods on protein content and amino acid composition of three forage legumes. Can J Plant Sci 59:707-712, 1979. Kemp J. G., Misener G.C. and Roach W.S.: Development of Emperical Formulae for Drying of Hay, Transactions of the ASAE 5(4):723-725, 1972. Kennedy W.K., Hesse W.H. and Johnson C.M.: The effect of herbicides on the drying rate of hay crops. Agronmy Journal 46:199-203, 1954. 83 Kepner R.A., Goss J.R. and Meyer J.H.: Evaluation of hay conditioning effects Agr Eng May, 1960. Klinner W.E.: Design and performance characteristics of an experimental crop conditioning system for difficult Climates. J Agric Engng Res 20:149-165, 1975. Kramer P.J.: Symposium: Responses of field crops to environmental factors, summary statements. Agronomy Journal 55:31-35, 1963. Leshem Y., Thaine R., Harris C.E. and Canaway R.J.: Water loss from cut grass with special reference to hay-making. Ann Appl Biol 72:89-104, 1972. Mears D.R. and Roberts W.J.: Methods of accelerating forage drying. Transactions of the ASAE 531-533, 1970. Menzies D.J. and O'Callaghan R.O.: The effect of temperature on the drying rate of grass. J agric Engng Res 16(3):213-222, 1971. Merva C.E.: Physioengjneering Principles, AVI Publishing Co., Inc., Westport, Connecticut, 1975. Michigan Agriculture Statistics (preliminary statistics) 1982. Meyer B.S. and Anderson D.B.: Plant Ph siolo D. VanNostrand Co. inc, New York, NY, 1939. Meyer and Anderson : (1939) O'Callaghan J.R. and Creig D.J.: Work on conservation problems at the University of Newcastle-Upon-Tyne., J Proc Instn agric Engrs 26:75-80, 1971. O'Callaghan J.R., Menzies D.J., and Baily P.H.: Digital Simulation of Agricultural Drier Performance. J Agr Engr 16(3):223, 1971. Pedersen T.T. and Buchele W.F.: Drying rate of alfalfa hay. Agricultural Engineering, Feb., 1960. Person N.K.Jr. and Sorenson J.W.Jr.: Drying hay with infrared radiation. Agricultural Engineering, 204-207, 226-227, 1962. 84 Person N.K.Jr. and Sorenson J.W.Jr.: Comparitive drying rates of selected forage crops. Transactions of the ASAE, 352-356, 1970. Pohl J. and Fechner M.: How a drying plant influences intensification in the development of forage production along industrial lines. Proceedings of the XII International Grassland Congress, Leipzig, 1407-1409, 1977. Poissingham: Surface wax structure in fresh and dried sultana grapes. Ann Botany 36(148):993, 1972. Priepke E.H. and Bruhn H.D.: Altering physical characteristics of alfalfa to increase the drying rate. Transactions of the ASAE, 1970. Rotz C.A., Sprott D.J. and Thomas J.W.: Interaction of mechanical and chemical conditioning of alfalfa. Transactions of the ASAE (in press), 1984. Schieferstein R.H. and Loomis W.E.: Wax deposits on leaf surfaces Plant Phys 31:240-247, 1956. Schonherr J.: Water permeability of isolated cuticular membranes: the effect of pH and cations on diffusion, hydrodynamic permeability and size of polar pores in the cutin matrix. Planta 128:113-126, 1976. Schonherr J. and Bukovac M.J.: Ion exchange properties of isolated tomato fruit cuticular membrane: Exchange capacity, nature of fixed charges and cation selectivity. Planta 109:73-93, 1973. Shepherd J.B., Wiseman H.G., Ely R.E., Melin C.G., Sweetman W.J., et al.: Experiments in harvesting and preserving alfalfa for dairy cattle feed. USDA Technical Bulletin No. 1079, 1954. Shepherd W.: Effects of a chemical desiccant on the speed of curing of hay. Australian Institute of agric Sci, Journal, 25:218-221, 1959. Shepperson G.: Ouick haymaking machinery and methods. Outlook on Agriculture, 1964. Slatyer R.O. and Taylor S.A.: Terminology in plant- and soil-water relations. Nature 187:922-924, 1960. 85 Stafelt M.G.: The water output of the guard cells of the stomata. Physiologia Plantarum 10:752-773, 1957. StetsonL.E., Ogden R.L. and Nelson 5.0.: Effects of radiofrequency electric fields on drying and carotene retention of chopped alfalfa. Transactions of the ASAE 407-410, 1933. Tullberg J.N. and Angus D.E.: Increasing the drying rate of lucerne by the use of chemicals. J Aust Inst Ag Sci Sept., 1972. Tullberg J.N. and Angus D.B.: The effect of potassium carbonate solution on the drying of lucerne. J Agric Sci, Camb. 91:551-556. 1978. Turner N.C.: Fusicoccin: a fungal toxin that opens stomata. Nature 223:1070-1071, 1969. USDA Agricultural Statistics, United States Department of Agriculture, 1982. Walker D.A. and Zelitch 1.: Some effects of metabolic inhibitors, temperature, and anaerobic conditions on the stomatal movement. Plant Phys 38:390-396, 1963. Whitney L.F., Agrawal H.M. and Livingston R.B.: Stomatal effects on high temperature, short-time drying of alfalfa leaves. Transactions of the ASAE 12(5):769-771, 1969. Wieghart M.: Acceleration of Forage Drying by Chemical Application at Cutting. (unpublished Thesis, Michigan State University, 1979 . Wieghart M., Thomas J.W. and Tesar M.B.: Hastening drying rate of cut alfalfa with Chemical treatment. J Animal Sci 51:1-9, 1980. Wolf D.D. and Carson E.W.: Respiration during Ddrying of alfalfa herbage. Crop Science 13:660-662, 1973. Wood J.G.M. and Parker J.: Respiration during the drying of hay. J agric Engng Res 16(3):179-191, 1971. Zelitch 1: Biochemical control of stomatal opening in leaves. Proc US Nat Sci 47:1423-1433, 1961. 86 Zelitch 1: Control of leaf stomata - their role in transpiration and photosynthesis. American Scientist 55(4):472-485, 1967. APPENDI X A Table A.1 86a CALCKHD, Program to Calculate the Drying Constant and Mass Transfer Coefficient _A 100: 110: 120: 130=C 143:C 150=C 160=C 170: 130: 190: 200210 210: 220:1 230: 240:2 250: 253=C 273: 280=C 293:3 300: 313=C 320=C 330:0 340=C 350: 362:13 37C: 350: 393=C 400=C 410:0 023=C 430=C 440=C 05C=C 460: 47;: 052: 090: 5C0: 518: 520: 530: 593:110 55:: 56C: 573:120 58:: 590: 620=C 613=C 620=C 530: 542: 650=C 563: $70=C 680=C 590: 700: 710=C 723: PROGRAM CALCRnccIxcuT.OUTpUT.TApE1.TARET.TAFEI=:~RUT.T4=EE.TA=59> connou T4420) DIRENSIOR AHC4400108MC4400) THIS PROGRAM CALCULATEs THE DRYING CONSTART. K. I THE "‘ss TRANSFER COEFFICIENT. H00 400 THEIR RESPECTIVE R2 VALUES FOR THREE OROuPsts or ROISTURE CONTENTS: START To 2 RC. 2 T0 1 RC AND 1 T0 :00 or TRIAL REUINO I 00 10 I=I.4OO ANCIIO=BHCCI):T(1)=J.0 CONTINUE PRINT 1 FORRAT (zen UHAT {:10 13.2) NAME =0RRAT 416.110) PRINT 3 IS THE DATA FILE NAME?) INPUT CONSTANTS FORMAT (IZQH ENTER ClofnolDHOEDFgNUH77 READ 'QCAoEHOAOHqBDflgNUM TAPE 1 - DATA TAPE. HITH TIVE IND HOTSTURE CONTENTS FDR THO SAMPLES AT EACH DATA COLLECTION POINT. READ (1013) ((T(J)9A”C(J)03HC(J)0J=10NUF)) FORHAT fCF120701X0F50301195503)) PRINT 00 T41)0T(2)9T4NUH) DO 130 IPASS=102 RUN THROUGH ANALYSIS FOR SIDE A0 THEN SIDE E IT - CCUNTER F03 NUHBER OF VALUES SUHHED ”H - COUNTER FOR PC RANGE (5-202-101-81 SUHK - VALUES TO BE AVERAGED TO DETERMINE K SUHHD = VALUES TC 8E AVERAGED TO DETERMINE HD Q o IT an 1 SUHK = SUHHD = 303 30 2C8 V = ZONU" IF (IP4550E002) GOTD 113 ZHC : AHCCH) THC : APCCH-l) SOTO 120 CONTINUE ZHC = BHCIH) Yflc : BMCIH-1) CONTINUE IT = IT 0 1 V0 = I CALCULATE THE SUM OF THE DRYING CONSTANTS FOD 4N INTERVtL SUPK O (C'ALOO((ZHC'EH7/(YHC-EH))7/(T(“)-T(H-1))) (THC-ZHC)I(T(H)-T(”-I)) SUHK : AHR : TC C‘LCULATE THE “ASS TRANSFER COEFFICIENT FIRST CALCULATE THE MOISTURE RATIO. AHRO AND THE CONCENTRATION OF HOISTUPE 11 THE PLANT IZHC'ADH)/t(17.67-ZHC)027052) SUHHD 0 (AHP/(CP-CA)) 3P = SUHPD = IF CUTSFF POINTS HAVE BEEN REACHED CHECK TO SEE 87 Table A.1 (continued) A 730=C 740: 750: 166: 770:300 780=C 790=c 800=C 813: $23: 8 0: 843:: 850=C 863:C 870: 880: 890: 900: 910:310 920:320 93:: 94::330 950: 960= 970: 98?: 990: 1000:000 1c10=c 1:20:C 1030:C ICCO: 1:5C= 1:60: 1:70: 1580: 1990:010 1100: 1110: 1120: 1130: 1140:200 1150:C 1160:: 1170: 1180: 1190: 1200: 1210: 1226=C 1230:C 1240:C 1250: 126C=610 1270: 1293: 1290:C 1300=C 1310=C 1320:C 1:30:20 IF ((ZHCALT020)AAVDACHH0E901D) SOTO 300 IF ((ZHC0LT01010ANDA(HH0E002II SOTO 400 SOTO 200 I CONTINUE CALCULATE AVERAGE VALUES (KOHD) FOR MC > 20: AKI = SUHK/IT AHDI 3 SUHHC/IT suaROUTIHE RSORO CALCULATES THE FIT OF THE THEORETICAL CURVE FRoH TH! AVERAGE R AND «O VALUES IF IIPASSoEOaI) CALL RSORDCAHCoAkloAHDIVRZAIRZEONOVHM03OUOCA0E“) IF (I'ASSAE002) CALL RSORD(BHCOAKI0AHDIORZAORZEONOOHH09DH0CA0E”) IF (IPASSoEOol) URITE (70310) HAPE IF CIPASS05002) IRITE (70329! FORMATCIH 0A801019H VALUES FOR SIDE A:) FORMATTIH .19H VALUES FOR SIDE 8:) HRITE (70330) AKIORZAoAHOIORZBOIT FORMATCIH 0CX02TH FOR "C05T02: K: 0‘50399‘05"923 0 9F6030/031X06H H0: 0F50303X05H82: 0F60303IO9HNUH pTS: 0I3) SUNK 3 SUHHD = 00: IT 3 0 MM : 2 SOTO 233 CONTINUE CALCULATE AVERAGE VALUES (KoHD) FOR 2 ) "C > I AKZ : SUHK/ZT ANOZ 3 SLFHO/IT IF (IpASS0E001’ CALL RSORDCAHCOAKZOAHDZOPZAORZE0N09”"9‘0"0C‘05“) IF (IPASSoEOoZ) CALL RSORO(BHCOAKZOAH020R2A0RZEONORHUORDH0CA0E") URITE (70919) AKZORZAOAHDZORZBOIT FORHAT (IH OEN02TH FOR IoOoLToNCALT020: K: 0F50309105H°23 9 0F6030/032105F ND: 0750303105552: 9‘60303X09HNUH FTS: 013) SUHK=SUMHD = 300 HM : 3 1T : S CONTINUE CALCULATE AVERAGE VALUES (KOHD) FOR MC < loC AK} : SUHK/IT AND! 3 SUHHD/IT IF (IPASSAEOAI) CALL RSGRDCAPCOAKSOAHDIORZAOPZCONOOPHOADFOC405“) IF (IPA550E002) CALL RSGRD¢8VC0AK3OAHDSORZAORCSONUAHHABDPOCAOEP) “H : 5 UNITE CUTPUT ON TA:E 7 URITE (70615) AK30°ZAAAPD309239IT FORMAT(IH 0EX925“ EUR :000LT0UC0LT0I033 “3 O‘C0PO3X05H92: 0 ¢F60301932105F "D: 0F50303X05592: oF6-303‘09HUU” cTS: 0130/) IF (I‘ASSAEcol) IRITE (302C) AKIO‘KZ0‘53OAHDIVAHDZIAH93 URITE K AND HO VALUES FOR SIDE A ON TAPE 8 HRITE K AND HO VALUES FOR SIDE 8 ON TAPE 9 FORMAT (6(F60401K71 88 Table A.1 (continued) 13‘03 IF (1“SSOEGOZI 1917: A902?) AKIQAKZQAK3AAFDIQAHDZOAHDB 1350:100 COATINLE 1360:' END 1370: SUBROUTINE RSORD(AHCoAKoAHDoREGRKoREORHoNOo'HoADHoCAoE“) 1380: COHHON TAAOO) 1390: OIHEASION AHCTAOJ)oAHAAOO).AHHO(AGC) 1000:C 1~10=C THIS SUBROUTINE CALCULATES A REGRESSION AND RETUPNS F-R or 1420=C THE ACTUAL DATA COMPARED TO THE THEORETICAL DATA 1030=C IOQO: IF (”MOEGO1, ISTART : 2 1.50: IF (HH.E0.S) ISTART : 2 1060: “AH(1) : AMHDIl) : ANCC1) 1Q10: ANCISTART-l) : ANHOAISTART-1) : AHCTTSTART-l) 1420: 00 100 J:ISTART9VC 1A90:C 1500=C CALCULATE TUE THEORETICAL CURVES. BASED ON AVERAGE K9 HO VALUES 1E10=C 1520: AHtJO = C(APCJ-11-EH)OEIP(-AKO(TCJ)-TIJ-1))))¢E” 1530: CP = (AHHD(J-1)OADH)/(CIT-GTOAFHD(J-I)) 0 27.52) 1540: AHHDIJ) : A“HDCJ'1) - «AND-(TAJT-TCJ-l)T-(c9-CA)) 1550: 100 CONTIUUE 1560: L : 0 1570: SUHA = suHK = suHH : 3.9 1580: ASUHK = ASLHH : SUHAZ : SUMKZ : SUHHZ : 0.: 159;: no 110 N : ISTARTq~O xeoo=c 1&10=C SUM THE VALUES ‘09 THE REGRESSION EQUATIONS 1620=C 1630: L = L 0 1 1640: SUNA : SUMA o AMEtN) 1650: SUMK : SLIK o AHTLJ 1660: SUFH : SUHH . AMHDAN) 1&70= ASUHK = ASUMK o (AHCCNIOAFCN)) 1590: ASUHH = ASUHH o (AHCAN)~AHHD(H)) 1&90: SU~A2 = SUHA2 o (=HC(~)-o2) 1750: SUHKZ = SU“K2 . (A‘CN)0'2) 1710: SUFHZ : SUHHZ o (:HHDLN)--2) 172C=110 CONTINUE 1730: R1K : ASUMK - CSJMAOSUFKIL) 1740: R1H : ASUHH - (SJHA-SUNH/L) 1750: RZK : (A$UHA2-ASUHAao2)IL)0(SUHK2-(SUHKo.2)/L))ooc.5 1760: RZH : c(SUHAz-CSUHA..2)/L)oASUMHZ-ASUHHc-2)IL))c¢f.5 1TTc=C 1780=C CALCULATE RcR 1790=C 1500: REORK : (R1KIR2K)--2 1610: REGRH : (R1H/R2H)--2 1E2C: ISTART : N0 e 1 1é3c= RETUR~ 1540: END # 89 Table A.2 BMDP3R, Program to Calculate a Two Exponential Best Fit Curve using BMDP (BioMedical Data Package) 18.:/ RCcLVH TITLE 15 '2 ixpouimTch - Fl‘A'. 1‘2:/ TLFUT UNIT : 55 . i‘2= VsfilAfiLES £9: 2. 1‘5: ‘CDVLT IS '(F1;.’olxo:5.')'. 2 ;- CALE : 32;. :1‘=/ Vt=iaEL; fl:V?S LR: T:M:,M‘:STUD§, 221: L5: = TI"EV”DISTJ°EA 23' I ”EE‘iSS CS‘ENEEKT IS "OISTUF . 25 : ‘L-T’bE: : 10 23?: g-,-I‘u"~*~'f~.'TiF?S ADE ?. 2" INSESENUEWT IS TIVE. 2’ = CC“VF°G£NCE = T. :1. 2: = ITL:LTIS’S = S‘ . 2=‘: F:LVZRS - é . 3';- :zt‘lT : f." 0 31'=/ FLF1”ETEF IKZTIAL = 3.39 -‘.19 2. ., -'.-_, 2. 2-- - .~. __ — . . -4. " """‘ ' :0 T 59 -00 o" o 33‘:/ PLZT 3:223L5L. 34 - ‘ZZE : 33.: . 39': V:‘I:ELE : TIME~VCICTU=€. 3.4 = '.o=“=L. 37.:1/ —..D c \U APPENDI X B 9O 0000000.00 00000.0.00 0000000..0 0000000.0. 0000000.0. b000000.h. 0000000.h. 0000000.0. 0000000.1. 00000.0.0. 0000000.0. 0000000... 0000000.0. h000000.0 0000000.0 0000000.h h000000.0 00000.0.1 0000000.0 0000000.0 0000000.0 0000000. 0000000.0 30000.0.00 0000000.00 0000000.h0 0000000.00 b000000.00 0000000.00 0000000.00 0000000.00 h0000.0..0 0000000.00 0000000.01 0000000.01 h000000.h1 0000000.01 0000000.01 0000000.01 00000.0.01 0000000..1 0000000.01 00000.0.01 h000000.00 0000000.h0 0000000.00 h000000.00 0509 0.0. 000. 050. 000. hhh. 000. 000. 000. 000. 00.. -NNHC FQ'OO FUWQN 0.0. assg F00 0000000'4-‘44ev- 000. 001. 00.. h... 00.. h... h... 00.. h... .0.. h... h... .0.. .0.. .0.. .0.. .0.. h... b... 00.. 00.. 00.. 01.. m muqm 010. 000. 000. 000. 010. 000. 000. 010.. .0... 000.. 000.. 000.. 0.0.. 000.. 000.. 000.0 .0..0 010.0 000.0 010.0 000.0 .00.0 100.0 0 41.0» ..0. 01.. 01.. 00.. 01.. ..0. 10.. 10.. ..0. 10.. 10.. 00.. 00.. 50.. 10.. 05.. 10.. 00.. 10.. 10.. 000. 010. 000. 000. < «cam macwus Ham you mammm zapsom am so mucmucoo ousumwoz mcu motwcfiumwa .0 OHAMHT 50000.0.10 0000000.00 0000000.00 00000.0..0 0000000.00 0000000.00 0000000.h0 h000000.00 0000000.00 0000000.10 0000000.00 h0000.0.00 0000000..0 0000000.00 0000000.0. h000000.0. 0000000.h. 0000000.0. 0000000.1. 00000.0.0. 0000000.0. 0000000... 00000.0... h000000.0 0000000.0 0000000.h h000000.0 00000.0.0 0000000.1 00000.0.0 0000000.. 00000.0. 0000000.0 00000.0.00 0000000.00 0000000.00 0000000.00 h000000.00 0000000.00 0000000.00 0000000.00 50000.0..0 0000000.00 0000000.01 0000000.01 h000000.h1 000mm%MHW1 01.. .0.. 000. 00.. 00.. 000. 000. 050. 100. 0.0. 010. 100. 150. N01. .00. 000. .00. 0.0.. .hh.0 «an. hon. «an. awn. ”on. own. a... one. was. omc. omv. .ov. .oc. .oc. m mvwm 000. 000. 010. >00. #00. 000. .00. .00. 010. 000. 0.1. 0.1. 001. 001. 001. 000. 000. 100. #00. 11B. 000. 000. 000. 00... 00... 000.. 011.. 000.. 100.. 010.. 00..0 001.0 000.0 0 41—0. 000. 000. 000. 000. 0.1. 010. 0.1. 100. 0.1. 0.1. 011. 011. .01. 001. < mufim 0000000.01 0000000.01 00000.0.01 0000000 .1 0000000.01 50000.0.01 5000000.00 0000000.00 0000000.00 h000000.00 50000.0.10 0000000.00 0000000.00 00000.0..0 0000000.00 0000000.00 0000000.h0 h000000.00 0000000.00 0000000.10 0000000.00 50000.0.00 0000000..0 0000000.00 0000000.0. b000000.0. 0000000.h. 0000000.0. 0000000.1. 00000.0.0. 0000000.0. 0000000... 00000.0... 5000000.0 0000000.0 0000000.b h000000.0 00000.0.0 0000000.1 50000.0.0 0000000.. 50000.0. 0000000.0 :3... 08.. 001. 001. .01. 001. .01. .01. 000. 000. 000. 000. 000. 000. 000. 0.0. .10. .10. 000. 00h. 105. 500. 000. 000. 010. .00. 000. 000. 000.. 0.... 00.. 000. 000. 1.1. 000. 000. .00. 0.0. .00. 000. .01. 000. 000. .50. 000. 'nflflfiflfluvcwuuwp .08. 0 00.0 m.m. .mc. mum. .m'. .m.. m.m. n.o. new. can. can. «.m. one. cow. own. mos. own. «.o. you. ooh. moo. aha» can. one. cum. mum. nno.. moo.. no... o.u.. o.«.. «a... one.. new.. m.o.. no... nmo.. one.“ no..u can." ohm." owe.“ .va.n wan.c .0... 1 00.0 . a<.ub 91 5000000.0. 050- 000. 0000000.50 010. 005. 0000000.00 010. 000. 0000000.0. 500. .00. 5000000.00 .00. 000. 0000000.00 000. 000. 0000000.1. 500. ..5. 0000005.00 005. 100. 50000...50 000. 000. 0000000.0. 100. 005. 0000000.10 105. 500. 000000..00 005. 050. 50000.0.0. 050. 000. 0000000.00 0.0. 050. 0000005.10 005. 015. 0000000... 000. 010. 50000.0.00 000. 500. 0000005.00 000. 005. ocoooo«.o. 001. 100. 0000000..0 100. 0.0. 0600000.00 050. 000. nnnnnnn.0 1.0. 110.. ooooooa.0. .50. .00.. OOOOOmo.00 050. 000. 5000000.0 100. 11... 0000000 0. 0.0.. 500.. 0000000.00 .00. 000. 0000000.5 050. 110.. 5000000.5. 000.. 500.. 000000...0 510. 000. 0000000.0 .05. .00.. OOOOOOO.5. 00... 01... 0000000.00 100.. 010. 0000000.1 0.0.. .00.. 0000000.0. 00... 0.0.. 0000005.0 .55.. 050.. 50000.0.0 000.. 550.. 0000000 1. 000.. 000.. 0000005.0 000.. 005.. .5wowowo.0 100.. 110.. 50000.0.0. 000.. 000.. 5000005.5 550.. 000.. 5000000.. 000.. 000.. cocoono.0. 001.. 001.. OOOOOOo.0 .50.0 010.. 50000.0. 000.0 010.0 0000000... 001.. .01.. nnnnnno.0 no..0 .00.0 ooooooo.o 050.0 010.0 nnnnnnn.o. 010.. 050.. 5oommmo.1 000.0 55..0 . 5000000.0 000.. 000.. oooooom.n 501.0 0.0.0 OOOOOOn.0 0.5.. 005.. cocoomo.0 050.0 001.0 0 41.3. nnnnnnn.5 000.. 000.. 0000000.. 005.0 .00.0 5000000.0 050.. 500.. 0606000. .00.0 000.0 50000.0.1 000.0 .o..0 ooooooo.o 0.1.0 050.0 50000.0.00 000. 0.0. OOOOOma.0 o...0 000.0 5000000.00 150. 000. oooooom.0 010.0 .00.0 0000000.50 150. 000. 0000000.0 000.0 0.1.0 1 41.05 nnnnnno.00 150. 010. 0000000. 051.0 110.0 5000000.00 000. 000. ooooooo.o 1.0.0 050.0 0000005.10 000. 010. 0000000 Om .00. .10. ooooomn.nm 000. 010. oooooom.o1 .00. .10. 0000000.00 000. 550. 0 41.05 0000000.51 000. .00. 50000.0..0 000. .00. 5wowwwm.w1 000. .00. OOOOOmn.On 000. .00. ooooooo.01 010. 000. 0000000.01 000. 550. 6600000.0m 000. «on. «000000.01 010. «an. nnnnnnm.51 000. 550. nnnnnno.01 500. «on. nnnnnnm.n1 Own. 010. 5000000.01 .00. 100. 0000000.01 000. 000. 50000.0.01 000. 000. ooooooo.01 500. 0.1. 50000.1.51 000. «O1. 0000000..1 000. .00. nnnnnno.m1 .00. 100. 0000001.m1 000. .01. 0000000301 000. 000. 0000000.01 000. 0.1. 0000001.01 000. .01. 0000000.00 000. 010. 50000.0.01 000. 011. 0000000.11 000. .01. 5000000.00 001. 000. 0000000..1 000. 011. 5000000.01 101. .01. OOOOOOn.5n 001. 001. nnnnno¢.01 000. 051. oooooo..01 101. 051. 0000000.0n 001. 001. nonnnnn.mn non. onv. nannnn...c mac. 5.m. 5ooooon.nn env. «o1. 5ouowou.on 000. 001. 500000..O1 001. .01. 50000.0.00 011. 001. 0000000.5n 111. 000. 50000.5.00 011. 011. 0000000.00 001. 001. 0000000.00 001. 001. 0000005 50 011. .01. 0000000..0 001. 001. 5000000.00 111. 051. 0000005.00 .01. 011. 50000.0..0 051. 001. 50000.0.00 051. 000. 50000.0.00 051. 001. 5000000.00 001. 100. ooooomo.un .91. can. 5ooooon.vn omq. mac. ooooooo.ou mac. nov. 0000000..0 000. .00. 0000001.00 5.0. 050. 0000000.50 0.0. 100. 50000.0..0 500. 500. nnnnnn1.00 000. 5.0. 5wowwwo.on 000. 100. 5000000.00 000. 010. 5000001..n 100. 050. 0360005.mn 000. 100. OOOOOOo.oa 000. 100. 0000000.00 «mm. 050. 0600000 10 000. 000. 02.. m mu.m < mu.m 02.. m mu.m < an.m 02.. m mefim < mc.m .cwncwucoov .0 magma 92 5555551 .5 115. 511. 5555555.5 555.. 555.. 5555551 55 115. 511. . 5555555.5 .55.5 555.. 5555551 51 555. 155. 5555555 55 .15. .5.. 55555.5 1 55..5 555.5 5555551.51 555. 551. 5555555.55 555. .5.. 555555..5 5.5.5 5.5.5 5555551.51 555. 551. 5555555 55 555. 555. 5555555.5 551.5 555.5 5555555.51 555. .15. 55555.5 15 555. .5.. 55555.5.. 555.5 555.5 5555555 51 1.1. 555. 5555555.55 555. 555. 5555555.5 555.5 515.5 0000000.11 1.1. 000. 0000000.50 000. 500. 50000.0.01 001. 050. 50000.0..0 500. 500. 50000.0.01 001. 000. 50000.0.00 0.0. 500. 5 40.05 5555555..1 551. 555. 5555555.51 5.5. 555. 555555..55 555. 515. 5555555 51 5.5. 515. 555555..55 555. 555. 55555.5.51 5.5. 515. 5555555.55 555. 55.. 555555. 15 555. 155. 5555555.51 555. 155. 5555555.55 555. 51.. 555555. 55 555. 515. 5555555.51 555. 155. 5555555.55 555. 51.. 555555..55 .55. 555. 5555555 11 555. 555. 5555555.55 555. 51.. 555555. .5 555. 5.5. 5555555 51 555. 555. 5555555.55 555. 55.. 555555..55 .55. 555. 5555555.51 555. 555. 5555555.15 555. 55.. 555555..55 155. ..5.. 5555555..1 115. 5.5. 55555.5.55 555. 51.. 555555. 55 5.5. 515.. 5555555.51 555. 555. 5555555.55 5... 55.. 555555..55 155. 55... 5555551.55 115. 5.5. 5555555..5 5... 55.. 555555..55 555. .5... 5555551.55 .55. 555. 55555.5.55 5... 55.. 555555..55 555. 555.. 5555551.55 .55. 5.5. 5555555.51 5... 55.. 555555. 15 5.5.. 155 . 5555551.55 555. 555. 5555555.51 15.. 51.. 5555555 55 515.. 555.. 5555555 55 555. 555. 5555555.51 15.. 5.5. 5555555 55 .5... 555.. 5555555.15 5.1. 555. 5555555.51 15.. 55.. 5555555 .5 15... .51.. 55555.5 55 5.1. .55. 5555555.51 .5.. 55.. 55555.5 55 155 . 555 . 5555555.55 551. .55. 5555555.51 .5.. 5.5. 55555.5 5. 555.. 555 . 5555555..5 511. 5.1. 5555555.51 15.. 555. 55555.5.5. 555.. 555.. 5555555.55 551. 511. 55555.5..1 .5.. 555. 5555555.5. 5.1.. 555.. 5555555.55 551. 551. 5555555 51 15.. 55.. 5555555 5. 555.. 555.. 555555..55 155. 551. 5555555.55 .5.. 555. 5555555.5. 555.. 5.5 . 555555..55 .55. .55. 5555555.55 55.. 555. 5555555 1. 555.. 555 . 555555..55 555. 555. 5555555.55 15.. 5.5. 5555555.5. 5.5.. 555 5 555555..55 555. 555. 5555555 55 .5.. 515. 5555555.5. 555.. 555.5 555555o.15 555. 555. 5555555 55 .5.. 515. 5555555... 555.. 55. 5 5555555.55 555. 555. 5555555.55 55.. 555. 5555555.5. 555.. 5.5.5 55555.5.55 .15. .55. 55555.5 55 55.. 555. 5555555 5 .55 5 555 5 5555555..5 555. 555. 5555555..5 55.. 555. 5555555 5 55. 5 555 5 5555555.5. 555. ..5. 5555555.55 55.. 555. 5555555 5 555 5 551 5 5555555.5. 555. 515. 5555555 55 55.. 555. 5555555.5 555.5 .55.5 5555555.5. 515. 555. 5555555 55 555. 555. 5555555.5 .51.5 115 5 5555555.5. 555.. 555. 5555555.55 55.. 555. 55555.5 1 1.5.5 555 5 5555555.5. 555.. ..5.. 5555555.55 555. 555. 5555555.5 555.5 555 5 55555.5 1. .5... 555.. 5555555.55 555. 555. 5555555 5 555.5 555 5 5555555.5. 555.. 55. . 55555.5 15 555. 555. 5555555.5 .55 5 555 5 5555555.5. 555.. 555.. 5555555.55 555. 555. 55555.5.. 55..5 55. 5 55555.5... 555 . 555.. 5555555 55 555. 5.1. 5555555 5 555.5 555.5 5555555 5. 551 . 155.. 55555.5 .5 555. 511. 5555555.5 555.. 551.. 5555555.5. 555. 551. 5555555.5 555.. 555.. 5555555 5. 555. 551. 5 .1.5. 55555.5.5 555.. 555.. 5555555.5. 555. 555. m 5 w 5 5 55.5 < 55.5 55.5 5 55.5 < 55.5 55.5 5 5. < 5. 55.5 ll. - AcmscwucouV .0 5.055 93 oo000m9.~. has. n.a. occooo~.. «om. moo.n cocoons.on nmm sum. 0000000... 000. .50.. 00000000 mum. «and 000000.20.” now 000. nannnno.o. ohm. an... swwmuo..vn was amp. 0000000.0 50... 000.. h00000h.00 00h 00h. swcuoom.o mnu.. ouv.. .. 4¢_a. swoomos.un no. n... 0000000.» .00.. "0.. 0000005..0 000 005. h000000.0 000.. 000.. 0000005.00 000 0.0. 0000000 0 000.. 000.. 0000000.00 000.. 000.. 000000h.00 500 000. 0000000 v 000.. h0..0 0000000.00 0.... 000.. 0000000.00 0.0 0h0. h0000...0 00..0 00'.0 0000000.00 N0... 00'.. 0000000.50 000 0.0. unannoo.u nuc.« .no.« nnnnnnu.on um... um... .mwow.o.ou maa .cm. swoomwo.. coo.« soo.« nannnoo.p« on... ocv.. sowww.c.mu 090.. noo.. 0000000.0 000.0 .00.0 0000000.h0 00... 0h'.. h0000.0.'0 .00.. 000.. pounce..pu an... up... annnnno.nn «a... moo.. soooo.o.on .o«.. no... nnnnnno.nu «a... ea... 0. a4—flh 000000'.00 O0... 00'.. 0000000..0 000.. Oh... oooooo..ou .o«.. o.m.. coo00mo.0u 'o«.. ..n.. nannnn..u. own . non." cooOOmo.o. m0n.. v.n.. 0000000.vv Q09. 000. 000000.... 000.0 000.0 h000000.0. 500.. 0.0.. 000000o.nv .qv. now. unannoo.o. ace.“ «cc.« panacea... ow... non.. 0000000.“. umc. «no. socoowo.o no..u seq.“ powwooo.o. mov.. .c'.. 0000000..' 00'. 0.0. 0000000 0 00..0 0.0.0 0000000.0. 000.. 000.. h00000h.0' 0.0. 000. 0000000.h 000.0 '00.0 0000000.'. ..0.. 0'0.. nonnnns.on .nm. nap. powwooo.n Ova." .oo.u nnnnnoo.n. «so.. now.. nonnnnp.on man. .no. nnnnnna.v non.“ v..." oooooo«.n. an... ovo.. 0000009.h0 000. .00. 0000000.0 000.0 0.0.0 0000000... 000.. 000.. 0000000.00 .00. 000. 000000..0 0.0.0 000.0 0000000.0. 5.0.. 000.. h000000.00 000. '00. 0000000.. .'h.0 00..0 00000.0.0 000.. 000.. 0000000.v0 000. 0'0. 0000000.0 000.0 00v.0 50000.0.0 00..0 v00.0 0000000.00. ..h. 000. 0000000.b 00..0 000.0 0000000.00 0'0. 050.. 0000000.0 '00.0 00..0 ooooooo..n nap. .oo.. o. 4<.u. annnnno.m mun.“ can.a h000000.00 500. '0... 0000000.' 00'.0 0h0.0 h000000.00 000. 00... 5000000.' 0'0.0 0".0 smmowwn.ou vvm. an«.. nannnnm.vm ace. a... nnnnnna.n mmm.u nov.n 0000000.h0 000. 000.. 0000000.00 00'. 00'. 0000000.0 000.0 000.0 0000000.00 0.0.. 000.. 0000000.00 00'. 0.'. 0000000.. v.0.0 '00.0 soooo.o.ma o.o.. ohn.. ooooOmo..m on'. use. .woooca. o.o.n "cw.“ h000001.00 000.. 000.. #000000.00 00v. 0". 0000000.0 .00.0 000.0 00000.. .00 0O. .. 00' .. h000000.0' 00' . 0" . oooo0m..vu no... ac... soooowo.ov onc. .mv. oooooov.'u ca... co... nnnnnom... one. o... m 4¢_u. ooooomo.vu «v... «a... nannno¢.oc a... new. 00000.0.00 0... .. 0.0.. 88000.0' 00'. 000. smmmwwn.nu no... a.n.. oooooo. we mow. mum. nnnnnnv.mm .mn. .mn. nnnnnno.nu no... can.. oooooo. n. o.m. cqm. nonnnn'.on .au. 0... oooooo..«n anu.. .nn.. .mcwm.p.ac can. ccm. oocoomw.sm .m«. .an. .ooowon.ua nu«.. nnm.. swcoo....v can. awn. OOOOOmv.wm men. .on. annnnon.v a.».« «mo.« .wowo.p.0v .pm. mom. swowwoc.mm «an. «ac. ooooomo.v opn.u oar.“ nnnnnnp.on .an. poo. .woocmc.cm man. mac. nannnn..n van.« no..n nnnnnnp.on n.u. hum. sowowmv.nm vvn. coy. ooooon..« can.“ has.“ nnnnnns.pn one. cow. nonnnov.um cvn. ocv. 05.9 0 000m < 000m 0508 0 mvfim < mc.w 0508 0 oc.m < 000m Aco=:.ucoov .m mfinmy 94 0000000.00 0000000.00 00000.0.v0 000000'.00 000000v.00 000000'..0 000000'.00 00000.'.0' 0000000.0' h000000.h' h000000.0' 0000000.0' 00000.0.vv 00000.0.0' 0000000.0' 0000000..' 0000000.0' 0000000.00 00000.0.00 50000.0.00 000000..00 000000..00 000000..'0 000000..00 h0000...00 0000000..0 h000000.00 0000000.00 0000000.00 00000.0.00 0000000.00 0000000.00 0000000.00 0000000.00 00000.0..0 00000.0.00 0000000.0. 0000000.0. 0000000.h. 0000000.0. 0000000.0. 00000.0.v. 000000h.0. h000000.0. 0000000... 000000h.0. 00000.0.0 0000000.0 0000000.h 0000000.0 0000000.0 00000.0.v 050B 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 0.0. 0.0. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. '00. 000. 000. 000. '00. .00. .00. .00. 0.0. 000. 000. 000. 000. 00'. 00'. 00'. .h'. .h'. 000. 000. 000. .00. 000. 005. ..0. 000. 0'0. 0'0.. 00... 000.. '0'.. '00.. 0 000m 0'.. 0'.. 0'.. 0'.. '0.. 000. 000. 000. 00.. 00.. 0.0. 0.0. 000. 000. '00. '00. 000. 000. '00. 0.0. 000. 000. 000. 000. 000. '00. 000. 000. 000. 000. 000. 000. 000. '0'. '0'. '0'. 00'. 00'. '00. 000. 000. 000. 000. 000. 0'0. 000. 000.. 00... 000.. 0".. 000.. < 000m 0000000.0 0000000.0 000000... 0000000.0 h000000.0' 0000000.0' 0000000." 0000000.0' 000000..0' 000000...' 000000..0' 000000..00 h0000...00 000000..00 h000000.00 0000000.00 0000000.'0 0000000.00 0000000.00 0000000.00 0000000.00 0000000.00 0000000.00 0000000.00 0000000.00 0000000.'0 h000000.00 0000000.00 00000.0..0 00000.0.00 0000000.0. h00000h.0. 0000000.h. 0000000.0. 000000h.0. 0000000.'. h000000.0. 0000000.0. 0000000... 0000000.0. 0000000.0 h000000.0 h000000.h h000000.0 0000000.0 0000000.' 000000..0 000%? c 0'0.. 000.0 000.0 000.0 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 000. 000. 000. 000. 0'0. 0'0. 0'0. 000. 0'0. 000. 000. 000. 0.0. 0.0. 0'0. 000. 000. 000. 00'. 00'. 000. 000. 000. 000. 00... 00'.. a mafia 000.. 000.0 .00.0 000.0 '. 41—0. .0.. .0.. .0.. 00.. .0.n 00.. 00.. 00.. 0b.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 0b.. 00.. 00.. 00.. v.0. 000. 000. 000. 000. 000. 0'0. 0'0. 000. 00'. 00'. .". 00'. '00. 000. 000. 000. 0v0.. < 000m 0000000.0 0000000.0 0000000..0 0000000.00 h000000.00 0000000.00 00000.0.00 00000.0.00 0000000.00 h00000h.'0 0000000.00 0000000.00 0000000..0 0000000.00 h000000.0' h000000.0' h000000.h' 0000000.0' 0000000.0v 0000000.v' 0000000.0v h000000.0' 0000000..' 000000'.00 000000'.00 00000.'.00 00000.'.v0 000000'.00 h000000.00 0000000..0 0000000.00 0000000.00 0000000.00 h000000.h0 0000000.00 h000000.00 000000h.0. h00000h.0. 000000h.h. 000000h.0. .0000000.0. 0000000.'. h000000.0. Ucfififi .00.. 0'..0 000.0 00.. 0b.. 00.. 00.. 00.. 00.. 00.. 00.. 00.. 0.0. 00.. 00.. 000. 0.0. 0.0. 000. 000. 000. 0'0. 000. 0'0. 0'0. 000. 000. 0'0. 000. 000. 000. 000. 000. 000. .00. 0'0. .00. 0'0. 000. 050. .0'. 0". 00'. .'0. 000. 000. 000. 0 000m 0'h.. '0..0 00..0 0. 40.0. 00.. 000. 000. 00.. 00.. 00.. 000. 000. 000. 000. 000. 0'0. 000. 000. .00. .00. 000. .00. .00. .00. .00. .00. 000. .'0. 000. .'0. .v0. 000. 00'. 000. 00'. 00'. 00'. 00'. 00'. 00'. 00'. 000. 000. 0.0. 000. 0.0. '.0. '00. < 000m AcmaafiuaouV .m magma 95 annamoo.n ..0.. o...v .wowooo.0n m... n... 00000m..0w 00o. mun. 0000000.. 000.0 000.0 0000000.00 00.. 00.. 0000000.0v 0.0. 000. 0000000.0 .90.0 0.0.0 0000000.?0 01.. 00.. 0000000.vv 0.0. 000. .0ooc.0.nn n... .n.. oooooon.nv 00.. men. 0000000.00 00.. 000. 000000..0' .00. 000. 0. .«.¢. . cocoomo..n nv.. 00.. .00000...c .nn. .0.. 0000000.00 0'.. 000. 000000..Ov .00. 000. nonnnn0.o. 00.. n... .ooow...mm .nn. 0... 0000'00.00 000. 0.0. 00000.0.00 0'.. 000. 000000..00 000. 00'. 0000000.0. 000. 000. 0000000.00 00.. 000. 000000. 00 0'0. 000. 0000000.0. 000. 0'0. 0000000.00 0'.. 000. 0000000 00 000. 0.0. ""'00.0. 000. 000. 0000000.'0 00.. 000. 0000000 00 000. 0". annom.o.0. nnv. can. 00000m0.nn 00.. a... nnnnnn0.vn non. mow. annoouo.n. gov. any. .0000.a.un 00.. a... .0000.0.nn man. 00'. ...0000.'. 00'. 00'. 0000000..0 00.. 000. 0000000.00 00'. 00'. ....woa.n. «an. ..o. oooooom.0« 0... con. .owwwom.0n ..c. «on. nnnnoom.n. ..0. can. nannnoo.o. .0.. 00m. .mwcwwm.m. mmv. 0.0. o....oa... can. «no. 00000mo.o. 0... man. 00000mm.ou omq. awn. cvvmwom.0. .c.. n00. .0000.o... .0.. man. nonnanm... n.v. can. umnoooa.m 00¢. a... .0000.o.0. .0.. man. 0000000.0. n.v. «on. oouoooa.o «on. a... .0000.o.m. .0.. u.n. nannnoo.m. .0.. 0.0. «9.0000 . 000.. 000. nnnnno..v. .0.. 00¢. nnnnnoo... onm. 0.0. nn«.nao.o 0'... cc... 000000.... .0.. now. .000000... mam. ave. ..«.cma.n on... ~.n.. nnnnnn.... a... 'mc. nonnnno... .on. new. .oooooo.v 000.. .vn.. manna..... a... ..c. .owmm.o... 0.0. 000. 0000000.. no... .0u.. .0000...0. 00.. nmm. .owow.o.0n can. 0... .wc.v00.n 0.... on... nannnom.m .0.. won. 0000000... 000. u... annom00.« on... now.“ 00000m0.o can. «no. .mmoow..o. emu. no.. 0mmm000.. n.n.n a...» 00000m0.. own. a... nonnnn.... 0... 00¢. annnooo. n00.m no... 000003..0 ..n. «no. nnnnnn..o. va.. ..0. .0000.0.m onq. mma. nnnnnn..m. man“ .00 . nannnom.v «on. .'... 000000.... ..0. «.0.. .. .<.¢. nannno..n a... ..v.. .000000... ..0.. no... nnnnnn... om0.. 0.... ooooomo... .00.. an... . . . .omoo.... ..... ....n ooooomo... cw... um... 0000000 00 0'. 00. 0000000.0 000.0 000.0 0000000 0. 0'0 . 000 . 0000000.00 00. 00.. 0000000.0 '00.. 00'.. 0000000..0 0'. 00.. 0000000.0 0.'.. 00'.. 0000000.00 0'. 0'.. 0. 41—00 0000000.0 000 . 000.. 00000.0.0' 0'.. '0.. 0000000.0 000.. 000.. 0000000.0v 00. 0'.. 0000000.0 000.. .00.. 0000000.0' 00.. 0'.. . . 0000000.' 000.. 000.. 00000mo.0v 00.. .0.. MMMMMM".MM ”mm. "H” 000000... ....u n...u 00000000v 00. .0. . 000000v.00 000. 00. .000000. a «00.0 000.0 ooooooo.wc av. n... ooooooc.vm can. am. 00000m0.. nov.u .av.u annnno..nc 00. v... .omooon.nm 00a. om. 0000000.0 000.. one.« .00000...v n... .n.. nannnnn..m can. .nn. .00000...q 0.. .n.. nnnnnnn..m 00.. men nnnnnn..0v av. .m.. nonnnnn.0m can. .nn m. .«.u. 000000.... a... n... 000000.... 00.. .nn. 000000.... nv.. 00.” .00000..ov 0.0. .nn. . . . 0000000 .0 on. 0.. 0000000.: .00. own. 0000000 .0 00. v0. 05.0 0 00.0 < 00.0 05.0 0 00.0 < 00.0 05.0 0 00.0 < 00.0 A0000.ucouv- .0 0.000 965 ......0... .... .... ........0. 0..” ...H 0000000.0 ..... ..0.. ......0... .... .0.. .00.....0. 00. 0... ......0.0. .0.. .... 0000...... 0... 0.. ....0.0 .. 0... .0.. 0000...... .0.. 00.. 0. ..... 000000.... 0... .... ....0..... 0... .... .......... .00. 0... ..0....... .0.. .0.. . . . . .......... ..0. ..0. ....0...0. ..o. ..0. .00.....0. ..0. ..0. 00000...0. ..n. ..0. ....0..... a... 0... 00000.. .. .0.. ..0. ....... .. 0.0. ..0. .......... .00.. ..0. ....0..... ..0. .0.. .......... ..0. 0... ..0....... ..0.. .00.. .......... 0... ..0. ....00.... .0.. ..0. .....0.... .00.. ..0.. 00.00..... .... .0.. .....0.... ..0. ..0. .....0..0. 0.... 0.... 0000...... .0.. .0.. ....0...0. ..o. .... .......... 0.... 0.... .00....... .... 00.. .......... .0.. 00.. .....0.... ..... .0... 0000...... 0... 0... 0000...... .... .0.. ....0..... ..... 00... ....0.. .. 0... .... ....00... 0.0.. ..... ....... .. ..... 0.... ..00...... .00.. .... ......... 0.... 0.... ....0...0. ..0.. ..0.. ....0...0. 0.0.. .00.. ....... . .00.. ..... ......... 0.... 0.... ......... ... . 00... ....0.... ..... .0... ......... ..... 000.. ....0.... .0... 0.... ...0..... 0.0.. .0. . 00000.... 00... 0.... ......... .0... ... . ...0.00.. ..... ..0.. 00000.... ..0.. 0.0.. 0000000.. ..0.. .00 . .....00. ..0.. ..... ........0 ..0.. 0.... .....00 0 ... . ... . ....... . ..... ..0.. .....00.0 .0... ..0.. ....000.. ..... ..... ..0..00.. ..... ..... .. 4.... .....00.. 0.... .0... .....00 . 00... 00... .....00 . 0.... 0.0.. ....000.. .0... ..... ....000.. 0.... ..... .....00. ..... ....0 .. 4.... ....00.... 0... o... .. 4.... .....0.... .... ...H ........0. ... . 00... mmmmnwu mu ”mun MM". .......... 00... 00... .......... ..... 0.... . . . .....0.... ...H 0..” 0000...... 0.... .0. . MMMMWMw.mm mum. “mm. ....... .. .0.. 000. .....0.... ..... ... . ......0... 0... .00. .....0. .. .... 000. .....0..0. 0.... ... . .0000.0... .... .00. ....0.. .. 0... .00. 000000.... ..0.. 0.0.. ......0... .0.. ..0. ....... 0. ..0. .00. .....0.... ..0.. .0... ...0..0... 0... .... ....0..... 0.0. ..0. ....00.... ... . 0.... 0000..0.0. .0.. 0... ....0..... ..0. 0... .......... 0.... 0.... ......0 .. n... .... ..0....... ..0. .00. .0......0. 0.... ..0.. ......0... .... 0... ....... 0. ..0. 0.0. .....0... ... . .0. . ......0... .... .00.. 0000...... 0... ..0. ......... 0.... ..... ....0.0... .... ..0.. 0000...... 0... ..0. ......... 0.0.. .0... ......0.0. .0.. 0.... ......0 .. .00. ..0. ........0 ..... .0... 0000..0.. ..0. .0... 000000.... ..0. 0.0. ........0 ..... ..... ......0.. ..0. ..0.. 0000.0.... .0.. ..0. 00000.... ..... .00.. .....00.. .0.. .0... ...1.0..0. 0.. .0. 000...... ..... .0... .....00.0 ..0.. 0.... .......... .... 00.. ....... . 0.... 0.... .....00 0 .0... .0... ....0..... .0.. .... 0000..... ..0 . ..0.. .....00 . ..... ..... ....... .. ... .0. .0..... .0. . .00 0 ....000.. .0... .0... 0 0. 0 .5.. . a... < m... .5.. . a... < a... .5.. . ... < .. {l‘ A000:.u:00v .0 0.00. 97 opphp'o.pu nan. .nv. nnnnnmo.m. paw. nan. vccm.nm..u nm.. ..n. «nunpvo.ou new. .mv. «auunmo.o. new. man. ammoonm.nn mo.. .nu. assumvo.nn own. hue. 00000m0.5. pan. mcu. ...oncm.nn no.. .nu. ammoovo.vn can. con. nnnnoco.o. ohm. non. cc'acmm..n cm.. «cu. aoon.vo.nn son. mam. 0000mc0.n. nan. non. «unpvom.o« ..u. ohm. «nusono.nn ..o. mum. ohssuvo.v. m.n. .mn. smoo.sm.m. 0.". can. nannonO..u one. can” ooooo.o.n. pan. p.n. poo.v~n.o. man. can. annnmno.0u woo. coo. ospspno.u. aha. nmn. nnnoOom.s. .mn. con. «nanvno.o. one. .00. nanomno... v.c. can. wcvmwom.m. can. .nn. «muo.no.o. no». haw nnnnnno.o. mac. n... camcoom.m. an“. «en. ooooo«o.n. .np. o.s. OOOOOno.o nun. .Om. unusoam.v. wen. nun. ...oono.o. cos. cps. poo.auo.o «on. who. «unaumm.n. ann. non.. ....cno.n. «as. can. soooouo.p ..p. mom. oooommm.u. awn. a... soo.vno.v. one. nan. «nunvuo.u awn. «no. ommmOOm... .mn. ave. coco.uo.n. 9.0. new. «unuauo.n ano.. «no.. cmmOnom.o .mo. ..5. «cva.uo.u. .oo. ano.. poc.m.o.v oou.. «o«.. c'cocom.m pow. coa. .....uo... «50.. on... anaco.o.n .no.. '09.. unusaom.c sno.. up... ooooc~o.o. .o... no«.. anonn.o.u o...« cm..n oooommm.a o.n.. an... opppp.o.m vo«.. o.n.. «unnu.o.. nus.“ ems." mammooo.n nou.. no... annom.o.o non.. on... nannooo. m.p.n omo.n “fiancee." nae." co..u ...on_o.p nvn.. opn.. ooomhoo.. vmc.~ ".0." nannn.o.o «as.. .m».. nnnnooo. vma.u oo..n aeoon.o.n noo.. «0a.. on 4¢_u» cvva..c.v own.“ p««.« ......o.n pso.« pvn.u cu 4¢.¢» ooooo.o.n onu.« oco.« ooooo.o.. occ.n o.u.n nannOOO. 0.0.9 '00.. omnmm.m.nn ohm. mo“. ..zwwm.mn mwn. hum. ccvm.nm..n man. paw. swoo.hm.vu .ha. sou” ooooo~m.0n a.n. cvu onhumpo.nu can. can. b“ 4