Changes between Version 4 and Version 5 of u/erica/GudonovMethodEuler

Timestamp:

04/23/13 13:05:24 (12 years ago)

Author:

Erica Kaminski

Comment:

—

u/erica/GudonovMethodEuler

@@ -21 +21 @@
 where F(i-1/2) is the numerical flux, which = F(u(i+1/2, 0)) (i.e. the physical flux evaluated using the solution to local RP)
 
-where delta_t is constrained by the delta_t <= delta_x/max[S(n)] <- note this allows a wave to travel a full delta_x in a time-step
+where delta_t is constrained by the delta_t <= delta_x/max[S(n)] <- here s(n) is a given wave speed. note this allows a wave to travel a full delta_x in a time-step
 
 we don't worry about wave-wave interaction based on linearity -- that is, assume wave interactions do not lead to wave accelerations
 
+we can specify the time constraint as follows, delta_t = cfl*delta_x/max[s(n)], if we let 0<cfl<=1. a cfl = 1 gives the most efficient time marching scheme (the largest time steps). max[s(n)] finds the max speed on the entire domain. 
+
 u(i, n+1) = average over the cell