| 4 | | |
| 5 | | == THE SIGNIFICANCE OF R_jet: Part 1== |
| 6 | | |
| 7 | | In the past few days we've been musing about the relationship between Rjet and the shape of the developing flow. Muse no more. Rjet matters!!!! It is a shape- and size-controlling variable of paramount significance (forcing us to make some strategic choices in this project)! |
| 8 | | Rjet really matters!!!! |
| 9 | | Rjet really matters!!!! |
| 10 | | Rjet really matters!!!! |
| 11 | | Rjet really matters!!!! |
| 12 | | Here's the originating problem.data file for one of the runs (the only variable that changes is Rjet) |
| 13 | | |
| 14 | | &ProblemData |
| 15 | | ! MODEL: jet 0deg (n=4e4, r=1000AU, T=100K) @ 200 km/s into 4e4 torus + AGB wind |
| 16 | | ! folder TapAGB45Rhn4e2v200namb4e4 |
| 17 | | ! |
| 18 | | ! BACKGROUND or “AMBIENT” SECTION. Values apply to origin |
| 19 | | tamb = 1d3 ! ambient temp, 1cu = 0.1K (100K=1000cu) |
| 20 | | namb = 4e4 ! ambient central density cm^-3. Usually 400 for 1/r^2 or tor$ |
| 21 | | stratified = t ! true = add a 1/r^2 background 'AGB stellar wind' |
| 22 | | torus = f ! true - add torus to the background |
| 23 | | torusalpha = 0.7 ! alpha and beta specify the geometry |
| 24 | | torusbeta = 10d0 ! see Frank & Mellema, 1994ApJ...430..800F |
| 25 | | rings = f ! true - add radial density modulations to AGB wind |
| 26 | | ! |
| 27 | | ! FLOW DESCRIPTION SECTION, values apply at origin at t=0 |
| 28 | | outflowType = 2 ! TYPE OF FLOW 1 cyl jet, 2 conical wind, 3 is clump |
| 29 | | njet = 4d2 ! flow density at launch zone, 1cu = 1cm^-3 |
| 30 | | Rjet = 0.5, 1.0. 2.0 ! flow radius at launch zone, 1cu = 500AU (outflowType=1 only) |
| 31 | | vjet = 2e7 ! flow velocity , 1cu = cm/s (100km/s=1e7cu) |
| 32 | | tjet = 1d3 ! flow temp, 1cu = 0.1K (100K=1000cu) |
| 33 | | tt = 0.0d0 ! flow accel time, 1cu = 8250y (0.02 = 165y) |
| 34 | | open_angle = 90d0 ! conical flow open angle (deg) |
| 35 | | tf = 45d0 ! conical flow Gaussian taper (deg) for njet and vjet; 0= disa$ |
| 36 | | sigma = 0d0 ! !toroidal.magnetic.energy / kinetic.energy, example 0.6 |
| 37 | | |
| 38 | | ! |
| 39 | | ! OTHER PARAMETERS |
| 40 | | lcooling = t ! radiative cooling? |
| 41 | | buff = 8 ! central refinement of a grid with a resolution 1/2 |
| 42 | | / |
| 43 | | |
| 44 | | |
| 45 | | |
| 46 | | Comments about the differences: |
| 47 | | |
| 48 | | The primary reason to make Rjet large is to better resolve its surface (which is the lower boundary of the flow) |
| 49 | | The primary reason to make Rjet small is to wipe out less of the background. This is surprisingly important. |
| 50 | | |
| 51 | | Comments about the flow size: |
| 52 | | |
| 53 | | 1. Increasing the surface area while holding njet fixed changes the mass and momentum of the flow. That is, njet is the density at the launching surface. Increasing Rjet increases the launch surface area and the total flow momentum. |
| 54 | | |
| 55 | | 2. Increasing Rjet means that the launch sphere covers more of the ambient medium. It is the innermost zone of the ambient medium where it density is largest, so flows starting at a large Rjet never see much of the density in the extrenal environment. This is especially important for tori (not used in the sims above). |
| 56 | | |
| 57 | | Both of these conditions (1 and 2) mean that flows launched with a large Rjet will penetrate significantly further. |
| 58 | | |
| 59 | | Comments about the flow shape: |
| 60 | | |
| 61 | | 1. About large-scale shapes: Much of the shaping of collimated flows occurs when they are very small. Launching them at Rjet = 0.5 or Rjet = 2 can have a profound impact on their overall shapes. |
| 62 | | |
| 63 | | 2. About features near the y axis: These change profoundly with the size of Rjet (that is, the details of the pixelated shape of the launch sphere), as we suspected that they would. Look at the attached figure near the y axis. Small Rjets are almost certain to produce dense artifacts near the y axis. However, there are still residual artifacts even when Rjet=2. |
| 64 | | |
| 65 | | 3. Look at the right panel of the figure and then at the others to its left. In the right panel you might imagine that the thin-shell instability has disturbed the surface. Not so fast. The shape of the wiggles along the boundary at 440 y is very clearly related to the wiggles at 110y -- all of which are flow artifacts. Thus the ripples at 440y are born of "original" artifacts. |
| 66 | | This problem becomes more severe when Rjet is smaller. But the artificial ripples in the large-scale flow are presisiten with any value of Rjet. And these artificial features are not confined to the region of the y axis. They're all along the leading edge of the flow! |
| 67 | | |
| 68 | | Conclusions: |
| 69 | | |
| 70 | | 1. For many types of flows (but probably not dense cylindrical jets or bullets) Rjet should be at least 2.0. Even then, of wiggles appear at the leading edge of the flow near the y-axis at early times then you can't trust the resulting structures as the flow edge grows. Any wiggles or ripples will be hydrodynamically amplified in time. The only way to control these is to run at higher spatial resolution. |
| 71 | | |
| 72 | | One way to remove some of the early wiggles that develop near the y axis is to launch the flow into the grid at 45 degrees. |
| 73 | | |
| 74 | | [This may be easier than it sounds. Martin already has a way of modulating the flow density of conical outflows around the boundary with a gaussian taper that he applies to a spherical wind. The form of the taper is |
| 75 | | {{{exp(-theta/user-half-width)^2}}}. That algorithm might be changed to {{{exp(-(theta-45deg)/user-half-width)^2}}}.If you want to sharpen the edges (by suppressing the Gaussian wings) then the taper can be {{{exp(-(theta-45deg)/user-half-width)^X}}} where X>2. |
| 76 | | |
| 77 | | The lay-down of the external medium prior to the start of the flow MUST somehow adjust itself for the size of the luanch sphere. After a moment's though I suspect that this can be done by scaling njet so that the density of the initial medium at the launch surface is held fixed as Rjet is changed. |
| 78 | | |
| 79 | | Note that njet ALWAYS describes the flow density at the launch surface. However, namb does not describe the density of the external medium at the launch surface. |
