rss
Update 8/10
Charge Exchange
Postprocessing done on the collisional ionization runs. Solar wind run is going well:
MHD
+----------------------+------------+--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+ | | | | Strong Field (4.5x10^1^ G) | Weak Field (4.5x10^-4^ G) | +======================+============+==============+============================================================================================================+================================================================================================================================================================================+ | | | No wind | Immediate pressure protections, but they slowly decrease. By the time the run ended (frame 10), none left. | Runs well. Slight acceleration, moving at max of ~5 km/s at 21 frames. Possibly from boundary condition? | | | No gravity +--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+ | | | Stellar wind | | | | Non-corotating frame +------------+--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+ | | | No wind | Pressure protections all over the place after a few frames, but they also disappear just before frame 10. | Runs well. Accelerates toward star, bending field lines. 20 frames. | | | Gravity +--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+ | | | Stellar wind | | | +----------------------+------------+--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+ | | | No wind | Pressure protections after a few steps, still many at frame 9.8 (furthest it got). | Runs well. Accelerates as expected, fixed at left & top because of boundary condition. 20 frames. | | | No gravity +--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+ | | | Stellar wind | | | | Corotating frame +------------+--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+ | | | No wind | Pressure protections after a few steps. Only 1 frame. | Runs well. Accelerates outward (I think as expected? Orbital velocity of planet is higher than Keplerian velocity of ambient). 4 frames (for some reason only ran on 4 cores). | | | Gravity +--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+ | | | Stellar wind | | | +----------------------+------------+--------------+------------------------------------------------------------------------------------------------------------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+
There still appears to be a problem with the stellar wind being injected. I'll check it without fields to make sure I can get it working properly. Also, the strong wind runs are looking suspiciously similar. I may run clean runs of everything for comparison purposes.
Run 1
Here are the last pressure protections from the simplest run:
q= 0.73E-04 -0.11E-01 0.37E-03 0.12E-01 0.39E+01 0.79E+00 -0.25E-01 -0.18E+01 0.00E+00 0.36E-04 Info%q= 0.73E-04 -0.11E-01 0.37E-03 0.12E-01 0.39E+01 0.79E+00 -0.25E-01 -0.18E+01 0.00E+00 0.36E-04 pos= 0.51875000E+02 0.11250000E+01 0.11125000E+02 i,j,k,mx,my,mz 20 45 25 20 53 20 Info%level 0 tnow= 0.10220561E+01 Temp is -17241.5756049894 q= 0.73E-04 -0.11E-01 0.37E-03 0.12E-01 0.39E+01 0.79E+00 -0.25E-01 -0.18E+01 0.00E+00 0.36E-04 Info%q= 0.73E-04 -0.11E-01 0.37E-03 0.12E-01 0.39E+01 0.79E+00 -0.25E-01 -0.18E+01 0.00E+00 0.36E-04 pos= 0.51875000E+02 0.11250000E+01 0.11125000E+02 i,j,k,mx,my,mz 0 45 25 20 53 20 Info%level 0 tnow= 0.10220561E+01 Temp is -17241.5756071120 Pressure protection encountered at 20 45 25 Pressure protection encountered at 0 45 25
Comments
No comments.