Acoustical applications of modulation beamforming technology

Abstract

A method for generating a spatial distribution of acoustic transmitting frequencies by a time modulated array (TMA) of transducers includes generating, by a tapped delay line, a plurality of pulsed sampling signals, wherein each pulsed sampling signal includes a series of frequency harmonics and successive signals of the plurality of pulsed sampling signals are separated by a predetermined delay time; mixing each of the plurality of pulsed sampling signals with a time-limited information signal wherein a plurality of mixer output signals is generated; bandpass filtering each of the plurality of mixer output signals; generating a first plurality of simultaneous TMA beams from the plurality of filtered and weighted output signals by driving a plurality of acoustic transducers in a spatial array of acoustic transducers, wherein each beam of the first plurality of simultaneous beams is associated with one of a plurality of transmitting carrier frequencies

Claims

1. A method for generating a spatial distribution of acoustic transmitting frequencies by a time modulated array (TMA) of transducers that produce simultaneous multiple beam steering in a non-linear propagation medium, comprising the steps of: generating, by a tapped delay line, a plurality of pulsed sampling signals, wherein each pulsed sampling signal includes a series of frequency harmonics and successive signals of the plurality of pulsed sampling signals are separated by a predetermined delay time; mixing each of the plurality of pulsed sampling signals with a time-limited information signal wherein a plurality of mixer output signals is generated; bandpass filtering each of the plurality of mixer output signals; generating a first plurality of simultaneous TMA beams from the plurality of filtered and weighted output signals by driving a plurality of acoustic transducers in a spatial array of acoustic transducers, wherein each beam of the first plurality of simultaneous beams is associated with one of a plurality of transmitting carrier frequencies, wherein a bearing-target doppler ambiguity, a uniform response in bearing, out-of-band spectra leakage interference, peak power level, and cross-beam reverberation are controlled in the first plurality of simultaneous beams, and a secondary plurality of beams is formed that includes a plurality of different primary carrier frequencies that intersect the first plurality of simultaneous TMA beams in a far-field, wherein a parametric sonar array is produced; receiving multiple echo frequencies of the first and second plurality of TMA simultaneous beams at a spectrum analyzer; and, processing the multiple echo frequencies of the first and second plurality of TMA simultaneous beams with the spectrum analyzer, wherein each spectral component of the multiple echo frequencies corresponds to a unique spatial bearing.

2. The method of claim 1, further comprising producing a spectral component with a phase slope that corresponds to a particular beam steering angle by time-modulating each array element in a spatial array.

3. The method of claim 1, where the spatial beamformer comprises two sub-beamformers, wherein one or both are TMA beamformers, whose output intersects in a specific region in the far-field with one or more common maximum response axes (MRAs).

4. The method of claim 1, wherein the spatial array is an N-point element linear array with constant element spacing d and a time-varying far-field pressure for a monochromatic signal is given as: $p\left(,t\right)={.Math.}_p{a}_p{.Math.}_m{W}_m\left\{\exp \left(j2m\left[\frac{\left(f_r+f_0+{\mathrm{pf}}_s\right)d\sin \left(\right)}{c}-\frac{p}{N}\right]\right)\exp \left(j2\left(f_r+f_0+{\mathrm{pf}}_s\right)t\right)\right\},$ wherein is a pointing angle, which is a direction of maximum response for a particular beam associated with a harmonic number, a.sub.p are weighting coefficients for each harmonic p, W.sub.m is a complex weighting for each element m, f.sub.r is a reference frequency, f.sub.s is a fundamental harmonic frequency, f.sub.o is a transmit frequency, p ranges from a first harmonic P.sub.o to a last harmonic P.sub.c, c is the speed of sound, and is +1 or 1 depending on whether a sum or difference frequency band is retained.

5. The method of claim 4, further comprising choosing a total delay time t equal to a sampling period, N=1/f.sub.s, wherein =d sin()/c, and the pointing angle is a function of the harmonic parameters (p) and (f.sub.s), the center frequency f.sub.c=f.sub.r+f.sub.o, and the array parameters (N) and (d).

6. The method of claim 5, wherein a sampling frequency f.sub.s is at least: f.sub.s=f.sub.D+6/T, wherein f.sub.D is the maximum doppler and T is the signal pulse length.

7. The method of claim 5, wherein out-of-band spectra leakage interference is controlled by setting 2f.sub.o(P.sub.cP.sub.o)f.sub.s, wherein (f.sub.s) is a sampling frequency that includes spreading due to source, target and receiver speeds.

8. The method of claim 7, wherein the pointing angle p=sin.sup.1[2p/(Nk.sub.pd)], wherein k.sub.p is a harmonic dependent wave-number that equals 2(f.sub.r+f.sub.o+pf.sub.s)/c.

9. The method of claim 4, wherein, when the information signal is s(t) has a complex Fourier spectrum S(f), the time-varying far-field pressure is expressed by: $p\left(,t\right)=\underset{p}{.Math.}{a}_p{\mathrm{df}}^{}S\left(f^{}\right)\left\{\underset{m}{.Math.}{W}_m\left\{\exp \left(j2m\left[\frac{\left(f^{}+f_r+f_0+{\mathrm{pf}}_s\right)d\sin \left(\right)}{c}-p/N\right]\right)\right\}\exp \left(j2\left(f^{}+f_r+f_0+{\mathrm{pf}}_s\right)t\right)\right\}.$

10. The method of claim 4, wherein the total peak power required for M tonals is (2M.sup.2P.sub.rms), wherein P.sub.rms is an rms power per tonal required for a specified source level.

11. The method of claim 4, wherein a transmit source level of the TMA SL.sub.TMA for the linear array with a constant input electrical power is equal to transmit source level of the linear array SL.sub.LIN reduced by the number of frequency coded beams according to
SL.sub.TMA=SL.sub.LIN10 log M, wherein SL.sub.TMA=20 log rms pressure on the MRA.

12. The method of claim 1, wherein the spatial array is an N-point element circular array with constant angular element spacing , and a time-varying far-field pressure for an arbitrary signal with complex Fourier spectrum S(f) is given as: $p\left(,t\right)=\underset{p}{.Math.}{a}_p{\mathrm{df}}^{}S\left(f^{}\right)\left\{\underset{m}{.Math.}{W}_m\exp \left(-j2\left[\left(f_p\left(R/c\right)\left(\cos \left(-m\right)-\cos \left(\right)+\mathrm{pm}/N\right)\right)\right]\right)\right\}\exp \left(j2f_{p}t\right),$ wherein is a pointing angle, which is a direction of maximum response for a particular beam associated with a harmonic number, a.sub.p are weighting coefficients for each harmonic p, W.sub.m is a complex weighting for each element m, f.sub.p=(f+f.sub.r+f.sub.o+pf.sub.s) where f.sub.r is a reference frequency, f.sub.s is a fundamental harmonic frequency, f.sub.o is a transmit frequency, p ranges from a first harmonic P.sub.o to a last harmonic P.sub.c, R is the radius of the circular array, c is the speed of sound, and is +1 or 1 depending on whether a sum or difference frequency band is retained.

13. The method of claim 12, further comprising additional phase weightings that enhance spatial filtering, wherein each element is phased back-to-a-line or to a circular arc for maximum signal echo-to-noise ratio.

14. The method of claim 1, further comprising selecting frequencies that decrease the coherence of the spatially interacting beams that transmit major lobes wherein one or more of a cross-beam masking clutter or a reverberation level is reduced.

15. The method of claim 1, wherein the spatial array is a planar array of N.sub.1-point elements with an inter-element d.sub.1-spacing in a first direction and N.sub.2-point elements with an inter-element d.sub.2-spacing in a second direction that is excited by complex signal s(t)=2u(t)cos .sub.otu(t) that represents a low-frequency amplitude modulation u (t) with carrier frequency .sub.o and whose complex Fourier spectrum is S()=U(+.sub.o)+U(.sub.o), wherein a time-varying far-field pressure for an arbitrary signal with complex Fourier spectrum S(f) is given as: $p\left(,,t\right)=\frac{1}{2}\underset{P_1,P_2}{.Math.}{a}_{p_1}{a}_{p_2}d{}^{}U\left({}^{}-{}_1{}_0\right)\left\{\exp \left(j\left(m_1-1\right)\left[\frac{{}_{12}{d}_1\cos \left(\right)}{c}-\frac{2p_1}{N_1}\right]\right)\right\}\left\{\exp \left(j\left(m_2-1\right)\left[\frac{{}_{12}{d}_2\cos \left(\right)}{c}-\frac{2p_2}{N_2}\right]\right)\right\}$ where and are the direction cosines of a far-field observation point with respect to the 1-direction and 2-direction, respectively, .sub.12=.sub.r.sub.1++.sub.2.sub.r.sub.2+p.sub.1.sub.s.sub.1+p.sub.2.sub.s.sub.2 is a frequency of a p.sub.1p.sub.2-harmonic, and c is the speed of sound, wherein the direction cosines of the individual beam MRA's are defined as:
cos =(2p.sub.1/N.sub.1)(c/d.sub.1.sub.12), and
cos =(2p.sub.2/N.sub.2)(c/d.sub.2.sub.12).

16. The method of claim 1, wherein the spatial array further comprises a dual beamformer that generates a set of collinear first difference primary frequencies, wherein the first beamformer includes a set of beams of a multi-element TMA array and a second beamformer produces a wide-beam that encompasses all major TMA transmission lobes for deep depth underwater applications, wherein a source level of the wide beam is equal to a peak source level of the TMA array.

17. The method of claim 16, wherein the set of collinear first difference primary frequencies {f.sub.i} is defined as $\begin{array}{c}\left\{f_i\right\}={.Math.\left\{\mathrm{Wide}\mathrm{Beam}\right\}-\left\{{\mathrm{TMA}}_i\right\}.Math.}_{\mathrm{abs}}\\ =\left\{16\mathrm{kHz}\right\}-(12.17\mathrm{kHz},13.\mathrm{kHz},13.83\mathrm{kHz},14.68\mathrm{kHz}\}\\ =\left\{3.83\mathrm{kHz},3.\mathrm{kHz},2.17\mathrm{kHz},1.34\mathrm{kHz}\right\}.\end{array}$

18. The method of claim 17, wherein the collinear first difference frequencies include transmission spectra whose echoes are input to a spectrum analyzer.

19. The method of claim 17, wherein the collinear first difference frequencies for the parametric array maintain a high efficiency in the demodulation process, wherein the frequency down-ratio is 3-to-10.

20. The method of claim 1, further comprising filtering each of the plurality of mixer output signals wherein one of a sum-frequency sideband or a difference frequency sideband is filtered from each of the plurality of mixer output signals.

21. The method of claim 1, further comprising complex weighting each of the plurality of mixer output signals.

22. The method of claim 1, wherein the bandpass filtering is one of a sum frequency band pass filtering or a difference frequency bandpass filtering.

23. The method of claim 1, wherein the different primary carrier frequencies of the plurality of different primary carrier frequencies are one of lower primary carrier frequencies or higher primary carrier frequencies.

24. An acoustic space-frequency diversity system, comprising: a tapped delay line includes a plurality of taps that generate a plurality of pulsed sampling signals, wherein each tap generates one of the plurality of pulsed sampling signals, each pulsed sampling signal includes a series of frequency harmonics, and successive signals of the plurality of pulsed sampling signals are separated by a predetermined delay time; a plurality of mixers, wherein each mixer of the plurality of mixers combines one of the plurality of pulsed sampling signals with a time limited information signal to generate one of a plurality of mixer output signals; and a spatial array of transducer elements, where each transducer element of the spatial array of transducer elements receives one of the plurality of mixer output signals, wherein the spatial array of transducer elements generates from the plurality of mixer output signals a plurality of simultaneous beams wherein each beam of the plurality of simultaneous beams is associated with one of a plurality of carrier frequencies, wherein simultaneous beams are subject to design constraints that control a bearing-target doppler ambiguity, a uniform response in bearing, out-of-band spectra leakage interference, peak power level, and cross-beam reverberation are controlled in the first plurality of simultaneous beams upon receive.

25. The system of claim 24, wherein spatial array of transducer elements includes 3-rectangular elements that generate simultaneous endosteal and periosteal bone healing beam patterns, wherein each rectangular element measures 25 mm3.4 mm with an inter-element spacing of 3.5 mm.

26. The system of claim 24, further comprising a set of bandpass filters, wherein each bandpass filter of the plurality of bandpass filters receives one of the plurality of mixer output signals and filters one of a sum-frequency side lobe or a difference frequency side lobe from the one of the plurality of mixer output signals.

27. The system of claim 24, wherein: the bearing-target Doppler ambiguity is controlled by a sampling frequency and an array geometry, the uniform response in bearing is controlled by one or more of the array geometry or harmonic separation choices, the out-of-band leakage is controlled by one or more of a bandpass filter or a frequency separation rule, peak power is controlled by a power amplifier design and a number of harmonic beams, and cross-beam reverberation is controlled by one or more of assigning different carrier frequencies to each beam or by spatial beam steering geometry.

28. The method of claim 1, wherein: the bearing-target Doppler ambiguity is controlled by a sampling frequency and an array geometry, the uniform response in bearing is controlled by one or more of the array geometry or harmonic separation choices, the out-of-band leakage is controlled by one or more of a bandpass filter or a frequency separation rule, peak power is controlled by a power amplifier design and a number of harmonic beams, and cross-beam reverberation is controlled by one or more of assigning different carrier frequencies to each beam or by spatial beam steering geometry.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIGS. 1A-B depict TMA principles in reducing cross-beam reverberation for multiple beam transmissions, according to an embodiment of the disclosure.

(2) FIG. 2 depicts an exemplary TMA configuration comprised of a linear array of N-equi-spaced elements, according to an embodiment of the disclosure.

(3) FIG. 3 depicts an exemplary TMA configuration comprised of a circular array of N-equi-spaced elements, according to an embodiment of the disclosure.

(4) FIG. 4 depicts an exemplary TMA configuration comprised of a planar array of N.sub.1N.sub.2-equi-spaced elements, according to an embodiment of the disclosure.

(5) FIG. 5 is a functional block diagram of a receiver for a TMA processor analyzer with omni receive beam output driving a multi-spectral line FFT receiver, according to an embodiment of the disclosure.

(6) FIG. 6 depicts exemplary BGS TMA transmissions, according to embodiments of the disclosure.

(7) FIG. 7 illustrates BGS configurations for: (A) normal operation; and (B) TMA operation, according to embodiments of the disclosure.

(8) FIG. 8 illustrates the application of TMA principles for underwater sonar operation to generate a parametric array, according to an embodiment of the disclosure.

DETAILED DESCRIPTION

(9) Embodiments of the disclosure are directed to a unique acoustic space-frequency diversity system that includes a time modulated array that generates a spatial distribution of transmitting frequencies that permit simultaneous multiple beam steering, where each beam is associated with a different carrier frequency. Each transmit beam can be steered to a different spatial direction in either the vertical or horizontal plane for maximum acoustic coverage and signal-to-masking ratio. The multiple transmitted frequencies can be received with a single omnidirectional piezoelectric hydrophone whose output is processed with a high-resolution fast Fourier transform (FFT) receiver, where each spectral component corresponds to a unique spatial bearing. Thus, by coding the spatial acoustic field in frequency, the receiver architecture complexity is minimized and most of the cross-beam clutter or reverberation masking level in the receiver display output can be significantly reduced.

(10) FIGS. 1A-B depict the basic principles behind the effectiveness of TMA transmissions for reducing the coherence of reverberation in multi-beam echoes, utilizing a common method for scanning called rotational directional transmission (RDT). RDT sequentially transmits a beam in a circular pattern at a constant frequency and narrow bandwidth, producing a coherent summation (i.e., a linear summation on the basis of amplitude) of reverberation in any one of the receive beams due to (1) the product of the transmit and receive major-lobes and (2) sidelobe interference from adjacent beams in a major lobe. This is shown in FIG. 1A, which shows a left beam, a center beam and a right beam of the same frequencies and phases. In FIGS. 1A-B, the vertical axis represents the signal strength or amplitude, the horizontal axis represents range, and the oblique axis represents frequency, where the bracket indicates the signal bandwidth. TMA operation tends to reduce this to an incoherent summation (i.e., a linear summation on the basis of energy), often referred to as spatial decoupling, by employing a separate frequency for each beam and spreading the set of transmit frequencies over a larger bandwidth. In addition, the spatial distribution of the transmitted frequencies will also tend to reduce the spatial variance of the background fluctuations, thereby improving the false alarm rate. This is shown in FIG. 1B, where the left, center and right beams now have different frequencies.

(11) Embodiments of the disclosure refine the frequency encoding procedure in the TMA beamforming process of operational systems to address the issues given above.

(12) First, the geometry of the linear, convex and planar arrays are presented in a series of signal representations that are clearly defined. Conducting these functional operations results in the general expression for the instantaneous spatial-temporal pressure field producing simultaneous beams as a function of frequency.

(13) Second, in practical TMA operation, the spatial distribution of the beams will depend upon the constraints of finite pulse length, finite array bandwidth, element directivity pattern, and target doppler. These, in turn, will impact on the following engineering and operational considerations: 1. Interbeam interference within the transmitting bandwidth; 2. Bearing-doppler ambiguity; and 3. Spectra leakage interference from outside the passband.

(14) Interbeam interference will result in a nonuniform signal response on the maximum response axes (MRA's) of a preformed beam (PFB) receiver. For efficient detection capability at all bearings, the spatial distribution of the maximum response axis should be (almost) constant, i.e., the separation frequency (f.sub.s) should be sufficiently large so that the sidelobes of each sinc function do not seriously overlap the adjacent sinc function, as expressed in EQ. (6), below.

(15) Bearing-doppler ambiguity is a consideration in those applications characterized by target and/or platform motion. Since normal TMA operation codes the spatial distribution of beams in frequency, if the frequency separation between beams is insufficient, a moving target at a given bearing may produce a doppler shifted echo frequency that lies within the frequency band assigned to an adjacent beam. Although the detection capability of sonar will not be meaningfully affected, the target will appear on the display at an erroneous bearing. This is referred to as a bearing-doppler ambiguity and must be eliminated to permit accurate target localization and tracking. Thus, in practice, the harmonic separation will have to be much greater than that indicated above to avoid bearing-doppler ambiguity due to the combined effects of target and sonar platform motion.

(16) For a moving source, target, and receiver, the frequency of the echo pulse (f.sub.e) is given by (Horton),
f.sub.e=f.sub.p[(cv.sub.t)(cv.sub.r)]/[(cv.sub.s)(cv.sub.t)]
where f.sub.p the transmitted frequency at the source, (c) is the velocity of the sound wave moving away from the source, (cv.sub.x) is the relative velocity of the sound wave removing source, (cv.sub.t) is the relative velocity of the transmitted pulse at the target, (cv.sub.t) is the relative velocity of the sound wave reflected from the moving target, and (cv.sub.r) is the relative velocity of the sound wave at the moving receiver.

(17) The total maximum doppler (f.sub.D) incurred is determined by computing the range of (f.sub.ef.sub.p) for the spread of source and receiver speeds and opening and closing target speeds. If it is assumed that v.sub.t=v.sub.t and v.sub.r=v.sub.s, (v.sub.t, v.sub.r)/c1, and c=4920 feet/sec, then the doppler shift f.sub.D is
f.sub.D=|f.sub.ef.sub.p|.sub.max=0.7f.sub.p|v.sub.sv.sub.t|.sub.max (in Hz).

(18) Combining previous results, the separation frequency f.sub.s required to eliminate bearing-doppler ambiguity and provide reasonably uniform responses in bearing is
f.sub.s=f.sub.D+6T=0.7f.sub.p|v.sub.sv.sub.t|.sub.max+6/(pulse length).

(19) Spectra (out-of-band) interference leakage can be reduced to insignificant levels (less than 50 dB re MRA) by judicious selection of the separation or sampling frequency f.sub.s and information frequency f.sub.o and the number of beams, namely, that
2f(P.sub.cP.sub.o)f.sub.s
where P.sub.c and P.sub.o are integers that denote the last and first harmonics, respectively, (P.sub.cP.sub.o) is the number of beams, (f.sub.s) is determined from the equation above, which includes the spreading due to source, target and receiver speeds.

(20) It has been shown that interbeam interference and bearing-doppler ambiguity can be eliminated by increasing the sampling sufficiently to cover the total doppler spread, including both platform and target. However, for a given piezoceramic element design, as the sampling frequency increases, the number of allowable frequency diversity beams decreases. Also, the information frequency f.sub.o must increase in accordance with the equation above, to avoid out-of-band leakage interference.

(21) Third, the transmit source level (SL.sub.LIN) on a maximum response axis (MRA) for a linear array with N elements with constant input electrical power will be reduced by the number of frequency coded beams, so that the transmit source level for the TMA is
SL.sub.TMA=SL.sub.LIN10 log M=[7.16+10 log P.sub.eN+10 log E.sub.ff+DI]10 log M
where the reference (for sonar operation) is to 1 microbar at a distance of 1 yard from the face of the transducer element, (P.sub.e) is the electrical power input per element, (E.sub.ff) is the transmitting frequency, (DI) is the directivity index of the transducer array, (M) is the number of beams formed, and each element contributes to the source level simultaneously. In an embodiment, SL=20 log rms pressure on the MRA.

(22) It is important to see how the number of beams impacts on the peak power rating of the power amplifier. For the case of M harmonics of peak amplitude (A), when the tonals are in phase and add coherently, the peak amplitude will be (MA) at the fundamental period (1/f.sub.s). This implies that if the rms power per tonal required to give a specified source level is P.sub.rms, then the total peak power required for (M) tonals is (2M.sup.2P.sub.rms). The design complexity and cost of the power amplifiers will increase with the number of tonals or TMA beams, i.e., whether it is a narrowband or wideband system.

(23) Formulation of Linear TMA

(24) According to an embodiment, simultaneous multiple beam frequency steering can be obtained by producing controlled frequency-varying phase shifts between adjacent transducer elements. A basic TMA acoustic system depicting the technique for a general array, in broad terms, can be considered to be comprised of the signal generator, which is an informational signal, the TMA beamformer, which includes a pulsed sampling signal, a tapped delay line, mixers, bandpass filters, and complex weightings, a power amplifier, considered normalized to 1 watt, and a transducer array.

(25) For a linear array, the formulation of a instantaneous pressure field p(, t) is based on FIG. 2 which depicts an exemplary TMA configuration comprised of a tapped delay line 202 that receives a pulse sampling signal 201, mixers 204 that combine the output of the tapped delay line with an information signal 203, sum frequency bandpass filters 205 and complex weightings 206 that act on the signals output from the mixers 204, and a linear array of N-equi-spaced transducer elements 207 with inter-element spacing d. Each element of this array 207 is fed from a tap on the delay line 202 which is excited by a series of frequency harmonics 201. According to embodiments, the transducer elements 207 map include power amplifiers. The requirement for simultaneous multiple beam frequency steering can be reduced to one of producing controlled frequency-varying phase shifts between adjacent transducer elements. The resulting pattern generated by this array consists of a set of distinct beams of different frequencies, pointing in different directions. The transmitting acoustic frequency and beam steering angle for each beam are uniquely determined by the harmonic number with the zeroth harmonic producing a beam whose maximum response axis (MRA) is perpendicular to the array. The linear TMA is analyzed as follows. According to embodiments, the tapped delay line 202, the mixers 204, bandpass filters 205 and complex weightings 206, along with the pulse sampling signal 201 and the information signal 203 make up the beam forming apparatus.

(26) Referring to FIG. 2, according to an embodiment, the harmonic generation function or pulse sampling signal E.sub.s(t) 201 is a finite Fourier series input to the tapped delay line 202 and is represented as
E.sub.s(t)=.sub.P.sub.0.sup.P.sup.ca.sub.p exp(j2(f.sub.r+pf.sub.s)t)(1)
where j={square root over (1)}, (f.sub.r) is a reference frequency, (f.sub.s) is the fundamental harmonic or separation frequency, and P.sub.o and P.sub.c are integers denoting the first and last harmonics, respectively. The weighting coefficients, a.sub.p, are included for generality. The fundamental frequency is chosen to satisfy the requirements of total doppler spread, including target and platform motion, and the effects of finite signal pulse length. The reference frequency is chosen to facilitate subsequent filtering operations. The total number of harmonics is (P.sub.cP.sub.o+1). If the time is referenced to tap m+1, the sampling signal output by the mth-tap is given by
.sub.mE.sub.s(t)=.sub.P.sub.0.sup.P.sup.ca.sub.p exp(j2(f.sub.r+pf.sub.s)(tm))(2)
where () is the delay between successive taps. Each E.sub.s(t) is then mixed or time-modulated by a mixer 204 with a time-limited information signal, s(t) 203, of the form
s(t)=p.sub.T(t)cos(2f.sub.ot),(3)
where f.sub.o is the transmit frequency p.sub.T(t) is a rectangular envelope of pulse length 2T defined by:

(27) $p_T\left(t\right)=\{\begin{array}{cc}1& \mathrm{if}0t2T\\ 0& \mathrm{otherwise}\end{array},$
whose Fourier transform is

(28) $\begin{array}{cc}=\frac{\sin \left(2\mathrm{fT}\right)}{f}\exp \left(-j2\mathrm{fT}\right).& \left(4\right)\end{array}$

(29) The output of the mth-mixer is

(30) $\begin{array}{cc}{E}_m(t\}=p_T\left(t\right)\cos \left(2f_{o}t\right){}_m{E}_s\left(t\right)=\frac{1}{2}{p}_T\left(t\right){.Math.}_{P_0}^{P_c}{a}_P\left\{\exp (j2\left(\left(f_r+f_0+{\mathrm{pf}}_s\right)t-m\left(f_r+p{f}_s\right)\right)+\exp \left(j2\left(\left(f_r-f_0+{\mathrm{pf}}_s\right)t\right)-m\left(f_r+p{f}_s\right)\right)\right\}.& \left(5\right)\end{array}$
Taking the Fourier transform of E.sub.m(t), the sum and difference spectra are

(31) $\begin{array}{cc}{E}_m\left(f\right)=\underset{P_0}{\stackrel{P_c}{.Math.}}{a}_p\exp \left(-j2\left(m\left(f_r+{\mathrm{pf}}_s\right)-\left(f_r+{\mathrm{pf}}_s\right)T\right)\right)\left\{\frac{\sin \left(2\left[f-\left(f_r+f_0+{\mathrm{pf}}_s\right)\right]T\right)}{\left[f-\left(f_r+f_0+{\mathrm{pf}}_s\right)\right]T}+\frac{\sin \left(2\left[f+\left(f_r+f_0+{\mathrm{pf}}_s\right)\right]T\right)}{\left[f+\left(f_r+f_0+{\mathrm{pf}}_s\right)\right]T}\right\}.& \left(6\right)\end{array}$

(32) It may be seen that the effect of a finite pulse length is to introduce a sin(x)/x (sinc function) amplitude distribution for each harmonic. Considering transmitting in the sum-frequency band, the lobes of the first sin(x)/x term in the curly brackets establish the degree of interbeam interference. The lobes of the second term reflect the spectra leakage interference from outside the passband. By applying the condition, 2f.sub.o>(P.sub.c+P.sub.o)f.sub.s, spectra leakage will be minimal.

(33) After filtering by bandpass filters 205, which filter one of the lobes in the curly bracket term of EQ. (6), the resultant signal waveform may be weighted by the shading coefficient W.sub.m 206 prior to driving the mth-array element and is given by
.sub.m(t)W.sub.mp.sub.T(t).sub.pa.sub.p{exp(j2[(f.sub.r+f.sub.0+pf.sub.s)tm(f.sub.r+pf.sub.s)])},(7)
where the index is +1 or 1 depending on whether the sum or difference frequency band, respectively, is retained. Note that the weighting by the complex shading coefficient W.sub.m 206 is optional and can be eliminated in other embodiments. Exemplary, non-limiting complex shading coefficients have unit magnitude and 0 phase. The approximation is due to the time-limited nature of the information signal. In the monochromatic case, T.fwdarw., and the approximation can be replaced by an equality. EQ. (7) may be further simplified by letting (f.sub.r) be an integral number of cycles, i.e., f.sub.r=q, where q is any integer and assuming T.fwdarw.,
.sub.m(t)=W.sub.m.sub.pa.sub.p{exp(j2[(f.sub.r+f.sub.0+pf.sub.s)tmpf.sub.s])}.(8)

(34) Continuing with FIG. 2, as a result of mixing the pulse sampling signal with the finite sinusoidal pulse at the signal generator output, sum and difference frequencies are produced at the mixer output. When the difference frequencies are filtered by the bandpass filters 205, the output of the bandpass filter 205 is quite similar in structure to the pulse sampling signal, but with several distinct differences. First, the carrier frequency of the zeroth harmonic is (f.sub.r+f.sub.o) instead of (f.sub.r). Second, the signal is attenuated at the leading and trailing edges. This is attributed to employing an ideal rectangular bandpass filter which corresponds to convolving the mixer output with a sin(x)/x impulse response. Only a part of the sinc function overlaps the signal at the leading and trailing edges during the convolution process, thus resulting in an attenuation of the signal at these points.

(35) The attenuation will be less with more realizable filter characteristics. In actuality, the signal does not end abruptly at t=2T, but continues to decay, eventually reaching a steady-state of zero amplitude. With the property of reciprocity between conjugate Fourier domains, a broadband filter will produce a more rapid decay than a narrowband filter. Finally, a primary difference is clearly shown by comparing EQS. (1) and (8): each spectral component has a phase slope defined by the harmonic number. It will be shown below that when .sub.m(t) is applied to a TMA planar array of spatial elements, every harmonic numberand thus, each spectral componentcorresponds to a particular beam steering angle.

(36) If each signal .sub.m(t) is applied to the mth-transducer element of the linear array 207, the time-varying sound pressure at an arbitrary point in the far-field is the vector sum of the signals {.sub.m(t)}, considering relative propagation time delays. The far-field TMA sound pressure p(, t) in the direction , where is the angle of the line to the far-field observation point 208 with respect to an array normal 209, may then be expressed as
p(,t)=.sub.m(tm),(9)
referred to the m=1 element where =d sin()/c and c is the velocity of sound.

(37) Carrying out the summation utilizing EQ. (8), the output of the linear multi-element transducer array 207 is

(38) $\begin{array}{cc}p\left(,t\right){.Math.}_p{a}_p{.Math.}_m{W}_m\left\{\exp \left(j2m\left[\frac{\left(f_r+f_0+{\mathrm{pf}}_s\right)d\sin \left(\right)}{c}-\frac{p}{N}\right]\right)\exp \left(j2\left(f_r+f_0+{\mathrm{pf}}_s\right)t\right)\right\},& \left(10\right)\end{array}$
where it can be seen that each beam in the summation is associated with one of the carrier frequencies.

(39) To improve the visualization and understanding of the TMA process, the above derivation of the beam pattern distribution expressed in EQ. (10) assumes a long signal pulse length to approximate the monochromatic case and negligible spectra leakage.

(40) For an information signal s(t) with complex Fourier spectrum S(f), the far-field time-varying pressure for a monochromatic signal given in equation (10) is weighted by the spectrum to give:

(41) $\begin{array}{cc}p\left(,t\right)=\underset{p}{.Math.}{a}_p{\mathrm{df}}^{}S\left(f^{}\right)\left\{\underset{m}{.Math.}{W}_m\left\{\exp \left(j2m\left[\frac{\left(f^{}+f_r+f_0+{\mathrm{pf}}_s\right)d\sin \left(\right)}{c}-\frac{p}{N}\right]\right)\right\}\exp \left(j2\left(f^{}+f_r+f_0+{\mathrm{pf}}_s\right)t\right)\right\},& \left(11\right)\end{array}$
where the total delay time N is chosen to be equal to the sampling period, N=1/f.sub.s.

(42) In a variation of the linear array according to an embodiment, the TMA beamformer for the linear array of N-equi-spaced transducer elements depicted in FIG. 2 can utilize the nonlinear properties of the acoustic medium by transmitting at higher power, sufficient to produce finite-amplitude waves in the far-field. This is further discussed below in for the section Underwater Mining and Communications.

(43) TMA Beam Steering Constraints

(44) The sum of EQ. (10) may be recognized as the beam pattern for a linear array of frequency (f.sub.r+f.sub.o+pf.sub.s) pointing in the direction
.sub.p=sin.sup.1[2pf.sub.s]/(k.sub.pd)](12)
where .sub.p is a pointing angle, which is a direction of maximum response for a particular beam associated with a harmonic number, (k.sub.p) is the harmonic dependent wave-number,
k.sub.p=2/.sub.p=2(f.sub.r+f.sub.o+pf.sub.s)/c.(13)

(45) By choosing the total delay time equal to the sampling period, N=1/f.sub.s, the pointing angle .sub.p is a function of the harmonic parameters (p) and (f.sub.s), the center frequency (f.sub.r+f.sub.o), and the array parameters (N) and (d); EQ. (12) simplifies to
.sub.p=sin.sup.1[2p/(Nk.sub.pd)].(14)

(46) For a given array geometry and fixed values of (f.sub.r, f.sub.o, f.sub.s), the range of (p) allowable in a TMA design according to an embodiment are those which satisfy the inequality,
[2p/(Nk.sub.pd)].sub.abs1.(15)

(47) According to an embodiment, the number of acceptable beams determines the spatial coverage and the required bandwidth of the transmitting system. Due to practical constraints of transducer bandwidth and the operational requirement that the harmonic beams have uniform beamwidth and sidelobe structure, the number of harmonics is TMA beamforming will be less than that permitted by EQ. (15).

(48) EQ. (15) expresses how the acoustic (or E/M) system parameters define the directions to which the beams are simultaneously steered. As the number of elements or element spacing increases, the MRAs of the harmonic beams will move closer to one another. For a fixed element spacing, increasing the center frequency (zeroth harmonic) will result in a narrower 3 dB beamwidth. This implies that as the separation frequency increases, the beamwidths will become narrower for positive harmonics and wider for negative harmonics, the latter accompanied with steering angles that diverge more widely from the zeroth harmonic. For separation frequencies much smaller than the center frequency, the harmonic wave-number will be practically constant, producing a symmetric distribution. If [2p/(Nk.sub.pd)].sub.abs1, then the beam steering angles will be multiples of [2p/(Nk.sub.pd)].sub.abs. If it is desired to steer a particular harmonic in a given direction, say .sub.p*, it can be done by defining the delay time .sub.p* such that
.sub.p*=k.sub.pd sin .sub.p*/(2pf.sub.s).(16)
Formulation of TMA Convex Array

(49) FIG. 3 depicts an exemplary TMA configuration comprised of a circular array of N-equi-spaced elements, according to an embodiment of the disclosure. In FIG. 3, elements that are similar to those in FIG. 2 have similar reference numbers. FIG. 3 differs from FIG. 2 in that the linear array of transducer elements 207 is replaced by a circular array of transducer elements 307 of radius (R) to the array center 311, where an angular separation of each transducer element is (). According to an embodiment, the complications introduced by arranging the transducer elements in a circular arc are geometrical; all other aspects of the TMA principle are the same as before. In an analysis, all times to the far-field points are referenced to a line through the apex of the circular sector which is perpendicular to the symmetry axis or the line bisecting the sector angle, as indicated in FIG. 3. The TMA pattern for a N-point convex array in the far-field may consequently be written for an arbitrary signal with spectrum S(f) as

(50) 0$\begin{array}{cc}p\left(,t\right)=\underset{p}{.Math.}{a}_p{\mathrm{df}}^{}S\left(f^{}\right)\{\underset{m}{.Math.}{W}_m\left\{\exp \left(-j2\left[\left(f_p\left(R/c\right)\left(\cos \left(-m\right)-\cos \left(\right)+\mathrm{pm}/N\right)\right)\right]\right)\right\}\exp \left(j2f_{p}t\right),& \left(17\right)\end{array}$
where f.sub.P=(f+f.sub.r+f.sub.o+pf.sub.s). The term in the square brackets is the conventional beam pattern steered to a direction defined by the harmonic number of the sampling signal. Note that the direction of each component beam cannot be determined directly, since the quantity in the exponent is not in general zero for all (m) when () corresponds to the MRA of a particular beam. Also, each harmonic represents a beam smeared out slightly due to the finite bandwidth of the signal s(t). As for linear arrays, this, as well as interbeam interference, is minimized for signals having long pulse lengths.

(51) According to an embodiment, additional phase weightings can be introduced to phase each element back-to-a-line or to a circular arc. This will act as a spatial filter to maximize the signal-to-noise ratio. Let this phasing circle 312 pass through the apex, A 313 of the circular transducer array 307, and be perpendicular to the symmetry axis. The phase weighting will be:
exp[j2(f.sub.p/c)[R(1cos )(1cos )],(18)
where () is the angular location of the mth-element with respect to the symmetry axis, () is the radius of the phasing circle, and =sin.sup.1(R/)sin . Note that ()(R), the radius of the circular transducer array, so that the center 310 of the phasing circle is behind the tapped delay line 302. The circular TMA formulation is completely characterized by EQS. (17) and (18).
Formulation of TMA Planar Array

(52) FIG. 4 depicts an exemplary TMA configuration comprised of a planar array of N.sub.1N.sub.2 spaced elements, according to an embodiment of the disclosure. In FIG. 4, elements that are similar to those in FIG. 2 have similar reference numbers. FIG. 4 depicts of a first tapped delay line 402 that receives a first pulse sampling signal 401, a second tapped delay line 412 that receives a second pulse sampling signal 411, mixers 404 that combine the outputs of the first and second tapped delay lines with an information signal 403, and a planar array 407 of N.sub.1 elements with spacing d.sub.1 in a first direction D1, and N.sub.2 elements with spacing d.sub.2 in a second direction D2. For convenience of illustration, the sum frequency bandpass filters and complex weightings are not shown. A transmission process includes two distinct linear modulation excitations, each applied simultaneously along one of the directions of the array, thereby providing a three-dimensional fan of beams coded with different frequencies. The TMA planar array system is briefly analyzed as follows.

(53) Following the method of analyzing the linear array, the pulse sampling signal E.sub.s(t) is again a finite harmonic series input to the tapped delay line, similar to that in EQ. (1),
E.sub.s.sub.1,2(t)=.sub.P.sub.0.sup.P.sup.ca.sub.p.sub.1,2 exp(j2(f.sub.r.sub.1,2+p.sub.1,2f.sub.s.sub.1,2)t),(19)
where the subscripts 1, 2 refer to either the first or second tapped delay line, and the output at the m.sub.1-tap in the D1 direction is;
m.sub.1E.sub.s.sub.1(t)=exp(j2f.sub.r.sub.1t).sub.p.sub.0.sup.p.sup.ca.sub.p.sub.1 exp(j2p.sub.1(f.sub.s.sub.1t(m.sub.11)/N.sub.1)).(20)

(54) For an information signal 303 of the form s(t)=2u(t)cos .sub.ot, where u(t) is the low-frequency amplitude modulation, the complex Fourier spectrum can be expressed as S()=U(+.sub.o)+U(.sub.o). Therefore, the output of the m.sub.1-mixer is

(55) $\begin{array}{cc}{E}_{m_1}\left(t\right)=\frac{1}{2}U\left({}^{}-{}_1{}_0\right)\exp \left(j\left({}^{}+{}_{r_1}\right)t\right)d{}^{}{.Math.}_{P_0}^{P_c}{a}_{p_1}\exp \left(j2p_1\left(f_{s_1}t-\left(m_1-1\right)/N_1\right)\right),& \left(21\right)\end{array}$
where the index, .sub.1, is +1 or 1 depending on whether the sum or difference frequency, respectively, is retained by the bandpass filter. The m.sub.1-mixer output, .sub.m.sub.1(t), is multiplied by the sampling signal at the m.sub.1-tap, which is similar to the excitation in EQ. (21), and is given by:
m.sub.2E.sub.s.sub.2(t)=exp(j2f.sub.r.sub.2t).sub.p.sub.0.sup.p.sup.ca.sub.p.sub.2 exp(j2p.sub.2(f.sub.s.sub.2t(m.sub.21)/N.sub.2)).(22)
The resultant signal is filtered again prior to driving the m.sub.1,2 array element, and is designated E.sub.m.sub.12(t).

(56) If each signal E.sub.m.sub.12(t) is applied to the m.sub.1,2-transducer element of the planar array 407, the time-varying sound pressure established in the far-field is the vector sum of the signals {E.sub.m.sub.12(t)}, considering relative propagation time delays (re elements m.sub.1, m.sub.2=1). The far-field TMA sound pressure for the TMA planar array 407 may consequently be written as

(57) $\begin{array}{cc}p\left(,,t\right)=\frac{1}{2}\underset{P_1,P_2}{.Math.}{a}_{p_1}{a}_{p_2}d{}^{}U\left({}^{}-{}_1{}_0\right)\left\{\exp \left(j\left(m_1-1\right)\left[\frac{{}_{12}{d}_1\cos \left(\right)}{c}-\frac{2p_1}{N_1}\right]\right)\right\}\left\{\exp \left(j\left(m_2-1\right)\left[\frac{{}_{12}{d}_2\cos \left(\right)}{c}-\frac{2p_2}{N_2}\right]\right)\right\}& \left(23\right)\end{array}$
where and are the direction cosines of the far-field observation point with respect to the 1-direction and 2-direction, respectively, and .sub.12=.sub.r.sub.2++.sub.2.sub.r.sub.2+p.sub.1.sub.s.sub.1+p.sub.2.sub.s.sub.2 is the frequency of the pipe-harmonic. The terms in the brackets are the conventional beam patterns, each steered to a direction defined by the harmonic number of the associated sampling signal. The direction cosines of the individual beam MRA's are defined in EQ. (24):
cos =(2p.sub.1/N.sub.1)(c/d.sub.1.sub.12),
cos =(2p.sub.2/N.sub.2)(c/d.sub.2.sub.12).(24)

(58) Each harmonic represents a beam smeared out slightly due to the finite bandwidth of s(t). This, as well as interbeam interference, is minimized for signals having long pulse lengths.

(59) TMA Receivers

(60) The linear, convex, and planar TMAs described above require a special receiver to process the echo returns. This is because space has now been coded into frequency and thus requires a fine-frequency detection and tracking receiver.

(61) FIG. 5 is a functional block diagram of a receiver TMA processor analyzer, implemented as a passive detector with a probability of detection (PD) above 90% and false alarm rate (FAR) below one per day. The TMA echo signatures, which is the RCV Beam Output, are processed by three contiguous modules: an input signal conditioning unit (ISCU) 510, a spectrum analyzer unit (SAU) 520, and a post Processing Unit (PPU) 530. The ISCU 510 is controlled by a system clock signal SYSCLK. The SAU 520 includes a spectrum analyzer, such an FFT processor 521, a band of tuned bandpass filters or a heterodyne correlator, and a magnitude spectral weighting unit 522 that takes the absolute magnitude of the output of the FFT processor 521. An exemplary post processing unit (PPU) 530 includes a cascade of proven processing technologies. The PPU 530 includes averaging logic and data memory for additional algorithmic processing 531 to enhance the SNR, and a unit 532 that performs automatic spectral line detection and tracking of the bearing-frequency, pattern recognition and ambient noise rejection. As shown, the output of the PPU is then sent to a Bone Fracture Healing Assessment Indicator (BFHAI).

(62) Applications of TMA Technology

(63) There are a number of applications of TMA technology that have not been considered to date (2020), such as medical therapeutic ultrasound, acoustic microscopy, and underwater mining, exploration, and communications.

(64) Bone Growth Stimulation (BGS)

(65) The TMA approach can promote Bone Growth Stimulation (BGS) for simultaneous endosteal and periosteal bone fracture healing which would require two different spatial beams with different steering angles as shown in FIG. 6.

(66) FIGS. 7A-B depict exemplary configurations for BGS where both endosteal and periosteal healing are simultaneously stimulated. This would not require an Acoustic Modal Converter (AMC) but an array comprised of three (3) rectangular strips to form a piston transducer. The power transmitted will be compatible with the beam nonuniformity ratio (BNR), spatial average-temporal average intensity (SATA), effective radiation area (ERA), and transmit frequency of the transducer proposed with the AMC.

(67) FIG. 7A is a block diagram for normal operation and uses a signal generator 71 to produce a low-intensity pulsed ultrasound (LIPUS) signal that includes a transmitting signal of 1 MHZ, a constant envelope, a pulse repetition rate of 1 kHz, a pulse width of 200 secs, and a duty cycle of 20%, that is sent to a power amplifier 73 via a signal conditioner 72 that provides an rms output power of 314 mW to a circular piston transducer 74 whose 6 dB bandwidth is about 400 kHz, i.e., from 0.8 to 1.2 MHz, SATA of 50 mW/cm.sup.2, acoustic rms power of 194 mW, efficiency 50%, with an ERA of 3.14 cm.sup.2.

(68) FIG. 7B is a block diagram for TMA system operation for the same signal parameters as in FIG. 7A that includes a 3-section TMA module 76 based on the general configuration shown in FIG. 1. Each TMA section of TMA module 76 in FIG. 5B includes a tapped delay line that receives a 50 KHz pulse sampling signal 77, a double-balanced mixer-suppressed carrier (DBM-SC) that combines the output of the delay line with an information signal 78, a bandpass filter (BPF), and a linear power amplifier (LPA) whose output is about 90-100 mW rms, and a transducer element 79. The delay of each step is 6.66667 secs. The information signal is 1 MHZ with a repetition rate of 1 KHz and a pulse width of 200 s. The information signal includes the 0.sup.th and 4.sup.th harmonics, but the 1.sup.st, 2.sup.nd and 3.sup.rd harmonics have been removed. The output of each delay step drives the RF input of the DBM-SC, where the intermediate frequency (IF) output is band-pass filtered and sent to the LPA. Each of the 3 LPA outputs drives a specific rectangular transducer slat, assuming an efficiency of 50%, that generates two unique acoustic beams simultaneously, one at bearing 0 for endosteal healing and the other at 31.5 for periosteal healing, as shown in FIG. 6. Each beam will have a SATA of about 52 mW/cm.sup.2, an acoustic rms power of about 135 mW, and an ERA 2.58 cm.sup.2. Note that when the TMA is employed, an absorptive AMC is not required to achieve periosteal healing as in FIG. 7A. Each transducer element measures about 25 mm3.4 mm, with a separation d of 3.5 mm.

(69) Acoustic Microscopy

(70) A TMA approach according to an embodiment can also be used for a Scanning Acoustic Microscope (SAM). A SAM device can provide unique information about the absorption, elastic properties, and density of tissue cells, which can improve the visualization and clinical assessment of consistent patterns of cellular malignancy.

(71) Acoustical bio-physics estimates a maximum SAM frequency of about 4-5 GHz, producing a wavelength of 347 nm. However, an issue with operating at this frequency is that the usual aluminum motor-driven spatial-positional-stepping system is fabricated by an aluminum extrusion process that produces a positional error of about 1 m (1 GHz). A TMA according to an embodiment can eliminate this effect by electronically producing the acoustic beams to cover the desired sector.

(72) Underwater Mining and Communications

(73) A TMA approach according to an embodiment can also be used for underwater communications and oceanographic mining and exploration. Accurate spatial communications in several specific bearings requires minimal sidelobe interference while mining and exploration requires special consideration of beams propagating to mineral beds at different ocean depths and bearings correlating to frequency dependent attenuation. For these applications the term bearing refers to either azimuthal or D/E (depth/elevation).

(74) FIG. 8 depicts two transmit beam configurations. The first uses TMA modulation beamforming technology with a 9-element linear array, such as that shown in FIG. 2, where each element is a square 2.292 inches transducer with an inter-element spacing of 2.4 inches and that produces a complex signal of 4 harmonics (1, 0, 1, 2) and simultaneously transmits the four beams covering a 3 dB 45 sector (20 to +25). This can be hull-mounted in the D/E plane on a surface ship. The transmit frequency at the 0.sup.th-harmonic is 13 kHz.

(75) The same 9-element linear array may include a second beamformer that utilizes complex element weightings to produce the wide-beam shown in FIG. 8 that covers the 4 narrow TMA beams and matches their peak source level. If the transmit power is sufficiently high, the nonlinear properties of the acoustic medium can be exploited to generate a set of first difference frequencies defined by the frequency of the wide-or-omnidirectional beam, such as f.sub.omni, minus each of the TMA frequencies f.sub.1, f.sub.0, f.sub.1, f.sub.2 shown.

(76) According to an embodiment, for the parametric array in FIG. 8 to maintain a high efficiency in the demodulation process, the frequency down-ratio should be about 3-to-10. In an exemplary embodiment, f.sub.omni will be about 15-19.0 kHz. Thus, very low frequency propagating signals, about 4.83 kHz, 4.00 kHz, 3.17 kHz and 2.34 kHz, can be generated with a relatively small high frequency array {TMA: 13.00 kHz; OMNI: 17.00 kHz}. The advantages of an active parametric sonar accrue from a wide bandwidth capability, high angular resolution with reduced sidelobes, and a nearly constant beamwidth over a very wide frequency range. Parametric sonar operational performance is usually better than a conventional linear system due to reduced propagation loss, reduced reverberation, fewer multipaths and an improved data rate. Thus, multi-dimensional arrays can be employed to trade-off design complexity for improved parametric performance.

(77) Unique Features of TMA Technology

(78) For Medical Applications

(79) A TMA approach according to an embodiment can be characterized by at least six (6) unique features for obtaining space-frequency diversity beamforming for various medical applications, namely: a. Inducing both periosteal and endosteal healing of a long-bone fracture at the same time; b. Enhancing the production of shear waves to promote maximum integrin response of the bone tissue extracellular matrix (ECM), in lieu of acoustic modal converters; c. Making small angular adjustments to the insonifying beams in the fracture channel and along the bone tissue to efficiently improve overall healing; d. Minimizing clutter interference in tissue parameter estimation; e. Coding tissue space in frequency to enhance parameter estimation with FFT processing on receive; and f. Scanning tissue regions of interest in the high MHz (25-200 MHz) and low GHz (1-5 GHz) frequency domains for Dermatological and Scanning Acoustic Microscope applications, respectively.

(80) From the foregoing, it will be appreciated by those skilled in the art that embodiments of the present disclosure can provide an effective method and apparatus that overcomes many of the limitations associated with the mechanical stimulation of biological materials.

(81) For Underwater Sonar Applications

(82) A TMA approach according to an embodiment can be characterized by at least eight (8) unique features for obtaining space-frequency diversity beamforming for underwater sonar applications, namely: a. Greater near-field coverage by eliminating the blind zone; b. More detection resistant due to shorter transmission time; c. Reduced cross-beam reverberation per receiving beam; d. Simultaneous bearing-detection information; e. Inherent reduction of mutual sonar interference; f. Lower input power requirement for given OMNI source level; g. Coding underwater volumes in frequency to enhance parameter estimation with FFT processing on receive; and h. Reduced hardware complexity for multi-beam formation.

(83) It will also be readily appreciated by one with ordinary skill in the art use a method and apparatus according to embodiments of the disclosure in a variety of sonar applications, such as multi-mode echo-ranging, the detection of deep moored mines, and submarine active sonar. It may also provide a multi-beam formation capability for circular and planar arrays.

(84) Although certain exemplary embodiments of the present disclosure have been specifically described herein, it will be apparent to those skilled in the art to which embodiments of the disclosure pertain that variations and modifications of the various embodiments shown and described herein may be made without departing from the spirit and scope of this disclosure.

REFERENCES

(85) U.S. Pat. No. 2,426,460 System for Locating a Radiated-Signal Reflector, Lewis H. M., Aug. 26, 1947. U.S. Pat. No. 2,852,772, Receiver Scanning System, Gitzdanner L. G., Sep. 16, 1958. U.S. Pat. No. 3,012,244, Receiver Scanning System, Langenwalter D. F. and Stevenson M., Dec. 5, 1961 1. Winder A. A. and Loda C. J., Introduction to Acoustical Space-Time Information Processing, Office of Naval Research, Washington, DC/ONR Report ACR-63, January 1963. 2. Horton J. W., Fundamentals of Sonar, U.S. Naval Institute, Annapolis, MD), 1957. 3. Urick R. J., Principles of Underwater Sound for Engineers, McGraw-Hill Book Co, New York, 1967. 4. Shelkunoff S. A., Theory of Antennas of Arbitrary Size and Shape, Proc. IRE, Vol. 29, 493-521, September 1941. 5. Shelkunoff S. A., A Mathematical Theory of Linear Arrays, BSTJ, Vol. 22, 80-107, January 1943. 6. Shnitkin H., Survey of Electronically Scanned Antennas, The Microwave Journal, Part 1:67-72, December 1960; Part 2: 57-64, January 1961. 7. Shanks H. E., A New Technique for Electronic Scanning, IRE Transactions on Antennas and Propagation, 162-166, March 1961. 8. Kummer W. H., Villeneuve A. T., Fong T. S., and Terrio F. G., Ultra-Low Sidelobes from Time-Modulated Arrays, IEEE Trans on Antennas and Propagation, 633-639, November 1963. 9. Davies D. E. N. and McCartney B. S., Cylindrical Arrays with Electronic Beam Steering, Proc. IEEE, Vol. 112, No. 3, March 1965. 10. Johnson, M. A., Phased-Array Beam Steering by Multiplex Sampling, Proc. IEEE, Vol. 56, No. 11, 1802-1811, November 1968. 11. Cole B. F. and Hanrahan J. J., Influence of Beamwidth and Multiple Transmissions on Reverberation-Limited Sonars, JASA, Vol. 43, No. 6, 1373-1377, 1968. 12. Eyring C., Christensen R., Raitt R., Reverberation in the Sea, JASA, Vol. 20, 462-475, 1948. 13. Westerfelt P. J., Parametric Acoustic Array, J. Acoust. Soc. Am., 35, 535-537, 1963. 14. Voglis G. M., A General Treatment of Modulation Scanning as Applied to Acoustic Linear Arrays, Ultrasonics, Part 1: Vol. 9, No. 3, 142-153, July 1971; Part 2: Vol. 10, No. 3, 103-113, May 1972. 15. Winder A. A., Sonar System Technology, IEEE Trans Sonics & Ultrasonics, SU-22, No. 5, 291-332, September 1975. 16. Haykin S., Multiple-Beam Sampler for Continuously Scanned Array Antennas, IRE Trans. on Antennas and Propagation, 526-528, July 1976. 17. Machado P., Li J, Blackman R., Liu J-B., Kepler C., Fang T., Muratore R., Winder J., Winder A. and Forsberg F., Comparison Between Clinically Available Low Intensity Pulsed Ultrasound (LIPUS) and a Novel BiModal Acoustic Signal System for Accelerating Fracture Healing, UFFC pending acceptance, Winter 2021. 18. Winder A. and Muratore R., Ultrasound Stimulation in Musculo-Skeletal Tissue Structures, WIPO International Publication No. WO 2020/154633, Jul. 30, 2020.