pdepe extract intermediate values

2 次查看（过去 30 天）

显示更早的评论

Guy Elisha 2022-6-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1749180-pdepe-extract-intermediate-values

编辑： Torsten 2022-7-12

采纳的回答： Bill Greene

Hi everyone,

I am solving a pdes using the pdepe function, but I have an issue.

Suppose you want to solve for u, and have your c, f, and s. However, s includes some function g(u)

Suppose that g is a sigmoid function that depends on the neighboring node.

For instance, g(i) = u(i) - sig(u(i-1) - threshold)

That way, the value of u at point i depends on neighboring nodes of previous time steps.

To overcome this issue, I would like to extract all my values of u at each time step into a function, calculate g and return back to the regular itteration. Basically calulate g at each time step, feed it into the giv equation, and proceed in solving.

Is there a way for me to extract all the solutions for u at each time step?

Thank you!

2 个评论
显示无隐藏无

Torsten 2022-6-28

编辑：Torsten 2022-6-28

Is there a way for me to extract all the solutions for u at each time step?

No, only the u-value(s) in the actual grid point x(i).

And the fact that you need the values of u in neighbouring nodes of previous time steps makes your problem a delay pde - even more complicated.

Guy Elisha 2022-6-28

Hi, thanks for the response!

That is what I thought.

Also, sorry for being not clear, I dont need it in the previous time step but of neighboring nodes.

In my s term, I have s = (u * g) where g is a sigmoid function (g = tanh(u(i+1)-threshold)), so that if u ahead of the current node is greater than some threshold, g=1, and if u ahead of the current node is smaller than some threshold, g=-1.

Hence, s depends on u ahead of the current node (u(i+1))

To check that, I need basically all values of u at each time step.

Do you have a suggestion of how I can go around it?

Thanks again!

请先登录，再进行评论。

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

采纳的回答

Bill Greene 2022-6-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1749180-pdepe-extract-intermediate-values#answer_995355

I have written a PDE solver, pde1dm, that has some similarities to pdepe and accepts the same input syntax as pdepe. Most input that is acceptable to pdepe will also work in pde1dm.

However pde1dm has some additional capabilities. If you pass the option vectorized='on' to pde1dm, your pde function will be passed a vector of x-values spanning your complete spatial domain along with u and DuDx matrices at these x values. This allows your pde function to, for example, evaluate an integral involving the solution over the complete spatial domain. (Note that if you use this vectorized option, your pde function is also expected to return c, f, and s matrices corresponding to all values in the x vector.) You may be able to use this capability for your particular application.

You can find more information and download pde1dm here.

19 个评论
显示 17更早的评论隐藏 17更早的评论

Guy Elisha 2022-6-29

Before I even tried using your function on my code, I used one of your examples, for the heat conduction problem. This is where I got the error warning:

Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 2.011891e-17

Main code:

addpath('github_repo/')

%function heatConduction

% transient heat conduction in a rod

n=12; % number of nodes in x

L=1; % length of the bar in m

x = linspace(0,L,n);

t = linspace(0,2000,30); % number of time points for output of results

m=0; % rectangular Cartesian coordinate system

u = pde1dm(m, @heatpde,@heatic,@heatbc,x,t);

figure; plot(t, u(:,end)); grid on;

xlabel('Time'); ylabel('Temperature');

title('Temperature at right end as a function of time');

figure; plot(x, u(end,:)); grid on;

xlabel('x'); ylabel('Temperature');

title('Temperature along the length at final time');

%end

function [c,f,s] = heatpde(x,t,u,DuDx)

rho=8940; % material density

cp=390; % material specific heat

k=385; % material thermal conductivity

c = rho*cp;

f = k*DuDx;

s = 0;

end

function [pl,ql,pr,qr] = heatbc(xl,Tl,xr,Tr,t)

T0=100; % temperature at left end, degrees C

pl = Tl-T0;

ql = 0;

pr = 0;

qr = 1;

end

function u0 = heatic(x)

T0=100; % temperature at left end, degrees C

if x==0

u0=T0;

else

u0=0;

end

Guy Elisha 2022-7-11

Thanks! Yes I saw that. Maybe my question was worded wrong.

In pdepe you can include variables and structs,

ex: sol = pdepe(0,@goveqns,@ICs,@BCs,xspan,tspan,a, b, c).

Is this an option in pde1dm?

Torsten 2022-7-11

编辑：Torsten 2022-7-12

在 MATLAB Online 中打开

Pass extra arguments directly to the functions where they are needed, e.g.

sol = pdepe(0,@(x,t,u,dudx)goveqns(x,t,u,dudx,a,b,c),@ICs,@BCs,xspan,tspan)

That's the usual way in "modern" MATLAB and is also supported by pde1dm:

https://de.mathworks.com/help/optim/ug/passing-extra-parameters.html

请先登录，再进行评论。

类别

Mathematics and Optimization Partial Differential Equation Toolbox Domain-Specific Modeling Structural Mechanics

在 Help Center 和 File Exchange 中查找有关 Structural Mechanics 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by