ESP Logo
 Elliott Sound Products Satellites and Subwoofers 

Satellites And Subwoofers

© 2004 - Robert C. White and Rod Elliott (ESP)
Page Created 26 August 2004
(Updated 17 Jan 2005)

HomeMain Index articlesArticles Index

Contents
1.0   Introduction

In recent years the satellite sub woofer type of loudspeaker system has become popular, this largely due to the advent of surround sound home theatre.  For this and also for music the sat+sub system has many things to recommend it, the use of one subwoofer is made possible because the range bellow 100Hz is, as indicated by most testing, non directional, thus allowing the subwoofer to be put in the best place for bass, and allowing the various other smaller speakers to be arranged in a way most suited to imaging and surround effect.  Another benefit is that the higher cut off frequency of the satellites potentially keeps the non linear distortion caused by large low bass cone excursions at a low level.

The low cost and compact size is however bought at a cost, this is that the sound pressure generated by such a system is limited by the use of a small driver down to a typical 80Hz crossover frequency.  This is exacerbated by the current fashion to denigrate the reflex enclosure, with people advised to block up those despised ports, and be transported to the sunlit uplands of sonic bliss thereby.

It is generally true that the transient response of vented enclosures is not as good as that from a (properly designed) sealed box, but we live in a world of compromise.  There is little point having a system with perfect transient response if the distortion at realistic listening levels is intolerable.  It is also pointless to try to convince someone who has a small listening area that they need large 3-way systems - especially if there is no physical room for such enclosures.  So, in this area of compromise, there has to be a solution that will provide a reasonable SPL with minimal excursion created distortion.

While sealed enclosures certainly have their place in the scheme of things, any attempt to use them in a 'full range' 2-way system will almost certainly cause excessive distortion at even a moderate SPL, whereas a solution is available that lowers distortion, raises the SPL limit and if carefully done will give very good performance.

Interestingly, any information about the alignment discussed in this article is very scarce on the Net, to the extent that it is almost non-existent.  While the author has recommended the method described for some time (mainly in forum sites), it is an alignment that the editor has also contemplated, but (until now) did not have the information to perform a theoretical design - any attempts had to be empirical.  That this is time consuming in the extreme is obvious, and it is probable that the Quasi-Butterworth, 5th order (QB5) alignment will see an increase in popularity in the future.

What follows is a discussion of the issues and methods for dealing with the most vexing problems, but first some discussion about the linear excursion limit is in order.


2.0   Description

Over part of its cone excursion the loudspeaker driver is to a good approximation a linear transducer, if its B×L (magnetic field strength × conductor length) product is constant, then the force generated by the voice coil is linearly proportional to the current in the voice coil, i.e. ...

F = ( B × L ) × I

Olsen, in his 1962 paper [1], identified two non linear elements that produce harmonic distortion as a function of cone displacement.  Typically as the voice coil moves out of the magnetic gap, the B×L product is no longer constant but starts to fall as the excursion increases, this can be described mathematically as ...

F(x) = a × x + bx²   and   F(x) = a × x + bx³

Where, F(x) is an expression relevant for force when B×L is falling, a & b are constants relevant to a particular driver, and x is the displacement.  The first of these expressions is relevant for the voice coil moving from its rest position asymmetrically, and the second when the motion is symmetrical about the centre position.

When excited by a sine wave the Fourier transform of these yields ...

x = A × cos ωt + B× cos 2ωt and x = A ×cos ωt + B× cos 3ωt

Where ω represents the angular frequency (equal to 2πf), x is the amplitude, ωt is the exciting frequency and A ,B are proportional constants, dependant upon signal level.

The second terms, B × cos 2ωt and B × cos 3ωt represent the second and third harmonics respectively.

For small displacements when BL is constant, the second term is vanishingly small, and the distortion produced comes from non linearities in the suspension and surround, and from air trapped in the magnetic gap.  The 2nd harmonic component is largely out of our control since its level is largely due to the design and construction of the driver.

The third harmonic is in our control however because it is directly related to the amount of signal we put in.  This where the linear excursion comes in....

The peak linear excursion limit is derived, unless it is specified by the manufacturer, by using the rule of thumb that it is the voice coil length, minus the magnetic gap height over two, this is true for the usual overhung voice coil.  Gander, [2] measured a wide selection of drivers and came up with the findings that calculating the linear excursion limit in this way is a very reliable indicator of the driver producing 3% third harmonic distortion, when driven to this limit.  In what follows this is taken as the maximum peak distortion that is acceptable.

As an example of what this means in practical terms, let us take a typical good quality satellite mid/woofer, the Vifa P17WJ, and put it in a sealed box.

Nominal impedance [ohm]8Air gap height [mm]6
Voice coil resistance [ohm]5.8Voice coil inductance [mH]0.55
Nominal power [W]40Eff. diaphragm Area [cm²]136
Short term max power [W]350Moving mass [g]14
Long term max power [W]150Magnet weight [g]415
Operating power [W]6.3Force factor [Bl]6.5
Sensitivity [dB]88VAS [litres]34.7
Frequency range [Hz]37-5000Qms1.55
Free air resonance [Hz]37Qes0.45
Voice coil diameter [mm]32Qts0.35
Voice coil height [mm]14
P17WJ-00-08 Parameters

In a box of 11.5 litres it will have a Q of 0.707 and an fs around 78Hz by the manufacturers specifications.  it has a linear cone excursion of around 4mm peak.  Using Small's expression for the acoustic power we can expect in an average room [3]

Pa = kp × Vd2 × f34

Where kp is the power output constant, it has a value of around 0.8 for a sealed box

Vd is the volume displacement (cone area × peak excursion) and f3 is the system's lower -3dB frequency.  This gives...

kp = 0.8 × (0.0136 × 0.004)2 × 784 = 0.088 Watts
Converting this to decibels we have...
107.5 - 10 × log10 (1 / 0.088) = 96.9dB

This is the peak sound level we can reproduce before exceeding the excursion limit.  Taking into account that "normal" program is said to have a 20dB dynamic range then our average level is around 77dB.  This falls far short of the standard 90dB with 20dB headroom requirement for a hi fi system.  Although it can play a lot louder, this will not be without producing large amounts of harsh sounding non linear distortion - the amount of distortion increases very rapidly beyond our excursion limit.

For instance the 110dB peak output causes a 4mm excursion at 189Hz.  At the 78Hz region, the driver needs an excursion of 16mm.

From the above it seems that the position is hopeless, crossing small satellites of this sort over to a single subwoofer causes large amounts of distortion.

There are of course several alternatives, these include three way satellites, stereo subwoofers crossing over at a higher frequency and so on, but all of these lack the basic simplicity and low cost of the two way satellite plus single subwoofer system.

Luckily there is a solution, this is in the form of the QB5 reflex alignment.

I know that immediately I say that word 'reflex', many people will throw their hands up in horror, make the sign of the cross, or other gestures that are best not considered.  Swimming against the tide of fashionable anti-reflex sentiment is difficult, but I nevertheless present the following for your consideration ...


2.1   The Basic Principles

Two people named Kreutz and Panzer derived the QB5 reflex alignments, and published their work in 1993, [4] the class II group of these have the characteristic that they optimise power handling by means of reducing the cone excursion needed for a given acoustic power output, and these are the most useful for our purpose.

In essence, what we will do is create a box that is specifically designed to have a peaked response at the low end.  This means that the box tuning (rather than the loudspeaker cone) is producing the larger proportion of the LF energy, so less power is needed to achieve a given SPL.  In order to prevent the system sounding like a boom-box, a 2nd order filter (whose Q is also obtained from the table) is placed before the power amp, tuned to the frequency obtained from the charts below.  The system as a whole will now have a vastly higher power output at the low frequency end, without the problem of excursion limiting and subsequent distortion.  On the down side (and like all things, there are compromises), the box will be larger than a conventionally tuned system, but this would be the case regardless if a larger woofer were to be used to try to maximise the excursion limited SPL.

The details of these derivations are as they say 'beyond the scope of this article', and for those interested the paper cited in the references is available from the AES.  Suffice to say however, that the benefits of filter assisted reflex boxes have been known and exploited for many years, but the derivation of these particular alignments was not achieved formally until Panzer came up with a solution to some very difficult mathematics.  Luckily for us this gave us three alignment tables, these are reproduced below.

[Reproduced from:- "Derivation of the Quasi-Butterworth 5 Alignments" , by J.J.M. Kreutz & J. Panzer, Journal of the Audio engineering Society, Vol. 42, No. 5, May 1994]

QtVas/Vb F3/Fskp Fb/Fs1/Qa T60Fa/Fs
Driver QtsBox Volume-3dB Freq.Power Ratio Box TuningFilter DampingSettling TimeFilter Freq.

The three tables stem from the value of a parameter used in the solutions of the equation set.  In practical terms, Group I is the optimum power handling vs. bass extension for drivers in the Qt range specified.  The Group II alignments are for maximum power handling, again within the Qt range covered.  Group III are included for completeness, as they represent the edge of usefulness and can be largely ignored.


2.2   Group I Alignments

QtVas/Vb F3/Fskp Fb/Fs1/Qa T60Fa/Fs
0.3241.9890.9353.911 0.9950.4674.7780.944
0.3182.1031.0006.436 1.000.5184.3151.000
0.3112.2391.0596.782 1.0080.5494.0741.055
0.3032.3961.1147.237 1.0180.5693.9331.110
0.2952.5741.1697.811 1.0300.5823.8441.164
0.2872.7731.2238.512 1.0460.5923.7851.219
0.2792.9911.2779.349 1.0630.5983.7461.274
0.2713.2291.33210.331 1.0830.6033.7171.329
0.2633.4861.38711.469 1.1040.6063.6971.385
0.2553.7621.44312.776 1.1270.6083.6821.441
0.2474.0561.49914.265 1.1520.6103.6711.498
0.2404.3681.55615.951 1.1780.6123.6631.555
0.2334.6981.61417.849 1.2050.6133.6561.613
0.2265.0461.67219.977 1.2340.6143.6511.671
0.2195.4121.73122.352 1.2630.6143.6471.730
0.2135.7951.79024.992 1.2930.6153.6431.789
0.2076.1961.84927.919 1.3240.6163.6411.848
0.2016.6141.90931.153 1.3560.6163.6391.908
0.1967.0501.96934.716 1.3890.6163.6371.968
0.1907.5032.03038.631 1.4220.6163.6352.029
0.1857.9732.09042.921 1.4560.6173.6342.090

2.3   Group II Alignments

QtVas/Vb F3/Fskp Fb/Fs1/Qa T60Fa/Fs
0.4450.5501.0006.436 1.0001.4144.3161.000
0.4250.7511.15210.152 1.1151.5273.9191.107
0.4150.8561.23112.587 1.1741.5583.8391.187
0.4050.9601.30315.176 1.2281.5763.7901.264
0.3941.0721.37318.090 1.2811.5883.7561.339
0.3841.1931.44421.426 1.3351.5953.7291.415
0.3731.3251.51525.264 1.3891.6013.7091.491
0.3621.4671.58829.678 1.4441.6043.6931.567
0.3521.6201.66134.742 1.5001.6073.6811.643
0.3421.7841.73640.528 1.5571.6093.6711.720
0.3311.9591.81147.115 1.6151.6113.6631.797
0.3222.1441.88754.580 1.6741.6123.6571.875
0.3122.3391.96363.009 1.7331.6133.6521.953
0.3032.5452.04072.486 1.7931.6143.6482.031
0.2942.7612.11883.102 1.8541.6153.6452.110
0.2862.9882.19694.952 1.9161.6153.6422.188
0.2783.2242.274108.134 1.9781.6163.6402.268
0.2703.4702.353122.749 2.041.6163.6382.347
0.2623.7262.432138.903 2.1031.6163.6362.426
0.2553.9932.511156.706 2.1661.6163.6352.506
0.2494.2692.590176.216 2.2291.6173.6342.586

2.4   Group III Alignments

QtVas/Vb F3/Fskp Fb/Fs1/Qa T60Fa/Fs
0.5140.5331.0294.976 1.0182.0203.8330.794
0.5170.5271.0305.136 1.0192.0083.8520.813
0.5200.5201.0315.295 1.021.9963.8710.832
0.5230.5131.0325.451 1.021.9873.8920.850
0.5260.5051.0325.603 1.0201.9783.9140.867
0.5300.4961.0325.749 1.0201.9703.9390.884
0.5340.4871.0315.887 1.0201.9633.9660.901
0.5380.4761.0306.016 1.0191.9573.9970.917
0.5430.4651.0286.133 1.0171.9514.0310.932
0.5490.4521.0246.234 1.0151.9474.0700.947
0.5550.4381.0206.319 1.0131.9424.1160.961
0.5620.4221.0156.383 1.0091.9384.1710.975
0.5690.4051.0086.423 1.0051.9354.2360.988
0.5770.3861.0006.434 1.0001.9324.3171.000
0.5870.3640.9906.251 0.9941.9294.4171.011
0.5980.3390.9776.035 0.9861.9274.5451.022
0.6100.3120.9625.780 0.9771.9244.7131.031
0.6250.2800.9435.482 0.9661.9224.9391.039
0.6430.2440.9215.132 0.9521.9215.2541.045
0.6640.2020.8944.719 0.9361.9195.7161.049
0.6910.1500.8614.224 0.9151.9186.4491.050

*   T60 is the time needed for a pulse input to settle by 60dB


2.5   Examples and Charts

Plugging our Vifa P17WJ driver into the relevant class II alignment, i.e. that for a Qt of 0.35, yields (after calculation) the following ...

Box volume21.4 litres
f361.5Hz
Box tuning frequency55.5 Hz
Fa (Aux. filter 12dB/Octave high pass)60.8Hz
Qa0.622
kp34.74

Note that in the above the power handling constant kp has gone up from 0.8 to an impressive 34.74, the practical effect of this is shown below...

Pa = 34.74 × (0.0136 × 0.004)2 × 61.54 = 1.47 Watts
dB SPL = 10 × log10(1.47) + 107.5 = 109.17dB

Our peak output before we reach our excursion limit has gone up by 12.3dB even with a lower f3, if we use the high pass filter from our subwoofer amplifier this should be set at the 61Hz frequency as a minimum, setting it higher results in an even greater increase in power handling.

Two of these speakers in a room can easily meet the 110dB SPL peak level that is the standard for domestic hi fi systems without exceeding the 4mm peak cone excursion limit, this results in a significant reduction in non linear distortion.  One must be careful not to exceed the long term power handling of the drivers - since there will be little loudspeaker distortion at even high levels, there will be no warning that the average electrical power is too high, so some common sense is required to prevent voicecoil damage.

Another fully worked example (similar speaker, but different parameters) will help understanding ...

Nominal impedance [ohm]6Air gap height [mm] -
Voice coil resistance [ohm]4.0Voice coil inductance [mH]0.4
Nominal power [W]50Eff. diaphragm Area [cm²]140
Short term max power [W]  -Moving mass [g]  -
Long term max power [W]  -Magnet weight [g]  -
Operating power [W]  -Force factor [Bl]4.1
Sensitivity [dB]87VAS [litres]49.4
Frequency range [Hz]  -Qms1.34
Free air resonance [Hz]35Qes0.58
Voice coil diameter [mm]  -Qts0.40
Voice coil height [mm]  -
P17WG-00-08 Parameters

The above data were obtained from WinISD, and some parameters are missing, or may be in error.  However, this was used for the simulations that follow, and the general idea is the same for any speaker.

DriverVifa P17WG-00-08
Qt0.40
fs35 Hz
Vas / Vb0.960Vb = 49.4 / 0.96Vb = 51.45 litres
f3 / fs1.303f3 = 35 × 1.303f3 = 45.6 Hz
fb / fs1.228fb = 35 × 1.228fb = 42.98 Hz
1 / Qa1.576Qa = 1 / 1.576Qa = 0.635
fa / fs1.264fa = 35 × 1.264fa = 44.24 Hz

The box in this case is a little over the normal volume for a critically aligned vented enclosure, and will have a response as shown below.  The normal reaction to this would be "Yuck!" because of the large peak of 2.45dB occurring just before rolloff.  This is where the electrical filter comes into play, and this filter should have a frequency of 44.24 Hz (fa) and a Q of 0.635 (Qa) to match this alignment.

fig 1
Figure 1 - Speaker Response

The above graph is the output from WinISD, with the box size and tuning frequencies adjusted according to the calculations above.  The actual f3 point of the speaker is 37.8Hz, but this will be raised by the electrical filter to the figure of 45.6Hz as obtained from the table.  The simplest filter to use for this is the equal component value Sallen-Key, and the schematic for this is shown below.  As noted above, the exact Q will not make a great deal of difference - it may look bad, but will be completely overshadowed by room nodes.  For the time being, we'll look at the system response with the correct Q, especially since it is quite easy to do.

fig 2
Figure 2 - Speaker (Red) and Combined (Green) Response

I used the Simetrix simulator to produce the above, and although it is not a perfect match for the loudspeaker response in Fig. 1, the result is very close.  The red trace is the uncorrected simulated loudspeaker response, while the green trace shows the result after correction with the auxiliary filter.  The simplest way to make the auxiliary filter is to use the 'equal component value Sallen-Key' filter topology.  This allows independent selection of frequency and Q, but the gain of the circuit will change.  The gain obtained from a low Q filter is small, and is unlikely to cause any problems in typical systems.

fig 3
Figure 3 - Auxiliary Correction Filter Circuit

Referring to the schematic above, you can see that R1=R2 and C1=C2.  The second half of the opamp would be used for the other channel.  Note that this circuit must be fed from a low impedance source (typically less than 1k).  The values needed are given by ...

f = 1 / (2π × R × C)

This is not affected by the gain setting (which alters the Q), so it is a flexible solution with no bad habits.  For the filter response we calculated above, first we shall select a convenient value for C - say 100nF.  The value for R is obtained by ...

R = 1 / (2π × C × f ) = 1 / ( 2π × 100E-9 × 44.24 ) = 35.98kΩ

It has to be said that the calculated resistance value is not going to be easy to get, but a 50k dual-gang pot will allow you to set the frequency wherever you like.  Alternatively, you could use a 33k and a 2.7k in series, or even just 33k with a still insignificant error.  Now that the frequency is decided, we need to calculate the Q.  This is not even a little bit painful, since the tables actually give the figure for damping (1/Qa), which simplifies things even more.

Q = 1 / (3 - K) or R3 = (2 - d) × R4

where K is the circuit gain (determined by the ratio of R3:R4), and d = 1/Q.  For minimum DC offset, the parallel combination of R3 and R4 should equal R1, but this is of no consequence if the power amplifier has a capacitor input as offset will be very low.  A value of 10k will be quite acceptable for our purposes.  From the table, we determined that d (1/Qa) is 1.576, so substituting ...

R3 = (2 - 1.576) × 10k = 0.424 × 10k = 4.24k

Alas, we have another non-standard value, but a 3.9k resistor in series with 330Ω gives 4.23k which will be quite acceptable.  In fact, you could just use the 3.9k resistor, since the error is very small indeed, and will be completely un-noticed in practice.  The values calculated and shown on the schematic were used in the simulation, and as is obvious, the error is small.

Using the same formula as before, we can calculate the acoustic power increase over a sealed box.  5mm linear excursion is assumed, from the specifications for the P17WG ...

Pa = kp × Vd2 × f34 Pa = 0.8 × (0.0136 × 0.005)2 × 624 = 34mW

and for the QB5 alignment we calculated ...

Pa = 15.176 × (0.0136 × 0.005)2 × 404 = 179mW

Not quite as impressive as the previous example, but note that we have a -3dB frequency of 40Hz instead of 62Hz, and still have a 7.2dB increase in effective acoustic power output.  Converting to dB SPL, we get ...

dB SPL = 10 × log10(0.179) + 107.5 = 100.03 dB

The is a substantial improvement over the 92.8dB SPL obtainable from a sealed box, plus we have over ½ octave more low frequency extension.  As before, if we were to raise the crossover frequency, there will be an equivalent increase in level before the excursion limit is reached.  As it stands, the arrangement just calculated would outperform the majority of systems using conventional vented boxes and the same driver.  The enclosure is only marginally larger than a 'critically aligned' system, but has an f3 of 40Hz rather than 34Hz.  Another plus is that there is a filter that will help reduce infrasonic excursions that produce no useful output, but can introduce considerable distortion because the cone is unloaded below the cutoff frequency.

As previously mentioned, the 2nd order high pass output from your sub amplifier can be used as the filter.  While it is true that this is usually a Butterworth type, with a Q of around 0.7, overall this does not make much difference at this low crossover frequency, as room modes and cancellation/reinforcement from walls predominate over any exact crossover characteristic.  Be aware that many 'plate' sub amps have a fixed high-pass filter, and only the low-pass section is variable.  Such units will almost certainly need the filter circuit described.

If the high-pass section of the crossover is variable this is also useful, as this allows an optimum setting for a smooth response to be found more easily.  It also increases power handling even more.  Referring to the original driver (the P17WJ), if kp remains the same and f3 is increased to 80Hz, the excursion limited peak output goes up to 113.7dB - this represents an increase of 16.8dB in peak output over a sealed box before excursion limit is reached.  One can use a peak input of 263 Watts in the case of the P17WJ, which has an unspecified 'short term' maximum power input rating of 350 Watts.

The method outlined here allows the maximum performance to be extracted from the single subwoofer plus two way satellite configuration, and in terms of 'bang for your buck' is hard to beat.

It does however use the now demonised reflex enclosure, but it must be emphasised that any negative effects from a port are more than compensated for by the very significant reduction in non linear distortion this technique confers.

In addition, it could be argued that the electrical filter brings the box back into alignment, and not only relieves the loudspeaker of undue stress and distortion, but by removing the low frequency component from the amplifier, the power actually needed could be reduced by perhaps 3dB or more.  While this may not sound like much, it is the difference between a 100W amp and a 50W amp, so the savings all round could be very significant.


Appendix

Curve fit solutions for the first two alignments are presented here.  These may be used to calculate the values, rather than using the lookup tables.

Class I Alignments

All of these have a correlation coefficient of at least 0.98, except for 1/Qa, with 0.91

Vas / Vb = 0.126 × (Qt^-2.4731)
F3 / Fs = -1.2576 - 1.9761 × ln(Qt)
kp = 221.419 - 1.401.7 × Qt + 2283.34 × Qt^2
1 / Qa = 0.4437 + 10.2652 × Qt^2 - 29.982 × Qt^3
Fb / Fs = 3.252 - 13.448 × Qt + 20.0398 × Qt^2
Fa / Fs = 4.4786 - 15.94 × Qt + 15.7596 × Qt^2

Class II Alignments

1/Qa has a correlation coefficient of 0.871

Vas / Vb = 16.883 - 69.941 × Qt + 75.2755 × Qt^2
F3 / Fs = 5.6721 - 15.168 × Qt + 10.6884 × Qt^2
kp = 9357.22 × e^(-15.987 × Qt)
1 / Qa = 1.3873 + 7.3126 × Qt^2 - 15.584 × Qt^3
Fb / Fs = 4.7413 - 12.556 × Qt + 9.4877 × Qt^2
Fa/Fs = -1.173 - 2.694 × ln(Qt)

References

[1]H.F. Olsen"Analysis of the effects of non linear elements upon the performance of back enclosed direct radiator loudspeaker mechanisms"   AES journal, Vol 10, pp.156-162, (Apr 1962)
[2]M.R. Gander"Moving coil loudspeaker topology as an indicator of linear excursion capability"   AES journal, Vol 29, No ½, (Jan./Feb 1981).
[3]R.H. Small"Vented box loudspeaker systems, Part II: Large signal analysis"   AES journal, Vol 21, No 6 (Jul/Aug 1973).
[4]J.J.M. Kreutz & J. Panzer"Derivation of the quasi Butterworth 5 alignments"   AES journal, Vol. 42, No 5 (May 1994)

 

HomeMain Index articlesArticles Index
Copyright Notice. This article, including but not limited to all text and diagrams, is the intellectual property of the reference authors quoted (in particular the AES), Robert C. White and Rod Elliott, and is © 2004. Reproduction or re-publication by any means whatsoever, whether electronic, mechanical or electro-mechanical, is strictly prohibited under International Copyright laws. The authors (Robert C White and Rod Elliott) grant the reader the right to use this information for personal use only, and further allow that one (1) copy may be made for reference. Commercial use is prohibited without express written authorisation from Robert C. White and Rod Elliott.
Page created and copyright © 26 August 2004./ Updated 17 Jan 05 - corrected errors in formulae.