SpectraSim Documentation
Parameter · 07 Run Settings

Highest frequency

acoustics.frequency.f_max Unit: Hz Step 07 · Run Settings

The highest frequency of the acoustic sweep. On auto it is the frequency at which the smallest inlet or outlet begins to carry its first non-plane mode: below it every port carries plane waves only.

Last updated

Theory

Besides the plane wave, a duct of radius aa carries transverse modes that propagate only above their cut-on frequency. For a circular duct with rigid walls the first one is the (1,0) mode, with kta=1.8412k_t a = 1.8412, the first zero of J1′J_1', equation (1). The auto value uses the smallest port, its equivalent radius A/π\sqrt{A/\pi} and the speed of sound of the medium: that of the Material step's gas for the uniform medium, or the solved temperature at that port for the medium from the CFD fields (Acoustic medium).

The computation remains valid above that frequency: the ports are modal, every mode whose cut-on is below fmaxf_{max} is included, plus Evanescent modes per port further modes, and the mode table on this step lists them. Two things change with fmaxf_{max}:

  • the mesh: the largest element is c/(fmaxNλ)c/(f_{max}N_\lambda), equation (2) (Elements per wavelength), so the number of unknowns grows roughly with fmax3f_{max}^3 where the wavelength, not the curvature, sets the size;
  • the insertion loss: above an outlet's first cut-on, the plane-wave radiation terminations are incomplete and the insertion loss is flagged there, unless the outlet uses the radiation-domain termination.

Equations

fmax=1.8412 c2πa,a=Aminπf_{max} = \frac{1.8412\,c}{2\pi a}, \qquad a = \sqrt{\frac{A_{min}}{\pi}}
(1)
hmax=cfmaxNλh_{max} = \frac{c}{f_{max} N_\lambda}
(2)
wheremeaningunit
fmaxf_{max}Highest frequency (auto)Hz
ccSpeed of sound (uniform medium, or at the smallest port)m/s
aaEquivalent radius of the smallest inlet or outletm
AminA_{min}Area of the smallest inlet or outlet (from the CAD)m²
hmaxh_{max}Largest element size of the acoustic meshm
NλN_\lambdaElements per wavelength-

Default and rule

Default
auto: first non-plane cut-on of the smallest portComputed by the programMeasured
Rule

Auto (equation 1). It requires the exported geometry (the port areas); until then the field reads "pending". For octave, 1/3-octave and custom spacing the field is locked and shows the upper edge of the last band or the last listed frequency.

Source

The cut-on of the (1,0) mode of a circular duct is classical duct acoustics. The program checks its numerical port modes against the zeros of Jm′J_m', computed numerically (tolerance 0.5 %; measured 0.01-0.17 % on the ports of case V2, a two-inlet single-chamber muffler with internal inlet and outlet tubes).

Cost measured on case V2 (cold air, 6 elements per wavelength, one worker, with an earlier sparse LU solver; the current solver is faster):

  • 1.0 kHz: 39.8k unknowns, 0.9 s per frequency;
  • 1.8 kHz: 144.6k unknowns, 11.6 s per frequency;
  • 2.6 kHz: 381k unknowns, 141 s per frequency.

Valid range > f_min; set by the band range or the list for octave, 1/3-octave and custom spacing

Values

Above fminf_{min}. Above the auto value, port modes beyond the plane wave are included automatically.

Cost

Doubling fmaxf_{max} multiplies the number of unknowns by up to 8 on a wavelength-limited mesh, and the solve time per frequency by more.