Help for solving differential equations

Ask a Question

Create a Post

Question:Help for solving differential equations

Posted: saaaaber 5 Product: Maple 18

September 01 2016

sbbbr-1.mw sbbbr-1.mw sbbbr-1.mw

S_gmi := 1-S_wmi-S_omi;

(1)

`τ_mafw` := -([K_f*k_rwf*sigma/(`μ_w`*B_w)])(diff(P_wm(x, t), x));

[-K_f(diff(P_wm(x, t), x))*(S_wfi(x, t))(diff(P_wm(x, t), x))*sigma(diff(P_wm(x, t), x))/(`μ_w`(diff(P_wm(x, t), x))*B_w(diff(P_wm(x, t), x)))]

[-K(diff(P_om(x, t), x))*(1-(S_gf(x, t))(diff(P_om(x, t), x))-(S_wfi(x, t))(diff(P_om(x, t), x)))*sigma(diff(P_om(x, t), x))/(`μ_o`(diff(P_om(x, t), x))*B_o(diff(P_om(x, t), x)))]

[-K(diff(P_gm(x, t), x))*(S_gf(x, t))(diff(P_gm(x, t), x))*sigma(diff(P_gm(x, t), x))/(`μ_g`(diff(P_gm(x, t), x))*B_g(diff(P_gm(x, t), x)))]

(2)

Eq1 := diff([K_f*k_rwf*(diff(P_wf(x, t), x))/(`μ_w`*B_w)], x)-q_w+`τ_mafw` = diff(`φ_f`*S_wf(x, t)/B_w, t);

[-K_f(diff(P_wm(x, t), x))*(S_wfi(x, t))(diff(P_wm(x, t), x))*sigma(diff(P_wm(x, t), x))/(`μ_w`(diff(P_wm(x, t), x))*B_w(diff(P_wm(x, t), x)))+K_f*(diff(S_wfi(x, t), x))*(diff(P_wf(x, t), x))/(`μ_w`*B_w)+K_f*S_wfi(x, t)*(diff(diff(P_wf(x, t), x), x))/(`μ_w`*B_w)]-q_w = `φ_f`*(diff(S_wf(x, t), t))/B_w

[-K(diff(P_om(x, t), x))*(1-(S_gf(x, t))(diff(P_om(x, t), x))-(S_wfi(x, t))(diff(P_om(x, t), x)))*sigma(diff(P_om(x, t), x))/(`μ_o`(diff(P_om(x, t), x))*B_o(diff(P_om(x, t), x)))+K_f*(-(diff(S_gf(x, t), x))-(diff(S_wfi(x, t), x)))*(diff(P_of(x, t), x))/(`μ_o`*B_o)+K_f*(1-S_gf(x, t)-S_wfi(x, t))*(diff(diff(P_of(x, t), x), x))/(`μ_o`*B_o)]-q_o = `φ_f`*(diff(S_of(x, t), t))*o/B_

[-K_f(diff(P_wm(x, t), x))*(S_wfi(x, t))(diff(P_wm(x, t), x))*sigma(diff(P_wm(x, t), x))/(`μ_w`(diff(P_wm(x, t), x))*B_w(diff(P_wm(x, t), x)))+4.8*K_*m*(S_wm(x, t)-.15)*(diff(P_wm(x, t), x))*(diff(S_wm(x, t), x))/(`μ_w`*B_w)+2.4*K_*m*(S_wm(x, t)-.15)^2*(diff(diff(P_wm(x, t), x), x))/(`μ_w`*B_w)]-q_w = `φ_m`*(diff(S_wm(x, t), t))/B_w

[-K(diff(P_om(x, t), x))*(1-(S_gf(x, t))(diff(P_om(x, t), x))-(S_wfi(x, t))(diff(P_om(x, t), x)))*sigma(diff(P_om(x, t), x))/(`μ_o`(diff(P_om(x, t), x))*B_o(diff(P_om(x, t), x)))+K_*m*k_rom*(diff(diff(P_om(x, t), x), x))/(`μ_o`*B_o)]-q_o = `φ_`*m*(diff(S_om(x, t), t))/B_o

(3)

IBC := P_gf(x, 0) = P_gfi, P_gm(x, 0) = P_gmi, P_of(x, 0) = P_ofi, P_om(x, 0) = P_omi, P_wf(x, 0) = P_wfi, P_wm(x, 0) = P_wmi, S_gf(x, 0) = S_gfi, S_gm(x, 0) = S_gmi, S_of(x, 0) = S_ofi, S_om(x, 0) = S_omi, S_wf(x, 0) = S_wfi, S_wm(x, 0) = S_wmi;

P_gf(x, 0) = P_gfi, P_gm(x, 0) = P_gmi, P_of(x, 0) = P_ofi, P_om(x, 0) = P_omi, P_wf(x, 0) = P_wfi, P_wm(x, 0) = P_wmi, S_gf(x, 0) = 1-S_wfi-S_ofi, S_gm(x, 0) = 1-S_wmi-S_omi, S_of(x, 0) = S_ofi, S_om(x, 0) = S_omi, S_wf(x, 0) = S_wfi, S_wm(x, 0) = S_wmi

(4)

vars := indets({Eq1, Eq10, Eq11, Eq12, Eq13, Eq2, Eq3, Eq4, Eq5, Eq6, Eq7, Eq8, Eq9}, function(identical(x, t))); pdsolve({Eq1, Eq10, Eq11, Eq12, Eq2, Eq3, Eq4, Eq5, Eq6, Eq7, Eq8, Eq9, IBC}, vars)

${P_gf(x, t), P_gm(x, t), P_of(x, t), P_om(x, t), P_wf(x, t), P_wm(x, t), S_gf(x, t), S_gm(x, t), S_of(x, t), S_om(x, t), S_wf(x, t), S_wfi(x, t), S_wm(x, t)}$

Error, (in pdsolve/sys/info) ambiguous input: the variables {`μ_g`, `μ_o`, `μ_w`, B_g, B_o, B_w, K_f, S_wfi} and the functions {`μ_g`(diff(P_gm(x, t), x)), `μ_o`(diff(P_om(x, t), x)), `μ_w`(diff(P_wm(x, t), x)), B_g(diff(P_gm(x, t), x)), B_o(diff(P_om(x, t), x)), B_w(diff(P_wm(x, t), x)), K_f(diff(P_wm(x, t), x)), S_wfi(x, t), (S_wfi(x, t))(diff(P_om(x, t), x)), (S_wfi(x, t))(diff(P_wm(x, t), x))} cannot both appear in the system

Download sbbbr-1.mw

help me for find the vars

122 views
More...

Contact Author

Subscribe

Generate PDF

Answer

Reply

Share

Branch

Create a new Post in reply

Ask a related Question

Flag

Answer this Question Add a Comment

Share via:

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post

Question: Help for solving differential equations

Question:Help for solving differential equations

122 views
More...

Contact Author

Subscribe

Generate PDF

Answer

Reply

Share

Branch

Create a new Post in reply

Ask a related Question

Flag

Preview

Preview

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

Save this setting as your default sorting preference?
Note: You can change your preference any time in your account settings
Don't show this again

From:
To:

Custom Message (optional):

Ask a Question

Create a Post

Question: Help for solving differential equations

Question:Help for solving differential equations

122 views More... Contact Author Subscribe Generate PDF Answer Reply Share Branch Create a new Post in reply Ask a related Question Flag

Preview

Preview

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

122 views
More...

Contact Author

Subscribe

Generate PDF

Answer

Reply

Share

Branch

Create a new Post in reply

Ask a related Question

Flag