Приложение

Исходная математическая модель:

#include <conio.h>

#include <stdio.h>

#include <math.h>

#define n 5

int main()

{ int i;

double L=0.0001,G=0.001,k=5e-5,F=1,dz=0.1,nu=0.142,Tk=3600,h,dt,V,Vb,Vv,m,xr,J;

double x1[n+2],x0[n+2]={0,0,0,0,0,0,0};

double y1[n+2],y0[n+2]={0,0,0,0,0,0,0};

clrscr();

x0[n+1]=0;

y0[0]=0.3;

V=F*dz;

Vb=V*nu;

Vv=V*(1-nu);

dt=10;

h=Tk/dt;

for(m=0;m<h;m++)

{

for(i=1;i<n+1;i++)

{ xr= y0[i]*18;

J=k*(xr-x0[i]);

x1[i]=x0[i]+dt*(L*(x0[i+1]-x0[i])+J)/Vb;

y1[i]=y0[i]+dt*(G*(y0[i-1]-y0[i])-J)/Vv;

}

for(i=1;i<n+1;i++)

{

x0[i]=x1[i];

y0[i]=y1[i];

printf("\n%d %lf %lf",i,x0[i],y0[i]);

} printf("\n\n");

getch();

}return 0;

}

Имитационная модель:

#include <conio.h>

#include <stdio.h>

#include <math.h>

# include <stdlib.h>

#define n 5

double xn(double m, double s, double min, double max)

{

int i;

double xr,x,sum;

m: for(i=0, sum=0;i<6;i++)

sum+=(double)rand()/RAND_MAX;

x=sqrt(2)*s*(sum-3)+m;

if(x<min||x>max) goto m;

return x;

}

int main()

{ int i,N=10000,l,n1,n2;

double L=0.0001,G,k=5e-5,F=1,dz=0.1,nu=0.142,Tk=3600,h,dt,V,Vb,Vv,m,xr,J,P1;

double x1[n+2],x0[n+2]={0,0,0,0,0,0,0};

double y1[n+2],y0[n+2]={0,0,0,0,0,0,0};

clrscr(); randomize ();

n1=0; n2=0;

for (l=0;l<N;l++)

{

G=xn(0.001,0.0001,0.0007,0.0013);

x0[n+1]=0;

y0[0]=0.3;

V=F*dz;

Vb=V*nu;

Vv=V*(1-nu);

dt=10;

h=Tk/dt;

for(m=0;m<h;m++)

{

for(i=1;i<n+1;i++)

{ xr= y0[i]*18;

J=k*(xr-x0[i]);

x1[i]=x0[i]+dt*(L*(x0[i+1]-x0[i])+J)/Vb;

y1[i]=y0[i]+dt*(G*(y0[i-1]-y0[i])-J)/Vv;

}

for(i=1;i<n+1;i++)

{

x0[i]=x1[i];

y0[i]=y1[i];

}

}

//printf("\n g=%g x1=%g yn=%g ",G,x0[1],y0[n]);getch();

if(y0[n]>0.05) n1++;

}

P1=(double)n1/N;

printf ("\n P1=%lf", P1);

getch();

return 0;

}

Список литературы

1. Кафаров В.В., Глебов М.Б. Математическое моделирование основных процессов химических производств: Учеб. пособие для вузов. – М.: Высш. шк., 1991. – 400с.: ил.

2. Александров И.А. Ректификационные и абсорбционные аппараты. Методы расчета и основы конструирования. – М.: Химия, 1978.

3. Дытнерский Ю.И. Основные процессы и аппараты химической технологии.

4. Плановский А.Н., Николаев П.И. Процессы и аппараты химической и нефтехимической технологии. – М.: Химия, 1972. – 496с.