create_data.m

雷诺方程的数值解法

function [Data]=create_Data()
h0=evalin('base','h0');
epsilong=evalin('base','epsilong');
R=evalin('base','R');
dx=evalin('base','dx');
dz=evalin('base','dz');
M=evalin('base','M');
N=evalin('base','N');
l=evalin('base','l');
dx1=dx/l;
dz1=dz/l;
for i=1:M
    for j=1:(N+1)
        x=(i-1)*dx;
        z=(j-1)*dz;
        h=h0*(1-epsilong*cos(x/R));
        h1=h/h0;
        a=3*h1^2*epsilong*l*sin(x/R)/R;
        A=h1^3/dx1^2+a/(2*dx1);
        B=h1^3/dx1^2-a/(2*dx1);
        C=-2*h1^3*(1/dx1^2+1/dz1^2);
        D=h1^3/dz1^2;
        E=h1^3/dz1^2;
        Data{i,j}=[x z h a A B C D E];
    end
end
end