Changes between Version 20 and Version 21 of GravoTurbulence


Ignore:
Timestamp:
11/08/11 16:19:00 (13 years ago)
Author:
Jonathan
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • GravoTurbulence

    v20 v21  
    3030
    3131
     32In all of these, there was initial pressure equilibrium (so different temperatues), although I was using a gamma=1.001 which assumes an isothermal gas.  Consider what happens when a gas at a density of 1 and a temperature of 100 K mixes with a gas at a density of 100 and a temperature of 1 K. 
     33
     34The new density is [[latex($\rho=100+1=101$)]]
     35
     36The new energy is [[latex($E=10^5+10^5=2\times 10^5$)]]
     37
     38The new pressure is [[latex($P=2\times 10^2$)]]
     39
     40And the new temperature is [[latex($T=200/101 \approx 2$)]]
     41
     42So the dense gas at 100 K is instantly heated to 2 K and the equation of state is violated
     43
     44
     45If instead the gas is at the same temperature everywhere, then initial differences in density will cause expansion
     46
     47[[Image(MultiPoleCollapse.png, width=400)]]
     48
     49[attachment:MultiPoleCollapse.gif movie]
     50
     51
     52
     53
     54
     55
     56
     57
     58
     59
     60
     61
     62
     63
     64
     65
     66
     67
     68
     69
     70
     71
     72
     73
     74
     75
     76
     77
     78
     79
     80
     81
     82An alternative is to use the IICooling and to start with a much larger domain at a lower density with a small perturbation that will cause it to self-consistently develop large density perturbations.
     83
     84
     85[[Image(MPCollBig.png, width=400)]]
     86
     87[attachment:MPCollBig.gif movie]
     88
     89
     90
     91
     92
     93
     94
     95
     96
     97
     98
     99
     100
     101
     102
     103
     104
     105
     106
     107
     108
     109
     110
     111
     112
     113
     114
     115
     116
     117
     118
     119
     120
     121
     122
     123
     124
     125
     126
     127
    32128
    33129== Velocity Perturbations ==
    34130
    35131Virial equilibrium requires [[latex($2T=V$)]] where [[latex($V=\frac{3GM^2}{5R}$)]]
     132
     133Or in the case of 2D [[latex($V=\frac{GM}{2}(\log(R)-\frac{1}{4})$)]]
     134
     135If the kinetic energy is put into solenoidal motions with a spectral index [[latex($\alpha = -1$)]] so that [[latex($\hat{v}(k)\propto k^{\alpha}$)]], then energy dissipates fairly quickly and the whole thing collapses
     136
     137
     138[[Image(MPCollVir.png, width=400)]]
     139
     140[attachment:MPCollVir.gif movie]
    36141
    37142