Changes between Version 1 and Version 2 of AstrobearRestarts

Timestamp:

10/24/14 11:47:38 (10 years ago)

Author:

Jonathan

Comment:

—

AstrobearRestarts

@@ -43 +43 @@
   END IF
 }}}
+
+
+=== How max speeds are calculated ===
+==== Elliptic max speed depends on the solver ====
+* Poisson
+{{{
+elliptic_maxspeed(level)=max(elliptic_maxspeed(level),sqrt(nDim*&
+                maxval(abs(Info%q(ir(1,1):ir(1,2),ir(2,1):ir(2,2),ir(3,1):ir(3,2),iPhiGas)-&
+                Info%q(il(1,1):il(1,2),il(2,1):il(2,2),il(3,1):il(3,2),iPhiGas)))))
+}}}
+==== Hydro max speed ====
+* MUSCL uses maximum speed from the riemann solvers - which is usually just the max of the left and right most waves
+{{{
+  max(ABS(SL),ABS(SR))
+}}}
+* sweep uses the GetMaxSpeed routine to find the fastest eigen speed on the grid after an update
+{{{
+   SUBROUTINE GetMaxSpeed(Info)
+      TYPE(InfoDef)  :: Info
+      REAL(KIND=qPrec), DIMENSION(:,:,:,:), POINTER :: q
+      REAL(KIND=qPREC) :: myMaxSpeed, fast_speed(3), A2, B2, B, dx, A
+      INTEGER :: i,j,k,level
+      level=Info%level
+      dx=levels(level)%dx
+      q=>Info%q
+      myMaxSpeed=0d0
+      DO i=1, Info%mX(1)
+         DO j=1, Info%mX(2)
+            DO k=1, Info%mX(3)
+               A=SoundSpeed(q(i,j,k,:))
+               A2=q(i,j,k,1)*A**2
+               IF (iE_rad /= 0) THEN
+                  A2=A2+4d0/9d0*q(i,j,k,iE_rad)*(1d0-exp(-q(i,j,k,iKappa)*dx))
+                  A=sqrt(A2/q(i,j,k,1))
+               END IF
+               IF (lMHD) THEN
+                  B2 = SUM(q(i,j,k,iBx:iBz)**2)
+                  B = A2+B2
+                  fast_speed(1:nDim) = sqrt(half*(B+sqrt(B**2-4d0*A2*q(i,j,k,iB(1:nDim))**2))/q(i,j,k,1))
+                  mymaxspeed= maxval(abs(q(i,j,k,imom(1:nDim)))/q(i,j,k,1)+fast_speed(1:nDim))
+               ELSE
+!                  mymaxspeed=maxval(abs(q(i,j,k,imom(1:nDim)))+A2)/q(i,j,k,1)                                                                                
+                  mymaxspeed=sqrt(sum(q(i,j,k,imom(1:nDim))**2))/q(i,j,k,1)+A
+               END IF
+               IF (mymaxspeed > maxspeed(level)) THEN
+                  maxspeed(level)=mymaxspeed
+                  maxspeed_location(level,:)=CellPos(Info, i, j, k)
+               END IF
+            END DO
+         END DO
+      END DO
+   END SUBROUTINE GetMaxSpeed
+}}}
+
+==== Particle max speed ====
+This just finds the fastest particle speed
+{{{
+ParticleMaxSpeed=max(ParticleMaxSpeed, sqrt(sum(vel(i,1:nDim)**2)))
+}}}