• 検索結果がありません。

described,inwhich,thesampleisfractionatedintotenortwentyfractionsandthe Thismethodiscomputersimulationwhichisperformedbyassumingthenumberof

N/A
N/A
Protected

Academic year: 2021

シェア "described,inwhich,thesampleisfractionatedintotenortwentyfractionsandthe Thismethodiscomputersimulationwhichisperformedbyassumingthenumberof"

Copied!
11
0
0

読み込み中.... (全文を見る)

全文

(1)

DistributionofDiffUsionCoefncientanditsConcentrationDependence

hy

Shin-ichiroIsHIDA*,NorioSENDA**,KozoNAKAMuRA*andKatsumasaKANEKo*

Synopsis

AnewmethodfordeterminingthedistributionofthediffusioncoefIicientandits

concentrationdependencebyana】ysis『ofthediffusioncurveoftheheterOgeneoussystemls

proposed・

Thismethodiscomputersimulationwhichisperformedbyassumingthenumberof componentandtheweightfractionofeachcomponentandbyselectingthediffUsioncoefIicient andtheconcentrationcoeflicientfrompseudorandomnumbers・

Fromtheresultsoftheapplicationtotheoretica11ysyntheticdiffusioncurves,itbecame Obviousthatthismethodhashighaccuracy.

Introduction

Whenthediffusioncoeflicientofunfractionatedpolymersampleinsolutionhasconcentration dependence,itisdifIiculttodeterminethedistributionofthediffUsionCoefficient,fUrthermore,

themolecularweightdistributionfromdiffusioneXperiment・

Sincetheconcentrationdependenceappearsastheskewnessofthediffusioncurvewhich isrelationShipofconcentrationgradientanddiffusiondistance,wepreviouslyproposedan analyticalmethodtoeliminatetheSkewnessandtoobtainthesynⅢnetricaldiffusioncurve whichshowsonlypolymolecu]arity.')ApplyingthelOgarithmicana1ysis2)andthesimu】taneous method3)forthissymmetrizeddiffusioncurve,thediffusioncoefIicientdistributioncanbe ObtainedtheoreticaHy・Bythesemethods,however,thesampleisfractionatedonlyintothree

orfourfractions・

Inthispaper,anewanalyticalmethodtoobtainthediffusioncoefIicientdistributionis described,inwhich,thesampleisfractionatedintotenortwentyfractionsandthe concentrationcoefIicientforeachdiffusioncoefficientshouldbeobtainedtheoretically.

Theoretical

SymmetrizationoftheDiffusionCurve

Indilutesolutionofhomogeneoussystem,there1ationShipbetweenthedilfusioncoefIicient Dandtheconcentrationccanberepresentedasfollows:

*Dept・oflndustrialChemistry

**NowatldemitsuKosanCo.

(2)

金沢大学工学部紀要7巻1号1973 D-Do(1+COA3°C/CO),

10

(1)

whereDoisthediffusioncoeflicientatinfinitedilution,coistheinitialconcentrationofthe solutionandbisthecharacteristicparameteroftheconcentrationdependence・

IftheaveragediffusioncoeflicientDmobtainedfromeq.(2)isused,eq.(1)canbe transformedintoeq.(3),sinceDmmeansweightaverageva1ue,

(2)

(3)

Dm-加2/2t,

Dm-Do(1+0.5COA).

Intheseequations,籾2isthesecondmomentofthediffusioncurveandノisthediffusiontime・

Now,inthesystemofheterOgeneityinmolecularweight,iftheconcentration,the diffusioncoeflicientandthecharacteristicparameteroftheconcentrationdependenceofeach componentaredesignatedascf,DdandAFz,respectively,D狐isrepresentedbyfollowing equation,

Dm-Zc`、‘(1+0.5cobf)/zc`. (4)

Ifthestandardizedconcentrationgradientandthestandardizeddiffusiondistanceare representedbyYandXtheweightaveragediffusioncoeHicientatanydistancecanbegiven

asfollows:

、=ZYkD`/■Yb,coh`=0.

Consideringtheconcentrationdependence,wehave DZZDK1+COルポ.c/CO)/ZY‘

DmZc`DK1+0.5coAF`)/Zc‘

(5)

(6)

Thisisrewrittenasfollows:

二m/肌十三鵲響 Dm

D瓜

(c/CO). (7)

Usingγdennedas

ZcilDc

、、=γ・DFV=γ・ Zci (8)

eq.(7)istransformedintothefollowingformuIa,

---.五m/皿十譜砦豐 ,mγ ,1 Zc`Di/zc‘ (c/CO), (9)

c/co-rM(.

 ̄。。

(10)

-10-

(3)

IsHIDA・SENDA・NAKAMuRA・KANEKo:DistributionofDiffiusionCoefficientandits11 ConcentrationDependence

Thefirsttermoftheright-handineq.(9)isD/Dmwhichisindependentoftheconcentration andthengivesadilfusioncurvesynmetryaboutc/c0-0.5.Consideringapairofconcentration

(c/CO)1and(c/CO)zwhichsatisfythecondition(c/CO),+(c/CO)2-1andapairofdiffusion coefIicient(D/D、),and(D/DmLwhichcorrespondtoeachConcentration,thefollowing

equationisgiven,

ZYtDzcohf (D/D、ルー(D/D、)ェ (11)

、m・ZYk (c/COルー(c/CO)1

Substitutingthisequationintoeq.(9),thevalueofthefirsttermoftheright-handis determinedcorrespondingtoeachconcentration・

Thistreatmentgivesthecurveof(、/D、)/γvs.c/coinwhichtheconentration

dependenceiseliminated・

Nowthefollowingequationcanbewritten,

I:(D/2M(./‘。)-1 (12)

Thenγ1sgwenasfollows:

●●

-H:(c/・・)。(Dハル(`/`.)ユ(、ハル‘(c/‘。).(13) (c/COルー(c/CO)1

F…h…I…hい…・ロ蓋会儲・mMMhM伽。、…whi…

nottheconcentrationdependencecanbeobtainedbyiterativemethodusingeqs.(14)and

(10),

(M)expI-l:(M)x`x}

Y= (14)

'二(DC/D)抑'-1:(M)jMY)肛

AsthediffusioncurvewithoutconcentrationdependencehastheadditiveprOperties,we getthefollowingequationsfromeq.(9),

`(器)一別

Y=

(zcjDD雑

急・‐三二=と),

(皿D`一弘)exp( (15)

(2元)牝(Zc`)雑

z…ユー('+MWC…(-糸.旦豐L)

一一 D|山

zwl+MM…D乳叩(-舟. z芸=L) (16)

-11-

(4)

金沢大学工学部記要7巻1号1973

12

ThediffusioncoefIicientatanydiffusiondistanceisgivenastheweightaverageandso nextequationisObtained,

ZY;iD`/己YH

(会ルー。‐ zQD`/Zc‘ (17)

TheweightaveragediffusioncoefIicientDw,iseXpressedas

Dw=(D、)‘。脆=。=zQD`/Zc`・

TheconcentrationaveragecoefIicientDH:isrepresentedbythefollowingformula,

(18)

鶚)…-器菫

昨( (19)

andtheaveragediffusioncoefIicientobtainedfromareaofthediffusioncurveisshownas follows:

昨( ZcfDz-斗も ZDI )。 (20)

ThevaluesofDw,D聡and、』Canbedeterminedbythefollowingrelationshipsfromthe analysisofthesymnetrizeddiifUsioncurveandc/covs.D/Dmcurve.

Dw=、、/γ,

、脆=、W/(D/Dm)c'@。=0.5,

(21)

(22)

DA=Dw/2元(Y;…)2. (23)

Thecharacteristicparameteroftheconcentrationdependencelsrelatedtoγbythenext

equation,

ノセザー2(γ-1)/CO・

Inheterogeneoussystem,ルristheaveragedefinedbythefollowingequation’

んヅーヱc‘D`ん‘/二cd、‘.

(24)

(25)

DeterminationoftheDistributionoftheDifTusionCoeflicient

Severalmethodshavebeenproposedfordeterminingthedistributionofthedijrusion coefHcientfromanalysisofthediffusioncurve,butinthesemethods,itmutsbeneedto assumeafUnction・Onthecontrary,thelOgarithmicanalysisandthesimultaneousmethodneed nosuchassumptionandcandirectlyestimatethedistributionfromthediffusioncurve・By thesemethods,thesampleishypotheticallyfractionatedinonly3or4fractions,butlOto20

-12-

(5)

IsHIDA・SENDA・NAKAMuRA・KANEKO:DistributionofDiffiusionCoefficientandits13 ConcentrationDependence

fractionsareobtamedbyournewtheoreticalmethoddescribedbelow・

Assummgthenlmberofcomponent〃(e,9.10or20),weightfractionofeachcomponent c‘(e、g、cユーc2-……..=cユo=0.1orcユーc2=………=c2o=0.05)andD‘/Dwvbyusing pseudorandomnumberS,thetheoreticallysyntheticdiffusioncurveswereObtainedbythe followingequation,

仏zy2慰云ノハ・叩(一命蒜)

AnaverageerrorisdefinedasfolloWS:

(26)

層-1/(Y-Y'm+(Y-W…+(Y-P胸雀十……+(Y-W…(27)

ThegroUpofD`/Dwwhichmakes5thesmallestisthenearesttothetruecomposition ofthesample・RepeatingthemodiHcation,agroupofD`/D”whosesyntheticdiffUsioncurve isverysimilartotheeXperimentaloneisObtained

Similarlyassumlng〃,c`,D`/Dwandco々`frompseudorandomnumbers,theC/c0-,/Dm curvescorrespondingtothosegroupsareObtainedbythefollowingequation,

ZqZq(D蝋ハ)雑(l+c伽./゜。)eXp(-- X2

(金)'一囚q(D`ハ)(,+・…)皿川禰)-払卸 2,`/、〃 (-,蒜D噸γ (28)

AnaverageerrorisdefinedasfoUows:

ど'-1/[(島)-(民)']:…+[(金)-(民)']:_翠十…+[(犬)-(金)']:…

(29)

SimilarlyselectingcoルガfrompseudorandomnumberforeachD`/Dwandchoosmgthe groUpofcoノセzwhichhasthesmallestvalueof§',thisgroUpisthenearesttothetruecohf ofthesample、Repeatingthemodification,agroupofcob`whosesyntheticc/co-D/Dm

curveisverysimilartothatofthesampleisobtamed・

Insuchway,thesampleisfractionatedtheoreticallyin〃fractions,eachofwhichhas arbitraryweightfractioncfandconcentrationcoefIicientcob`.

TheoreticalExample

Theaccuracyofthenewanalyticalmethodischeckedbythreetheoreticalexamples・

ComputerFACOM230-25/35inKANAZAWAUniversityisusedforcalculation・

EachexampleiscomposedoflOfractionsasshowninTablesl,2and3.Thedistribution functionofexamplelisresembledtoGral6n,s・Integraldistributionofexample2islinear.

-13-

(6)

14金沢大学工学部紀要7巻1号1973

ExampIe3hasthedistributionassameasthatofPVAc-acetonesystem.

TablelComponentsoftheoreticalexamplel

Dz Cf

cohi

0.06 0.1 4.07

0.10 0.1 2.55

0.13 01 1.91

0.16 0.1 1.36

0.20 0.1 1.02

0.25 0.1 0.83

0.31 0.1 0.61

0.42 0.1 0.39

0.60 0.1 0.20

1.0 0.1 0.07

Table2Componentsoftheoreticalexample2

Di

Ct

coAFd

1 0.1

10 2 0.1

9 3 0.1

8 4 0.1

7 5 0.1

6 6 0.1

5 7 0.1

4 8 0.1

3 9 0.1

2 lO ql l

Table3Componentsoftheoreticalexample3

Dz

Cf

cohi 1.73 0.1 4.07

2.13 0.1 2.55

2.44 01 1.91

2.84 0.1 1.36

3.27 0.1 1.02

3.72 0.1 0.83

4.44 01 0.61

5.85 0.1 0.39

8.96 01 0.20

13.97 0.1 0.07

Figs、1,2and3showthesyntheticdiffusioncurvesobtainedfromabovetablesandthe symmetrizedone.

0.5

0.4

門0.3

0.2

0.1

-4-3-2-1

OX

1234

Fig.10riginalandsymmetrizeddiffusioncurveofexamplel

-14-

ロワ■

二一|』一

(7)

IsHIDA・SENDA・NAKAMuRA・KANEKo:DistributionofDiffiusionCoefficientandits 15 ConcentrationDependence

0.5

0.4

0.3 門

0.2

0.1

-4-3-2-1

OX

1234

Fig.20riginalandsymmetrizeddiffusioncurveofexample2

0.5

0.4

0.3 門

0.2

0.1

-4-3-2-1 OX 1234

Fig.30iginalandsymmetrizeddiffusioncurveofexample3

-15-

■■■ ̄

ロ‐-し‐  ̄

ログ 一一一一一

(8)

金沢大学工学部紀要7巻1号1973

16

c/coandD/DmofthesediffusioncurvesareshowninFi9.4,5 Therelationshipbetween

and6,respectively.

層ロベロ

/ ノ ノ ノ 1.5

、】

0 0.5 //

c/CD

1.0

ノノ

//

Fig.4c/covs.D/DmCurveforexamplel //

//

冨口へ。 //

1.0 //

/ 0.5

0 0.5

c/CO

1.0

Fig.5c/covs.D/DmCurveforexample2

-16-

(9)

IsHIDA・SENDA・NAKAMuRA・KANEKo:DistributionofDiffiusionCoefficientandits 17 ConcentrationDependence

1.5

ミロヘロ 1.0

0.5

0

00.5LO C/、

Fig.6c/covs.D/DmCurveforexample3

1.0 1.0

宮。》旨二目、国【飼団。]pH 5

口◎垣ロー官芭ロ『飼団①]日

0 0

00.51.0 DiffusionCoefficient

Fi9.7Distributionofdiffusioncoefficientof examplel

(circle:analyzedbythenewmethod,

solidline:theoretical)

0510 DiffUsionCOefficient

Fi9.8Distributionofdiffusioncofficientof example2

(circle:analyzedbythenewmethod,

solidline:theoretical)

-17-

(10)

金沢大学工学部紀要7巻1号1973

18

1.0

〃。

ロ○垣ゴロ眉]の『ロ[飼肖図&pH

CD 0.5

0 510

DiffUsionCoefficient 0 15

Fig.9Distributionofdiffusioncoefficientofexample3 (circle:analyzedbythenewmethod,solidline:theoretical)

ItisobviousfromFigs、7,8and9thatthedistributionsofthediIfusioncoefIicient caIculatedbyournewmethodaregoodagreementwiththegivendistributions・

TheconcentrationcoefIicientversusthediffusioncoefIicientareplottedmFi9.10,11 and12.

壁(ざ ●●

::鴎亀、

0.

0 0.5

DiffusionCoefficient

1.0

Fig.10Aplotofconcentrationcoefficientvs・diffusion coefficientforexamplel

(circle:analyzedbythenewmethod,solidline

:theoretical)

-18-

(11)

lsHIDA・SENDA・NAKAMuRA・KANEKo:DistributionofDiffiusionCoefficientandits

ConcentrationDependence 19

10

曇CQ

0

0510 DiffUsionCoefficient

FigllAplotofconcentrationcoefficientvs・Diffusion coefficientfOrexample2

(circle:analyzedbythenewmethod,solid line:theoretical)

二。。 稗

go 0 gQo

0510l5 DiffUsionCoefficient

Fi9.12Aplotofconcentrationcoefficientvs・diffusioncoefficient forexample3

(circle:analyzedbythenewmethod,solidline:theoretical)

Fromtheseresults,itisshownthatthenewmethodhasthehighaccuracy.

Acknowledgement

TheauthorsareverygratefultothestaffoftheComputerCenterofKanazawa University.

Reference

l)K・Kaneko,S・Ishida,T・Mizuno,Y・KuraandY・Tamura,J・Tamura,J・Chem・SOC・Japan,

PureChemSec.,91,36(1970).

2)K・Kaneko,J・Chem・SOC・Japan,Ind・Che1,.sec.,45,61(1942).

3)Y・Nakayama,MDissertation,KanazawaUniversity(1971).

(RecievedSept、20,1972)

-19-

参照

関連したドキュメント

Construction of Lyapunov functionals is usually used for investigation of stability of hereditary systems which are described by functional differential equations or Volterra

The problem is modelled by the Stefan problem with a modified Gibbs-Thomson law, which includes the anisotropic mean curvature corresponding to a surface energy that depends on

Gosset polytopes are described in section 3: since they have very large groups of symmetry, the technique used in [ST] for regular polytopes is convenient

Instead an elementary random occurrence will be denoted by the variable (though unpredictable) element x of the (now Cartesian) sample space, and a general random variable will

External morphologies of three major edible crustaceans, prawns, crabs, and squillas, are described and compared. Additionally, an example of summary of observation results

Here we purpose, firstly, to establish analogous results for collocation with respect to Chebyshev nodes of first kind (and to compare them with the results of [7]) and, secondly,

Lang, The generalized Hardy operators with kernel and variable integral limits in Banach function spaces, J.. Sinnamon, Mapping properties of integral averaging operators,

Figure 2. The biofilm system is described by three phases: the actual biofilm V 2 , the concentration boundary layer V 1 , and the bulk liquid, which is described in the model by