http://rxnkin.usc.es/index.php?action=history&feed=atom&title=Statistical_Thermodynamics_and_Trasnport_Phenomena
Statistical Thermodynamics and Trasnport Phenomena - Revision history
2026-04-19T10:03:04Z
Revision history for this page on the wiki
MediaWiki 1.32.0
http://rxnkin.usc.es/index.php?title=Statistical_Thermodynamics_and_Trasnport_Phenomena&diff=831&oldid=prev
Wikiadmin: Blanked the page
2020-07-21T14:10:19Z
<p>Blanked the page</p>
<a href="http://rxnkin.usc.es/index.php?title=Statistical_Thermodynamics_and_Trasnport_Phenomena&diff=831&oldid=830">Show changes</a>
Wikiadmin
http://rxnkin.usc.es/index.php?title=Statistical_Thermodynamics_and_Trasnport_Phenomena&diff=830&oldid=prev
Wikiadmin at 14:06, 21 July 2020
2020-07-21T14:06:51Z
<p></p>
<table class="diff diff-contentalign-left" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">Revision as of 14:06, 21 July 2020</td>
</tr><tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l1" >Line 1:</td>
<td colspan="2" class="diff-lineno">Line 1:</td></tr>
<tr><td colspan="2"> </td><td class='diff-marker'>+</td><td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div><ins style="font-weight: bold; text-decoration: none;"><html lang="en"></ins></div></td></tr>
<tr><td colspan="2"> </td><td class='diff-marker'>+</td><td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div><ins style="font-weight: bold; text-decoration: none;"></ins></div></td></tr>
<tr><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div><div id="glowscript" class="glowscript"></div></td><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div><div id="glowscript" class="glowscript"></div></td></tr>
<tr><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></div></td><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></div></td></tr>
<tr><td colspan="2" class="diff-lineno" id="mw-diff-left-l524" >Line 524:</td>
<td colspan="2" class="diff-lineno">Line 526:</td></tr>
<tr><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>//--><!]]></script></div></td><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>//--><!]]></script></div></td></tr>
<tr><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></div></div></td><td class='diff-marker'> </td><td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></div></div></td></tr>
<tr><td colspan="2"> </td><td class='diff-marker'>+</td><td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div><ins style="font-weight: bold; text-decoration: none;"></ins></div></td></tr>
<tr><td colspan="2"> </td><td class='diff-marker'>+</td><td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div><ins style="font-weight: bold; text-decoration: none;"></html></ins></div></td></tr>
</table>
Wikiadmin
http://rxnkin.usc.es/index.php?title=Statistical_Thermodynamics_and_Trasnport_Phenomena&diff=829&oldid=prev
Wikiadmin: Created page with "<div id="glowscript" class="glowscript"> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <link type="text/css" href="https://s3.amazonaws.com/glowscript/cs..."
2020-07-21T14:05:25Z
<p>Created page with "<div id="glowscript" class="glowscript"> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <link type="text/css" href="https://s3.amazonaws.com/glowscript/cs..."</p>
<p><b>New page</b></p><div><div id="glowscript" class="glowscript"><br />
<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><br />
<link type="text/css" href="https://s3.amazonaws.com/glowscript/css/redmond/2.1/jquery-ui.custom.css" rel="stylesheet" /><br />
<link type="text/css" href="https://s3.amazonaws.com/glowscript/css/ide.css" rel="stylesheet" /><br />
<script type="text/javascript" src="https://s3.amazonaws.com/glowscript/lib/jquery/2.1/jquery.min.js"></script><br />
<script type="text/javascript" src="https://s3.amazonaws.com/glowscript/lib/jquery/2.1/jquery-ui.custom.min.js"></script><br />
<script type="text/javascript" src="https://s3.amazonaws.com/glowscript/package/glow.3.0.min.js"></script><br />
<script type="text/javascript" src="https://s3.amazonaws.com/glowscript/package/RSrun.3.0.min.js"></script><br />
<script type="text/javascript"><!--//--><![CDATA[//><!--<br />
<br />
// START JAVASCRIPT<br />
;(function() {;<br />
var RS_modules = {};<br />
RS_modules.pythonize = {};<br />
<br />
(function(){<br />
function strings() {<br />
var string_funcs, exclude, name;<br />
string_funcs = set("capitalize strip lstrip rstrip islower isupper isspace lower upper swapcase center count endswith startswith find rfind index rindex format join ljust rjust partition rpartition replace split rsplit splitlines zfill".split(" "));<br />
if (!arguments.length) {<br />
exclude = (function(){<br />
var s = RS_set();<br />
s.jsset.add("split");<br />
s.jsset.add("replace");<br />
return s;<br />
})();<br />
} else if (arguments[0]) {<br />
exclude = Array.prototype.slice.call(arguments);<br />
} else {<br />
exclude = null;<br />
}<br />
if (exclude) {<br />
string_funcs = string_funcs.difference(set(exclude));<br />
}<br />
var RS_Iter0 = string_funcs;<br />
RS_Iter0 = ((typeof RS_Iter0[Symbol.iterator] === "function") ? (RS_Iter0 instanceof Map ? RS_Iter0.keys() : RS_Iter0) : Object.keys(RS_Iter0));<br />
for (var RS_Index0 of RS_Iter0) {<br />
name = RS_Index0;<br />
(RS_expr_temp = String.prototype)[(typeof name === "number" && name < 0) ? RS_expr_temp.length + name : name] = (RS_expr_temp = RS_str.prototype)[(typeof name === "number" && name < 0) ? RS_expr_temp.length + name : name];<br />
}<br />
};<br />
if (!strings.__module__) Object.defineProperties(strings, {<br />
__module__ : {value: "pythonize"}<br />
});<br />
<br />
RS_modules.pythonize.strings = strings;<br />
})();<br />
async function __main__() {<br />
"use strict";<br />
<br />
var version, display, vector, print, arange, __name__, type, scene, RS_ls, win, Natoms, NA, L, V, gray, mass, Ratom, k, T, dt, animation, s, d, r, boxbottom, boxtop, vert1, vert2, vert3, vert4, Atoms, p, apos, mfp, pavg, x, y, z, theta, phi, px, py, pz, i, deltav, nhisto, histo, xmax, ymax, gg, theory, dv, v, accum, vdist, nc, gmfp, cnts, ti, tiref, gti, t, ncoll, aposold, dist, v_square, vi, P, hitlist, j, ptot, posi, posj, vj, vrel, a, rrel, dx, dy, value, alpha, deltat, mtot, pcmi, pcmj, ij, sigma, gmfp0, MFP, Zco, z0, loc, nsteps;<br />
version = RS_list_decorate([ "3.0", "glowscript" ]);<br />
Array.prototype['+'] = function(r) {return this.concat(r)}<br />
Array.prototype['*'] = function(r) {return __array_times_number(this, r)}<br />
window.__GSlang = "vpython";<br />
display = canvas;<br />
vector = vec;<br />
print = GSprint;<br />
arange = range;<br />
__name__ = "__main__";<br />
type = pytype;<br />
scene = canvas();<br />
var strings = RS_modules.pythonize.strings;<br />
<br />
strings();<br />
"5";<br />
win = 500;<br />
"7";<br />
Natoms = 200;<br />
"8";<br />
NA = 6.022e23;<br />
"10";<br />
L = 1;<br />
"11";<br />
V = GS_power(L, 3);<br />
"12";<br />
gray = (await color.gray(.7));<br />
"13";<br />
mass = .004["/"](NA);<br />
"14";<br />
Ratom = .03;<br />
"15";<br />
k = 1.38064852e-23;<br />
"16";<br />
T = float(input("Temperature of the simulation(K): "));<br />
"17";<br />
dt = 1e-5;<br />
"19";<br />
animation = RS_interpolate_kwargs.call(this, canvas, [RS_desugar_kwargs({width: win, height: win, align: "left"})]);<br />
"20";<br />
animation.range = L;<br />
"21";<br />
animation.title = "Simulation with "["+"](str(Natoms))["+"](" \"hard-sphere\" gas particles at T = ")["+"](str(T))["+"](" K");<br />
"32";<br />
s = "<b>Theoretical and simulated speed distributions.<\/b>\n Initially all atoms have the same speed, \n but collisions change the speeds of the \n colliding atoms. \n \n One of the atoms is marked and leaves a \n trail so you can follow its path.\n \n While speeds are realistic, the collision\n frequency and mean free path are not\n";<br />
"33";<br />
animation.caption = s;<br />
"35";<br />
d = L["/"](2)["+"](Ratom);<br />
"36";<br />
r = .005;<br />
"37";<br />
boxbottom = RS_interpolate_kwargs.call(this, curve, [RS_desugar_kwargs({color: gray, radius: r})]);<br />
"38";<br />
boxbottom.append(RS_list_decorate([ vector(d["-u"](), d["-u"](), d["-u"]()), vector(d["-u"](), d["-u"](), d), vector(d, d["-u"](), d), vector(d, d["-u"](), d["-u"]()), vector(d["-u"](), d["-u"](), d["-u"]()) ]));<br />
"39";<br />
boxtop = RS_interpolate_kwargs.call(this, curve, [RS_desugar_kwargs({color: gray, radius: r})]);<br />
"40";<br />
boxtop.append(RS_list_decorate([ vector(d["-u"](), d, d["-u"]()), vector(d["-u"](), d, d), vector(d, d, d), vector(d, d, d["-u"]()), vector(d["-u"](), d, d["-u"]()) ]));<br />
"41";<br />
vert1 = RS_interpolate_kwargs.call(this, curve, [RS_desugar_kwargs({color: gray, radius: r})]);<br />
"42";<br />
vert2 = RS_interpolate_kwargs.call(this, curve, [RS_desugar_kwargs({color: gray, radius: r})]);<br />
"43";<br />
vert3 = RS_interpolate_kwargs.call(this, curve, [RS_desugar_kwargs({color: gray, radius: r})]);<br />
"44";<br />
vert4 = RS_interpolate_kwargs.call(this, curve, [RS_desugar_kwargs({color: gray, radius: r})]);<br />
"45";<br />
vert1.append(RS_list_decorate([ vector(d["-u"](), d["-u"](), d["-u"]()), vector(d["-u"](), d, d["-u"]()) ]));<br />
"46";<br />
vert2.append(RS_list_decorate([ vector(d["-u"](), d["-u"](), d), vector(d["-u"](), d, d) ]));<br />
"47";<br />
vert3.append(RS_list_decorate([ vector(d, d["-u"](), d), vector(d, d, d) ]));<br />
"48";<br />
vert4.append(RS_list_decorate([ vector(d, d["-u"](), d["-u"]()), vector(d, d, d["-u"]()) ]));<br />
"50";<br />
Atoms = RS_list_decorate([]);<br />
"51";<br />
p = RS_list_decorate([]);<br />
"52";<br />
apos = RS_list_decorate([]);<br />
"53";<br />
mfp = RS_list_decorate([ 0 ])["*"](Natoms);<br />
"54";<br />
pavg = sqrt(2["*"](mass)["*"](1.5)["*"](k)["*"](T));<br />
"56";<br />
for (var RS_Index1 = 0; RS_Index1["<"](Natoms); RS_Index1++) {<br />
i = RS_Index1;<br />
"57";<br />
x = L["*"](random())["-"](L["/"](2));<br />
"58";<br />
y = L["*"](random())["-"](L["/"](2));<br />
"59";<br />
z = L["*"](random())["-"](L["/"](2));<br />
"60";<br />
if ((i === 0 || typeof i === "object" && RS_equals(i, 0))) {<br />
"61";<br />
Atoms.append(RS_interpolate_kwargs.call(this, sphere, [RS_desugar_kwargs({pos: vector(x, y, z), radius: Ratom, color: color.cyan, make_trail: true, retain: 100, trail_radius: .3["*"](Ratom)})]));<br />
"62";<br />
} else {<br />
Atoms.append(RS_interpolate_kwargs.call(this, sphere, [RS_desugar_kwargs({pos: vector(x, y, z), radius: Ratom, color: gray})]));<br />
}<br />
"63";<br />
apos.append(vec(x, y, z));<br />
"64";<br />
theta = pi["*"](random());<br />
"65";<br />
phi = 2["*"](pi)["*"](random());<br />
"66";<br />
px = pavg["*"](sin(theta))["*"](cos(phi));<br />
"67";<br />
py = pavg["*"](sin(theta))["*"](sin(phi));<br />
"68";<br />
pz = pavg["*"](cos(theta));<br />
"69";<br />
p.append(vector(px, py, pz));<br />
}<br />
"71";<br />
deltav = 100;<br />
"73";<br />
async function barx(v) {<br />
"74";<br />
return int(v["/"](deltav));<br />
};<br />
if (!barx.__argnames__) Object.defineProperties(barx, {<br />
__argnames__ : {value: ["v"]},<br />
__module__ : {value: null}<br />
});<br />
<br />
"77";<br />
nhisto = int(5e3["/"](deltav));<br />
"79";<br />
histo = RS_list_decorate([]);<br />
"80";<br />
for (var RS_Index2 = 0; RS_Index2["<"](nhisto); RS_Index2++) {<br />
i = RS_Index2;<br />
histo.append(0);<br />
}<br />
"81";<br />
histo[RS_bound_index((await barx(pavg["/"](mass))), histo)] = Natoms;<br />
"83";<br />
xmax = sqrt(2["*"](k)["*"](T)["/"](mass))["+"](2.5["*"](sqrt(k["*"](T)["/"](mass))));<br />
"84";<br />
ymax = Natoms["*"](4)["*"](GS_power((mass["/"](2["*"](pi)["*"](k)["*"](T))), .5))["/"](exp(1))["*"](deltav)["*"](1.2);<br />
"85";<br />
gg = RS_interpolate_kwargs.call(this, graph, [RS_desugar_kwargs({width: win, height: .4["*"](win), xmax: xmax, align: "left", xtitle: "speed, m/s", ytitle: "Number of atoms", ymax: ymax})]);<br />
"88";<br />
theory = RS_interpolate_kwargs.call(this, gcurve, [RS_desugar_kwargs({color: color.blue})]);<br />
"89";<br />
dv = 10;<br />
"90";<br />
var RS_Iter3 = range(0, 5001["+"](dv), dv);<br />
RS_Iter3 = ((typeof RS_Iter3[Symbol.iterator] === "function") ? (RS_Iter3 instanceof Map ? RS_Iter3.keys() : RS_Iter3) : Object.keys(RS_Iter3));<br />
for (var RS_Index3 of RS_Iter3) {<br />
v = RS_Index3;<br />
"91";<br />
(await theory.plot(v, deltav["/"](dv)["*"](Natoms)["*"](4)["*"](pi)["*"](GS_power((mass["/"](2["*"](pi)["*"](k)["*"](T))), 1.5))["*"](exp(.5["-u"]()["*"](mass)["*"](GS_power(v, 2))["/"](k["*"](T))))["*"](GS_power(v, 2))["*"](dv)));<br />
}<br />
"93";<br />
accum = RS_list_decorate([]);<br />
"95";<br />
var RS_Iter4 = range(int(5e3["/"](deltav)));<br />
RS_Iter4 = ((typeof RS_Iter4[Symbol.iterator] === "function") ? (RS_Iter4 instanceof Map ? RS_Iter4.keys() : RS_Iter4) : Object.keys(RS_Iter4));<br />
for (var RS_Index4 of RS_Iter4) {<br />
i = RS_Index4;<br />
accum.append(RS_list_decorate([ deltav["*"](i["+"](.5)), 0 ]));<br />
}<br />
"96";<br />
vdist = RS_interpolate_kwargs.call(this, gvbars, [RS_desugar_kwargs({color: color.red, delta: deltav})]);<br />
"98";<br />
async function interchange(v1, v2) {<br />
var RS_ls, barx1, barx2;<br />
"99";<br />
barx1 = (await barx(v1));<br />
"100";<br />
barx2 = (await barx(v2));<br />
"101";<br />
if ((barx1 === barx2 || typeof barx1 === "object" && RS_equals(barx1, barx2))) {<br />
return;<br />
}<br />
"102";<br />
if (barx1[">="](len(histo)) || barx2[">="](len(histo))) {<br />
return;<br />
}<br />
"103";<br />
histo[(typeof barx1 === "number" && barx1["<"](0)) ? histo.length["+"](barx1) : barx1]=histo[(typeof barx1 === "number" && barx1["<"](0)) ? histo.length["+"](barx1) : barx1]["-"](1);<br />
"104";<br />
histo[(typeof barx2 === "number" && barx2["<"](0)) ? histo.length["+"](barx2) : barx2]=histo[(typeof barx2 === "number" && barx2["<"](0)) ? histo.length["+"](barx2) : barx2]["+"](1);<br />
};<br />
if (!interchange.__argnames__) Object.defineProperties(interchange, {<br />
__argnames__ : {value: ["v1", "v2"]},<br />
__module__ : {value: null}<br />
});<br />
<br />
"106";<br />
async function checkCollisions() {<br />
var RS_ls, hitlist, r2, ai, aj, dr, j, i;<br />
"107";<br />
hitlist = RS_list_decorate([]);<br />
"108";<br />
r2 = 2["*"](Ratom);<br />
"109";<br />
r2=r2["*"](r2);<br />
"111";<br />
for (var RS_Index5 = 0; RS_Index5["<"](Natoms); RS_Index5++) {<br />
i = RS_Index5;<br />
"112";<br />
ai = apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i];<br />
"113";<br />
for (var RS_Index6 = 0; RS_Index6["<"](i); RS_Index6++) {<br />
j = RS_Index6;<br />
"114";<br />
aj = apos[(typeof j === "number" && j["<"](0)) ? apos.length["+"](j) : j];<br />
"115";<br />
dr = ai["-"](aj);<br />
"117";<br />
if (mag2(dr)["<"](r2)) {<br />
hitlist.append(RS_list_decorate([ i, j ]));<br />
}<br />
}<br />
}<br />
"119";<br />
return hitlist;<br />
};<br />
if (!checkCollisions.__module__) Object.defineProperties(checkCollisions, {<br />
__module__ : {value: null}<br />
});<br />
<br />
"121";<br />
nhisto = 0;<br />
"122";<br />
z = 0;<br />
"123";<br />
nc = 0;<br />
"124";<br />
gmfp = RS_list_decorate([ 0 ])["*"](Natoms);<br />
"125";<br />
cnts = RS_list_decorate([ 0 ])["*"](Natoms);<br />
"126";<br />
ti = RS_list_decorate([ 0 ])["*"](Natoms);<br />
"127";<br />
tiref = RS_list_decorate([ 0 ])["*"](Natoms);<br />
"128";<br />
gti = RS_list_decorate([ 0 ])["*"](Natoms);<br />
"129";<br />
t = 0;<br />
"130";<br />
ncoll = 0;<br />
"132";<br />
for (var RS_Index7 = 0; RS_Index7["<"](5e4); RS_Index7++) {<br />
nsteps = RS_Index7;<br />
"133";<br />
(await rate(300));<br />
"135";<br />
var RS_Iter8 = range(len(accum));<br />
RS_Iter8 = ((typeof RS_Iter8[Symbol.iterator] === "function") ? (RS_Iter8 instanceof Map ? RS_Iter8.keys() : RS_Iter8) : Object.keys(RS_Iter8));<br />
for (var RS_Index8 of RS_Iter8) {<br />
i = RS_Index8;<br />
accum[(typeof i === "number" && i["<"](0)) ? accum.length["+"](i) : i][1] = nhisto["*"](accum[(typeof i === "number" && i["<"](0)) ? accum.length["+"](i) : i][1])["+"](histo[(typeof i === "number" && i["<"](0)) ? histo.length["+"](i) : i])["/"](nhisto["+"](1));<br />
}<br />
"136";<br />
if (RS_equals(nhisto["%"](10), 0)) {<br />
"137";<br />
vdist.data = accum;<br />
}<br />
"138";<br />
nhisto=nhisto["+"](1);<br />
"141";<br />
t = t["+"](dt);<br />
"142";<br />
for (var RS_Index9 = 0; RS_Index9["<"](Natoms); RS_Index9++) {<br />
i = RS_Index9;<br />
"143";<br />
aposold = apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i];<br />
"144";<br />
Atoms[(typeof i === "number" && i["<"](0)) ? Atoms.length["+"](i) : i].pos = apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i] = apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i]["+"](p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i]["/"](mass)["*"](dt));<br />
"145";<br />
ti[(typeof i === "number" && i["<"](0)) ? ti.length["+"](i) : i] = t["-"](tiref[(typeof i === "number" && i["<"](0)) ? tiref.length["+"](i) : i]);<br />
"147";<br />
dist = GS_power((GS_power((apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i].x["-"](aposold.x)), 2)["+"](GS_power((apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i].y["-"](aposold.y)), 2))["+"](GS_power((apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i].z["-"](aposold.z)), 2))), .5);<br />
"148";<br />
mfp[(typeof i === "number" && i["<"](0)) ? mfp.length["+"](i) : i] = mfp[(typeof i === "number" && i["<"](0)) ? mfp.length["+"](i) : i]["+"](dist);<br />
}<br />
"150";<br />
v_square = 0;<br />
"151";<br />
for (var RS_Index10 = 0; RS_Index10["<"](Natoms); RS_Index10++) {<br />
i = RS_Index10;<br />
"152";<br />
vi = p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i]["/"](mass);<br />
"153";<br />
v_square = v_square["+"](vi.x["*"](vi.x))["+"](vi.y["*"](vi.y))["+"](vi.z["*"](vi.z));<br />
}<br />
"155";<br />
P = mass["*"](v_square)["/"](3)["/"](V);<br />
"159";<br />
hitlist = (await checkCollisions());<br />
"164";<br />
var RS_Iter11 = hitlist;<br />
RS_Iter11 = ((typeof RS_Iter11[Symbol.iterator] === "function") ? (RS_Iter11 instanceof Map ? RS_Iter11.keys() : RS_Iter11) : Object.keys(RS_Iter11));<br />
for (var RS_Index11 of RS_Iter11) {<br />
ij = RS_Index11;<br />
"165";<br />
i = ij[0];<br />
"166";<br />
j = ij[1];<br />
"167";<br />
gmfp[(typeof i === "number" && i["<"](0)) ? gmfp.length["+"](i) : i]=gmfp[(typeof i === "number" && i["<"](0)) ? gmfp.length["+"](i) : i]["+"](mfp[(typeof i === "number" && i["<"](0)) ? mfp.length["+"](i) : i]);<br />
"168";<br />
gmfp[(typeof j === "number" && j["<"](0)) ? gmfp.length["+"](j) : j]=gmfp[(typeof j === "number" && j["<"](0)) ? gmfp.length["+"](j) : j]["+"](mfp[(typeof j === "number" && j["<"](0)) ? mfp.length["+"](j) : j]);<br />
"169";<br />
gti[(typeof i === "number" && i["<"](0)) ? gti.length["+"](i) : i]=gti[(typeof i === "number" && i["<"](0)) ? gti.length["+"](i) : i]["+"](ti[(typeof i === "number" && i["<"](0)) ? ti.length["+"](i) : i]);<br />
"170";<br />
gti[(typeof j === "number" && j["<"](0)) ? gti.length["+"](j) : j]=gti[(typeof j === "number" && j["<"](0)) ? gti.length["+"](j) : j]["+"](ti[(typeof j === "number" && j["<"](0)) ? ti.length["+"](j) : j]);<br />
"171";<br />
cnts[(typeof i === "number" && i["<"](0)) ? cnts.length["+"](i) : i]=cnts[(typeof i === "number" && i["<"](0)) ? cnts.length["+"](i) : i]["+"](1);<br />
"172";<br />
cnts[(typeof j === "number" && j["<"](0)) ? cnts.length["+"](j) : j]=cnts[(typeof j === "number" && j["<"](0)) ? cnts.length["+"](j) : j]["+"](1);<br />
"173";<br />
mfp[(typeof i === "number" && i["<"](0)) ? mfp.length["+"](i) : i] = 0;<br />
"174";<br />
mfp[(typeof j === "number" && j["<"](0)) ? mfp.length["+"](j) : j] = 0;<br />
"175";<br />
tiref[(typeof i === "number" && i["<"](0)) ? tiref.length["+"](i) : i] = t;<br />
"176";<br />
tiref[(typeof j === "number" && j["<"](0)) ? tiref.length["+"](j) : j] = t;<br />
"177";<br />
ptot = p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i]["+"](p[(typeof j === "number" && j["<"](0)) ? p.length["+"](j) : j]);<br />
"178";<br />
posi = apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i];<br />
"179";<br />
posj = apos[(typeof j === "number" && j["<"](0)) ? apos.length["+"](j) : j];<br />
"180";<br />
vi = p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i]["/"](mass);<br />
"181";<br />
vj = p[(typeof j === "number" && j["<"](0)) ? p.length["+"](j) : j]["/"](mass);<br />
"182";<br />
vrel = vj["-"](vi);<br />
"183";<br />
a = vrel.mag2;<br />
"184";<br />
if ((a === 0 || typeof a === "object" && RS_equals(a, 0))) {<br />
continue;<br />
}<br />
"185";<br />
rrel = posi["-"](posj);<br />
"186";<br />
if (rrel.mag["<"](Ratom)) {<br />
"188";<br />
continue;<br />
}<br />
"191";<br />
dx = dot(rrel, vrel.hat);<br />
"192";<br />
dy = cross(rrel, vrel.hat).mag;<br />
"195";<br />
value = dy["/"](2["*"](Ratom));<br />
"196";<br />
if (value[">"](1)) {<br />
value = 1;<br />
}<br />
"197";<br />
if (value["<"](1["-u"]())) {<br />
value = 1["-u"]();<br />
}<br />
"198";<br />
alpha = asin(value);<br />
"199";<br />
d = 2["*"](Ratom)["*"](cos(alpha))["-"](dx);<br />
"200";<br />
deltat = d["/"](vrel.mag);<br />
"204";<br />
posi = posi["-"](vi["*"](deltat));<br />
"205";<br />
posj = posj["-"](vj["*"](deltat));<br />
"207";<br />
mtot = 2["*"](mass);<br />
"208";<br />
pcmi = p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i]["-"](ptot["*"](mass)["/"](mtot));<br />
"209";<br />
pcmj = p[(typeof j === "number" && j["<"](0)) ? p.length["+"](j) : j]["-"](ptot["*"](mass)["/"](mtot));<br />
"210";<br />
rrel = norm(rrel);<br />
"211";<br />
pcmi = pcmi["-"](2["*"]((await pcmi.dot(rrel)))["*"](rrel));<br />
"212";<br />
pcmj = pcmj["-"](2["*"]((await pcmj.dot(rrel)))["*"](rrel));<br />
"213";<br />
p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i] = pcmi["+"](ptot["*"](mass)["/"](mtot));<br />
"214";<br />
p[(typeof j === "number" && j["<"](0)) ? p.length["+"](j) : j] = pcmj["+"](ptot["*"](mass)["/"](mtot));<br />
"215";<br />
apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i] = posi["+"](p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i]["/"](mass)["*"](deltat));<br />
"216";<br />
apos[(typeof j === "number" && j["<"](0)) ? apos.length["+"](j) : j] = posj["+"](p[(typeof j === "number" && j["<"](0)) ? p.length["+"](j) : j]["/"](mass)["*"](deltat));<br />
"217";<br />
(await interchange(vi.mag, p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].mag["/"](mass)));<br />
"218";<br />
(await interchange(vj.mag, p[(typeof j === "number" && j["<"](0)) ? p.length["+"](j) : j].mag["/"](mass)));<br />
}<br />
"219";<br />
ncoll=ncoll["+"](len(hitlist));<br />
"220";<br />
sigma = pi["*"](GS_power((2["*"](Ratom)), 2));<br />
"221";<br />
gmfp0 = k["*"](T)["/"](P)["/"](sigma["*"](GS_power(2, .5)));<br />
"222";<br />
MFP = 0;<br />
"223";<br />
Zco = 0;<br />
"224";<br />
for (var RS_Index12 = 0; RS_Index12["<"](Natoms); RS_Index12++) {<br />
i = RS_Index12;<br />
"225";<br />
if (cnts[(typeof i === "number" && i["<"](0)) ? cnts.length["+"](i) : i][">"](0)) {<br />
"226";<br />
MFP=MFP["+"](gmfp[(typeof i === "number" && i["<"](0)) ? gmfp.length["+"](i) : i]["/"](cnts[(typeof i === "number" && i["<"](0)) ? cnts.length["+"](i) : i]));<br />
"227";<br />
Zco=Zco["+"](cnts[(typeof i === "number" && i["<"](0)) ? cnts.length["+"](i) : i]["/"](gti[(typeof i === "number" && i["<"](0)) ? gti.length["+"](i) : i]));<br />
}<br />
}<br />
"229";<br />
z0 = Natoms["/"](V)["*"](sigma)["*"](sqrt(2))["*"](sqrt(8["*"](k)["*"](T)["/"](pi)["/"](mass)));<br />
"230";<br />
if ((nsteps === 0 || typeof nsteps === "object" && RS_equals(nsteps, 0))) {<br />
print(" Time Coll. freq. Mean free path");<br />
}<br />
"231";<br />
if (RS_equals(nsteps["%"](500), 0)) {<br />
print(t, ncoll["/"](t)["/"](Natoms)["*"](2), z0, MFP["/"](Natoms), gmfp0);<br />
}<br />
"233";<br />
for (var RS_Index13 = 0; RS_Index13["<"](Natoms); RS_Index13++) {<br />
i = RS_Index13;<br />
"234";<br />
loc = apos[(typeof i === "number" && i["<"](0)) ? apos.length["+"](i) : i];<br />
"235";<br />
if (abs(loc.x)[">"](L["/"](2))) {<br />
"236";<br />
if (loc.x["<"](0)) {<br />
p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].x = abs(p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].x);<br />
} else {<br />
p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].x = abs(p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].x)["-u"]();<br />
}<br />
"239";<br />
}<br />
if (abs(loc.y)[">"](L["/"](2))) {<br />
"240";<br />
if (loc.y["<"](0)) {<br />
p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].y = abs(p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].y);<br />
} else {<br />
p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].y = abs(p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].y)["-u"]();<br />
}<br />
"243";<br />
}<br />
if (abs(loc.z)[">"](L["/"](2))) {<br />
"244";<br />
if (loc.z["<"](0)) {<br />
p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].z = abs(p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].z);<br />
} else {<br />
p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].z = abs(p[(typeof i === "number" && i["<"](0)) ? p.length["+"](i) : i].z)["-u"]();<br />
}<br />
}<br />
}<br />
}<br />
};<br />
if (!__main__.__module__) Object.defineProperties(__main__, {<br />
__module__ : {value: null}<br />
});<br />
<br />
;$(function(){ window.__context = { glowscript_container: $("#glowscript").removeAttr("id") }; __main__() })})()<br />
// END JAVASCRIPT<br />
<br />
//--><!]]></script><br />
</div></div>
Wikiadmin