I want to simulate the following programme on pdepe solver with third type inlet boundary of Heaviside nature with t0 = 140 days and a zero flux ourlet boundary

2 次查看（过去 30 天）

显示更早的评论

Thomas TJOCK-MBAGA 2022-4-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1688744-i-want-to-simulate-the-following-programme-on-pdepe-solver-with-third-type-inlet-boundary-of-heavisi

评论： Sam Chak 2022-4-5

采纳的回答： Torsten

NumericBharati2018concPel.m

fonction [p1, q1, pr, qr]= pdex1bc (x1, u1, xr, ur) P1 =1; q1 = 1; pr = 0.71*(1+0.2*1)*Ur; qr = 1;

7 个评论
显示 5更早的评论隐藏 5更早的评论

Thomas TJOCK-MBAGA 2022-4-5

can't get hold of the pdf version. When I click on the link in any browser even in search sites I just see the HTML version

Sam Chak 2022-4-5

Step 1: Click View PDF.

Step 2: Click the Download button to save the PDF.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Torsten 2022-4-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1688744-i-want-to-simulate-the-following-programme-on-pdepe-solver-with-third-type-inlet-boundary-of-heavisi#answer_935149

编辑：Torsten 2022-4-5

在 MATLAB Online 中打开

pl = ul - C0 * ( 1 + exp(-2*k*t0) )/( 1 + exp(2*k*(t-t0)) );
ql = 0.0;

I suggest you plot the function

f(t) = C0 * ( 1 + exp(-2*k*t0) )/( 1 + exp(2*k*(t-t0)) )

to deduce a suitable value for k.

Are you sure about your setting for f in pdex1pde ?

Note that the flow velocity will be

v = 0.71*(1+0.2*x)

, but that this setting will generate an artificial source term

s = -0.71*0.2*u

in your model.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Torsten 2022-4-5

@Thomas TJOCK-MBAGA Comment moved here:

Yes I'm sure about the expression of f. I have computed many numericals solution with pdepde solver usingba first type inlet boundary with this expression of f and other in the samedi form. They match well with analytical solution using GITT ans other analytical solutions. Now i wanted to plot the same problem but using third type boundary and using also a pulse boundary (Heaviside source).

Torsten 2022-4-5

编辑：Torsten 2022-4-5

在 MATLAB Online 中打开

Ok. Then in your code variables, which boundary condition at the inlet do you want to set ?

0.6*dC/dx = 0.71*(C-f(t))

with f(t) defined as above ?

请先登录，再进行评论。

更多回答（1 个）

Thomas TJOCK-MBAGA 2022-4-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1688744-i-want-to-simulate-the-following-programme-on-pdepe-solver-with-third-type-inlet-boundary-of-heavisi#answer_935245

Yes this is the boundary condition in the inlet that i wanted to set but written as follow -0.6*dC/dx + 0.71*C = C0 And un the ourlet dC/dx = 0.

Now i also wanted to set thr inlet boundary written like this: 0.6*dC/dx + 0.71*C = C0 for 0 < t =< t0 and 0.6*dC/dx + 0.71*C = 0 for t >t0 That is the definition of the pulse boundary.

5 个评论
显示 3更早的评论隐藏 3更早的评论

Torsten 2022-4-5

在 MATLAB Online 中打开

If the diffusive flux is given by

f = 0.6*(1+0.2*x)^(1.5)*DuDx-0.71*(1+0.2*x)*u;

the settings for the left boundary are

pl = 0.71 * C0 * ( 1 + exp(-k*t0) )/( 1 + exp(k*(t-t0)) );
ql = 1.0;

(graph the function f(t) = C0 * ( 1 + exp(-k*t0) )/( 1 + exp(k*(t-t0)) ) with your values for C0 and t0 to get a suitable k)

and for the right boundary

pr = 0.71*1.2*ur;
qr = 1.0;

Sam Chak 2022-4-5

@Torsten

Got to forgive him. If the View PDF and Download buttons can be missed, then it is not surprising that the "Comment" field was overlooked too. However, it is believed that the replies were genuine.

请先登录，再进行评论。

请先登录，再回答此问题。