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

1 visualizzazione (ultimi 30 giorni)
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 Commenti
Thomas TJOCK-MBAGA
Thomas TJOCK-MBAGA il 5 Apr 2022
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

Accedi per commentare.

Risposta accettata

Torsten
Torsten il 5 Apr 2022
Modificato: Torsten il 5 Apr 2022
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 Commenti
Torsten
Torsten il 5 Apr 2022
@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
Torsten il 5 Apr 2022
Modificato: Torsten il 5 Apr 2022
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 ?

Accedi per commentare.

Più risposte (1)

Thomas TJOCK-MBAGA
Thomas TJOCK-MBAGA il 5 Apr 2022

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 Commenti
Torsten
Torsten il 5 Apr 2022
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
Sam Chak il 5 Apr 2022
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.

Accedi per commentare.

Categorie

Scopri di più su Mathematics in Help Center e File Exchange

Prodotti

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by