rss
pnStudy: Jet velocity Vs Jet Radius
In the pnStudy, the velocity of the jet is set according to the radius:
!======== J E T=========: IF (outflowType == collimated) then q(i,j,k,itracer4)=1d0 fact=1d0 !10 mar 2014 !fact=exp(-(x**2+z**2)/jet_width**2) ! b 4 10 mar 2014 qjet(1)=njet/nScale*fact v=vjet/velscale*(/0d0,1d0,0d0/)/Rjet*fact*timef !10 mar 2014 !v=vjet/velscale*(/0d0,y,0d0/)/Rjet*fact*timef !b 4 10 mar 2014 qjet(imom(1:nDim))=v(1:nDim)*qjet(1)*& !ramp up velocity mom_flow !5 may 2014, time dependent mom flux requested by bbalick. qjet(iE)=qjet(1)*tjet/TempScale*gamma7
This makes the Jet velocity is not vJet (which is given in problem.data) when Rjet !=0.
- Rjet=2d0, vJet-2e7 as in Problem.data
outflowType = 1 ! TYPE OF FLOW 1 cyl jet, 2 conical wind, 3 is clump njet = 4d4 ! flow density at launch zone, 1cu = 1cm^-3 Rjet = 2d0 ! flow radius at launch zone, 1cu = 500AU vjet = 2e7 ! flow velocity , 1cu = cm/s (100km/s=1e7cu) tjet = 1000d0 ! flow temp, 1cu = 0.1K (100K=1000cu) tt = 0.0d0 ! flow accel time, 1cu = 8250y (0.02 = 165y) open_angle = 00d0 ! conical flow open angle (deg) tf = 15d0 ! conical flow Gaussian taper (deg) for njet and vjet; 0= disable sigma = 0d0 ! !toroidal.magnetic.energy / kinetic.energy, example 0.6
Here's the file of Scales.data
TIMESCALE = 260347122628.507 , LSCALE = 7.479899800000000E+015, MSCALE = 2.099937121547526E+026, RSCALE = 5.017864740000001E-022, VELSCALE = 28730.4876830661 , PSCALE = 4.141950900000000E-013, NSCALE = 300.000000000000 , BSCALE = 2.281431350619136E-006, TEMPSCALE = 10.0000000000000 , SCALEGRAV = 2.269614763656989E-006
Velocity plot shows the jet velocity is 100 km/s instead of 200 km/s. And this can also be calculated from the distance the jet travels during 600 yrs:
| t=0 | t=660 y |
 | 
|
- Rjet=1d0, vJet-2e7 as in Problem.data
outflowType = 1 ! TYPE OF FLOW 1 cyl jet, 2 conical wind, 3 is clump njet = 4d4 ! flow density at launch zone, 1cu = 1cm^-3 Rjet = 1d0 ! flow radius at launch zone, 1cu = 500AU vjet = 2e7 ! flow velocity , 1cu = cm/s (100km/s=1e7cu) tjet = 1000d0 ! flow temp, 1cu = 0.1K (100K=1000cu) tt = 0.0d0 ! flow accel time, 1cu = 8250y (0.02 = 165y) open_angle = 00d0 ! conical flow open angle (deg) tf = 15d0 ! conical flow Gaussian taper (deg) for njet and vjet; 0= disable sigma = 0d0 ! !toroidal.magnetic.energy / kinetic.energy, example 0.6
Scales.data are same
Velocity plot shows the jet velocity is 200 km/s. And this can also be calculated from the distance the jet travels during 600 yrs:
| t=0 | t=660 y |
 | 
|
- Rjet=4d0, vJet-2e7 as in Problem.data
outflowType = 1 ! TYPE OF FLOW 1 cyl jet, 2 conical wind, 3 is clump njet = 4d4 ! flow density at launch zone, 1cu = 1cm^-3 Rjet = 4d0 ! flow radius at launch zone, 1cu = 500AU vjet = 2e7 ! flow velocity , 1cu = cm/s (100km/s=1e7cu) tjet = 1000d0 ! flow temp, 1cu = 0.1K (100K=1000cu) tt = 0.0d0 ! flow accel time, 1cu = 8250y (0.02 = 165y) open_angle = 00d0 ! conical flow open angle (deg) tf = 15d0 ! conical flow Gaussian taper (deg) for njet and vjet; 0= disable sigma = 0d0 ! !toroidal.magnetic.energy / kinetic.energy, example 0.6
Scales.data are same
Velocity plot shows the jet velocity is 50 km/s instead of 200 km/s. And this can also be calculated from the distance the jet travels during 600 yrs:
| t=0 | t=660 y |
 | 
|
Comments
No comments.