Diffusion of Heptane in Polyethylene Vinyl Acetate: Modelisation and Experimentation

Similar documents
J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

Ethanol Diffusion in Polyethylene Vinyl Acetate Modelling and Experimentation

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

Solution in semi infinite diffusion couples (error function analysis)

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

ELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION

Numerical Simulation of the Dispersion of a Plume of Exhaust Gases from Diesel and Petrol Engine Vehicles

Density Matrix Description of NMR BCMB/CHEM 8190

Polymerization Technology Laboratory Course

Density Matrix Description of NMR BCMB/CHEM 8190

Notes on the stability of dynamic systems and the use of Eigen Values.

On One Analytic Method of. Constructing Program Controls

( ) () we define the interaction representation by the unitary transformation () = ()

Study of a CE Mechanism in Cyclic Chronopotentiometry with Spherical Electrodes

2.1 Constitutive Theory

Nonequilibrium models for a multi component reactive distillation column

Motion in Two Dimensions

FI 3103 Quantum Physics

Transcription: Messenger RNA, mrna, is produced and transported to Ribosomes

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

P R = P 0. The system is shown on the next figure:

Anisotropic Behaviors and Its Application on Sheet Metal Stamping Processes

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

Chapter Lagrangian Interpolation

A Paper presentation on. Department of Hydrology, Indian Institute of Technology, Roorkee

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

TSS = SST + SSE An orthogonal partition of the total SS

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Let s treat the problem of the response of a system to an applied external force. Again,

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Response of MDOF systems

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2)

Linear Response Theory: The connection between QFT and experiments

Robustness Experiments with Two Variance Components

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION

Multi-Fuel and Mixed-Mode IC Engine Combustion Simulation with a Detailed Chemistry Based Progress Variable Library Approach

[Link to MIT-Lab 6P.1 goes here.] After completing the lab, fill in the following blanks: Numerical. Simulation s Calculations

Math 128b Project. Jude Yuen

Numerical simulation of a solar chimney power plant in the southern region of Iran

Chapter 5. Circuit Theorems

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

FTCS Solution to the Heat Equation

On computing differential transform of nonlinear non-autonomous functions and its applications

A General Lagrangian Approach to Simulate Pollutant Dispersion in Atmosphere for Low-wind Condition

Real time processing with low cost uncooled plane array IR camera-application to flash nondestructive

PHYS 1443 Section 001 Lecture #4

Graduate Macroeconomics 2 Problem set 5. - Solutions

Molecular Dynamics Simulation Study forgtransport Properties of Diatomic Liquids

T q (heat generation) Figure 0-1: Slab heated with constant source 2 = q k

Mechanics Physics 151

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

Spurious oscillations and conservation errors in interface-capturing schemes

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

First-order piecewise-linear dynamic circuits

Mechanics Physics 151

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

Lecture 9: Dynamic Properties

II. Light is a Ray (Geometrical Optics)

Implementation of Quantized State Systems in MATLAB/Simulink

MOLDI : A Fluid Permeation Model to Calculate the Annulus Composition in Flexible Pipes

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

MONTE CARLO ALGORITHM FOR CLASPING SEARCH AND NEUTRON LEAKAGE

Cubic Bezier Homotopy Function for Solving Exponential Equations

EVALUATION OF FORCE COEFFICIENTS FOR A 2-D ANGLE SECTION USING REALIZABLE k-ε TURBULENCE MODEL

Gear System Time-varying Reliability Analysis Based on Elastomer Dynamics

Variants of Pegasos. December 11, 2009

A Simulation Based Optimal Control System For Water Resources

Relative Intermolecular Orientation Probed via Molecular Heat Transport

NUMERICAL SIMULATION AND EXPERIMENTAL INVESTIGATION FOR INDOOR AIR ENVIRONMENT IN AN OFFICE ROOM

PHYS 705: Classical Mechanics. Canonical Transformation

2/20/2013. EE 101 Midterm 2 Review

Existence of Time Periodic Solutions for the Ginzburg-Landau Equations. model of superconductivity

ANALYSIS OF THE DUAL PHASE LAG BIO-HEAT TRANSFER EQUATION WITH CONSTANT AND TIME-DEPENDENT HEAT FLUX CONDITIONS ON SKIN SURFACE

Pendulum Dynamics. = Ft tangential direction (2) radial direction (1)

Scattering at an Interface: Oblique Incidence

Chapter 6 Conservation Equations for Multiphase- Multicomponent Flow Through Porous Media

EFFECT OF HEAT FLUX RATIO FROM BOTH SIDE-WALLS ON THERMAL- FLUID FLOW IN CHANNEL

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Volatility Interpolation

Part II CONTINUOUS TIME STOCHASTIC PROCESSES

5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

Computing Relevance, Similarity: The Vector Space Model

Transient Response in Electric Circuits

THEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that

Advanced Machine Learning & Perception

On elements with index of the form 2 a 3 b in a parametric family of biquadratic elds

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

FRACTIONAL OPTICAL SOLITARY WAVE SOLUTIONS OF THE HIGHER-ORDER NONLINEAR SCHRÖDINGER EQUATION

Comparison of Differences between Power Means 1

Bayesian Inference of the GARCH model with Rational Errors

Transcription:

IOSR Journal of Appled hemsry (IOSR-JA) e-issn: 78-5736.Volume 7, Issue 6 Ver. I. (Jun. 4), PP 8-86 Dffuson of Hepane n Polyehylene Vnyl Aceae: odelsaon and Expermenaon Rachd Aman *, Façal oubarak, hammed El Koual, ohammed Talb, Abdlhak El Brouz, SamaYousf, ona Bouamran / Research Laboraory of Analycal hemsry and Physcal hemsry of aerals Scences Faculy Ben 'sk, orocco / Research Laboraory of Physcal hemsry and Elecrochemsry Faculy of Scence and Technology ohammeda. orocco Absrac: Ths arcle descrbes he mass ransfer n food packagng, where he plasc s brough no conac wh food, ceran ransfers mplemened accordng o a complex process. The mehod we used n hs work, we coupled modelng wh expermenaon n order o undersand he packagng plasc ouch nerface ransfers / food. The numercal model akes no accoun he expermenal sudy. The knecs of dffuson ransfer and dffusvy of lqud smulaor were deermned. The analycal and numercal model has been developed whch ams o gve beer nformaon on he concenraon of lqud smulaor nsde he package (PEVA) and smulae a few hours ransfers lasng n realy several monhs. Keywords: Packagng maeral, Dffuson, odelsaon, Fne dfference mehod. I. Inroducon Food conac maerals play an essenal role n he food conservaon and proecon. They have also ncessan markeng funcon and have o be recoverable. When he polymer s n conac wh he packagng, some ransfers ake place, hs conac nfluence on packaged producs as well as on packagng maerals. The goal of hs arcle s o undersand hese ransfers, by couplng modelng wh expermenaon and onsruc a mahemacal model and a numercal model able o descrbe he dffuson process. II. aerals and mehods The sudy he lqud dffuson no a polymer smulaor s based on he followng smplfyng assumpons [3]: The dsrbuon s n accordance wh fck s laws. The dffuson coeffcen s ndependen of concenraon. The dffuson n he sphere s radale. The polymer was sphercal n shape wh a consan radus, as he amoun of ehanol s very small. The chemcal does no evaporae on surface. A. Analycal processng[] [] [5] When he dffuson s radal, Fck s second law expressng he dffuson equaon under ransen condons s n he general form... D r. ( eq) r r r Where s he concenraon a me and a a dsance r from he cener of he sphere. When he dffusvy D s consan, he dffuson equaon akes he form of: D.. ( eq) r r r Analycal soluons can be obaned when he dffusvy s consan. Problems wh a concenraon-dependen dffusvy need numercal mehods. By pung: U. r( eq3) Equaon () becomes: U U D. ( eq4) r 8 Page

Dffuson of Hepane n Polyehylene Vnyl Aceae: odelsaon and Expermenaon Equaon (4) s smlar o he equaon obaned for dffuson n one dmenson hrough he plane shee. ase sudy: The nal dsrbuon n he sphere s consan:. In hs case: Inal Boundary r R U (eq5) r U r R U R ( eq 6) The soluon of he equaon of dffuson becomes: n r, R n r n² ²..sn exp Deq7. r n n R R² Ths soluon n hs form canno be used for he cener of he sphere (r ). As he lm of, when Equaon (7) for he cener of he sphere s reduced o: r,. n n² ² R² n.exp D eq8 B- Analycal model The oal amoun of dffusng subsance gong no or leavng he sphere s gven by negrang Fck s frs law accordng o me. By consderng: And: We can oban: 6 n..exp D ( eq) n R Anoher expresson of he soluon of he equaon of dffuson n he sphere (equaon 5) s gven by: r, R n R r n R r. erfc erfc eq r ( D) D The knecs for he maer ranspored s: 6 D nr D. erfc 3 ( eq) R D R ase of shor mes: Equaon (8) s very useful for shor mes because can be reduced o: 6 D R sn x x n 6 4 R 3 ( eq3) 3. D. r x rr. d( eq9) - Numercal model - Fne dfference mehod-[4] [5] Analycal soluons can be obaned when he dffusvy s consan. Problems wh a concenraondependen dffusvy need numercal mehods. In hs case, he problem mus be solved by usng he numercal Fne dfference mehod. ase sudy: The amoun of he remanng ehanol whn he sphere. The prncple of he mehod s as follows: he sphere s dvded no N sphercal membranes of consan hckness r and he maer balance s evaluaed whn each membrane durng he ncremen of me. 83 Page

Dffuson of Hepane n Polyehylene Vnyl Aceae: odelsaon and Expermenaon The amoun of chemcal locaed n he sphere a me s obaned by negrang wh respec o space he concenraon of he chemcal a hs me: R 4 r. r,. dr( eq4) Ths expresson can be rewren usng fne dfference mehod. 3 9 3 4 r. n n n n n. n ( eq5 4 8 8 ) D- Expermenal procedure The maeral used s polyehylene vnyl aceae (also known as PEVA) s he copolymer of ehylene and vnyl aceae. And our produc smulaor used s hepane or n-hepane s he sragh-chan alkane wh he chemcal formula H 3 (H ) 5 H 3 or 7H6 and hs densy s:.6795 g ml. Hepane s a colourless lqud and odor s perolc. onacng: The conacng sample of polyehylene vnyl aceae s carred ou wh hepane a 5. Durng he conac, we measured he specfc mass of he sphere each me o sudy he evoluon of he mass sphere. III. Resuls And Dscussons The percenage of hepane mass varaon nsde our plasc sphere (fg.) s gven by he followng equaon: m m m ( eq6) m Effec of The dffuson coeffcen s gven by hs relaon: En cm²/s:. R D 6 6 ( eq7) D 4,4. -6 cm²/s Fgure - shows he varaon of he amoun of hepane n he PEVA wh me, we noe ha he hepane n polyehylene vnyl aceae mass ncreases wh he conac me. Unl equlbrum reached. Fgure 3- show he amoun of hepane n he maeral accordng o smulaon me for each soluon : analycal soluon, numercal soluon and experence. We noce from hs fgure ha hree soluons gve he same varaon n he amoun of hepane absorbed as a funcon of me. So we concluded ha he model s valdaed. Fgure 4 - show he profle of he concenraon of hepane afer every mn. a- The hepane concenraon whn he sphere s low, n conras o surfaces. b- Ths concenraon profle s amed o gve good nformaon on hepane nsde PEVA, s for each pon n our sample we can easly deermne s concenraon. IV. Fgures Fg.- Schema of he crcular cross secon of a sphere of radus R. 84 Page

asse m en % Dffuson of Hepane n Polyehylene Vnyl Aceae: odelsaon and Expermenaon 9 8 7 6 5 4 3 4 6 8 emps (mn) Fg.- Varaon of he amoun of hepane n he PEVA Seres Le profl de la dffuson d hepane dans le PEVA 3 x -3 Hsore de la masse dans le maérau.5.5 Soluon numérque Expérence Soluon analyque...3.4.6.7.8.9 Temps de smulaon (h) Fg.3- The amoun of hepane n he maeral as a funcon of me for he smulaon of hree soluons. Fg.4- The profle of he concenraon of hepane afer every mn 85 Page

Dffuson of Hepane n Polyehylene Vnyl Aceae: odelsaon and Expermenaon V. oncluson Through hs work, we conrbued o he sudy and developmen of new mehods couplng expermens wh modelng o undersand he behavor of plasc packagng n conac wh he food producs. The sudy was conduced by weghed followng he evoluon of he mass ransferred over me. The polyehylene vnyl aceae conac wh hepane (consdered smulang agen) a a emperaure of 5 showed ha he amoun ofhepane n he polymer ncreases wh me. odel valdaon was made by comparng he heorecal resuls wh expermenal resuls. The resulng profle gves beer nformaon on concenraons of hepane nsde he package polyehylene vnyl aceae. Ths model we allow a few hours o smulae mass ransfer n realy lasng several monhs. References []. Tocc, E.; Gugluzza, A.; De Lorenzo, L.; acchone,.; De Luca, G.; Drol, E. J. embr. Sc. 8, 33, 36-37. []. Hölck, O.; Heuchel,.; Böhnng,.; Hofmann, D. J. Polym. Sc., Par B: Polym. Phys. 8, 46, 59-7.[3]: Vergnaud, «Dffuson of a chemcal hrough he lqud locaed n a polymer: odelng and expermens», polymer esng 3 (994) 35-45 [3]. ark Geoghegan, 3, «Weng a polymer surfaces and nerfaces» Prog. Polym. Sc. 8,6 3 [4]. ROSA I.D., VERGNAUD J.., 997, «Transfer of conamnan no sold from a bole made of blayer polymers wh a recycled and a vrgn layer - effec of he hcknesses of hese polymer layers», Journal of Appled Polymer Scence, Vol. 66, 9-3. 86 Page