Coupled Nonlinear Differential Equations

  • Thread starter Thomas_W
  • Start date
3
0
Hey,

I need your help to solve the following set of coupled differential equations numerically.

dn(t,z)/dt=I^5(t,z)+I(t,z)*n(t,z)

dI(t,z)/dz=I^5(t,z)-α(n(t,z))*I(t,z)

where I(t,0)=I0*exp(-4ln2(t/Δt)^2) and n(t,0)=0 and n(-certrain time,z)=0. Some constant parameters I did not show here.
α(n(t,z)) is just a parameter which depends linear on n(t,z)

I tried to solve it with Mathematica and NDsolve.. since it works fine if I just solve the first equation without the z dependence. But Mathematica seems to be unhappy with the boundary conditions.

Thanks!

Thomas
 

haruspex

Science Advisor
Homework Helper
Insights Author
Gold Member
2018 Award
31,964
4,741
There doesn't seem to be enough information. There's nothing to set how n varies with z or I with t. You said it was ok if you just took the first equation and treated z as constant, but then you have two unknown dependent variables, n and I, and only 1 equation.
 
3
0
Thanks for your fast response.

Yeah I know, but what I solved was
dn(t)/dt=I^5(t)+I(t)*n(t)
with a time dependent I(t).

n depends on z only due to the fact that I is depending on z, there is no direct dependence. dI(t,z)/dt I only can give for the case I(t,0).. so I just adapt these coupled equations from a paper.

Thank you

Thomas
 
81
1
Hey,

I need your help to solve the following set of coupled differential equations numerically.

dn(t,z)/dt=I^5(t,z)+I(t,z)*n(t,z)

dI(t,z)/dz=I^5(t,z)-α(n(t,z))*I(t,z)

where I(t,0)=I0*exp(-4ln2(t/Δt)^2) and n(t,0)=0 and n(-certrain time,z)=0. Some constant parameters I did not show here.
α(n(t,z)) is just a parameter which depends linear on n(t,z)

I tried to solve it with Mathematica and NDsolve.. since it works fine if I just solve the first equation without the z dependence. But Mathematica seems to be unhappy with the boundary conditions.

Thanks!

Thomas
This is why mathematica is not solving thsi.

"where I(t,0)=I0*exp(-4ln2(t/Δt)^2)"

are you sure this is so, what if Δt=0? then what?
 
3
0
This is why mathematica is not solving thsi.

"where I(t,0)=I0*exp(-4ln2(t/Δt)^2)"

are you sure this is so, what if Δt=0? then what?
Δt is not zero. so I is just a gaussian distribution around Δt
 

hunt_mat

Homework Helper
1,718
16
Can't the method of characteristics be used here?
 

Related Threads for: Coupled Nonlinear Differential Equations

Replies
1
Views
6K
Replies
13
Views
2K
Replies
2
Views
1K
Replies
1
Views
3K
  • Last Post
Replies
5
Views
2K
  • Last Post
Replies
2
Views
3K
  • Last Post
Replies
10
Views
3K
Replies
0
Views
2K

Hot Threads

Top