Question

只有几个问题，我希望有人能抽出时间回答:)。

如果我们有COUPLED模型示例怎么办：n个独立变量的系统X和n非线性偏微分方程PDEf（X，PDEf（X））相对于取决于X的时间，PDEf（X）（偏微分方程取决于变量X）。你能提一些建议吗？这是一个例子：

假设c是输出或期望的变量。假设r是自变量。偏微分方程看起来像：

∂c/∂t= D * 1 / r +∂c/∂r+ 2（D *（∂^ 2 c）/（∂r^ 2）） d =常数 r = 0：0.1：Rp-Matlab语法，如何在Modelica中表示相同（我使用积分器，但没有工作）？

这是一个代码（不起作用）：

model PDEtest
/* Boundary conditions
1. delta(c)/delta(r)=0  for  r=0
2. delta(c)/delta(r)=-j*d for  r=Rp*/
parameter Real Rp=88*1e-3; // length
parameter Real initialConc=1000;
parameter Real Dp=1e-14;
parameter Integer np=10; // num. of points
Real cp[np](start=fill(initialConc,np));
Modelica.Blocks.Continuous.Integrator r(k=1); // independent x1
Real j;
protected 
parameter Real dr=Rp/np;
parameter Real ts= 0.01; // for using when loop (sample(0,ts) )
algorithm 
j:=sin(time); // this should be indepedent variable like x2
r.u:=dr;
while r.y<=Rp loop
for i in 2:np-1 loop
der(cp[i]):=2*Dp/r.y+(cp[i]-cp[i-1])/dr+2*(Dp*(cp[i+1]-2*cp[i]+cp[i-1])/dr^2);
end for;
if r.y==Rp then
cp[np]:=-j*Dp;
end if;
cp[1]:=if time >=0 then initialConc else initialConc;
end while;
annotation (uses(Modelica(version="3.2")));
end PDEtest;

以下是更多问题：

此代码在OpenModelica 1.8.1中不起作用，也不适用于Dymola 2013demo。我们怎样才能有变量c的连续函数，而不是函数数组？
我们可以将数组cp的值放在combiTable中吗？怎么样？
如果相反“算法”保持“方程式”代码不能成功检查。为什么？在OpenModelica中，错误是：无法展平模型：S。
有没有简化的方法来使用耦合的一组方程（PDE）？我知道Modelica中的PDE库，但我认为它们很复杂。我想写一个函数来解决PDE并在“主模型”中调用这些函数，这样函数的输出就是“c”的连续函数。我不知道如何处理函数数组。
如果我们像在Matlab中“说话”那样，你能给我一些如何理解Modelica语言的建议吗？例如：自变量r的值，我们可以在Matlab中说明，比如r = 0：TimeStep：Rp ...如何在Modelica中做同样的事情？请解释一下“方程式”部分是如何工作的，是否与Matlab有相似之处，是否有必要的财务方法？干杯：）

Answer 1

很难回答你的问题，因为你假设Modelica~Matlab，但事实并非如此。所以我不会评论你的代码，因为它确实是错的。让我给你一个burger equation的示例模型。也许你可以用它作为起点。

model burgereqn
 Real u[N+2](start=u0);
 parameter Real h = 1/(N+1);
 parameter Integer N = 10;
 parameter Real v = 234;
 parameter Real Pi = 3.14159265358979;
 parameter Real u0[N+2]={((sin(2*Pi*x[i]))+0.5*sin(Pi*x[i])) for i in 1:N+2};
 parameter Real x[N+2] = { h*i for i in 1:N+2};
equation
 der(u[1]) = 0;
 for i in 2:N+1 loop
    der(u[i]) = - ((u[i+1]^2-u[i-1]^2)/(4*(x[i+1]-x[i-1])))
                + (v/(x[i+1]-x[i-1])^2)*(u[i+1]-2*u[i]+u[i+1]);
 end for;
 der(u[N+2]) = 0;
end burgereqn;

您的进一步问题：

cp是一个连续变量，数组代表每个离散点。
据我所知，为什么你应该这样做您想要的解决方案变量。
你应该尝试使用几乎总是方程式部分算法部分通常用于函数中。我很漂亮确定你可以用方程表示你的欲望行为。
我不知道那个图书馆，但是pde上的难点就是离散化和解决它自我。你可能会遇到问题虽然通常使用modelica工具解决pde Modelica工具没有针对pdes的专门求解算法。
请考虑进一步参考该问题。你可以从Modelica.org开始。

在modelica中使用PDE的耦合系统

1 个答案: