Bergh-Tijdeman Solution: Extended Kutin-Svete Formulation

For calculating the frequency response of a line-cavity system with N tubes and cavities in series.

% Implement the Bergh-Tijdeman Calculations for N tubes and N volumes. This

% version DOES take into account the improvements made by Kutin and Svete,

% and extends Kutin and Svete's formulation to be compatible with

% multi-tubed systems.

% ===================== RELEVANT PUBLICATIONS =============================

% Bergh, H., and Tijdeman, H. Theoretical and Experimental Results for the

% Dynamic Response of Pressure Measuring Systems.

% https://vegvesen.brage.unit.no/vegvesen-xmlui/bitstream/handle/11250/

% 193925/NLR%20TR%20F.238.pdf?sequence=1.

%

% Kutin, J., and Svete, A. “On the Theory of the Frequency Response of Gas

% and Liquid Pressure Measurement Systems with Connecting Tubes.”

% Measurement science & technology, Vol. 29, No. 12, 2018, p. 125108.

% https://doi.org/10.1088/1361-6501/aae884.

% =========================================================================

% ======================== INPUTS and OUTPUTS =============================

% INPUTS:

% - Tubes: 2D Matrix of the lengths and Radii of the "thin tubes"

% portion of the Bergh-Tijdeman calculations. Expecting to be in the

% form of [2xN] although dimensions will be checked.

% Example: Tubes = [1xN Tube Radii; ...

% 1xN Tube Lengths];

%

% - Vols: 2D Matrix of the lengths and radii of the sensing volumes

% prescribed by Bergh-Tijdeman. Must be the same length as "Tubes" so

% fill any tube discontinuities with a Radius and Length of 0.

% Example: Vols = [1xN Vol Radii; ...

% 1xN Vol Lengths];

%

% - f: Frequency in Hz. Can either be 1x2 vector to indicate lower and

% upper limits, or a custom 1D vector of frequencies

%

% - med: The gas being used for the tube/cavity system. Input as string

% to get Room temperature properties for select gases, or as

% structure with necessary gas properties.

% >>>>> SEE MORE IN "Gas Properties Import" INSTRUCTIONS BELOW <<<<<<

%

% - doplots: Decide whether or not generate plots of the calculated

% pressure ratios.

%

% OUTPUTS:

% - G(iω): The Complex Frequency Response Function (Pout/Pin). The

% magnitude can be calculated from abs(G) and the phase lag (aka the

% argument of the complex frequency) can be calculated using angle(G)

% ========================================================================

% ======================== REQUIRED FILES =================================

% - besselzero.m

% - For calculating the Nth zero of the Bessel Functions (used in the

% calculation of nv)

% ======================== OPTIONAL FILES =================================

% - AirProperties_kPaK.m

% - If wanting to automatically obtain all gas properties for Air

% using just the temperature and pressure

% - AirPropTable2.m

% - For obtaining air properties below 0 degrees Celcius using just

% temperature and pressure

% - DrawTubeVolSystem.m

% - For generating a scaled drawing of the described geometry

% =========================================================================