From 94067499f2552297bbcddfbd155728f30d585b85 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Sun, 11 Jul 2021 06:24:27 +0200 Subject: [PATCH] Includes calculation of the Clausing Factor from an analytic expression (obtained from the Etienne Staub Master thesis) and an approximation. --- .../@MOTSimulator/calculateClausingFactor.m | 9 +++++++++ 1 file changed, 9 insertions(+) create mode 100644 MOT Capture Process Simulation/@MOTSimulator/calculateClausingFactor.m diff --git a/MOT Capture Process Simulation/@MOTSimulator/calculateClausingFactor.m b/MOT Capture Process Simulation/@MOTSimulator/calculateClausingFactor.m new file mode 100644 index 0000000..0c48d56 --- /dev/null +++ b/MOT Capture Process Simulation/@MOTSimulator/calculateClausingFactor.m @@ -0,0 +1,9 @@ +function ret = calculateClausingFactor(this) + + ClausingFactorApproximation = (8 * this.NozzleRadius) / (3 * this.NozzleLength); + + alpha = 0.5 - (3*this.Beta^2)^-1 * ((1 - (2*this.Beta^3) + ((2*this.Beta^2) - 1) * sqrt(1+this.Beta^2)) / (sqrt(1+this.Beta^2) - (this.Beta^2 * asinh((this.Beta^2)^-1)))); + ClausingFactorAnalytic = 1 + (2/3) * (1 - (2 * alpha)) * (this.Beta - sqrt(1 - this.Beta^2)) + (2/3) * (1 + alpha) * this.Beta^(-2) * (1 - sqrt(1 + this.Beta^2)); + + ret = ClausingFactorAnalytic; +end \ No newline at end of file