How to Define Temperature in a MATLAB Mass Transfer Simulation?

[Meaningless]

Hi Guys,

I'm trying to do a simple first-order, 1d (x-direction) simulation of mass transfer with the finite element toolbox in MATLAB for the following closed set of equations:

(Please note: "X" is not to be confused with the lowercase spatial dimension "x")

However, every time I cannot define the temperature (at least with exponential and erfc) without an error. Does anyone know if there is a proper syntax for this?

I appreciate your help.

 

[jedishrfu]

I assume t is an array of numbers, if so then you need to the ".*"

https://help.scilab.org/doc/5.5.2/en_US/m2sci_elementwise_multiplication.html

as this is a very common Matlab error and especially because you made it more difficult by not saying what error it was and/or showing your code.

 

[Meaningless]

This is the syntax I directly entered into the toolbox (I have also directly replaced "P" with the T-dependent expression). I also set x = X for the calculation. I am pretty sure one or two of my variables is not coded properly, particularly the temperature, T.

X' = (K/TX_t)*((10^(A-(B/TX_t + C)))^2)/X_t

TX_t = T0 + C*(t*exp(-X^2/t) + (X/2)*erfc(X/t))

 

[jedishrfu]

Okay, I'd start with building up the expression until it fails so then you have an even simpler failure to show us.

Perhaps:

t=(1:0.1:11)           // create t array of (1,1.1,1.2,1.3 ... 10)
x=(0:0.1:10)         // create X array of (0,0.1,0.2,0.3 ... 10)      same number of elements as t

y = erfc(x./t)          // notice the use of ./ because here you want to do an element by element /