Update 3/5: Wind Tunnel Runs 006AB&D, as well as reviewing Run 005B

Work in the past two weeks has been focused on testing an AMR scheme that only forces maximal refinement within a radius of the particle. Currently the module is working fine, except that it cannot have the box be relaxed first before introducing the particle.

Finished Runs

WT Run 006A

Same setup with Run 004A. No accretion, highest resolution 0.5CU, softening length (spline) is resolved by 8 pixels. Radius of highest refinement is set to 40CU, which is more than twice as large as the Bondi radius 18.7CU. Significant run time improvement (order of 8) & low storage requirement (~300M per frame vs. ~4G per frame w/ normal AMR).

WT Run 006B

Similar setup with Run 006A, except resolution is improved by a factor of two to 0.25CU, and softening length is resolved by 16 pixels. To compensate, radius of highest refinement is set to 20CU, still larger than Bondi radius. Similar run time as the previous run.

Log density movie with mesh grids:

WT Run 006C

Similar to 006B, but with radius of refinement set at 40CU. The runtime becomes too long (longer than the normal AMR), and was left unfinished.

WT Run 006D

Similar to 006A, low resolution, but larger radius of highest refinement - 80CU. The run time is less than 2 days, similar to the AMR case. But the result looks unreliable.

Log density movie with mesh grids and velocity vectors:

Comparison with Run 004A

More one Run 005B

Zoomed-in log density movie: Zoomed-in Mach number and velocity vector movie;

Attachments (7)

Comments

1. Baowei Liu -- 7 years ago

Hi Bo, you may try put the SetErrorFlag code in SUBROUTINE ProblemSetErrFlag(Info) inside this IF statement:

    IF (lPlaceParticle) THEN
    END IF