Error using pde.TransientThermalResults (line 56) Solution does not correspond to time-dependent PDE.

Question

Ben Sheard 2017-9-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/355555-error-using-pde-transientthermalresults-line-56-solution-does-not-correspond-to-time-dependent-pde

评论： Ben Sheard 2017-9-13

I have a transient thermal script that produces this error

Error using pde.TransientThermalResults (line 56)
Solution does not correspond to time-dependent PDE.
Error in pde.ThermalModel/solve (line 119)
    R = pde.TransientThermalResults(self,u,tlist);
Error in thermalmodel4 (line 328)
result = solve(thermalModelT,tlist);

I believe it may be a problem with my q func?

qFunc = @(region,state) (loadcurrentcurve(state.time)/Vcell).*...
                        ((OCVcurve(SoCcurve(state.time))-Vcurve(SoCcurve(state.time)))+...
                        ((state.u).*dScurve(SoCcurve(state.time))/F));
internalHeatSource(thermalModelT,qFunc);

My q function also uses the curve fitting toolbox

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Alan Weiss 2017-9-7

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/355555-error-using-pde-transientthermalresults-line-56-solution-does-not-correspond-to-time-dependent-pde#answer_280627

How are you including Vcell and F in the function qFunc? Are you passing them in a nested function? Perhaps you need to refresh yourself on how to pass extra parameters. I assume that you defined your other functions such as SoCcurve and OCVcurve correctly.

Alan Weiss

MATLAB mathematical toolbox documentation

11 个评论
显示 9更早的评论隐藏 9更早的评论

Ben Sheard 2017-9-11

编辑：Ben Sheard 2017-9-11

在 MATLAB Online 中打开

for example, a section my qfunc

    function q = qfunc(region,state)
  Di = 18 *(10^-3);      % m,             cell diameter
  R = Di/2;              % m,             cell radius
  He = 65.00 *(10^-3);   % m,             cell height
  Vcell = He * pi*(R)^2; % m^3,           simple volume of cell
  F = 96486;             %C.mol^-1,       Faraday constant
  n = 1;                 % electrons transfered for this chemistry
  nF = n.*F;
  %%Current and Voltage
  currentfunc = @loadcurrentcurve;
  voltagefunc = @loadvoltcurve;
  current = currentfunc(state.time);
  voltage = voltagefunc(state.time);
  %%SoC
  SoCfunc = @loadSoCcurve;
  SoC = SoCfunc(state.time);
      if     SoC > 100
             SoC = 100;
      elseif SoC < 0
             SoC = 0;
      end
  %%OCV
  OCVfunc = @OCVcurve;
  OCV = OCVfunc(SoC);
  %%dS
  dSfunc = load('dScurvefit.mat');
  dS = dSfunc(OCV);
  %%Temperature & Y    
  T = state.u;
  y = region.y;
  %%q, final equation
  q =          y.*...
               (current./Vcell).*...
               ((T.*dS./nF)+...
               voltage-OCV);
  fprintf('cfor time %s completed y point %s',state.time,region.y);
  fprintf('\n')
  end

and my "loadcurrentcurve" function

function I = loadcurrentcurve(time)
%%-- DEFINE FILE --
% define the file with the load data
%the load data must have curent in the cell, cell voltage, cell SoC
%(0-100%), and time in seconds, this is the data that will be simulated in
%the advanced model
loadfilename = 'load.xlsx';  % title of OCV data file in same folder
loadsheet = 1;               % OCV data file sheet number with data on
loadMaxrow = 2000;           % last row with data on, in OCV file
loadMinrow = 2;              % first row with data on, in OCV file
loadIcol = 'A';              % colum with current, in amps
loadtimecol = 'D';           % colum with time
%%-- DEFINE RANGES --
loadMaxrow = num2str(loadMaxrow);                                  % convert row value to string
loadMinrow = num2str(loadMinrow);                                  % convert row value to string
loadIrange = ([loadIcol,loadMinrow,':',loadIcol,loadMaxrow]);      % load current data range
loadtrange = ([loadtimecol,loadMinrow,':',loadtimecol,loadMaxrow]);% load current data range
%%-- READ FILE --
Iload = xlsread(loadfilename,loadsheet,loadIrange);     % read load current from file and assign to variable load current
tlist = xlsread(loadfilename,loadsheet,loadtrange);     % read load time from file and assign to variable tlist
%%-- CREATE FIT --
Fitfunc = @createFit;
CurrentFit = Fitfunc(tlist, Iload);             % fit curve to data time/current data
%%-- FINAL VALUE --
I = CurrentFit(time);

as you can see, this is incredibly slow, and actually the line

"dSfunc = load('dScurvefit.mat'); dS = dSfunc(OCV);"

doesn't work, but this is a different problem. I tried this as the dS curve fit was a particular slow function as it needs to be derived every time.

Hoping you can help, I'm pulling my hair out!

Ben Sheard 2017-9-13

Hi Alan, I'm glad too, now onto my next problem!

Thank you again for your help.

Maybe a section on how to pass extra parameters into the PDE toolbox would be a good idea, especially because the toolbox stops you passing them as arguments and specifically have to be passed as parameters. Maybe a page on the error code or something just saying a link to how to pass parameters with a specific example. (Like the one you typed above)

Thank you

Ben

Ben Sheard 2017-9-13

Hi Alan,

Hoping you can point me in the right direction with the next step of my code.

Thank you again

Ben

请先登录，再进行评论。

Answer 2

Svetlana Pease 2017-9-7

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/355555-error-using-pde-transientthermalresults-line-56-solution-does-not-correspond-to-time-dependent-pde#answer_280662

编辑：Svetlana Pease 2017-9-7

在 MATLAB Online 中打开

Hi Ben, I suggest checking carefully which times qFunc is using vs. the times you provide to the solve function as the tlist argument. Your definition of qFunc does not directly depend on time, as in, for example,

qFunc = @(region,state)state.time*...

My suspicion is, and I can be mistaken here, that qFunc either does not depend on time at all or there is a mismatch between the times used for qFunc and solution times you pass to solve as the tlist argument.

Regards,

Svetlana Pease

Technical Writer, MathWorks Documentation Group

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Ben Sheard 2017-9-8

在 MATLAB Online 中打开

"mismatch between the times used for qFunc and solution times you pass to solve as the tlist argument."

i suspect this is the reason too, especially with the error

Error using pde.TransientThermalResults (line 86)
The number of input solution times is not consistent with the results.
Error in pde.ThermalModel/solve (line 119)
    R = pde.TransientThermalResults(self,u,tlist);
Error in thermalmodel4 (line 332)
result = solve(thermalModelT,tlist);

I just need to figure out why...

No the qfunc doesn't directly depend on time, but the SoC and V both do, therefore require time as an Input. And then qfunc directly depends on these.

I don't see why they wouldn't have the same times, I'll use the debug process Alan suggested and spend some more time on this today

请先登录，再进行评论。