How to solve a system of PDAEs with eigenvalue

  • Context: Graduate 
  • Thread starter Thread starter ktsharp
  • Start date Start date
  • Tags Tags
    Eigenvalue System
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
2 replies · 2K views
ktsharp
Messages
8
Reaction score
0
I have the following system of partial differential algebraic equations:

[ tex ] \frac{1}{H_p}\frac{\partial H_p}{\partial t} = - \frac{\partial W_p}{\partial x} - \frac{f1(H_p,c_p,W_p)}{H_p}, [ \tex ]
[tex]\frac{1}{H_p}\frac{\partial}{\partial t}(H_p c_p} = - \frac{\partial}{\partial x}(W_p c_p) - \frac{f2(H_p,W_p,c_p)}{H_p}, [\tex]<br /> [tex]0 = f2(H_p,c_p,W_p) + f3(H_p,c_p,W_p). [\tex]<br /> <br /> with the following conditions:<br /> Hp(x,0) = 5<br /> Wp(x,0) = s1(x)<br /> cp(x,0) = s2(x)<br /> <br /> Hp(0,t) = s3(t)<br /> Wp(0,t) = W0<br /> Wp(L,t) = 0<br /> d(cp)/dx (L,t) = 0<br /> <br /> How can I solve this numerically?[/tex][/tex]
 
Physics news on Phys.org
type "/tex" not "\tex" (or "$$") :wink:
ktsharp said:
I have the following system of partial differential algebraic equations:

[tex]\frac{1}{H_p}\frac{\partial H_p}{\partial t} = - \frac{\partial W_p}{\partial x} - \frac{f1(H_p,c_p,W_p)}{H_p},[/tex]
[tex]\frac{1}{H_p}\frac{\partial}{\partial t}(H_p c_p) = - \frac{\partial}{\partial x}(W_p c_p) - \frac{f2(H_p,W_p,c_p)}{H_p},[/tex]
[tex]0 = f2(H_p,c_p,W_p) + f3(H_p,c_p,W_p).[/tex]
 
  • Like
Likes   Reactions: 1 person
[UPDATED]

I have the following system of partial differential algebraic equations:

[tex]\frac{1}{H}\frac{\partial H}{\partial t} = - \frac{\partial W}{\partial x} - \frac{f_1(H,c,W)}{H},[/tex]
[tex]\frac{1}{H}\frac{\partial}{\partial t}(H c) = - \frac{\partial}{\partial x}(W c) - \frac{f_2(H,W,c)}{H},[/tex]
[tex]0 = f_2(H,c,W) + f_3(H,c,W).[/tex]

with the following conditions:
H(x,0) = 5
W(x,0) = s1(x)
c(x,0) = s2(x)

H(0,t) = s3(t)
W(0,t) = W0
W(L,t) = 0
d(c)/dx (L,t) = 0

How can I solve this numerically?