ncombustion.m

来自「LiScNLS is a Matlab application for the 」· M 代码 · 共 9 行

function [Nu,Npu,Npup]=Ncombustion(x,u,up)
D=5.5*1.e8;Yo1=0.21;theta=2.6;T0=0.118;M=6;
Yo=-u+Yo1*exp(M*(x-1));
Yf=-u+1-exp(M*(x-1));
E=exp(-theta./(T0+u));
Nu=D*Yo.*Yf.*E;
Yop=-1;Yfp=-1;Ep=theta./(T0+u).^2.*E;
Npu=D*(-Yf-Yo+Yf.*Yo.*(theta./(T0+u).^2)).*E;
Npup=0*Nu;