Second ODE, initial conditions are zeros at infinity

Click For Summary
SUMMARY

The discussion focuses on solving a second-order ordinary differential equation (ODE) for the temperature profile in the interstellar medium during a phase transition. The equation is given as \(\frac{d^2T}{dx^2} = \frac{f(T)}{T^2} - \frac{1}{T}\) with boundary conditions \(T(x=-\infty) = 1\) and \(\frac{dT}{dx}(x=-\infty) = 0\). The user seeks numerical solutions and initial conditions, while a suggested approach involves multiplying the equation by \(\frac{dT}{dx}\) and using an iterative scheme for recursion equations to integrate from negative infinity to a specified \(x\).

PREREQUISITES
  • Understanding of second-order ordinary differential equations (ODEs)
  • Familiarity with numerical integration techniques
  • Knowledge of boundary value problems in differential equations
  • Experience with Taylor series expansions in mathematical analysis
NEXT STEPS
  • Research numerical methods for solving boundary value problems, such as the shooting method
  • Explore iterative schemes for solving differential equations, focusing on recursion relations
  • Learn about Taylor series expansions and their application in approximating solutions to ODEs
  • Investigate the properties of well-behaved functions in the context of differential equations
USEFUL FOR

Mathematicians, physicists, and engineers working on thermal dynamics in astrophysics, particularly those interested in solving complex differential equations in non-linear systems.

omyojj
Messages
32
Reaction score
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.
 
Physics news on Phys.org


One approach is to multiply by dT/dx this eqaution to get:

1/2 d/dx (dT/dx)^2 = (f(T)/T^2-1/T)dT/dx = (f(T)-T)/T^2 dT/dx

Now perhaps use an iterating scheme of recursion equations, and integrate from minus infinity to some arbitrary x.

Hope I helped somehow.
 


hmm..can you elaborate on the iterating scheme of recursion eq. or share some links I can refer to?
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
2K
Replies
5
Views
2K
  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 0 ·
Replies
0
Views
3K
  • · Replies 7 ·
Replies
7
Views
5K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
Replies
10
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K