Changes between Version 46 and Version 47 of 1DPulsedJets
- Timestamp:
- 03/30/12 12:45:47 (13 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
1DPulsedJets
v46 v47 225 225 226 226 227 where v1 is the ambient velocity, and M is the ambient mach number. Remember that the mach number M = v1/c where c is the ambient sound speed, and c = sqrt(gamma*p1/rho1) where p1 is ambient pressure, and rho1is ambient density.228 229 The post-shock density and pressure ( rho2 and p2respectively) can be found by using mass flux and momentum flux conservation across the shock:227 where v1 is the ambient velocity, and M is the ambient mach number. Remember that the mach number M = v1/c where c is the ambient sound speed, and [[latex($c = \sqrt{\frac{\gamma P1}{\rho1}}$)]] where P1 is ambient pressure, and [[latex($\rho1$)]] is ambient density. 228 229 The post-shock density and pressure ([[latex($\rho2$)]] and P_2_ respectively) can be found by using mass flux and momentum flux conservation across the shock: 230 230 231 231 These post-shock values become the boundary conditions for the fluid equations in the cooling region: 232 232 233 where rho, v, and p represent the density, velocity, and pressure in the cooling region as functions of x, and lambdais the cooling rate.233 where [[latex($\rho$)]], v, and P represent the density, velocity, and pressure in the cooling region as functions of x, and [[latex($\Lambda$)]] is the cooling rate. 234 234 235 235 [[BR]] 236 236 === Initial Parameters === 237 237 n1 = 60 particles/cc 238 238 239 v1 = 10^7^ cm/s 240 239 241 T1 = 10^4^ K 240 242 241 243 For analytic cooling, 242 [[latex($\Lambda = n^{2}*\alpha*T^{\beta}$)]] 244 245 [[latex($\Lambda = n^{2} \alpha T^{\beta}$)]] 246 243 247 where n and T and the number density and temperature in the post-shock region respectively. 244 248 245 249 For these simulations, 246 250 247 beta = 2 248 alpha = 1.2345 * 10^-34^ erg*cm^3^/s/K^2^ 251 [[latex($\beta$)]] = 2 252 253 [[latex($\alpha$)]] = 1.23786 * 10^-34^ erg*cm^3^/s/K^2^ 249 254 250 255 cell length = 10^18^ cm 256 251 257 problem domain = 400 cells 258 252 259 final time ~ 4000 years 253 260