Version 10 (modified by 12 years ago) ( diff ) | ,
---|

### Tracking the mass that falls onto the BE sphere

I would like to calculate the mass that falls onto a shell of dr=0.1RBe at r=Rbe from the ambient over the course of a simulation. There may be a few ways to do this, but easiest may be what came to mind today: find M(t) of the sphere with r=1.1Rbe, using Visit's query option for total mass in box.

#### The method

First, I had to determine the right sum to use in Visit. I see that the 'query' option had 2 sums: weighted and not. Additionally, I found that the sum query is for the current plot. If that plot is a slice of the full simulation, the query will return only a sum of variables from that slice's mesh. Using my knowledge of the BE sphere's total mass from my simulation output —- that is, 151 solar mass, I calculated the mass of the ambient quite easily, given the ambient is of uniform density:

Putting this into computational units, I discovered that **weighted variable sum** was the correct query to get the total rho/cell (and hence mass) in the grid.

The plan was to calculate M(t, 1.1Rbe) using the following formula:

### Attachments (8)

- M(t)_with_critBE.png (17.8 KB ) - added by 12 years ago.
- t0.png (26.7 KB ) - added by 12 years ago.
- t.276.png (28.5 KB ) - added by 12 years ago.
- t0.567.png (28.1 KB ) - added by 12 years ago.
- ambTimeQuery.png (33.3 KB ) - added by 12 years ago.
- rhoTimeQueryAmb.gif (1.3 MB ) - added by 12 years ago.
- timeQueryRho.png (34.3 KB ) - added by 12 years ago.
- asymMesh.png (71.5 KB ) - added by 12 years ago.

Download all attachments as: .zip

**Note:**See TracWiki for help on using the wiki.