File tree Expand file tree Collapse file tree 1 file changed +9
-9
lines changed
PySDM/physics/homogeneous_liquid_nucleation_rate Expand file tree Collapse file tree 1 file changed +9
-9
lines changed Original file line number Diff line number Diff line change @@ -13,20 +13,20 @@ def __init__(self, _):
13
13
14
14
@staticmethod
15
15
def j_liq_homo (const , T , S , e_s ):
16
- m1 = const .Mv / const .N_A # kg per molecule
17
- v1 = m1 / const .rho_w # m3 per molecule
18
- N1 = (S * e_s ) / (m1 * const .Rv * T ) # molecules per m3
16
+ mass_per_moleculue = const .Mv / const .N_A
17
+ volume_per_molecule = mass_per_moleculue / const .rho_w
18
+ N1 = (S * e_s ) / (mass_per_moleculue * const .Rv * T ) # molecules per m3
19
19
return (
20
- ((2 * const .sgm_w ) / (np .pi * m1 )) ** (1 / 2 )
21
- * (v1 * N1 ** 2 / S )
20
+ ((2 * const .sgm_w ) / (np .pi * mass_per_moleculue )) ** (1 / 2 )
21
+ * (volume_per_molecule * N1 ** 2 / S )
22
22
* np .exp (
23
- (- 16 * np .pi * v1 ** 2 * const .sgm_w ** 3 )
23
+ (- 16 * np .pi * volume_per_molecule ** 2 * const .sgm_w ** 3 )
24
24
/ (3 * const .k_B ** 3 * T ** 3 * np .log (S ) ** 2 )
25
25
)
26
26
)
27
27
28
28
@staticmethod
29
29
def r_liq_homo (const , T , S ):
30
- m1 = const .Mv / const .N_A # kg per molecule
31
- v1 = m1 / const .rho_w # m3 per molecule
32
- return (2 * const .sgm_w * v1 ) / (const .k_B * T * np .log (S ))
30
+ mass_per_moleculue = const .Mv / const .N_A
31
+ volume_per_molecule = mass_per_moleculue / const .rho_w
32
+ return (2 * const .sgm_w * volume_per_molecule ) / (const .k_B * T * np .log (S ))
You can’t perform that action at this time.
0 commit comments