Changes between Version 1 and Version 2 of u/adebrech/Matlab/ParkerWind
- Timestamp:
- 09/08/17 10:48:52 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/adebrech/Matlab/ParkerWind
v1 v2 57 57 First 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,,. 58 58 59 [[Image( attachment:ParkerApprox.jpg:Matlab)]]59 [[Image(wiki:u/adebrech/Matlab:ParkerApprox.jpg)]] 60 60 61 [[Image( Matlab/ParkerWrong.jpg)]]61 [[Image(wiki:u/adebrech/Matlab:ParkerWrong.jpg)]] 62 62 63 63 Looking 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,,. … … 76 76 Correct & incorrect plots: 77 77 78 [[Image( parkerwind.jpg)]]78 [[Image(wiki:u/adebrech/Matlab:parkerwind.jpg)]] 79 79 80 [[Image( ParkerWrongApprox.jpg)]]80 [[Image(wiki:u/adebrech/Matlab:ParkerWrongApprox.jpg)]] 81 81 82 82 == Numerical Integration of DE ==