Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Solar Modelling

  1. May 2, 2004 #1
    I am trying to model a star using the old shooting method technique. I've encountered problems after my very first pass of integration. I know that I need to use a multi-dimensional Newton-Raphson technique and i think that i've not done this correctly. Could anybody please try and explain this technique.
  2. jcsd
  3. Jun 7, 2004 #2


    User Avatar
    Science Advisor

    Start with the equations of stellar structure and the equation of hydrostatic equilibrium, apply the boundary conditions you know: core temperature and pressure, and T=0 and P=0 at r=Rsun. The XYZ composition and the Saha equation gives you the abundance of free ions and the opacity. Use the chain rule to change the independent variable from radius fraction to mass fraction. Begin by assuming constant density. Solve the equation set to obtain a revised density function, and repeat until rho(M) converges. Runga Kutta might give you the power fraction enclosed from M=0 to M=Mi. If there's a name for doing this, I never heard what it was.

    Jerry Abbott
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook