Changes between Version 46 and Version 47 of 1DPulsedJets


Ignore:
Timestamp:
03/30/12 12:45:47 (13 years ago)
Author:
ehansen
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • 1DPulsedJets

    v46 v47  
    225225
    226226
    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 rho1 is ambient density.
    228 
    229 The post-shock density and pressure (rho2 and p2 respectively) can be found by using mass flux and momentum flux conservation across the shock:
     227where 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
     229The 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:
    230230
    231231These post-shock values become the boundary conditions for the fluid equations in the cooling region:
    232232
    233 where rho, v, and p represent the density, velocity, and pressure in the cooling region as functions of x, and lambda is the cooling rate.
     233where [[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.
    234234
    235235[[BR]]
    236236=== Initial Parameters ===
    237237n1 = 60 particles/cc
     238
    238239v1 = 10^7^ cm/s
     240
    239241T1 = 10^4^ K
    240242
    241243For analytic cooling,
    242 [[latex($\Lambda = n^{2}*\alpha*T^{\beta}$)]]
     244
     245[[latex($\Lambda = n^{2} \alpha T^{\beta}$)]]
     246
    243247where n and T and the number density and temperature in the post-shock region respectively.
    244248
    245249For these simulations,
    246250
    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^
    249254
    250255cell length = 10^18^ cm
     256
    251257problem domain = 400 cells
     258
    252259final time ~ 4000 years
    253260