omyojj
- 32
- 0
second ODE, initial conditions are zeros at infinity!
I want to know the temperature profile of phase transition layer in the interstellar medium.
For stationary solution, the dimensionless differential equation I ended up with is
\frac{d^2T}{dx^2} = \frac{f(T)}{T^2} - \frac{1}{T}
where f(T) is some complicated but well-behaved function.
Boundary conditions are
T(x=-\infty) = 1 ,
\frac{dT}{dx}(x=-\infty) = 0,
However, f(T=1) = 1, one obtains \frac{d^2T}{dx^2}(x=-\infty) = 0
How do I solve this numerically? where should I start the integration? and what should be the initial condition?
Do I need to Taylor expand the differential equation?
Thank you for your attention.
I want to know the temperature profile of phase transition layer in the interstellar medium.
For stationary solution, the dimensionless differential equation I ended up with is
\frac{d^2T}{dx^2} = \frac{f(T)}{T^2} - \frac{1}{T}
where f(T) is some complicated but well-behaved function.
Boundary conditions are
T(x=-\infty) = 1 ,
\frac{dT}{dx}(x=-\infty) = 0,
However, f(T=1) = 1, one obtains \frac{d^2T}{dx^2}(x=-\infty) = 0
How do I solve this numerically? where should I start the integration? and what should be the initial condition?
Do I need to Taylor expand the differential equation?
Thank you for your attention.