Changes between Version 1 and Version 2 of u/adebrech/Matlab/ParkerWind


Ignore:
Timestamp:
09/08/17 10:48:52 (7 years ago)
Author:
adebrech
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/adebrech/Matlab/ParkerWind

    v1 v2  
    5757First attempted to solve equation in Mathematica - solution appears correct up to the sonic radius, but the incorrect solution (v -> 0) appears for r > r,,s,,. Next attempted to plot approximations for r << r,,s,, and r >> r,,s,, in Matlab, but as would be expected, the approximations fail near r = r,,s,,. Numerically solving the equation in Matlab gives the opposite case of Mathematica - correct solution for r > r,,s,,, but not for r < r,,s,,.
    5858
    59 [[Image(attachment:ParkerApprox.jpg:Matlab)]]
     59[[Image(wiki:u/adebrech/Matlab:ParkerApprox.jpg)]]
    6060
    61 [[Image(Matlab/ParkerWrong.jpg)]]
     61[[Image(wiki:u/adebrech/Matlab:ParkerWrong.jpg)]]
    6262
    6363Looking at Jonathan's code for calculating the Parker wind to find the correct approximation, try to implement it in a different way to test understanding. Works for r > r,,s,,, but not r < r,,s,,.
     
    7676Correct & incorrect plots:
    7777
    78 [[Image(parkerwind.jpg)]]
     78[[Image(wiki:u/adebrech/Matlab:parkerwind.jpg)]]
    7979
    80 [[Image(ParkerWrongApprox.jpg)]]
     80[[Image(wiki:u/adebrech/Matlab:ParkerWrongApprox.jpg)]]
    8181
    8282== Numerical Integration of DE ==