Time-dependent thermal conductivity

2 visualizaciones (últimos 30 días)
Giulia Ulpiani
Giulia Ulpiani el 28 de Ag. de 2019
Comentada: Giulia Ulpiani el 29 de Ag. de 2019
Dear Matlab community,
I am using the PDE toolbox to solve a transient thermal model. Among thermal properties, thermal conductivity changes over time. I used the function handle:
k = @(~,state) kin-(kin-kfin)/(duration/state.time);
But when I run the code, it prompts that: "Coefficient evaluation function, "@(~,state) kin-(kin-kfin)/(duration/state.time)", was requested to
calculate coefficients at 2500 locations so should have returned a matrix with 2500 columns. Instead it returned a matrix with 1
columns.
How can I solve this issue?
Thank you
  1 comentario
Bjorn Gustavsson
Bjorn Gustavsson el 28 de Ag. de 2019
If your thermal conductivity varies equally over the entire region independently of temperature or spatial coordinate just use repmat of your eqpression to the desired size, or multiply it with ones(size(x)) or whatever spatial coordinate the pde-solver calls your function k with.
HTH

Iniciar sesión para comentar.

Respuesta aceptada

Ravi Kumar
Ravi Kumar el 28 de Ag. de 2019
Function should return values of k at all points requested by the solver, points are in the first argument. Update your funciton to:
k = @(location,state) ones(size(location.x))*kin-(kin-kfin)/(duration/state.time);
  1 comentario
Giulia Ulpiani
Giulia Ulpiani el 29 de Ag. de 2019
Thank you very much. That solved the problem.

Iniciar sesión para comentar.

Más respuestas (0)

Community Treasure Hunt

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

Start Hunting!

Translated by