DRYING CHARACTERISTICS OF SWEET ORANGE PEEL

Similar documents
A Mathematical Modeling Study of Hot Air Drying for Some Agricultural Products

The Role Analysis of Superheated Steam Injection to Improve Performance in Thin Heavy Oil Reservoir

Marquette, Frontenac, St. Croix, La Crescent Training Trial

None. Tim Pratt & Oliver Coe. Ronald Duyvestijn. Chris Durnford

Biomass Dryers. Presented By GenCap Technologies Inc. 9/9/2010

Dissecting Proprietary Stormwater Treatment BMPs to Develop Practical Solutions Unbiased Research and Case Studies

Under-Flame Fire Tests on Timber Decks

EXPERIMENTAL INVESTIGATION OF THE EFFECT OF DRYING AIR RELATIVE HUMIDITY IN DRYING KINETICS OF KIWI

SA-SAMS. Guideline for using Physical Resources and Infrastructure. Register Module. September 2009

CMG GardenNotes #631 Tree Placement: Right Plant, Right Place

775 Combination Carbon Monoxide/Combustible Gas Detector. 775 Combination Carbon Monoxide/Combustible Gas Detector

CSI-Thermal Multifamily/Commercial. SRCC OG-100 Incentive Calculator Guide

Minimization of Fire Risks in Waste-to-Energy and Waste Processing Plants Early warning fire detection and automatic extinguishing

Dissecting Proprietary Stormwater Treatment BMPs to Develop Practical Solutions Unbiased Research and Case Studies

Cassava Chip Drying by Using a Small-Scale Hot-Air Microwave Oven

Solar Dryer Construction and Use Manual Version

CITY OF TUCSON, ARIZONA DEPARTMENT OF TRANSPORTATION ENGINEERING DIVISION ACTIVE PRACTICES GUIDELINES

An Experimental Study on Drying of Pistacia Terebinthus in a Fixed Bed Dryer

IEC Functional Safety Assessment. GE OIL & GAS India Private Limited

Influence of temperature and drying air velocity on the kinetics of convective drying in the food industry

manual Polymount Plate Cleaner installation protocol Setting / maintenance compressed air + pump pg. 16 Maintenance table settings pg.

The Right Door for Hope, Recovery and Wellness. Adopted 5/13/02. This procedure shall apply to The Right Door for Hope, Recovery and Wellness.

Study of the Drying Kinetics of Lemon Grass

Event Permit. (Operational Permit 2009 Fire Code, Section 105) Permit No.: Activity location: Date(s) of Event: Times: Contact Person: Telephone:

Cleaning Glazed Porcelain and Ceramic Standard Cleaning Guidelines for Glazed Products

ENGINEERING AIR HANDLING PRODUCT. January 2011 CATALOGUE MSA (A)

CMG GardenNotes #615 Pruning Mature Shade Trees

Study of the Drying Kinetics of Baccaurea angulata Merr. (Belimbing Dayak) Fruit

Drying Kinetics of Chilli Pepper in a force Convection Indirect Solar Drying

Model : SP601-15, 15T, 20, 20T, 30, 30T

CMG GardenNotes #101 IPM and Plant Health Care. Gardening and the Environment. Integrated Pest Management, IPM

Air conditioning can account for up to 50 per cent of your summer electricity bill.

CC-280E Thermoelectric Cigar Cooler

SOUTH METRO FIRE RESCUE AUTHORITY LIFE SAFETY BUREAU SUBMITTAL REQUIREMENTS

CHAPTER SEVENTEEN GUARDRAIL ITEMS

DX100 Field Maintenance Guide

Assessment of the Rules on Heater De-Rating

Design of a Geothermal Energy Dryer for Tea Withering and Drying in Wayang Windu Geothermal Field

Mathematical modelling of thin-layer drying of carrot

Safe Work Method Statement (SWMS) Hand held petrol HEDGE TRIMMER operation

Instruction manual. Version PRETREATmaker

Apartment Manager Fire Safety Manual

Signage 76A, Rani Jhansi Road, Jhandewalan, New Delhi , Mob. :

University of Arkansas for Medical Sciences. Part IV - Fire Prevention & Life Safety Management Plan FY18

CHAPTER SEVEN LANDSCAPE AND IRRIGATION ITEMS

Istruction Manual OP2

AWARENESS ON ELECTRICAL SAFETY. Safety Symposium & Exposition August; ITC Sonar, Kolkata

HOT WORK PERMIT PROGRAM SAMPLE POLICY

YCC205 RFP for Landscape Renewal December 7, Condo Corporation YCC 205 RFP for Landscape Renewal

Begonia boliviensis Copacabana F1

Policy 1.1-a has been omitted for brevity

Design and Fabrication of Experimental Dryer for Studying Agricultural Products

CITY OF MERCER ISLAND DEVELOPMENT SERVICES GROUP

Sensors Monitoring Systems. Crank-Case Protection System (CPS)

OTDRs Bi-Directional Testing

Evaluation and Selection of Thin-layer Models for Drying Kinetics of Apricot (cv.nasiry)

Guidelines Alarm Systems - Remote Services. Document: G1-23-EN Publication: May

PROTECTION OF KITCHEN COOKING EQUIPMENT

Perfekt Glide-24 Perfekt Glide-30 Perfekt Glide-36 Built-in Range Hood Instruction Manual

Mathematical modeling of drying characteristics of chilli in hot air oven and fluidized bed dryers

Range Hood Instruction Manual

GMG GardenNotes #713 Block Style Layout in Raised Bed Vegetable Gardens. Block Style Garden Layout

Sitesh Kumar Pandey1, Richa Soni2, Abhishek Shrivastav3 1,3

Learning Objectives. At the end of this unit, the student will be able to:

7th Edition Global Health & Sustainability

A New Approach for Modeling of Hot Air-microwave Thin Layer Drying of Soybean

Touchless Substation Project. Eric Rosenberger, PE Megan Coyle

PUBLIC WORKS DEPARTMENT. Charles D. Herbertson, Public Works Director/City Engineer

Prepared By: Lisa McAvoy 16/Nov/17 dd/mmm/yy. Approved By: Veronica Harris-McAllister 17/Nov/17 dd/mmm/yy

STORMWATER DESIGN MANUAL

Installation and Operating Instructions

Electric Cooktop Owners Manual

CMG GardenNotes #723 Growing Vegetables in a Hobby Greenhouse

How and Why Mission-critical Cooling Systems Differ from Common Air Conditioners

ISOBUS Diagnostics Functional Overview

PHANTASY ELECTRIC FIRE (TF-1561ERA) FEATURING THE WHISPERGLOW FLAME EFFECT

Electrical devices may only be mounted and connected by electrically skilled persons.

BTC6750GL/ BTC9750GL Cooker Hood

BHC605 60cm touch control ceramic hob

ELIGIBLE MEASURES (For Program fiscal year July 1, 2017 to June 30, 2018)

Hot Work Permit Program

Automatic Fire Detectors

ph REGULATOR INSTALLATION AND MAINTENANCE MANUAL

Characteristics of sunflower seed drying and microwave energy consumption

2 and 4 slice Accents toasters with Pause and Check (certain models only) Please read and keep these instructions

Guidance for Sanitary Appliances for the declaration of Performance (DoP) and CE-marking according to the Construction Products Regulation 305/2011

GO Rail Network Electrification. Transit Project Assessment Process. Public Meetings Round #3 June July 2017

Operating Instructions

locatewater Sensor 4. Configuring the locatewater sensor on a sensorprobe unit

BHI770SS 77 cm Flexi Zone Induction hob with slider control

CMG GardenNotes #264 Irrigation Equipment

Max pressure for your bathing pleasure

Experimental Analysis of Open, Simple and Modified Greenhouse Dryers for Drying Potato Flakes under Forced Convection

I308. Central Park Precinct

Avigilon Control Center 5 System Integration Guide

City of Bethlehem Northside and Southside

DO NOT PRINT THIS PAGE

CMG GardenNotes #113 Diagnosing Root and Soil Disorders On Landscape Trees

Range Hood Instruction Manual

International Journal of Modern Trends in Engineering and Research e-issn No.: , Date: 2-4 July, 2015

Transcription:

DRYIG CHARACTERISTICS OF SWEET ORAGE PEEL Sankalpa, K.B., Ramachandra, C.T., Udaykumar idni, Sharanaguda Hiregudar, Beladhadi, R.V & Dinesha, B.L. Department f Prcessing and Fd Engineering, Cllege f Agricultural Engineering, University f Agricultural Sciences, Raichur 584 104, Karnataka, India. Crrespnding authrs email: sankalpabgwda.7@gmail.cm ABSTRACT The aim f this wrk is t investigate the drying characteristics f sweet range peel in three different drying methds i.e., slar tunnel drying, ht air drying and dehumidified air drying. The drying study shwed that the times taken fr drying f sweet range peel in three different drying methd. The drying data were fitted t 5 thin-layer drying mdels, viz.,ewtn, Page, Hendersn-Pabis, Lgarithmic and Midilli-Kucuk mdels The perfrmances f these mdels were cmpared using the determinatin f cefficient (R ), reduced chi-square (χ ) and rt mean square errr (RMSE) between the bserved and predicted misture ratis. Amng the five tested mdel Midilli-Kucuk mdel described the best fit t the experimental data. KEYWORDS: Dehumidified air Drying, Ht air drying, Drying characteristics, Drying mdels and Sweet range peel. ITRODUCTIO Sweet range (Citrus sinensis Osbeck) belngs t sub family Aurantideae which is categrised under family Rutaceae (Milind and Dev, 01). In sweet range 16-0% is cmpsed f peel and it is main by-prduct btained frm sweet range fruit juice prcessing industries. It is the highly perishable and if nt prcessed the peel becme waste and in turn may becme a pssible surce f envirnmental pllutin. By-prduct recvery frm fruit wastes can imprve the verall ecnmics f the prcessing units and the prblem f envirnmental pllutin als can be reduced cnsiderably (Kumar et al., 011). Arund the wrld, 3100 tnnes f range and ther citrus fruits are annually prcessed ut f which the estimated annual waste is 15.6 millin MT (Djilas, 009). These are mainly used fr animal feeds, due t their high fibre cntent; they culd represent an interesting surce f dietary fibre (Larrauri et al., 1997). Since the fruit byprduct has high misture cntent it need efficient drying methd t cnvert wet sweet range peel int dry frm. Slar tunnel drying is an imprved methd ver the pen yard sun drying since the farmer results in shrter drying time and safer prduct. wadays ht air drying is quite cmmn in many industries due t easing f peratin, lw investment and peratin cst cmpare t ther advanced technique. Dehumidified air drying is ne f the imprved drying methds which dry the sample at relatively less temperature and humidity. Garau et al. (006) cnducted experiments n drying kinetics, mdelling and functinal prperties f range skin, Silva et al. (011) prpsed cmbined apprach t ptimize the drying prcess f flavnid-rich leaves using thin layer mdelling, Pnkham et al. (01) cnducted experiment n thin layer mdelling f cmbined farinfrared radiatin and air drying f a ring shaped pineapple, Galvez et al. (010) determined the effective misture diffusivity and mathematical mdelling f the drying curves fr the live-waste cake. In all the stated studies they cnsidered thin layer drying mdels viz., ewtn, Page, Hendersn and Pabis, lgarithmic and Midilli Kucuk mdels. There are n studies shwing drying characteristics f sweet range peel. S in present study effrts were undertaken t describe the drying characteristics f sweet range peel with respect t three drying (slar tunnel drying, ht air drying and dehumidified air drying) by cnsidering the abve said mdels. MATERIAL AD METHODS Sample preparatin Fresh sweet range (cv. Sathgudi) were selected unifrmly accrding t maturity, clur, size and freshness. The fruits were washed, weighed and peeled by using knife. Then peels were cut int small pieces (0mm 10mm), after that it was washed using ht water at 90 C fr 5 min. After that peel surface misture was remved by using muslin clth then it was fllwed by three drying methds with bed thickens f 5mm. Drying prcess and mathematical mdeling Sweet range peel was dried using three drying methds till cnstant weight is btained. In ht air drying, ht air drying at 60 ±1 C was selected accrding t results f Garau et al. (007) in case f dehumidified air drying 45 ±1 C at 15 ±1% RH was selected which is less than ht air drying temperature. In slar tunnel drying, the average temperature recrded was 5 C when average ambient temperature was 37 C. The mathematical mdels namely, ewtn, Page, Hendersn-Pabis, Lgarithmic and Midilli- Kucuk mdels were selected n their ability t best fit the experimental data. The mdels are; ewtn mdel: = exp (-Kθ)...(1) 171

Drying traits f sweet range peel Page mdel: = exp (-Kθ n)...() Hendersn- Pabis mdel: = a*exp (-Kθ)...(3) Lgarithmic mdel: = a*exp ( Kθ) + c...(4) Midilli-Kucuk mdel: = bθ + a*exp ( Kθ n )...(5) = Misture rati M-M e The misture rati () is dented by...(6) M-Me M e = equilibrium misture cntent, (% d.b.) M = misture cntent at any time, θ (% d.b.) M = initial misture cntent (% d.b.) K, n, a, b and c = drying rate cnstants θ = drying time (min) Drying rate was calculated accrding t methdlgy explained by Chakraverty 1981. Amunt f misture remved Drying Rate...(7) Ttal bne dry weight f sample in gm Time taken 100 The drying parameters f all the mdels were estimated by using Matlab versin 7.0 sftware. The fit quality f the prpsed mdels n the experimental data was evaluated using linear regressin analysis using curve fitting tl in MATLAB. The statistical parameters standard square errr (SSE) and rt mean square errr (RMSE) were calculated emplying the fllwing equatins. 100 P i0 - i=0 RMSE = df 1 SSE =( - ) P p...(8)...(9) p...(10) i=1 ( p ) i1 z...(11) = bserved misture rati P = predicted misture rati df = degrees f freedm =. f data pints z =. f cnstants RESULTS & DISCUSSIO Drying behaviur f sweet range peel Irrespective f drying methd, an increasing trend was bserved in the reductin f misture cntent in beginning f the drying. As the drying prceeded, the lss f misture in sweet range peel decreased with drying time. The reductin in misture cntent, drying rate and misture rati with respect t drying time is given in Table 1. TABLE 1: Final Misture cntent, average drying rate and drying time fr drying f sweet range peel in three drying methds Sl.. Dehumidified air drying Ht air drying slar tunnel drying 1 Drying time, h 3.5 5 1 Final Misture cntent, % d.b. 6.5 8.09 8.73 3 Average drying rate, g f water/min. per 100 gm f bne dry materials 0.19 0.14 0.06 17

Fig. 1 shws the influence f drying air temperature n drying rate with respect t drying time. It can be seen that, the drying prcess mainly cnsisted f three drying perids i.e., heating up, cnstant rate and falling rate perid. Ht air drying at temperature f 60±1 C shwed nly the falling rate perid, which was due t mderate temperature f drying. In ht air drying, the drying rate started frm 0.33 t 0.01 g f water/min. Per 100 g f bne dry materials at 60 ±1 C. Cnstant rate f drying was nt bserved in the drying perid, this might be due t thin layer arrangement and t rapid heating f peel, similar bservatin was bserved by Giri and Prasad (007) during micrwave and ht air drying f mushrm. Drying rate [gm f water/ min 100 gm f bne dry material] 0.6 0.5 0.4 0.3 0. 0.1 0 0 1 3 4 5 6 7 8 9 10 11 1 13 Average drying time (h) FIGURE 1: Drying rate curve fr different drying methds with respect t drying time Abbreviatin: STD, slar tunnel drying; HD, ht air drying; DD and dehumidified drying DD HD ST D In dehumidified air drying prcess at temperature f 45 ±1 C and 15 ±1% RH, nly falling rate perid was detected and the drying rate was bserved frm 0.51 t 0.0 g f water/min. per 100g f bne dry materials, due t lwer RH f drying air. Drying rate curves shwed a fast increase at the beginning f the prcess ( i.e., 0.51 g f water/min. per 100 gm f bne dry materials) due t rapid sample heating and a subsequent decrease f the drying rate ( i.e., 0.0 g f water/min. per 100 gm f bne dry materials). The results f present research are in agreement with previus studies dealing with drying f range peel in micrwave drying (Ghanem et al., 01) and similar result was als bserved fr drying f nin in dehumidified air dryer (Guda et al., 014). TABLE : Cnstants f drying mdels Sl.. Mdel Cnstants Drying methd STD TD DD 1 ewtn k 0.73 0.6084 1.071 Page k 0.375 0.499 1.017 3 Hendersn- Pabis 4 Lgarithmic 5 Midilli- Kucuk n 1.094 1.95 1.406 a 1.05 1.18 1.038 k 0.797 0.634 1.104 a 1.045 1.18 1.09 c -0.0338 0.1661-0.0638 k 0.548 0.4501 0.947 a 0.993 0.990 0.998 b -0.0006-0.0051 0.0019 k 0.344 0.485 1.03 n 1.089 1.58 1.48 The drying rate f sweet range peel varied frm 0.19 g f water/min. per 100 g f bne dry materials in the initial stage f drying t 0.01 g f water/min. per 100 gm f bne dry materials in final stage f drying in slar tunnel drying. In this drying, the drying rate was mainly dependent n varying drying temperature. Here, the temperature varied accrding t the climatic cnditin. Frm the data btained during investigatin, it was bserved that, the cnstant rate perid f drying was absent during the entire perid f drying and the drying tk place under the falling rate perid in all drying methds. Similar type f results were bserved in the case f range peel dried in ht air dryer (Garau et al., 007), micrwave drying f three citrus peel (Ghanem et al., 01), cabinet drying f pmegranate peel (gueira et 173

Drying traits f sweet range peel al., 01) and fr thin layer drying fr apple slices (Mabruk et al., 01). Mathematical mdelling f drying f sweet range peel in different drying methds The drying data btained during thin layer drying was fitted int five different drying mdels viz., ewtn, Page, Hendersn-Pabis, lgarithmic and Midilli-Kucuk mdels. The values f misture rati were determined fr all drying methds. Drying cnstants and estimated values f statistical parameters f all five mdels were given in Tables and 3, respectively. TABLE 3: Estimated values f statistical parameters f ewtn, Page, Hendersn-Pabis, Lgarithmic and Midilli-Kucuk mdels used fr fur drying methds Sl.. Parameter Methd Mdel ewtn Page Hendersn- Pabis Lgarithmic Midilli- Kucuk 1 R HD 0.9803 0.996 0.9831 0.9960 0.997 STD 0.9968 0.9993 0.9976 0.9990 0.9940 DD 0.983 0.999 0.9841 0.9898 0.9993 STD 0.0064 0.0015 0.0049 0.001 0.001 SSE HD 0.000 0.0041 0.0179 0.004 0.003 DD 0.0168 0.0008 0.015 0.0097 0.0007 STD 0.9968 0.9993 0.9975 0.9989 0.9930 3 Adjusted-R HD 0.9803 0.9957 0.9810 0.9949 0.9960 DD 0.983 0.9991 0.9814 0.9858 0.9987 STD 0.0164 0.0080 0.0147 0.0097 0.0077 4 RMSE HD 0.0481 0.05 0.0473 0.045 0.015 DD 0.0490 0.0113 0.0503 0.0440 0.0131 STD -14.553-3.4384-1.0647 1.8173-0.8350 5 P (%) HD -48.3636-1.930-6.7916 9.354-4.350 DD -44.097 18.344-39.9801 30.8935 9.8384 STD.679E-004 61E-005.145E-004 9.47E-005 5.88E-005 6 χ HD.3301E-003 5.143E-005 7.0799E-003 6.090E-004 4.70E-004 DD.3466E-003 1.496E-004.4664E-003.0908E-003 1.654E-004 The Midilli-Kucuk mdel successfully described the relatinship between misture rati and drying time with the highest R value and lwer χ, P, SSE and RMSE values. The experimental and predicted drying curves f dried sweet range peel are given in Fig.. Similar result was fund fr drying f nin (Guda et al. 014). Present result is als similar t the value (R =0.99) f Garavand et al. (011) fr thin layer drying f tmat and they als reprted Midilli-Kucuk mdel as a gd estimatin mdel fr drying prcess. The results btained are als in line with the value f Mirzaee et al. (011), fr thin layer drying f aprict and they reprted Midilli-Kucuk is the best mdel with highest R value f 0.999. 1.1 Misture rati [(M-Me)/M-Me)] 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0. DD actual HD actual STD actual DD predicted HD predicted STD predicted 0.1 0 0 1 3 4 5 6 7 8 9 10 11 1 13 Drying time (h) FIGURE : Experimental and midilli-kucuk mdel predicted misture rati fr different drying methds (Abbreviatin: STD, slar tunnel drying; HD, ht air drying; and DD dehumidified drying) 174

COCLUSIO Amng three drying techniques dehumidified air drying required less drying time (3.5 h) with lwest misture cntent (6.5 % d.b.). Experimental results shw that dehumidified drying is the best methd t dry the fresh sweet range peel. Amng the five drying mdels tested namely, ewtn, Page, Hendersn-Pabis, Lgarithmic and Midilli-Kucuk, the Midilli-Kucuk mdel described the best fit f the experimental data with higher R value and lwest SSE, χ, RMSE and P values. REFERECES Chakraerty, A. (1981) Pst harvest technlgy f cereals, pulses and ilseeds. 3 rd editin. Oxfrd and IBH publishing C. Pvt. Ltd. ew Delhi. Djilas, S. (009) By-prducts f fruits prcessing as a surce f phytchemicals- A review. Chem. Ind. Chem. Eng. Q. 15 (4): 191-0. Galvez, A.V., Miranda, M., Diaz, L.P., Lpez, L., Rdriguez, K. and Scala, K.D. (010) Effective misture diffusivity determinatin and mathematical mdelling f the drying curves f the live-waste cake. Biresurce Technl. 101, 765-770. Garau, M.C., Simal, S., Femenia, A. & Rssell, C. (006) Drying f range skin: drying kinetics mdelling and functinal prperties. J. Fd Engg. 75, 88 95. Garau, M.C., Simal, S., Rssell, C. and Femenia, A. (007) Effect f air-drying temperature n physicchemical prperties f dietary fibre and antixidant capacity f range ( Citrus aurantium v. Canneta) byprducts. J. Fd Chem. 104 (3), 1014-104. Garavand, A.T., Rafieea, S. and Keyhani, A. (011) Mathematical mdeling f thin layer drying kinetics f tmat influence f air dryer cnditins. Int. T. J. Eng. (), 147-160. Ghanem,., Mihubi, D., Kechau,. and Mihubi,. B. (01) Micrwave dehydratin f three citrus peel cultivars: Effect n water and il retentin capacities, clur, shrinkage and ttal phenls cntent. J. Ind. Crps Prd. 40, 167-177. Giri, S.K. and Prasad, S. (007) Drying kinetics and rehydratin characteristics f micrwave-vacuum and cnvective ht-air dried mushrms. J Fd Eng. 78, 51-51. Guda, G.P., Ramachandra, C.T. and idni, U. (014) Dehydratin f Onins with Different Drying Methds. Curr. Trends Tech. Sci. 3, 10-16. Kumar, A.K., arayani, A., Subanthini, D. and Jayakumar, M. (011) Antimicrbial activity and phytchemical analysis f citrus fruit peels-utilizatin f fruit waste. J. Sci. Technl. 3 (6), 5414-541. Larrauri, Rupbrez, P., Brav, L. and Calixt, S.F. (1997) High dietary fibre peels: assciated pwders frm range and lime plyphenls and antixidant capacity, Fd Res. Int. 9 (8), 757-76. Mabruk, S.B., Benali, E. and Queslati, H. (01) Experimental study and numerical mdelling f drying characteristics f apple slices, Fd Biprd. Prcess. 94 (4), 719-78. Milind, P. and Dev, C. (01) Orange: Range f benefits: A review. Int. Res. J. Pharm. 3 (7): 59-63. Mirzaee, E., Rafiee, S. and Keyhani, A. (011) Evaluatin and selectin f thin-layer mdels fr drying kinetics f aprict (cv. asiry). Agric Eng Int: CIGR J. 1 (), 111-116. gueira, R.I., Paim, D.R.F., Crnej, F.E.P., Marian, E. S., Barret, A.S. and Freitas, S.P. (01) Drying kinetics f pmegranate (Punica granatum) peels, Pster presented in internatinal cnference n agricultural engineering held at valencia, Spain during July 8-1, paper. P1765. Pnkham, K., Mees,., Spnrnnarit, S. and Siriamrnpun, S. (01) Mdeling f cmbined far - infrared radiatin and air drying f a ring shapedpineapple with/withut shrinkage. Fd Biprd. prcess. 90, 155 164. Silva, E.M., Silva, J.S., Pena, R.S. and Rgez, H. (011) A cmbined apprach t ptimize the drying prcess f flavnid-rich leaves ( Inga edulis) using experimental design and mathematical mdelling. Fd Biprd prcess. 89, 39 46. 175