Removed the multiplication of the full Clausing Factor since the correct flux coming out of the oven is obtained by accounting for the clipping by the aperture which is done by numerically integrating the angular distribution till the exit divergence angle that is limited by the aperture to obtain what is called the "reduced" Clausing Factor. This is calculated in the initialVelocitySampling and multiplied to this flux to get the "reduced" flux.

This commit is contained in:
Karthik 2021-07-11 06:28:50 +02:00
parent 94067499f2
commit 04b395b6ad

View File

@ -3,6 +3,8 @@ function ret = calculateFreeMolecularRegimeFlux(this)
%See Expected atomic flux section in Barbiero
Dy164VapourPressure = 133.322*exp(11.4103-2.3785e+04./(-219.4821+this.OvenTemperatureinKelvin)).*100; % Vapor Pressure Dysprosium for the given oven temperature
Dy164DensityinOven = Dy164VapourPressure/(Helper.PhysicsConstants.BoltzmannConstant*this.OvenTemperatureinKelvin);
ClausingFactorApproximation = (8 * this.NozzleRadius) / (3 * this.NozzleLength);
ret = Helper.PhysicsConstants.Dy164IsotopicAbundance * 1/4 * Dy164DensityinOven * this.AverageVelocity * ClausingFactorApproximation * pi * this.NozzleRadius.^2;
end
ret = Helper.PhysicsConstants.Dy164IsotopicAbundance * 1/4 * Dy164DensityinOven * this.AverageVelocity * pi * this.NozzleRadius.^2;
% Needs to be multiplied with the "Clausing Factor" which here would be
% the probability not for the full solid angle but the angle subtended
% by the aperture of the oven at its mouth.
end