非线性系数型PDE求解问题求助

jiyi0612

mph文件：

非线性系数型PDE求解问题求助

悬赏6个币吧，我就剩这几个了，略表心意。搞了个把月了，没有搞定。多谢了，多多指点下！

已知：数学模型为一个对流-扩散方程，如下：

采用COMSOL中的系数型偏微分方程（或者经典对流扩散方程）求解。

求解：稳态条件下的压强分布，也就是u的分布。

稳态时求解变量 u 以及非线性系数 da、c、β 和 f 分别为：

上式中 ，方程的系数和源项为关于 u、x、y 的非线性项，也就是c、β 和 f 为u、x、y 的函数。其中Q为x,y,u的函数。Λx为常数，即下表中A，Λy为常数0，详见下表：

参数名称	数值		变量名称	变量表达式
L	1.00E-03[m]	区域长度	A	(6*miu*Ux*L/pa)/h0^2	Λx常数
alfa	0.01[rad]	倾斜角度	Hx	（-L）*cos(alfa)/h0+1	常数
h0	5.00E-09[m]	初始高度	D0	pa*h0/miu*sqrt(2*R*T0)	特征雷诺数
miu	2.08E-05	粘度
Ux	25[m/s]	气流速度
pa	1.01E+05pa	初始大气压
R	287.03	气体常数
T0	293[K]	初始温度

方程系数项中的Q和H比较复杂，Q为（x,y,u）的函数，H为(x,y)的函数，我详细解释下：

            

Hsurf	0 .05*(y>0&&y<0.15)+0*(y>0.15&&y<0.65)+0.05*(y>0.65&&y<0.8)
H	(L-x)*(sin(alfa))+Hsurf/h0
D	D0*u*H
Qcon	D/6
Qp	1/6*D+1.0162+0.40134*log(1+1.2477/D)
Q	Qp/Qcon

求解域：为1mm长、0.8mm宽的长方形区域。（附件mph文件）

材料属性：气体为空气。

边界条件：Dirichlet 边界条件方法定义求解域压力场边界，定义空气轴承无量纲的压力边界初始条件：u(0，Y)=1，u(1，Y)=1，u(X，0)=1，u(X，B/L)=1。

问题描述word：文档

mph文件：

jiyi0612

我把帖子修改了下，更加详细了，请各位高手指点下迷津。
就是那个悬赏在哪里设置哦。没有找到。还是说要到一定级别才能设置悬赏呢？

jiyi0612

问题还没有解决呢，还望高人指点呀。

青梅煮酒28

这个论坛没高人指点的话，去中仿社区再发帖试试

liuzongfahjx

我也想知道如何设置，有高人指点吗？