3D Mach stem update

These are the base parameters that I have used for previous simulations. The ambient temperature could be lowered if we want higher Mach numbers for the clumps.

Ambient/Wind Clumps
density (1/cc) 1e3 5e5
velocity (km/s) 27.2644 0, 10
temperature (K) 5e3 10

So the density contrast is 500, and the clumps are initially in pressure equilibrium with the ambient.

We are in the reference frame of the slow clump, which is why the ambient is moving and the slow clump is given velocity = 0. Thus the slow clump has a Mach number of approximately 3.66, and the fast clump is at approximately M = 5.

The total simulation time ~ 50 years.

The clump separation is 3 clump radii.

I have completed 2 runs: the base run described above, and a run with twice the clump densities. Thus the higher density run has a nclump/nambient ratio of 1000. Below are the emission maps from the 2 runs at 4 different inclination angles: 0, 30, 60, and 90 degrees.

movie

We do not have to continue with runs that have different density contrasts. Perhaps it would be better to study how clump separation and velocity affects the emission. I have a grid of runs below that we might try:

Clump Separation (rclump) Fast Clump M
3 5
3 10
6 5
6 10

So this is just 4 runs to start. I would want to leave the slow clump alone, so it would remain at M = 3.66. The exact numbers can change, but I was basically thinking of doing what I already have + investigating a faster clump + investigating a clump that is farther away.

Also, do we want to move the clump so far away that it does not travel through the slow clump's wake? Or perhaps this is desirable, and we want to see an even closer approach?

Attachments (2)

Download all attachments as: .zip

Comments

No comments.