THE GOLDEN BALLROOM/BUNKER

====================================================================== FRCMFD v4.0 — TIER-2 MANIFOLD SWEEP ====================================================================== Start time: 2026-06-20 11:57:09 ====================================================================== Mounted at /content/drive [LEDGER] /content/TIER2_SWEEP_20260620_115712/FRCMFD_TIER2_SUMMARY.csv initialized [STEP 1] Output directory: /content/TIER2_SWEEP_20260620_115712 ====================================================================== GRID REFINEMENT PHASE ====================================================================== [RUN] GRID_G64_dx04_S1000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 64x64, dx=0.4, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 283.910943, Steps: 1000 ====================================================================== Step 100/1000: H_drift=1.275464e-03, R=5.854, R/L=0.46 Step 200/1000: H_drift=3.902280e-03, R=5.878, R/L=0.46 Step 300/1000: H_drift=5.639228e-03, R=5.970, R/L=0.47 Step 400/1000: H_drift=5.346345e-03, R=6.076, R/L=0.47 Step 500/1000: H_drift=3.855076e-03, R=6.197, R/L=0.48 Step 600/1000: H_drift=2.807438e-03, R=6.348, R/L=0.50 Step 700/1000: H_drift=2.968603e-03, R=6.562, R/L=0.51 Step 800/1000: H_drift=3.785675e-03, R=6.813, R/L=0.53 Step 900/1000: H_drift=4.304887e-03, R=7.133, R/L=0.56 Final H_drift = 0.004305, R_final = 7.133 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G64_dx04_S1000 complete [RUN] GRID_G128_dx02_S1000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 290.064733, Steps: 1000 ====================================================================== Step 100/1000: H_drift=3.916639e-03, R=5.863, R/L=0.46 Step 200/1000: H_drift=5.321095e-03, R=5.892, R/L=0.46 Step 300/1000: H_drift=2.715310e-03, R=5.967, R/L=0.47 Step 400/1000: H_drift=3.674965e-03, R=6.055, R/L=0.47 Step 500/1000: H_drift=4.111494e-03, R=6.179, R/L=0.48 Step 600/1000: H_drift=3.262380e-03, R=6.359, R/L=0.50 Step 700/1000: H_drift=3.491728e-03, R=6.556, R/L=0.51 Step 800/1000: H_drift=3.931548e-03, R=6.825, R/L=0.53 Step 900/1000: H_drift=3.461845e-03, R=7.142, R/L=0.56 Final H_drift = 0.003462, R_final = 7.142 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G128_dx02_S1000 complete [RUN] GRID_G128_dx04_S1000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.4, L=51.2, R_max=25.6 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 297.917229, Steps: 1000 ====================================================================== Step 100/1000: H_drift=1.179828e-03, R=5.844, R/L=0.23 Step 200/1000: H_drift=3.602807e-03, R=5.908, R/L=0.23 Step 300/1000: H_drift=5.188394e-03, R=5.964, R/L=0.23 Step 400/1000: H_drift=4.892908e-03, R=6.058, R/L=0.24 Step 500/1000: H_drift=3.507888e-03, R=6.194, R/L=0.24 Step 600/1000: H_drift=2.558644e-03, R=6.349, R/L=0.25 Step 700/1000: H_drift=2.730859e-03, R=6.564, R/L=0.26 Step 800/1000: H_drift=3.484885e-03, R=6.817, R/L=0.27 Step 900/1000: H_drift=3.936570e-03, R=7.160, R/L=0.28 Final H_drift = 0.003937, R_final = 7.160 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G128_dx04_S1000 complete [RUN] GRID_G256_dx02_S1000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 256x256, dx=0.2, L=51.2, R_max=25.6 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 302.363978, Steps: 1000 ====================================================================== Step 100/1000: H_drift=3.545015e-03, R=5.858, R/L=0.23 Step 200/1000: H_drift=4.769093e-03, R=5.892, R/L=0.23 Step 300/1000: H_drift=2.470043e-03, R=5.964, R/L=0.23 Step 400/1000: H_drift=3.371096e-03, R=6.060, R/L=0.24 Step 500/1000: H_drift=3.652701e-03, R=6.193, R/L=0.24 Step 600/1000: H_drift=2.949405e-03, R=6.357, R/L=0.25 Step 700/1000: H_drift=3.200302e-03, R=6.567, R/L=0.26 Step 800/1000: H_drift=3.502582e-03, R=6.825, R/L=0.27 Step 900/1000: H_drift=3.115730e-03, R=7.143, R/L=0.28 Final H_drift = 0.003116, R_final = 7.143 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G256_dx02_S1000 complete [RUN] GRID_G64_dx04_S5000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 64x64, dx=0.4, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 283.910943, Steps: 5000 ====================================================================== Step 500/5000: H_drift=3.855076e-03, R=6.197, R/L=0.48 Step 1000/5000: H_drift=4.164080e-03, R=7.691, R/L=0.60 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=3.899958e-03, R=10.243, R/L=0.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=3.513088e-03, R=-1.008, R/L=-0.08 Step 2500/5000: H_drift=3.414713e-03, R=9.443, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=3.597348e-03, R=10.805, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=3.751383e-03, R=12.307, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=3.685533e-03, R=14.964, R/L=1.17 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=3.526571e-03, R=0.000, R/L=0.00 Final H_drift = 0.003569, R_final = 2.054 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G64_dx04_S5000 complete [RUN] GRID_G128_dx02_S5000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 290.064733, Steps: 5000 ====================================================================== Step 500/5000: H_drift=4.111494e-03, R=6.179, R/L=0.48 Step 1000/5000: H_drift=3.401371e-03, R=7.693, R/L=0.60 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=3.602152e-03, R=10.191, R/L=0.80 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=3.676460e-03, R=-0.702, R/L=-0.05 Step 2500/5000: H_drift=3.462290e-03, R=9.412, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=3.680619e-03, R=10.781, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=3.492845e-03, R=12.247, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=3.597672e-03, R=14.983, R/L=1.17 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=3.588153e-03, R=0.000, R/L=0.00 Final H_drift = 0.003639, R_final = 2.049 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G128_dx02_S5000 complete [RUN] GRID_G128_dx04_S5000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.4, L=51.2, R_max=25.6 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 297.917229, Steps: 5000 ====================================================================== Step 500/5000: H_drift=3.507888e-03, R=6.194, R/L=0.24 Step 1000/5000: H_drift=3.778514e-03, R=7.539, R/L=0.29 Step 1500/5000: H_drift=3.559497e-03, R=12.663, R/L=0.49 Step 2000/5000: H_drift=3.226903e-03, R=199.499, R/L=7.79 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=3.136544e-03, R=10.436, R/L=0.41 Step 3000/5000: H_drift=3.294412e-03, R=12.516, R/L=0.49 Step 3500/5000: H_drift=3.437385e-03, R=14.434, R/L=0.56 Step 4000/5000: H_drift=3.383600e-03, R=17.253, R/L=0.67 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4500/5000: H_drift=3.226287e-03, R=25.728, R/L=1.00 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.003227, R_final = 2.041 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G128_dx04_S5000 complete [RUN] GRID_G256_dx02_S5000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 256x256, dx=0.2, L=51.2, R_max=25.6 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 302.363978, Steps: 5000 ====================================================================== Step 500/5000: H_drift=3.652701e-03, R=6.193, R/L=0.24 Step 1000/5000: H_drift=3.113524e-03, R=7.543, R/L=0.29 Step 1500/5000: H_drift=3.221663e-03, R=12.658, R/L=0.49 Step 2000/5000: H_drift=3.321561e-03, R=146.045, R/L=5.70 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=3.141207e-03, R=10.442, R/L=0.41 Step 3000/5000: H_drift=3.298810e-03, R=12.487, R/L=0.49 Step 3500/5000: H_drift=3.185606e-03, R=14.438, R/L=0.56 Step 4000/5000: H_drift=3.227647e-03, R=17.264, R/L=0.67 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4500/5000: H_drift=3.241948e-03, R=25.805, R/L=1.01 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.003252, R_final = 2.042 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G256_dx02_S5000 complete [RUN] GRID_G64_dx04_S10000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 64x64, dx=0.4, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 283.910943, Steps: 10000 ====================================================================== Step 1000/10000: H_drift=4.164080e-03, R=7.691, R/L=0.60 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/10000: H_drift=3.513088e-03, R=-1.008, R/L=-0.08 Step 3000/10000: H_drift=3.597348e-03, R=10.805, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/10000: H_drift=3.685533e-03, R=14.964, R/L=1.17 ⚠️ BOUNDARY WARNING (R > 80%) Step 5000/10000: H_drift=3.485282e-03, R=2.282, R/L=0.18 Step 6000/10000: H_drift=3.654776e-03, R=4.228, R/L=0.33 Step 7000/10000: H_drift=3.462538e-03, R=6.931, R/L=0.54 Step 8000/10000: H_drift=3.656968e-03, R=10.322, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 9000/10000: H_drift=3.657886e-03, R=13.413, R/L=1.05 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.003491, R_final = 15.922 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G64_dx04_S10000 complete [RUN] GRID_G128_dx02_S10000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 290.064733, Steps: 10000 ====================================================================== Step 1000/10000: H_drift=3.401371e-03, R=7.693, R/L=0.60 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/10000: H_drift=3.676460e-03, R=-0.702, R/L=-0.05 Step 3000/10000: H_drift=3.680619e-03, R=10.781, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/10000: H_drift=3.597672e-03, R=14.983, R/L=1.17 ⚠️ BOUNDARY WARNING (R > 80%) Step 5000/10000: H_drift=3.515215e-03, R=2.285, R/L=0.18 Step 6000/10000: H_drift=3.501072e-03, R=4.223, R/L=0.33 Step 7000/10000: H_drift=3.483482e-03, R=6.925, R/L=0.54 Step 8000/10000: H_drift=3.664732e-03, R=10.278, R/L=0.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 9000/10000: H_drift=3.624227e-03, R=13.284, R/L=1.04 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.003563, R_final = 15.929 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G128_dx02_S10000 complete [RUN] GRID_G128_dx04_S10000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.4, L=51.2, R_max=25.6 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 297.917229, Steps: 10000 ====================================================================== Step 1000/10000: H_drift=3.778514e-03, R=7.539, R/L=0.29 Step 2000/10000: H_drift=3.226903e-03, R=199.499, R/L=7.79 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/10000: H_drift=3.294412e-03, R=12.516, R/L=0.49 Step 4000/10000: H_drift=3.383600e-03, R=17.253, R/L=0.67 ⚠️ BOUNDARY CAUTION (R > 60%) Step 5000/10000: H_drift=3.161573e-03, R=2.267, R/L=0.09 Step 6000/10000: H_drift=3.325463e-03, R=4.230, R/L=0.17 Step 7000/10000: H_drift=3.220102e-03, R=7.151, R/L=0.28 Step 8000/10000: H_drift=3.209462e-03, R=11.684, R/L=0.46 Step 9000/10000: H_drift=3.220940e-03, R=15.398, R/L=0.60 ⚠️ BOUNDARY CAUTION (R > 60%) Final H_drift = 0.003172, R_final = 18.751 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G128_dx04_S10000 complete [RUN] GRID_G256_dx02_S10000 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 256x256, dx=0.2, L=51.2, R_max=25.6 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 302.363978, Steps: 10000 ====================================================================== Step 1000/10000: H_drift=3.113524e-03, R=7.543, R/L=0.29 Step 2000/10000: H_drift=3.321561e-03, R=146.045, R/L=5.70 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/10000: H_drift=3.298810e-03, R=12.487, R/L=0.49 Step 4000/10000: H_drift=3.227647e-03, R=17.264, R/L=0.67 ⚠️ BOUNDARY CAUTION (R > 60%) Step 5000/10000: H_drift=3.170011e-03, R=2.284, R/L=0.09 Step 6000/10000: H_drift=3.183239e-03, R=4.212, R/L=0.16 Step 7000/10000: H_drift=3.221101e-03, R=7.141, R/L=0.28 Step 8000/10000: H_drift=3.230643e-03, R=11.687, R/L=0.46 Step 9000/10000: H_drift=3.209379e-03, R=15.405, R/L=0.60 ⚠️ BOUNDARY CAUTION (R > 60%) Final H_drift = 0.003218, R_final = 18.698 ====================================================================== Simulation complete ====================================================================== ✅ GRID_G256_dx02_S10000 complete ====================================================================== MANIFOLD SWEEP PHASE (45 NODES) ====================================================================== [RUN 1/45] T2_K0.00_B0.2_G0.1 κ=0.0, β=0.2, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.2, γ=0.1 Initial H = 137.117847, Steps: 5000 ====================================================================== Step 500/5000: H_drift=8.693586e-03, R=6.014, R/L=0.47 Step 1000/5000: H_drift=7.189141e-03, R=6.599, R/L=0.52 Step 1500/5000: H_drift=7.618167e-03, R=8.093, R/L=0.63 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=7.785326e-03, R=9.640, R/L=0.75 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2500/5000: H_drift=7.341888e-03, R=11.996, R/L=0.94 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=7.804424e-03, R=75.507, R/L=5.90 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=7.398889e-03, R=7.358, R/L=0.57 Step 4000/5000: H_drift=7.609733e-03, R=9.469, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4500/5000: H_drift=7.579589e-03, R=10.685, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.007685, R_final = 11.488 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B0.2_G0.1 complete [RUN 2/45] T2_K0.00_B0.2_G0.2 κ=0.0, β=0.2, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.2, γ=0.2 Initial H = 176.604882, Steps: 5000 ====================================================================== Step 500/5000: H_drift=6.750968e-03, R=6.178, R/L=0.48 Step 1000/5000: H_drift=5.583742e-03, R=7.406, R/L=0.58 Step 1500/5000: H_drift=5.915367e-03, R=9.356, R/L=0.73 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=6.041437e-03, R=11.341, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=5.692142e-03, R=19.253, R/L=1.50 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=6.047386e-03, R=5.151, R/L=0.40 Step 3500/5000: H_drift=5.732170e-03, R=8.845, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4000/5000: H_drift=5.895612e-03, R=10.720, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=5.868063e-03, R=12.377, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.005948, R_final = 14.541 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B0.2_G0.2 complete [RUN 3/45] T2_K0.00_B0.2_G0.4 κ=0.0, β=0.2, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.2, γ=0.4 Initial H = 255.578953, Steps: 5000 ====================================================================== Step 500/5000: H_drift=4.666285e-03, R=6.527, R/L=0.51 Step 1000/5000: H_drift=3.859487e-03, R=8.757, R/L=0.68 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=4.085197e-03, R=10.904, R/L=0.85 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=4.166554e-03, R=16.420, R/L=1.28 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=3.921279e-03, R=2.954, R/L=0.23 Step 3000/5000: H_drift=4.165370e-03, R=8.687, R/L=0.68 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3500/5000: H_drift=3.936965e-03, R=11.393, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=4.035729e-03, R=61.389, R/L=4.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=4.018518e-03, R=1.739, R/L=0.14 Final H_drift = 0.004083, R_final = 2.394 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B0.2_G0.4 complete [RUN 4/45] T2_K0.00_B0.5_G0.1 κ=0.0, β=0.5, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.1 Initial H = 250.577697, Steps: 5000 ====================================================================== Step 500/5000: H_drift=4.758627e-03, R=6.020, R/L=0.47 Step 1000/5000: H_drift=3.936405e-03, R=6.674, R/L=0.52 Step 1500/5000: H_drift=4.170431e-03, R=8.837, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=4.259136e-03, R=16.782, R/L=1.31 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=4.013068e-03, R=9.412, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=4.264719e-03, R=10.601, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=4.045569e-03, R=11.456, R/L=0.90 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=4.167429e-03, R=12.354, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=4.158012e-03, R=14.392, R/L=1.12 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.004219, R_final = 0.000 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B0.5_G0.1 complete [RUN 5/45] T2_K0.00_B0.5_G0.2 κ=0.0, β=0.5, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.2 Initial H = 290.064733, Steps: 5000 ====================================================================== Step 500/5000: H_drift=4.111494e-03, R=6.179, R/L=0.48 Step 1000/5000: H_drift=3.401371e-03, R=7.693, R/L=0.60 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=3.602152e-03, R=10.191, R/L=0.80 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=3.676460e-03, R=-0.702, R/L=-0.05 Step 2500/5000: H_drift=3.462290e-03, R=9.412, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=3.680619e-03, R=10.781, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=3.492845e-03, R=12.247, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=3.597672e-03, R=14.983, R/L=1.17 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=3.588153e-03, R=0.000, R/L=0.00 Final H_drift = 0.003639, R_final = 2.049 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B0.5_G0.2 complete [RUN 6/45] T2_K0.00_B0.5_G0.4 κ=0.0, β=0.5, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=0.5, γ=0.4 Initial H = 369.038803, Steps: 5000 ====================================================================== Step 500/5000: H_drift=3.232517e-03, R=6.537, R/L=0.51 Step 1000/5000: H_drift=2.673828e-03, R=9.002, R/L=0.70 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=2.828919e-03, R=12.665, R/L=0.99 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=2.884271e-03, R=6.861, R/L=0.54 Step 2500/5000: H_drift=2.716196e-03, R=10.007, R/L=0.78 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=2.886849e-03, R=12.312, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.734732e-03, R=16.669, R/L=1.30 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=2.816194e-03, R=1.802, R/L=0.14 Step 4500/5000: H_drift=2.808123e-03, R=2.815, R/L=0.22 Final H_drift = 0.002845, R_final = 3.625 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B0.5_G0.4 complete [RUN 7/45] T2_K0.00_B1.0_G0.1 κ=0.0, β=1.0, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=1.0, γ=0.1 Initial H = 439.677449, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.713246e-03, R=6.030, R/L=0.47 Step 1000/5000: H_drift=2.244852e-03, R=6.908, R/L=0.54 Step 1500/5000: H_drift=2.376335e-03, R=26.378, R/L=2.06 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=2.424189e-03, R=10.622, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=2.283615e-03, R=11.366, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=2.430812e-03, R=12.119, R/L=0.95 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.310680e-03, R=13.623, R/L=1.06 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=2.381127e-03, R=0.000, R/L=0.00 Step 4500/5000: H_drift=2.371550e-03, R=3.159, R/L=0.25 Final H_drift = 0.002402, R_final = 4.260 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B1.0_G0.1 complete [RUN 8/45] T2_K0.00_B1.0_G0.2 κ=0.0, β=1.0, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=1.0, γ=0.2 Initial H = 479.164484, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.490009e-03, R=6.202, R/L=0.48 Step 1000/5000: H_drift=2.060056e-03, R=8.200, R/L=0.64 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=2.179704e-03, R=3.462, R/L=0.27 Step 2000/5000: H_drift=2.222816e-03, R=10.359, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=2.094588e-03, R=11.428, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=2.230820e-03, R=12.868, R/L=1.01 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.119909e-03, R=17.168, R/L=1.34 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=2.182473e-03, R=2.266, R/L=0.18 Step 4500/5000: H_drift=2.172033e-03, R=3.673, R/L=0.29 Final H_drift = 0.002200, R_final = 4.653 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B1.0_G0.2 complete [RUN 9/45] T2_K0.00_B1.0_G0.4 κ=0.0, β=1.0, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.0, β=1.0, γ=0.4 Initial H = 558.138555, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.138194e-03, R=6.556, R/L=0.51 Step 1000/5000: H_drift=1.768431e-03, R=9.546, R/L=0.75 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=1.869384e-03, R=7.086, R/L=0.55 Step 2000/5000: H_drift=1.905981e-03, R=10.348, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.797163e-03, R=12.307, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.912838e-03, R=16.806, R/L=1.31 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.815408e-03, R=2.087, R/L=0.16 Step 4000/5000: H_drift=1.867317e-03, R=3.359, R/L=0.26 Step 4500/5000: H_drift=1.857979e-03, R=4.841, R/L=0.38 Final H_drift = 0.001879, R_final = 6.571 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.00_B1.0_G0.4 complete [RUN 10/45] T2_K0.05_B0.2_G0.1 κ=0.05, β=0.2, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=0.2, γ=0.1 Initial H = 143.585378, Steps: 5000 ====================================================================== Step 500/5000: H_drift=5.838538e-03, R=6.025, R/L=0.47 Step 1000/5000: H_drift=5.068112e-05, R=6.631, R/L=0.52 Step 1500/5000: H_drift=3.189795e-03, R=8.173, R/L=0.64 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=4.633787e-03, R=9.730, R/L=0.76 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2500/5000: H_drift=5.423148e-03, R=12.135, R/L=0.95 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=5.408038e-03, R=261.967, R/L=20.47 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=7.296211e-03, R=7.394, R/L=0.58 Step 4000/5000: H_drift=9.947969e-03, R=9.450, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4500/5000: H_drift=1.377779e-02, R=10.680, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.016911, R_final = 11.508 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B0.2_G0.1 complete [RUN 11/45] T2_K0.05_B0.2_G0.2 κ=0.05, β=0.2, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=0.2, γ=0.2 Initial H = 183.072413, Steps: 5000 ====================================================================== Step 500/5000: H_drift=4.620244e-03, R=6.188, R/L=0.48 Step 1000/5000: H_drift=3.513277e-04, R=7.480, R/L=0.58 Step 1500/5000: H_drift=1.390269e-03, R=9.398, R/L=0.73 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=1.793127e-03, R=11.413, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=2.039824e-03, R=19.801, R/L=1.55 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=2.147326e-03, R=5.260, R/L=0.41 Step 3500/5000: H_drift=4.089192e-03, R=8.876, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4000/5000: H_drift=6.590788e-03, R=10.753, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=9.595525e-03, R=12.452, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.011650, R_final = 14.668 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B0.2_G0.2 complete [RUN 12/45] T2_K0.05_B0.2_G0.4 κ=0.05, β=0.2, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=0.2, γ=0.4 Initial H = 262.046484, Steps: 5000 ====================================================================== Step 500/5000: H_drift=3.282569e-03, R=6.532, R/L=0.51 Step 1000/5000: H_drift=7.242694e-04, R=8.782, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=2.971141e-04, R=10.942, R/L=0.85 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=7.838404e-04, R=16.584, R/L=1.30 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.050916e-03, R=3.038, R/L=0.24 Step 3000/5000: H_drift=1.105798e-03, R=8.713, R/L=0.68 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3500/5000: H_drift=1.532154e-05, R=11.442, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=8.686179e-04, R=18.684, R/L=1.46 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=1.534953e-03, R=1.754, R/L=0.14 Final H_drift = 0.001731, R_final = 2.398 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B0.2_G0.4 complete [RUN 13/45] T2_K0.05_B0.5_G0.1 κ=0.05, β=0.5, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=0.5, γ=0.1 Initial H = 257.045229, Steps: 5000 ====================================================================== Step 500/5000: H_drift=3.288520e-03, R=6.031, R/L=0.47 Step 1000/5000: H_drift=2.515171e-04, R=6.720, R/L=0.52 Step 1500/5000: H_drift=8.966133e-04, R=8.939, R/L=0.70 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=9.738042e-04, R=17.391, R/L=1.36 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=9.399374e-04, R=9.383, R/L=0.73 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=9.247541e-04, R=10.572, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.306946e-03, R=11.439, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=4.027482e-03, R=12.368, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=5.891116e-03, R=14.485, R/L=1.13 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.006935, R_final = 0.000 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B0.5_G0.1 complete [RUN 14/45] T2_K0.05_B0.5_G0.2 κ=0.05, β=0.5, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=0.5, γ=0.2 Initial H = 296.532264, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.875458e-03, R=6.197, R/L=0.48 Step 1000/5000: H_drift=4.446990e-04, R=7.742, R/L=0.60 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=1.645409e-04, R=10.278, R/L=0.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=1.315824e-04, R=-0.217, R/L=-0.02 Step 2500/5000: H_drift=3.468034e-04, R=9.404, R/L=0.73 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=3.925761e-04, R=10.777, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=7.088245e-04, R=12.258, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=1.797733e-03, R=15.083, R/L=1.18 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=2.702999e-03, R=0.143, R/L=0.01 Final H_drift = 0.002997, R_final = 2.049 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B0.5_G0.2 complete [RUN 15/45] T2_K0.05_B0.5_G0.4 κ=0.05, β=0.5, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=0.5, γ=0.4 Initial H = 375.506334, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.308215e-03, R=6.547, R/L=0.51 Step 1000/5000: H_drift=6.669195e-04, R=9.037, R/L=0.71 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=6.745591e-04, R=12.722, R/L=0.99 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=1.353413e-03, R=6.885, R/L=0.54 Step 2500/5000: H_drift=1.768355e-03, R=10.003, R/L=0.78 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=1.975951e-03, R=12.328, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.533166e-03, R=16.853, R/L=1.32 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=1.447988e-03, R=1.853, R/L=0.14 Step 4500/5000: H_drift=1.596937e-03, R=2.867, R/L=0.22 Final H_drift = 0.001888, R_final = 3.665 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B0.5_G0.4 complete [RUN 16/45] T2_K0.05_B1.0_G0.1 κ=0.05, β=1.0, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=1.0, γ=0.1 Initial H = 446.144980, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.920245e-03, R=6.037, R/L=0.47 Step 1000/5000: H_drift=3.990864e-04, R=6.969, R/L=0.54 Step 1500/5000: H_drift=2.456781e-04, R=29.136, R/L=2.28 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=7.491348e-04, R=10.582, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.088620e-03, R=11.329, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.211469e-03, R=12.108, R/L=0.95 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=6.534512e-04, R=13.668, R/L=1.07 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=3.366259e-04, R=0.000, R/L=0.00 Step 4500/5000: H_drift=3.357640e-04, R=3.202, R/L=0.25 Final H_drift = 0.000614, R_final = 4.260 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B1.0_G0.1 complete [RUN 17/45] T2_K0.05_B1.0_G0.2 κ=0.05, β=1.0, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=1.0, γ=0.2 Initial H = 485.632015, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.778807e-03, R=6.212, R/L=0.49 Step 1000/5000: H_drift=4.916184e-04, R=8.257, R/L=0.65 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=5.381667e-04, R=3.628, R/L=0.28 Step 2000/5000: H_drift=1.153893e-03, R=10.334, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.536522e-03, R=11.409, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.710831e-03, R=12.885, R/L=1.01 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.372792e-03, R=17.495, R/L=1.37 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=1.366430e-03, R=2.317, R/L=0.18 Step 4500/5000: H_drift=1.619723e-03, R=3.678, R/L=0.29 Final H_drift = 0.001970, R_final = 4.663 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B1.0_G0.2 complete [RUN 18/45] T2_K0.05_B1.0_G0.4 κ=0.05, β=1.0, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.05, β=1.0, γ=0.4 Initial H = 564.606086, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.554189e-03, R=6.569, R/L=0.51 Step 1000/5000: H_drift=6.143684e-04, R=9.580, R/L=0.75 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=9.203815e-04, R=7.097, R/L=0.55 Step 2000/5000: H_drift=1.654725e-03, R=10.336, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=2.065822e-03, R=12.324, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=2.314154e-03, R=16.917, R/L=1.32 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.269919e-03, R=2.136, R/L=0.17 Step 4000/5000: H_drift=2.581915e-03, R=3.385, R/L=0.26 Step 4500/5000: H_drift=3.016666e-03, R=4.873, R/L=0.38 Final H_drift = 0.003385, R_final = 6.612 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.05_B1.0_G0.4 complete [RUN 19/45] T2_K0.10_B0.2_G0.1 κ=0.1, β=0.2, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=0.2, γ=0.1 Initial H = 150.052909, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.430950e-03, R=6.035, R/L=0.47 Step 1000/5000: H_drift=8.809401e-03, R=6.671, R/L=0.52 Step 1500/5000: H_drift=1.575511e-02, R=8.248, R/L=0.64 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=1.848298e-02, R=9.816, R/L=0.77 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2500/5000: H_drift=1.960238e-02, R=12.266, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=2.092595e-02, R=-194.572, R/L=-15.20 Step 3500/5000: H_drift=2.619302e-02, R=7.421, R/L=0.58 Step 4000/5000: H_drift=3.456733e-02, R=9.428, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4500/5000: H_drift=4.616283e-02, R=10.673, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.056600, R_final = 11.516 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B0.2_G0.1 complete [RUN 20/45] T2_K0.10_B0.2_G0.2 κ=0.1, β=0.2, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=0.2, γ=0.2 Initial H = 189.539944, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.017899e-03, R=6.204, R/L=0.48 Step 1000/5000: H_drift=6.090057e-03, R=7.536, R/L=0.59 Step 1500/5000: H_drift=1.009215e-02, R=9.439, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=1.100252e-02, R=11.474, R/L=0.90 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.164045e-02, R=20.400, R/L=1.59 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.357310e-02, R=5.355, R/L=0.42 Step 3500/5000: H_drift=1.935075e-02, R=8.880, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4000/5000: H_drift=2.740882e-02, R=10.785, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=3.664121e-02, R=12.514, R/L=0.98 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.043310, R_final = 14.795 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B0.2_G0.2 complete [RUN 21/45] T2_K0.10_B0.2_G0.4 κ=0.1, β=0.2, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=0.2, γ=0.4 Initial H = 268.514015, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.550385e-03, R=6.544, R/L=0.51 Step 1000/5000: H_drift=3.208006e-03, R=8.810, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=4.381365e-03, R=10.980, R/L=0.86 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=3.656677e-03, R=16.737, R/L=1.31 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=3.571593e-03, R=3.174, R/L=0.25 Step 3000/5000: H_drift=4.895783e-03, R=8.728, R/L=0.68 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3500/5000: H_drift=8.142948e-03, R=11.472, R/L=0.90 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=1.100041e-02, R=17.272, R/L=1.35 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=1.285288e-02, R=1.742, R/L=0.14 Final H_drift = 0.013566, R_final = 2.367 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B0.2_G0.4 complete [RUN 22/45] T2_K0.10_B0.5_G0.1 κ=0.1, β=0.5, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=0.5, γ=0.1 Initial H = 263.512760, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.446804e-03, R=6.039, R/L=0.47 Step 1000/5000: H_drift=4.368065e-03, R=6.767, R/L=0.53 Step 1500/5000: H_drift=7.028687e-03, R=9.035, R/L=0.71 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=7.269951e-03, R=18.017, R/L=1.41 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=7.424486e-03, R=9.346, R/L=0.73 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=8.814435e-03, R=10.530, R/L=0.82 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.311385e-02, R=11.388, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=1.878151e-02, R=12.380, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=2.457044e-02, R=14.567, R/L=1.14 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.027987, R_final = 0.000 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B0.5_G0.1 complete [RUN 23/45] T2_K0.10_B0.5_G0.2 κ=0.1, β=0.5, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=0.5, γ=0.2 Initial H = 302.999795, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.315533e-03, R=6.205, R/L=0.48 Step 1000/5000: H_drift=3.281625e-03, R=7.788, R/L=0.61 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=4.789234e-03, R=10.352, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=4.361849e-03, R=0.281, R/L=0.02 Step 2500/5000: H_drift=4.288955e-03, R=9.389, R/L=0.73 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=5.531096e-03, R=10.770, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=8.937874e-03, R=12.269, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=1.249755e-02, R=15.173, R/L=1.19 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=1.510272e-02, R=0.470, R/L=0.04 Final H_drift = 0.015995, R_final = 1.983 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B0.5_G0.2 complete [RUN 24/45] T2_K0.10_B0.5_G0.4 κ=0.1, β=0.5, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=0.5, γ=0.4 Initial H = 381.973866, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.130424e-03, R=6.554, R/L=0.51 Step 1000/5000: H_drift=1.881506e-03, R=9.062, R/L=0.71 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=2.057199e-03, R=12.772, R/L=1.00 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=9.159955e-04, R=6.900, R/L=0.54 Step 2500/5000: H_drift=4.630251e-04, R=10.009, R/L=0.78 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=9.426361e-04, R=12.349, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.220132e-03, R=17.070, R/L=1.33 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=2.683994e-03, R=1.859, R/L=0.15 Step 4500/5000: H_drift=2.447754e-03, R=2.888, R/L=0.23 Final H_drift = 0.002139, R_final = 3.705 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B0.5_G0.4 complete [RUN 25/45] T2_K0.10_B1.0_G0.1 κ=0.1, β=1.0, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=1.0, γ=0.1 Initial H = 452.612511, Steps: 5000 ====================================================================== Step 500/5000: H_drift=9.019463e-04, R=6.040, R/L=0.47 Step 1000/5000: H_drift=1.951340e-03, R=7.028, R/L=0.55 Step 1500/5000: H_drift=2.416142e-03, R=34.123, R/L=2.67 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=1.557847e-03, R=10.542, R/L=0.82 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.205975e-03, R=11.292, R/L=0.88 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.830547e-03, R=12.095, R/L=0.94 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=3.498598e-03, R=13.704, R/L=1.07 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=4.554487e-03, R=0.000, R/L=0.00 Step 4500/5000: H_drift=4.557709e-03, R=3.240, R/L=0.25 Final H_drift = 0.004061, R_final = 4.272 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B1.0_G0.1 complete [RUN 26/45] T2_K0.10_B1.0_G0.2 κ=0.1, β=1.0, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=1.0, γ=0.2 Initial H = 492.099546, Steps: 5000 ====================================================================== Step 500/5000: H_drift=8.637050e-04, R=6.223, R/L=0.49 Step 1000/5000: H_drift=1.509142e-03, R=8.306, R/L=0.65 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=1.547765e-03, R=3.791, R/L=0.30 Step 2000/5000: H_drift=4.841176e-04, R=10.312, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=3.096106e-05, R=11.397, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=3.800240e-04, R=12.894, R/L=1.01 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.327942e-03, R=0.000, R/L=0.00 Step 4000/5000: H_drift=1.488657e-03, R=2.355, R/L=0.18 Step 4500/5000: H_drift=9.744061e-04, R=3.698, R/L=0.29 Final H_drift = 0.000530, R_final = 4.686 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B1.0_G0.2 complete [RUN 27/45] T2_K0.10_B1.0_G0.4 κ=0.1, β=1.0, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.1, β=1.0, γ=0.4 Initial H = 571.073617, Steps: 5000 ====================================================================== Step 500/5000: H_drift=8.005625e-04, R=6.589, R/L=0.51 Step 1000/5000: H_drift=8.633742e-04, R=9.613, R/L=0.75 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=3.483264e-04, R=7.103, R/L=0.55 Step 2000/5000: H_drift=9.403297e-04, R=10.321, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.529042e-03, R=12.340, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.629653e-03, R=17.064, R/L=1.33 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.564633e-03, R=2.181, R/L=0.17 Step 4000/5000: H_drift=2.133648e-03, R=3.410, R/L=0.27 Step 4500/5000: H_drift=2.812036e-03, R=4.904, R/L=0.38 Final H_drift = 0.003217, R_final = 6.664 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.10_B1.0_G0.4 complete [RUN 28/45] T2_K0.20_B0.2_G0.1 κ=0.2, β=0.2, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=0.2, γ=0.1 Initial H = 162.987971, Steps: 5000 ====================================================================== Step 500/5000: H_drift=5.703306e-03, R=6.045, R/L=0.47 Step 1000/5000: H_drift=2.938644e-02, R=6.738, R/L=0.53 Step 1500/5000: H_drift=4.331935e-02, R=8.388, R/L=0.66 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=4.710627e-02, R=9.973, R/L=0.78 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2500/5000: H_drift=4.908052e-02, R=12.503, R/L=0.98 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=5.612323e-02, R=-35.752, R/L=-2.79 Step 3500/5000: H_drift=7.489208e-02, R=7.436, R/L=0.58 Step 4000/5000: H_drift=1.080366e-01, R=9.367, R/L=0.73 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4500/5000: H_drift=1.610239e-01, R=10.621, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.220811, R_final = 11.495 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B0.2_G0.1 complete [RUN 29/45] T2_K0.20_B0.2_G0.2 κ=0.2, β=0.2, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=0.2, γ=0.2 Initial H = 202.475006, Steps: 5000 ====================================================================== Step 500/5000: H_drift=4.362412e-03, R=6.221, R/L=0.49 Step 1000/5000: H_drift=2.159452e-02, R=7.635, R/L=0.60 Step 1500/5000: H_drift=2.994271e-02, R=9.516, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=3.171448e-02, R=11.587, R/L=0.91 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=3.500788e-02, R=21.761, R/L=1.70 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=4.554249e-02, R=5.503, R/L=0.43 Step 3500/5000: H_drift=6.847488e-02, R=8.888, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4000/5000: H_drift=1.038811e-01, R=10.814, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=1.491762e-01, R=12.604, R/L=0.98 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.187166, R_final = 14.960 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B0.2_G0.2 complete [RUN 30/45] T2_K0.20_B0.2_G0.4 κ=0.2, β=0.2, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=0.2, γ=0.4 Initial H = 281.449077, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.824269e-03, R=6.556, R/L=0.51 Step 1000/5000: H_drift=1.294676e-02, R=8.860, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=1.568704e-02, R=11.041, R/L=0.86 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=1.521097e-02, R=17.052, R/L=1.33 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.809257e-02, R=3.365, R/L=0.26 Step 3000/5000: H_drift=2.690623e-02, R=8.764, R/L=0.68 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3500/5000: H_drift=4.064613e-02, R=11.543, R/L=0.90 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=5.316925e-02, R=16.870, R/L=1.32 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=6.086515e-02, R=1.553, R/L=0.12 Final H_drift = 0.063969, R_final = 2.082 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B0.2_G0.4 complete [RUN 31/45] T2_K0.20_B0.5_G0.1 κ=0.2, β=0.5, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=0.5, γ=0.1 Initial H = 276.447822, Steps: 5000 ====================================================================== Step 500/5000: H_drift=3.204211e-03, R=6.045, R/L=0.47 Step 1000/5000: H_drift=1.575451e-02, R=6.847, R/L=0.53 Step 1500/5000: H_drift=2.135612e-02, R=9.213, R/L=0.72 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=2.198434e-02, R=19.586, R/L=1.53 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=2.437583e-02, R=9.268, R/L=0.72 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=3.301512e-02, R=10.451, R/L=0.82 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=5.119821e-02, R=11.356, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=7.763328e-02, R=12.362, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=1.079869e-01, R=14.601, R/L=1.14 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 0.129137, R_final = 0.000 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B0.5_G0.1 complete [RUN 32/45] T2_K0.20_B0.5_G0.2 κ=0.2, β=0.5, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=0.5, γ=0.2 Initial H = 315.934857, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.660265e-03, R=6.225, R/L=0.49 Step 1000/5000: H_drift=1.255450e-02, R=7.879, R/L=0.62 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=1.586958e-02, R=10.476, R/L=0.82 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=1.564276e-02, R=1.014, R/L=0.08 Step 2500/5000: H_drift=1.806175e-02, R=9.345, R/L=0.73 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=2.631495e-02, R=10.741, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=4.114217e-02, R=12.299, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=5.765702e-02, R=15.310, R/L=1.20 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=6.990366e-02, R=0.513, R/L=0.04 Final H_drift = 0.074095, R_final = 1.566 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B0.5_G0.2 complete [RUN 33/45] T2_K0.20_B0.5_G0.4 κ=0.2, β=0.5, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=0.5, γ=0.4 Initial H = 394.908928, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.908969e-03, R=6.571, R/L=0.51 Step 1000/5000: H_drift=8.314682e-03, R=9.122, R/L=0.71 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=8.925051e-03, R=12.858, R/L=1.00 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=7.586200e-03, R=6.928, R/L=0.54 Step 2500/5000: H_drift=9.084584e-03, R=10.018, R/L=0.78 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=1.382345e-02, R=12.409, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.932104e-02, R=17.593, R/L=1.37 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=2.162059e-02, R=1.764, R/L=0.14 Step 4500/5000: H_drift=2.132251e-02, R=2.889, R/L=0.23 Final H_drift = 0.021882, R_final = 3.775 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B0.5_G0.4 complete [RUN 34/45] T2_K0.20_B1.0_G0.1 κ=0.2, β=1.0, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=1.0, γ=0.1 Initial H = 465.547573, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.748983e-03, R=6.055, R/L=0.47 Step 1000/5000: H_drift=7.900615e-03, R=7.156, R/L=0.56 Step 1500/5000: H_drift=8.997700e-03, R=51.070, R/L=3.99 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=7.959507e-03, R=10.469, R/L=0.82 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=9.360847e-03, R=11.240, R/L=0.88 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.442746e-02, R=12.068, R/L=0.94 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.177254e-02, R=13.763, R/L=1.08 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=2.661996e-02, R=0.000, R/L=0.00 Step 4500/5000: H_drift=2.660677e-02, R=3.248, R/L=0.25 Final H_drift = 0.025265, R_final = 4.357 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B1.0_G0.1 complete [RUN 35/45] T2_K0.20_B1.0_G0.2 κ=0.2, β=1.0, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=1.0, γ=0.2 Initial H = 505.034609, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.525560e-03, R=6.238, R/L=0.49 Step 1000/5000: H_drift=6.596580e-03, R=8.403, R/L=0.66 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=6.831590e-03, R=4.044, R/L=0.32 Step 2000/5000: H_drift=5.464896e-03, R=10.262, R/L=0.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=6.499719e-03, R=11.370, R/L=0.89 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.014994e-02, R=12.920, R/L=1.01 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.418815e-02, R=0.000, R/L=0.00 Step 4000/5000: H_drift=1.509041e-02, R=2.277, R/L=0.18 Step 4500/5000: H_drift=1.375298e-02, R=3.770, R/L=0.29 Final H_drift = 0.014011, R_final = 4.788 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B1.0_G0.2 complete [RUN 36/45] T2_K0.20_B1.0_G0.4 κ=0.2, β=1.0, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.2, β=1.0, γ=0.4 Initial H = 584.008679, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.175916e-03, R=6.604, R/L=0.52 Step 1000/5000: H_drift=4.653211e-03, R=9.678, R/L=0.76 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=3.755606e-03, R=7.124, R/L=0.56 Step 2000/5000: H_drift=1.949927e-03, R=10.296, R/L=0.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=2.178323e-03, R=12.365, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=3.381316e-03, R=17.414, R/L=1.36 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=3.722694e-03, R=2.254, R/L=0.18 Step 4000/5000: H_drift=2.521544e-03, R=3.441, R/L=0.27 Step 4500/5000: H_drift=2.087347e-03, R=4.968, R/L=0.39 Final H_drift = 0.002711, R_final = 6.754 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.20_B1.0_G0.4 complete [RUN 37/45] T2_K0.40_B0.2_G0.1 κ=0.4, β=0.2, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=0.2, γ=0.1 Initial H = 188.858096, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.554422e-02, R=6.080, R/L=0.48 Step 1000/5000: H_drift=7.578117e-02, R=6.876, R/L=0.54 Step 1500/5000: H_drift=9.675931e-02, R=8.626, R/L=0.67 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=9.660154e-02, R=10.249, R/L=0.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=1.025330e-01, R=12.974, R/L=1.01 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.353844e-01, R=-9.485, R/L=-0.74 Step 3500/5000: H_drift=2.267575e-01, R=7.412, R/L=0.58 Step 4000/5000: H_drift=4.570218e-01, R=9.208, R/L=0.72 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4500/5000: H_drift=1.050085e+00, R=10.435, R/L=0.82 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 2.164488, R_final = 11.275 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B0.2_G0.1 complete [RUN 38/45] T2_K0.40_B0.2_G0.2 κ=0.4, β=0.2, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=0.2, γ=0.2 Initial H = 228.345131, Steps: 5000 ====================================================================== Step 500/5000: H_drift=2.054345e-02, R=6.251, R/L=0.49 Step 1000/5000: H_drift=5.807421e-02, R=7.802, R/L=0.61 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=7.173409e-02, R=9.651, R/L=0.75 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=7.487556e-02, R=11.789, R/L=0.92 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=9.277119e-02, R=26.927, R/L=2.10 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=1.495088e-01, R=5.708, R/L=0.45 Step 3500/5000: H_drift=2.890534e-01, R=8.869, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 4000/5000: H_drift=5.909196e-01, R=10.814, R/L=0.84 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=1.183540e+00, R=12.560, R/L=0.98 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 2.037220, R_final = 14.535 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B0.2_G0.2 complete [RUN 39/45] T2_K0.40_B0.2_G0.4 κ=0.4, β=0.2, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=0.2, γ=0.4 Initial H = 307.319202, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.443705e-02, R=6.600, R/L=0.52 Step 1000/5000: H_drift=3.714112e-02, R=8.953, R/L=0.70 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=4.254045e-02, R=11.149, R/L=0.87 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=4.672382e-02, R=17.987, R/L=1.41 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=7.047614e-02, R=3.692, R/L=0.29 Step 3000/5000: H_drift=1.299475e-01, R=8.818, R/L=0.69 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3500/5000: H_drift=2.294396e-01, R=11.623, R/L=0.91 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=3.459403e-01, R=16.821, R/L=1.31 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=4.529755e-01, R=0.741, R/L=0.06 Final H_drift = 0.532002, R_final = 0.748 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B0.2_G0.4 complete [RUN 40/45] T2_K0.40_B0.5_G0.1 κ=0.4, β=0.5, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=0.5, γ=0.1 Initial H = 302.317947, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.553226e-02, R=6.084, R/L=0.48 Step 1000/5000: H_drift=4.362972e-02, R=7.010, R/L=0.55 Step 1500/5000: H_drift=5.295182e-02, R=9.489, R/L=0.74 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2000/5000: H_drift=5.524084e-02, R=25.215, R/L=1.97 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=7.192139e-02, R=9.122, R/L=0.71 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=1.235157e-01, R=10.283, R/L=0.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.467877e-01, R=11.188, R/L=0.87 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=5.051028e-01, R=12.196, R/L=0.95 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=9.985787e-01, R=14.037, R/L=1.10 ⚠️ BOUNDARY WARNING (R > 80%) Final H_drift = 1.707893, R_final = 15.755 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B0.5_G0.1 complete [RUN 41/45] T2_K0.40_B0.5_G0.2 κ=0.4, β=0.5, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=0.5, γ=0.2 Initial H = 341.804982, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.335667e-02, R=6.260, R/L=0.49 Step 1000/5000: H_drift=3.571975e-02, R=8.042, R/L=0.63 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=4.170827e-02, R=10.658, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=4.481480e-02, R=2.060, R/L=0.16 Step 2500/5000: H_drift=6.467375e-02, R=9.248, R/L=0.72 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=1.195597e-01, R=10.671, R/L=0.83 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=2.302200e-01, R=12.291, R/L=0.96 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=3.999636e-01, R=15.308, R/L=1.20 ⚠️ BOUNDARY WARNING (R > 80%) Step 4500/5000: H_drift=6.082213e-01, R=0.000, R/L=0.00 Final H_drift = 0.805239, R_final = 0.405 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B0.5_G0.2 complete [RUN 42/45] T2_K0.40_B0.5_G0.4 κ=0.4, β=0.5, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=0.5, γ=0.4 Initial H = 420.779052, Steps: 5000 ====================================================================== Step 500/5000: H_drift=1.025863e-02, R=6.613, R/L=0.52 Step 1000/5000: H_drift=2.489940e-02, R=9.221, R/L=0.72 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=2.657670e-02, R=13.097, R/L=1.02 ⚠️ BOUNDARY WARNING (R > 80%) Step 2000/5000: H_drift=2.926647e-02, R=6.955, R/L=0.54 Step 2500/5000: H_drift=4.733524e-02, R=10.012, R/L=0.78 ⚠️ BOUNDARY CAUTION (R > 60%) Step 3000/5000: H_drift=8.457320e-02, R=12.487, R/L=0.98 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.275472e-01, R=0.000, R/L=0.00 Step 4000/5000: H_drift=1.514655e-01, R=1.017, R/L=0.08 Step 4500/5000: H_drift=1.512250e-01, R=1.268, R/L=0.10 Final H_drift = 0.145774, R_final = 3.439 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B0.5_G0.4 complete [RUN 43/45] T2_K0.40_B1.0_G0.1 κ=0.4, β=1.0, γ=0.1 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=1.0, γ=0.1 Initial H = 491.417698, Steps: 5000 ====================================================================== Step 500/5000: H_drift=9.128372e-03, R=6.091, R/L=0.48 Step 1000/5000: H_drift=2.328009e-02, R=7.502, R/L=0.59 Step 1500/5000: H_drift=2.546283e-02, R=-112.909, R/L=-8.82 Step 2000/5000: H_drift=2.751491e-02, R=10.325, R/L=0.81 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=4.342578e-02, R=11.127, R/L=0.87 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=8.221192e-02, R=12.000, R/L=0.94 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.436435e-01, R=13.682, R/L=1.07 ⚠️ BOUNDARY WARNING (R > 80%) Step 4000/5000: H_drift=2.018074e-01, R=0.000, R/L=0.00 Step 4500/5000: H_drift=2.229858e-01, R=0.746, R/L=0.06 Final H_drift = 0.214494, R_final = 6.569 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B1.0_G0.1 complete [RUN 44/45] T2_K0.40_B1.0_G0.2 κ=0.4, β=1.0, γ=0.2 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=1.0, γ=0.2 Initial H = 530.904733, Steps: 5000 ====================================================================== Step 500/5000: H_drift=8.212887e-03, R=6.274, R/L=0.49 Step 1000/5000: H_drift=1.989946e-02, R=8.584, R/L=0.67 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=2.068571e-02, R=4.447, R/L=0.35 Step 2000/5000: H_drift=2.248437e-02, R=10.170, R/L=0.79 ⚠️ BOUNDARY CAUTION (R > 60%) Step 2500/5000: H_drift=3.694728e-02, R=11.306, R/L=0.88 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=6.691448e-02, R=12.942, R/L=1.01 ⚠️ BOUNDARY WARNING (R > 80%) Step 3500/5000: H_drift=1.002958e-01, R=0.000, R/L=0.00 Step 4000/5000: H_drift=1.123036e-01, R=1.013, R/L=0.08 Step 4500/5000: H_drift=1.013953e-01, R=4.348, R/L=0.34 Final H_drift = 0.096071, R_final = 6.119 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B1.0_G0.2 complete [RUN 45/45] T2_K0.40_B1.0_G0.4 κ=0.4, β=1.0, γ=0.4 ====================================================================== FRCMFD v4.0 — Tier-2 Run ====================================================================== Grid: 128x128, dx=0.2, L=25.6, R_max=12.8 Params: κ=0.4, β=1.0, γ=0.4 Initial H = 609.878804, Steps: 5000 ====================================================================== Step 500/5000: H_drift=6.755928e-03, R=6.642, R/L=0.52 Step 1000/5000: H_drift=1.476898e-02, R=9.802, R/L=0.77 ⚠️ BOUNDARY CAUTION (R > 60%) Step 1500/5000: H_drift=1.358304e-02, R=7.137, R/L=0.56 Step 2000/5000: H_drift=1.438585e-02, R=10.262, R/L=0.80 ⚠️ BOUNDARY WARNING (R > 80%) Step 2500/5000: H_drift=2.385914e-02, R=12.396, R/L=0.97 ⚠️ BOUNDARY WARNING (R > 80%) Step 3000/5000: H_drift=3.621418e-02, R=0.000, R/L=0.00 Step 3500/5000: H_drift=3.946448e-02, R=1.756, R/L=0.14 Step 4000/5000: H_drift=3.438302e-02, R=3.583, R/L=0.28 Step 4500/5000: H_drift=3.938422e-02, R=5.250, R/L=0.41 Final H_drift = 0.051203, R_final = 7.047 ====================================================================== Simulation complete ====================================================================== ✅ T2_K0.40_B1.0_G0.4 complete ====================================================================== GRID INDEPENDENCE RATIO (GIR) ANALYSIS ====================================================================== κ=0.0, β=0.5, γ=0.2 dx1=0.2, dx2=0.4 GIR_H=0.3818, GIR_eta=0.0027 Grid independent: ❌ κ=0.0, β=0.5, γ=0.2 dx1=0.2, dx2=0.4 GIR_H=0.0190, GIR_eta=0.0007 Grid independent: ✅ κ=0.0, β=0.5, γ=0.2 dx1=0.2, dx2=0.4 GIR_H=0.0852, GIR_eta=0.1496 Grid independent: ❌ ✅ GIR analysis saved to /content/TIER2_SWEEP_20260620_115712/FRCMFD_GIR_ANALYSIS.csv ====================================================================== GENERATING HEATMAPS ====================================================================== ✅ Saved: /content/TIER2_SWEEP_20260620_115712/FRCMFD_TIER2_HEATMAPS.png [STEP 2] Master ZIP created: /content/TIER2_SWEEP_20260620_115712.zip [STEP 3] Backed up to Drive: /content/drive/MyDrive/FRCMFD_Data/TIER2_SWEEP_20260620_115712 [STEP 4] Download triggered [STEP 5] Verification: ✅ output_dir ✅ ledger ================================================== ✓ Colab workspace saved ✓ Google Drive backup saved ✓ Download package created ================================================== === FINAL STATUS REPORT === OUTPUT DIRECTORY: /content/TIER2_SWEEP_20260620_115712 GOOGLE DRIVE BACKUP: True MASTER ZIP: /content/TIER2_SWEEP_20260620_115712.zip FILE COUNT: 3 STATUS: SUCCESS =========================== ====================================================================== ✅ TIER-2 SUITE COMPLETE ====================================================================== Total runs: 57 Ledger saved to: /content/TIER2_SWEEP_20260620_115712/FRCMFD_TIER2_SUMMARY.csv Output directory: /content/TIER2_SWEEP_20260620_115712 ====================================================================== ====================================================================== ✅ TIER-2 EXECUTION COMPLETE ====================================================================== Results: 57 nodes processed Ledger: /content/TIER2_SWEEP_20260620_115712/FRCMFD_TIER2_SUMMARY.csv ====================================================================== #!/usr/bin/env python3 """ FRCMFD v4.0 — TIER-2 MANIFOLD SWEEP (FULLY INTEGRATED) ======================================================== Purpose: Map the (κ, β, γ) parameter space over 45 nodes with grid refinement, boundary monitoring, and CM tracking. Protocol: Layer 1-3 Data Preservation (Steps 1-6) Key Features: - Grid refinement (N=64,128,256) with fixed domain - 45-node manifold sweep (κ×β×γ) - Boundary monitoring (3-level flag) - Grid Independence Ratio (GIR) analysis - Center-of-mass tracking (v_CM) - Confinement diagnostics (R0, Rf, Rmax, eta, v_R) - Full data preservation (ZIP, Drive, Local) """ import os import sys import json import csv import glob import zipfile import shutil import numpy as np import matplotlib.pyplot as plt import pandas as pd from datetime import datetime import time # Try to import Colab-specific modules try: from google.colab import drive, files IN_COLAB = True except ImportError: IN_COLAB = False print("=" * 70) print("FRCMFD v4.0 — TIER-2 MANIFOLD SWEEP") print("=" * 70) print(f"Start time: {datetime.now().strftime('%Y-%m-%d %H:%M:%S')}") print("=" * 70) # ============================================================ # PART 1: THE SOLVER (Tier-2 with CM tracking) # ============================================================ class FRCMFDSolver: """Tier-2 solver with complete Hamiltonian and CM tracking.""" def __init__(self, params): self.params = params self.N = params.get('N', 64) self.dx = params.get('dx', 0.4) self.dt = params.get('dt', 0.001) self.L = self.N * self.dx self.R_max = self.L / 2.0 self.c_S = params.get('c_S', 1.0) self.c_Psi = params.get('c_Psi', 1.0) self.beta = params.get('beta', 0.5) self.gamma = params.get('gamma', 0.2) self.m2 = params.get('m2', 0.1) self.kappa = params.get('kappa', 0.0) self.s_max = params.get('s_max', 2.0) # Grid self.x = np.linspace(-self.N*self.dx/2, self.N*self.dx/2, self.N) self.y = np.linspace(-self.N*self.dx/2, self.N*self.dx/2, self.N) self.X, self.Y = np.meshgrid(self.x, self.y, indexing='ij') self.R2 = self.X**2 + self.Y**2 self.history = { 'step': [], 't': [], 'energy': [], 'norm_psi': [], 'norm_S': [], 'plateau_radius': [], 'saturation_penetration': [], 'cm_x': [], 'cm_y': [] } self.output_dir = None self.run_id = None def initialize_fields(self, amplitude=2.0, width=8.0, psi_amplitude=1.0, psi_width=5.0): S = amplitude * np.exp(-self.R2 / (2 * width**2)) Psi = psi_amplitude * np.exp(-self.R2 / (2 * psi_width**2)) * np.exp(1j * np.arctan2(self.Y, self.X)) dS_dt = np.zeros_like(S) dPsi_dt = np.zeros_like(Psi, dtype=np.complex128) return S, Psi, dS_dt, dPsi_dt def compute_laplacian(self, F): return (np.roll(F, -1, axis=0) + np.roll(F, 1, axis=0) + np.roll(F, -1, axis=1) + np.roll(F, 1, axis=1) - 4*F) / self.dx**2 def compute_gradients(self, F): grad_x = (np.roll(F, -1, axis=1) - np.roll(F, 1, axis=1)) / (2 * self.dx) grad_y = (np.roll(F, -1, axis=0) - np.roll(F, 1, axis=0)) / (2 * self.dx) return grad_x, grad_y def compute_hamiltonian(self, S, Psi, dS_dt, dPsi_dt): grad_Psi_x, grad_Psi_y = self.compute_gradients(Psi) grad_Psi_sq = np.abs(grad_Psi_x)**2 + np.abs(grad_Psi_y)**2 psi_sq = np.abs(Psi)**2 grad_S_x, grad_S_y = self.compute_gradients(S) grad_S_sq = grad_S_x**2 + grad_S_y**2 H_density = ( 0.5 * np.abs(dPsi_dt)**2 + 0.5 * self.c_Psi**2 * grad_Psi_sq + 0.5 * self.m2 * psi_sq + self.kappa * S * psi_sq + 0.5 * dS_dt**2 + 0.5 * self.c_S**2 * grad_S_sq + 0.5 * self.beta * S**2 + 0.25 * self.gamma * S**4 ) return np.sum(H_density) * self.dx**2 def compute_rhs(self, S, Psi, dS_dt, dPsi_dt): lap_S = self.compute_laplacian(S) lap_Psi = self.compute_laplacian(Psi) psi_sq = np.abs(Psi)**2 d2S = self.c_S**2 * lap_S - self.beta * S - self.gamma * S**3 - self.kappa * psi_sq d2Psi = self.c_Psi**2 * lap_Psi - self.m2 * Psi - self.kappa * S * Psi return d2S, d2Psi def apply_saturation(self, S, dS_dt): mask = S > self.s_max S[mask] = self.s_max dS_dt[mask] = 0.0 return S, dS_dt def compute_plateau_radius(self, S): center = S[self.N//2, self.N//2] threshold = 0.5 * center r = np.sqrt(self.R2) mask = S > threshold if np.any(mask): return np.sum(r[mask] * S[mask]) / np.sum(S[mask]) return 0.0 def compute_center_of_mass(self, Psi): psi_sq = np.abs(Psi)**2 total = np.sum(psi_sq) if total == 0: return 0.0, 0.0 cx = np.sum(self.X * psi_sq) / total cy = np.sum(self.Y * psi_sq) / total return cx, cy def rk4_step(self, S, Psi, dS_dt, dPsi_dt): dt = self.dt d2S1, d2Psi1 = self.compute_rhs(S, Psi, dS_dt, dPsi_dt) k1_S, k1_dS, k1_Psi, k1_dPsi = dS_dt, d2S1, dPsi_dt, d2Psi1 S2 = S + 0.5*dt*k1_S dS2 = dS_dt + 0.5*dt*k1_dS Psi2 = Psi + 0.5*dt*k1_Psi dPsi2 = dPsi_dt + 0.5*dt*k1_dPsi d2S2, d2Psi2 = self.compute_rhs(S2, Psi2, dS2, dPsi2) k2_S, k2_dS, k2_Psi, k2_dPsi = dS2, d2S2, dPsi2, d2Psi2 S3 = S + 0.5*dt*k2_S dS3 = dS_dt + 0.5*dt*k2_dS Psi3 = Psi + 0.5*dt*k2_Psi dPsi3 = dPsi_dt + 0.5*dt*k2_dPsi d2S3, d2Psi3 = self.compute_rhs(S3, Psi3, dS3, dPsi3) k3_S, k3_dS, k3_Psi, k3_dPsi = dS3, d2S3, dPsi3, d2Psi3 S4 = S + dt*k3_S dS4 = dS_dt + dt*k3_dS Psi4 = Psi + dt*k3_Psi dPsi4 = dPsi_dt + dt*k3_dPsi d2S4, d2Psi4 = self.compute_rhs(S4, Psi4, dS4, dPsi4) k4_S, k4_dS, k4_Psi, k4_dPsi = dS4, d2S4, dPsi4, d2Psi4 S_new = S + (dt/6)*(k1_S + 2*k2_S + 2*k3_S + k4_S) dS_new = dS_dt + (dt/6)*(k1_dS + 2*k2_dS + 2*k3_dS + k4_dS) Psi_new = Psi + (dt/6)*(k1_Psi + 2*k2_Psi + 2*k3_Psi + k4_Psi) dPsi_new = dPsi_dt + (dt/6)*(k1_dPsi + 2*k2_dPsi + 2*k3_dPsi + k4_dPsi) S_new, dS_new = self.apply_saturation(S_new, dS_new) return S_new, Psi_new, dS_new, dPsi_new def run_simulation(self, steps, save_interval=100, verbose=True): S, Psi, dS_dt, dPsi_dt = self.initialize_fields() H0 = self.compute_hamiltonian(S, Psi, dS_dt, dPsi_dt) if verbose: print("=" * 70) print("FRCMFD v4.0 — Tier-2 Run") print("=" * 70) print(f"Grid: {self.N}x{self.N}, dx={self.dx}, L={self.L:.1f}, R_max={self.R_max:.1f}") print(f"Params: κ={self.kappa}, β={self.beta}, γ={self.gamma}") print(f"Initial H = {H0:.6f}, Steps: {steps}") print("=" * 70) for step in range(steps): S, Psi, dS_dt, dPsi_dt = self.rk4_step(S, Psi, dS_dt, dPsi_dt) if step % save_interval == 0: t = step * self.dt H = self.compute_hamiltonian(S, Psi, dS_dt, dPsi_dt) H_drift = abs(H - H0) / abs(H0) if abs(H0) > 0 else 0 psi_norm = np.sum(np.abs(Psi)**2) * self.dx**2 S_norm = np.sum(S**2) * self.dx**2 plateau_R = self.compute_plateau_radius(S) sat_pen = np.max(np.maximum(0, S - self.s_max)) cx, cy = self.compute_center_of_mass(Psi) self.history['step'].append(step) self.history['t'].append(t) self.history['energy'].append(H) self.history['norm_psi'].append(psi_norm) self.history['norm_S'].append(S_norm) self.history['plateau_radius'].append(plateau_R) self.history['saturation_penetration'].append(sat_pen) self.history['cm_x'].append(cx) self.history['cm_y'].append(cy) # Boundary monitoring R_over_L = plateau_R / self.R_max if self.R_max > 0 else 0 if verbose and step % (steps // 10) == 0 and step > 0: boundary_str = "" if R_over_L > 0.8: boundary_str = " ⚠️ BOUNDARY WARNING (R > 80%)" elif R_over_L > 0.6: boundary_str = " ⚠️ BOUNDARY CAUTION (R > 60%)" print(f" Step {step:6d}/{steps}: H_drift={H_drift:.6e}, R={plateau_R:.3f}, R/L={R_over_L:.2f}{boundary_str}") if verbose: H_final = self.history['energy'][-1] H_drift = abs(H_final - H0) / abs(H0) R_final = self.history['plateau_radius'][-1] if self.history['plateau_radius'] else 0 print(f"\nFinal H_drift = {H_drift:.6f}, R_final = {R_final:.3f}") print("=" * 70) print("Simulation complete") print("=" * 70) return self.history def analyze_results(self, steps, L, R_max): """Tier-2 analysis with boundary flags and CM tracking.""" if len(self.history['energy']) < 2: return None E0 = self.history['energy'][0] Ef = self.history['energy'][-1] energy_drift = abs(Ef - E0) / abs(E0) if abs(E0) > 0 else 0 R = np.array(self.history['plateau_radius']) R0 = R[0] if len(R) > 0 else 0 Rf = R[-1] if len(R) > 0 else 0 Rmax = np.max(R) if len(R) > 0 else 0 # Center of mass drift cm_x = np.array(self.history['cm_x']) cm_y = np.array(self.history['cm_y']) if len(cm_x) > 1: v_cm_x = (cm_x[-1] - cm_x[0]) / (steps * self.dt) if steps * self.dt > 0 else 0 v_cm_y = (cm_y[-1] - cm_y[0]) / (steps * self.dt) if steps * self.dt > 0 else 0 v_cm = np.sqrt(v_cm_x**2 + v_cm_y**2) else: v_cm = 0.0 t_final = steps * self.dt v_R = (Rf - R0) / t_final if t_final > 0 else 0 eta = Rf / R0 if R0 > 0 else 0 R_over_L = Rf / R_max if R_max > 0 else 0 # Boundary flag: 0=safe, 1=caution, 2=contaminated if R_over_L > 0.8: boundary_flag = 2 elif R_over_L > 0.6: boundary_flag = 1 else: boundary_flag = 0 sat_pen = max(self.history['saturation_penetration']) if self.history['saturation_penetration'] else 0 saturation_flag = 1 if sat_pen > 0 else 0 # Confinement classification (preliminary) if boundary_flag == 0: if eta < 1.1 and v_R < 0.01: confinement_class = "CANDIDATE_CONFINEMENT" elif v_R > 0.05: confinement_class = "DISPERSIVE" else: confinement_class = "MARGINAL" else: confinement_class = "BOUNDARY_CONTAMINATED" return { 'energy_drift': energy_drift, 'R0': R0, 'Rf': Rf, 'Rmax': Rmax, 'v_R': v_R, 'eta': eta, 'R_over_L': R_over_L, 'boundary_flag': boundary_flag, 'saturation_flag': saturation_flag, 'v_cm': v_cm, 'confinement_class': confinement_class, 'H0': E0, 'Hf': Ef } # ============================================================ # PART 2: TIER-2 SWEEP EXECUTOR # ============================================================ class Tier2Sweep: def __init__(self): # Fixed numerical parameters self.base_params = { 'dt': 0.001, 'c_S': 1.0, 'c_Psi': 1.0, 'm2': 0.1, 's_max': 2.0 } # Grid configurations self.grids = [ {'N': 64, 'dx': 0.4, 'label': 'G64_dx04'}, {'N': 128, 'dx': 0.2, 'label': 'G128_dx02'}, {'N': 128, 'dx': 0.4, 'label': 'G128_dx04'}, {'N': 256, 'dx': 0.2, 'label': 'G256_dx02'} ] # Physical parameter sweeps self.kappa_values = [0.0, 0.05, 0.10, 0.20, 0.40] self.beta_values = [0.2, 0.5, 1.0] self.gamma_values = [0.1, 0.2, 0.4] # Step counts self.step_values = [1000, 5000, 10000] # Output self.timestamp = datetime.now().strftime("%Y%m%d_%H%M%S") self.run_id = f"TIER2_SWEEP_{self.timestamp}" self.output_dir = os.path.join(os.getcwd(), self.run_id) os.makedirs(self.output_dir, exist_ok=True) self.ledger_file = os.path.join(self.output_dir, "FRCMFD_TIER2_SUMMARY.csv") self._init_ledger() self.results = [] def _init_ledger(self): with open(self.ledger_file, 'w', newline='') as f: writer = csv.writer(f) writer.writerow([ "run_id", "timestamp", "N", "dx", "dt", "L", "R_max", "kappa", "beta", "gamma", "steps", "H0", "Hf", "delta_H", "R0", "Rf", "Rmax", "v_R", "eta", "R_over_L", "boundary_flag", "saturation_flag", "v_cm", "confinement_class", "solver_version" ]) print(f"[LEDGER] {self.ledger_file} initialized") def _append_ledger(self, row): with open(self.ledger_file, 'a', newline='') as f: writer = csv.writer(f) writer.writerow(row) def _get_solver_version(self): return "TIER2_v1.0" def run_grid_refinement(self): """Run grid refinement test.""" print("\n" + "=" * 70) print("GRID REFINEMENT PHASE") print("=" * 70) ref_kappa = 0.0 ref_beta = 0.5 ref_gamma = 0.2 for steps in self.step_values: for grid in self.grids: run_id = f"GRID_{grid['label']}_S{steps}" print(f"\n[RUN] {run_id}") params = self.base_params.copy() params.update(grid) params['kappa'] = ref_kappa params['beta'] = ref_beta params['gamma'] = ref_gamma solver = FRCMFDSolver(params) history = solver.run_simulation(steps=steps, save_interval=100, verbose=True) analysis = solver.analyze_results(steps, solver.L, solver.R_max) if analysis: self.results.append({ 'run_id': run_id, 'N': grid['N'], 'dx': grid['dx'], 'L': solver.L, 'R_max': solver.R_max, 'kappa': ref_kappa, 'beta': ref_beta, 'gamma': ref_gamma, 'steps': steps, 'analysis': analysis }) self._append_ledger([ run_id, datetime.now().isoformat(), grid['N'], grid['dx'], params['dt'], solver.L, solver.R_max, ref_kappa, ref_beta, ref_gamma, steps, analysis['H0'], analysis['Hf'], analysis['energy_drift'], analysis['R0'], analysis['Rf'], analysis['Rmax'], analysis['v_R'], analysis['eta'], analysis['R_over_L'], analysis['boundary_flag'], analysis['saturation_flag'], analysis['v_cm'], analysis['confinement_class'], self._get_solver_version() ]) print(f" ✅ {run_id} complete") else: print(f" ❌ {run_id} failed") def run_manifold_sweep(self): """Run the full 45-node manifold sweep.""" print("\n" + "=" * 70) print("MANIFOLD SWEEP PHASE (45 NODES)") print("=" * 70) grid = {'N': 128, 'dx': 0.2} steps = 5000 total_nodes = len(self.kappa_values) * len(self.beta_values) * len(self.gamma_values) node_count = 0 for kappa in self.kappa_values: for beta in self.beta_values: for gamma in self.gamma_values: node_count += 1 run_id = f"T2_K{kappa:.2f}_B{beta:.1f}_G{gamma:.1f}" print(f"\n[RUN {node_count}/{total_nodes}] {run_id}") print(f" κ={kappa}, β={beta}, γ={gamma}") params = self.base_params.copy() params.update(grid) params['kappa'] = kappa params['beta'] = beta params['gamma'] = gamma solver = FRCMFDSolver(params) history = solver.run_simulation(steps=steps, save_interval=100, verbose=True) analysis = solver.analyze_results(steps, solver.L, solver.R_max) if analysis: self.results.append({ 'run_id': run_id, 'N': grid['N'], 'dx': grid['dx'], 'L': solver.L, 'R_max': solver.R_max, 'kappa': kappa, 'beta': beta, 'gamma': gamma, 'steps': steps, 'analysis': analysis }) self._append_ledger([ run_id, datetime.now().isoformat(), grid['N'], grid['dx'], params['dt'], solver.L, solver.R_max, kappa, beta, gamma, steps, analysis['H0'], analysis['Hf'], analysis['energy_drift'], analysis['R0'], analysis['Rf'], analysis['Rmax'], analysis['v_R'], analysis['eta'], analysis['R_over_L'], analysis['boundary_flag'], analysis['saturation_flag'], analysis['v_cm'], analysis['confinement_class'], self._get_solver_version() ]) print(f" ✅ {run_id} complete") else: print(f" ❌ {run_id} failed") def compute_gir(self): """Compute Grid Independence Ratio for paired runs.""" print("\n" + "=" * 70) print("GRID INDEPENDENCE RATIO (GIR) ANALYSIS") print("=" * 70) gir_results = [] groups = {} for r in self.results: key = (r['kappa'], r['beta'], r['gamma'], r['steps']) if key not in groups: groups[key] = [] groups[key].append(r) for key, runs in groups.items(): if len(runs) < 2: continue runs.sort(key=lambda x: x['dx']) for i in range(len(runs) - 1): r1 = runs[i] r2 = runs[i+1] if r1['dx'] == r2['dx']: continue delta_H1 = r1['analysis']['energy_drift'] delta_H2 = r2['analysis']['energy_drift'] eta1 = r1['analysis']['eta'] eta2 = r2['analysis']['eta'] GIR_H = abs(delta_H1 - delta_H2) / abs(delta_H1) if abs(delta_H1) > 0 else 0 GIR_eta = abs(eta1 - eta2) / abs(eta1) if abs(eta1) > 0 else 0 grid_independent = 1 if (GIR_H < 0.05 and GIR_eta < 0.05) else 0 gir_results.append({ 'kappa': r1['kappa'], 'beta': r1['beta'], 'gamma': r1['gamma'], 'steps': r1['steps'], 'dx1': r1['dx'], 'dx2': r2['dx'], 'GIR_H': GIR_H, 'GIR_eta': GIR_eta, 'grid_independent': grid_independent }) print(f" κ={r1['kappa']}, β={r1['beta']}, γ={r1['gamma']}") print(f" dx1={r1['dx']}, dx2={r2['dx']}") print(f" GIR_H={GIR_H:.4f}, GIR_eta={GIR_eta:.4f}") print(f" Grid independent: {'✅' if grid_independent else '❌'}") if gir_results: gir_df = pd.DataFrame(gir_results) gir_path = os.path.join(self.output_dir, "FRCMFD_GIR_ANALYSIS.csv") gir_df.to_csv(gir_path, index=False) print(f"\n✅ GIR analysis saved to {gir_path}") return gir_results def generate_heatmaps(self): """Generate heatmaps from the manifold sweep results.""" if not self.results: print("No results to plot") return print("\n" + "=" * 70) print("GENERATING HEATMAPS") print("=" * 70) manifold_runs = [r for r in self.results if r['kappa'] in self.kappa_values] if not manifold_runs: print("No manifold runs found") return df = pd.DataFrame([{ 'kappa': r['kappa'], 'beta': r['beta'], 'gamma': r['gamma'], 'delta_H': r['analysis']['energy_drift'], 'v_R': r['analysis']['v_R'], 'eta': r['analysis']['eta'], 'Rf': r['analysis']['Rf'], 'Rmax': r['analysis']['Rmax'], 'boundary_flag': r['analysis']['boundary_flag'], 'confinement_class': r['analysis']['confinement_class'], 'v_cm': r['analysis']['v_cm'] } for r in manifold_runs]) fig, axes = plt.subplots(2, 3, figsize=(18, 12)) fig.suptitle('FRCMFD v4.0 — Tier-2 Manifold Mapping', fontsize=16, fontweight='bold') # Plot 1: delta_H vs kappa ax = axes[0, 0] for gamma in sorted(df['gamma'].unique()): subset = df[df['gamma'] == gamma] grouped = subset.groupby(['kappa', 'beta'])['delta_H'].mean().reset_index() for beta in sorted(df['beta'].unique()): beta_subset = grouped[grouped['beta'] == beta] ax.plot(beta_subset['kappa'], beta_subset['delta_H'], 'o-', label=f'γ={gamma}, β={beta}', linewidth=1.5, markersize=4) ax.set_xlabel('Kappa (κ)') ax.set_ylabel('Energy Drift (ΔH)') ax.set_title('Energy Drift vs κ') ax.legend(loc='best', fontsize=7) ax.grid(True, alpha=0.3) # Plot 2: v_R vs kappa ax = axes[0, 1] for gamma in sorted(df['gamma'].unique()): subset = df[df['gamma'] == gamma] grouped = subset.groupby(['kappa', 'beta'])['v_R'].mean().reset_index() for beta in sorted(df['beta'].unique()): beta_subset = grouped[grouped['beta'] == beta] ax.plot(beta_subset['kappa'], beta_subset['v_R'], 's-', label=f'γ={gamma}, β={beta}', linewidth=1.5, markersize=4) ax.axhline(y=0, color='black', linestyle='--', alpha=0.5) ax.set_xlabel('Kappa (κ)') ax.set_ylabel('Spreading Rate (v_R)') ax.set_title('Spreading Rate vs κ') ax.legend(loc='best', fontsize=7) ax.grid(True, alpha=0.3) # Plot 3: eta vs kappa ax = axes[0, 2] for gamma in sorted(df['gamma'].unique()): subset = df[df['gamma'] == gamma] grouped = subset.groupby(['kappa', 'beta'])['eta'].mean().reset_index() for beta in sorted(df['beta'].unique()): beta_subset = grouped[grouped['beta'] == beta] ax.plot(beta_subset['kappa'], beta_subset['eta'], '^-', label=f'γ={gamma}, β={beta}', linewidth=1.5, markersize=4) ax.axhline(y=1, color='black', linestyle='--', alpha=0.5) ax.set_xlabel('Kappa (κ)') ax.set_ylabel('Dimensionless Spreading (η)') ax.set_title('Dimensionless Spreading vs κ') ax.legend(loc='best', fontsize=7) ax.grid(True, alpha=0.3) # Plot 4: Confinement class distribution ax = axes[1, 0] class_counts = df['confinement_class'].value_counts() colors = {'CANDIDATE_CONFINEMENT': 'green', 'MARGINAL': 'yellow', 'DISPERSIVE': 'blue', 'BOUNDARY_CONTAMINATED': 'red'} ax.bar(class_counts.index, class_counts.values, color=[colors.get(c, 'gray') for c in class_counts.index]) ax.set_xlabel('Confinement Class') ax.set_ylabel('Count') ax.set_title('Confinement Class Distribution') ax.tick_params(axis='x', rotation=45) ax.grid(True, alpha=0.3) # Plot 5: v_R vs boundary flag ax = axes[1, 1] for flag in sorted(df['boundary_flag'].unique()): subset = df[df['boundary_flag'] == flag] ax.scatter([flag]*len(subset), subset['v_R'], alpha=0.6, label=f'Flag {flag}') ax.set_xlabel('Boundary Flag') ax.set_ylabel('Spreading Rate (v_R)') ax.set_title('v_R vs Boundary Flag') ax.legend() ax.grid(True, alpha=0.3) # Plot 6: v_cm distribution ax = axes[1, 2] ax.hist(df['v_cm'], bins=20, color='purple', alpha=0.7) ax.set_xlabel('Center of Mass Velocity (v_cm)') ax.set_ylabel('Frequency') ax.set_title('CM Velocity Distribution') ax.grid(True, alpha=0.3) plt.tight_layout() plot_path = os.path.join(self.output_dir, 'FRCMFD_TIER2_HEATMAPS.png') plt.savefig(plot_path, dpi=150) print(f"✅ Saved: {plot_path}") plt.show() def create_master_zip(self): """STEP 2: Create master ZIP archive.""" zip_name = f"{self.run_id}.zip" zip_path = os.path.join(os.getcwd(), zip_name) with zipfile.ZipFile(zip_path, 'w', zipfile.ZIP_DEFLATED) as zf: for root, dirs, files in os.walk(self.output_dir): for file in files: file_path = os.path.join(root, file) arcname = os.path.relpath(file_path, os.path.dirname(self.output_dir)) zf.write(file_path, arcname) print(f"[STEP 2] Master ZIP created: {zip_path}") return zip_path def backup_to_drive(self, zip_path): """STEP 3: Backup to Google Drive.""" if not IN_COLAB: print("[STEP 3] Drive backup skipped (not in Colab)") return False try: drive_base = "/content/drive/MyDrive/FRCMFD_Data" os.makedirs(drive_base, exist_ok=True) drive_dir = os.path.join(drive_base, self.run_id) shutil.copytree(self.output_dir, drive_dir, dirs_exist_ok=True) drive_zip = os.path.join(drive_base, os.path.basename(zip_path)) shutil.copy2(zip_path, drive_zip) print(f"[STEP 3] Backed up to Drive: {drive_dir}") return True except Exception as e: print(f"[STEP 3] Drive backup failed: {e}") return False def download_to_local(self, zip_path): """STEP 4: Download to local machine.""" try: if IN_COLAB: files.download(zip_path) print("[STEP 4] Download triggered") return True else: print("[STEP 4] Not in Colab — manual download required") return False except Exception as e: print(f"[STEP 4] Download failed: {e}") return False def verify_files(self): """STEP 5: Verify files exist.""" checks = { 'output_dir': os.path.exists(self.output_dir), 'ledger': os.path.exists(self.ledger_file) } print("[STEP 5] Verification:") for check, status in checks.items(): print(f" {'✅' if status else '❌'} {check}") return all(checks.values()) def final_status_report(self, zip_path, drive_success, verify_success): """STEP 6: Final status report.""" print("\n" + "=" * 50) print("✓ Colab workspace saved") print("✓ Google Drive backup saved" if drive_success else "⚠️ Drive backup skipped") print("✓ Download package created") print("=" * 50) print("\n=== FINAL STATUS REPORT ===") print(f"OUTPUT DIRECTORY: {os.path.abspath(self.output_dir)}") print(f"GOOGLE DRIVE BACKUP: {drive_success}") print(f"MASTER ZIP: {zip_path}") print(f"FILE COUNT: {len(os.listdir(self.output_dir))}") print(f"STATUS: {'SUCCESS' if verify_success else 'FAILURE'}") print("===========================") def run_suite(self): """Execute the full Tier-2 suite.""" print(f"\n[STEP 1] Output directory: {self.output_dir}") # Part 1: Grid refinement self.run_grid_refinement() # Part 2: Manifold sweep (45 nodes) self.run_manifold_sweep() # Part 3: GIR analysis gir = self.compute_gir() # Part 4: Heatmaps self.generate_heatmaps() # Part 5: Data preservation zip_path = self.create_master_zip() drive_success = self.backup_to_drive(zip_path) self.download_to_local(zip_path) verify_success = self.verify_files() self.final_status_report(zip_path, drive_success, verify_success) print("\n" + "=" * 70) print("✅ TIER-2 SUITE COMPLETE") print("=" * 70) print(f"Total runs: {len(self.results)}") print(f"Ledger saved to: {self.ledger_file}") print(f"Output directory: {self.output_dir}") print("=" * 70) return self.results, gir # ============================================================ # PART 3: EXECUTE # ============================================================ if __name__ == "__main__": if IN_COLAB: drive.mount('/content/drive', force_remount=True) sweep = Tier2Sweep() results, gir = sweep.run_suite() print("\n" + "=" * 70) print("✅ TIER-2 EXECUTION COMPLETE") print("=" * 70) print(f"Results: {len(results)} nodes processed") print(f"Ledger: {sweep.ledger_file}") print("=" * 70)