CN102510548A - Output adjustment method for sound energy obtained by sound wave interaction in nonlinear medium - Google Patents

Output adjustment method for sound energy obtained by sound wave interaction in nonlinear medium Download PDF

Info

Publication number
CN102510548A
CN102510548A CN2011103161251A CN201110316125A CN102510548A CN 102510548 A CN102510548 A CN 102510548A CN 2011103161251 A CN2011103161251 A CN 2011103161251A CN 201110316125 A CN201110316125 A CN 201110316125A CN 102510548 A CN102510548 A CN 102510548A
Authority
CN
China
Prior art keywords
wave
amplitude
sound
omega
weak signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2011103161251A
Other languages
Chinese (zh)
Other versions
CN102510548B (en
Inventor
杨德森
兰朝凤
时胜国
李思纯
洪连进
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN 201110316125 priority Critical patent/CN102510548B/en
Publication of CN102510548A publication Critical patent/CN102510548A/en
Application granted granted Critical
Publication of CN102510548B publication Critical patent/CN102510548B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

The invention provides an output adjustment method for sound energy obtained by sound wave interaction in a nonlinear medium. The output adjustment method comprises the following steps of: (a) performing nonlinear interaction on a pump wave with the frequency of w3 and a weak signal wave with the frequency of omega1, and generating a resonance wave with the frequency of omega2; (b) respectively calculating amplitude values B1(x), B2(x) and B3(x) of the three waves at the displacement x after interaction according to the frequencies omega3, omega1 and omega2 of the pump wave, the weak signal wave and the resonance wave; and (c) adjusting the output energy of the three waves according to the change characteristics of the amplitude values B1(x), B2(x) and B3(x) of the pump wave, the weak signal wave and the resonance wave. According to the output adjustment method, the energy of the generated sound wave is changed according to a pulse rule on the basis of a basic theory of interaction of sound waves in optics and hydroacoustics; therefore, the output energy of each wave can be adjusted according to actual requirements.

Description

The output adjusting method of acoustic energy after sound wave interacts in the nonlinear dielectric
Technical field
The present invention relates to sound wave nonlinear interaction field, relate in particular to the interactional output adjusting method of nonlinear parameter between the sound wave.
Background technology
The nonlinear interaction that proposes utilization sound and sound since nineteen sixty Vista Wei Erte has formed since the parametric array notion, and the development of nonlinear acoustics is more and more faster, uses more and more extensivelyr, and its research and application obtained a lot of new progresses.The nonlinear interaction of underwater acoustic wave can be thought to have formed variable element reflector or variable element receiver in interaction zone.The notion of variable element interaction process is born in the study of radio, when its essence is electric capacity or the variation of inductance generating period property in certain oscillating circuit, just can make ultra-weak electronic signal occur amplifying or weakening.Take place in the process of nonlinear interaction between the sound wave, the parametric interaction problem is occupied extremely important status.In the process of acoustic wave energy transmission, what play the role of a nucleus is to induce diffusion process, makes weak signal ripple and pump wave interaction produce each order harmonics simultaneously; Show as the unsteadiness of applicator quantum splitting; The vibration of ripple, the transfer of spectrum energy etc., but observe law of conservation of energy in the interaction process.The interactional condition of nonlinear parameter takes place being medium between the sound wave, to have one of strong non-linear or interactional sound wave be high-power sound wave, can exciting media non-linear, this sound wave can low-frequency sound wave also can be a low-frequency sound wave.Existing a lot of scholar's research the parametric interaction problem of sound wave, Fenlon utilizes fourier progression expanding method, has provided the interact form of expression of each frequency content of back of a plurality of large amplitude simple harmonic quantity sound waves.Be the basis with the Fenlon theory in the prior art, studied in the air dielectric good amplitude wave to the inhibition of sound wave, also relevant for of the research of good amplitude wave to the scale-up problem of smooth sea.Utilize two hyperacoustic nonlinear interactions theories in the spectral factorization method research aqueous medium in addition, the variation of low frequency wave energy has been discussed.Do not generate ω after interacting about the sound wave variable element 2Three row magnitudes of acoustic waves are with the variation characteristic of distance and frequency, the acoustic wave energy energy transfer process with propagation distance and frequency change after the new sound field.
Summary of the invention
The object of the present invention is to provide the interact output adjusting method of back acoustic energy of sound wave in a kind of nonlinear dielectric of the output energy that can regulate each train wave according to actual needs.
The objective of the invention is to realize like this:
May further comprise the steps:
(a) frequency is ω 3Pump ripple and frequency be ω 1Weak signal ripple generation nonlinear interaction, the generation frequency is ω 2Resonance wave;
(b) according to said pump ripple, weak signal ripple resonant wave frequency ω 3, ω 1And ω 2, the amplitude B at three train wave displacement x places after calculating interacts respectively 1(x), B 2(x) and B 3(x);
(c) according to the pump wave amplitude B that obtains 3(x), weak signal wave amplitude B 1(x) and resonance wave amplitude B 2(x) variation characteristics realize the output energy adjustment to three train waves.
The present invention can also comprise:
1, said resonance wave is and frequently resonance wave or difference frequency resonance wave, i.e.
Figure BDA0000099677330000021
perhaps
Figure BDA0000099677330000022
2, the amplitude B at said calculating three train wave displacement x places 1(x), B 2(x) and B 3(x) method is:
Step (b1) is calculated and is obtained at the displacement x place sound wave vibration velocity v (x) according to pump ripple and weak signal ripple interactional Burgers (Burgers) equation in nonlinear dielectric;
The sound wave vibration velocity v (x) that step (b2) obtains according to step (b1) calculates the wave amplitude equation that obtains after three row sound waves interact;
Three train waves after step (b3) is obtained to interact by the wave amplitude Equation for Calculating are at the amplitude B at displacement x place 1(x), B 2(x) and B 3(x).
3, said Burgers (Burgers) equation is:
Wherein, v is a sound wave vibration velocity,
Figure BDA0000099677330000024
B/A is the nonlinear parameter of medium, is the ratio of quadratic term coefficient and linear coefficient in the state equation taylor series expansion, and it is the basic parameter of nonlinear acoustics, c 0Be the static velocity of sound, ρ 0Be the density of nonlinear dielectric,
Figure BDA0000099677330000025
Be time delay, x is a measuring distance.B is the medium coefficient of viscosity,
Figure BDA0000099677330000026
ζ is for cutting the coefficient of viscosity, and η is the body coefficient of viscosity,
Figure BDA0000099677330000027
Be temperature conductivity coefficient, c v, c pBe electric capacity specific heat and voltage specific heat;
The sine wave that sound source is sent at the x=0 place, sound wave vibration velocity v (x) is expressed as at the displacement x place to calculate acquisition:
Figure BDA0000099677330000028
Wherein, A 1(x), A 2(x), A 3(x) be the complex amplitude of three row sound waves, c is a constant.
4, nonlinear dielectric is a water.
5, the wave amplitude equation after three row sound waves interact:
Figure BDA0000099677330000032
Wherein,
Figure BDA0000099677330000034
is the phase mismatch factor.If
Figure BDA0000099677330000035
then three ripples is phase matched; Be equivalent to three phonon conservations of momentum,
Figure BDA0000099677330000036
is dissipation factor.
6, when the resonance wave that produces be and during the frequency resonance wave; Promptly when
Figure BDA0000099677330000037
, the complex amplitude of three train waves is expressed as the form of real amplitude and phase place:
Wherein, B 1(x), B 2(x), B 3(x) and
Figure BDA0000099677330000039
Be frequencies omega 1, ω 2, ω 3The real number amplitude and the phase constant of sound wave;
Dissipation factor The pump intensity of wave does not change initial condition because generate resonance wave
Figure BDA00000996773300000311
Figure BDA00000996773300000312
The time, trying to achieve the back frequency of sound wave that interacts is ω 1, ω 2Amplitude do
Figure BDA00000996773300000313
7, when the resonance wave that produces is the difference frequency resonance wave; Promptly when , the complex amplitude of three train waves is expressed as the form of real amplitude and phase place:
Figure BDA00000996773300000315
Wherein, B 1(x), B 2(x), B 3(x) and
Figure BDA0000099677330000041
Be frequencies omega 1, ω 2, ω 3The real number amplitude and the phase constant of sound wave;
During dissipation factor
Figure BDA0000099677330000042
initial condition
Figure BDA0000099677330000043
Figure BDA0000099677330000044
, the amplitude of trying to achieve the back three row sound waves that interact is:
Figure BDA0000099677330000045
Figure BDA0000099677330000046
Figure BDA0000099677330000047
Figure BDA0000099677330000048
is the Jacobi elliptic function;
Figure BDA0000099677330000049
wherein,
Figure BDA00000996773300000410
8, according to the back pump wave amplitude B that interacts 3(x), weak signal wave amplitude B 1(x) resonant wave amplitude B 2(x), regulate the output energy of weak signal ripple or resonance wave with the variation characteristics of propagation distance.
9, according to the back weak signal wave amplitude B that interacts 1(x) resonant wave amplitude B 2(x) with pump wave frequency ω 3The variation characteristics, regulate the output energy of weak signal ripple or resonance wave;
Perhaps, according to the back weak signal wave amplitude B that interacts 1(x)) resonant wave amplitude B 2(x) with weak signal wave frequency ω 1The variation characteristics, regulate the output energy of weak signal ripple or resonance wave;
Perhaps, according to the back weak signal wave amplitude B that interacts 1(x) resonant wave amplitude B 2(x) with pump wave amplitude B 3(0) variation characteristics, the output energy of adjusting weak signal ripple or resonance wave.
Method of the present invention in conjunction with the interactional basic principle of sound wave in optics and the marine acoustics, is considered the dissipation effect that high-power sound wave is propagated, sound wave ω in the aqueous medium in nonlinear dielectric 1And ω 3The second order nonlinear effect that (weak signal ripple and pump ripple) causes because of nonlinear interaction generates and frequency sound wave ω 2Or difference frequency sound wave ω 2, utilized runge kutta method (Runger-Kutta) numerical analysis the sound wave variable element back that interacts generate ω 2Resonance wave after three row magnitudes of acoustic waves with the variation characteristic of distance and frequency, acoustic wave energy energy transfer process with propagation distance and frequency change; Confirm that the energy variation that generates sound wave presents the pulsation rule, can regulate the output energy of each train wave according to actual needs.
Description of drawings
Fig. 1 is that magnitudes of acoustic waves is with the propagation distance change curve;
Fig. 2 is that acoustic wave energy is with the propagation distance change curve;
Fig. 3 is that acoustic wave energy is with pump wave frequency change curve;
Fig. 4 is weak signal wave amplitude (ω 1Ripple) with pump wave amplitude and frequency variation curve;
Fig. 5 is weak signal wave energy (ω 1Ripple) with pump wave amplitude and frequency variation curve;
Fig. 6 is that magnitudes of acoustic waves is with the propagation distance change curve; Wherein
(a) be pump wave frequency ω 3During for 2kHz, weak signal wave amplitude (ω 1Ripple), resonance wave amplitude (ω 2Ripple), pump wave amplitude (ω 3Ripple) with the propagation distance change curve;
(b) be pump wave frequency ω 3During for 5kHz, weak signal wave amplitude (ω 1Ripple), resonance wave amplitude (ω 2Ripple), pump wave amplitude (ω 3Ripple) with the propagation distance change curve;
(c) pump wave frequency ω 3During for 10kHz, weak signal wave amplitude (ω 1Ripple), resonance wave amplitude (ω 2Ripple), pump wave amplitude (ω 3Ripple) with the propagation distance change curve;
Fig. 7 is different pump wave frequency ω 3The time, resonance wave energy (ω 2Ripple) with the propagation distance change curve;
Fig. 8 is different pump wave amplitude (ω 3Ripple) time, weak signal wave energy (ω 1Ripple) with the propagation distance change curve.
Fig. 9 is the flow chart of the interactional output adjusting method of nonlinear parameter between the sound wave.
Embodiment
To combine accompanying drawing and embodiment that technical scheme of the present invention is explained in more detail below.
The present invention is the basis with Burgers (Burgers) equation; Draw the expression formula of the amplitude of each sound field after weak signal ripple and the pump wave interaction; Utilize the emulation of Runger-Kutta method to obtain amplitude or the energy of each sound field change curve with propagation distance; The non-linear variable element interaction of the sound wave back energy transfer process that obtains directly perceived is utilized said transfer, realizes the output energy adjustment to three train waves.
The variable element interaction phenomenon of sound wave is explained with the quantum mechanics language, promptly can regard high frequency pump ripple ω as 3Phonon splits into two low frequency omega 1, ω 2Phonon perhaps splits into the more process of low frequency phonon.The essence of three-wave interaction is the non-linear variable element of sound wave to have taken place interact.In aqueous medium, as pump ripple ω 3With weak signal sound wave ω 1Nonlinear interaction takes place generate ω 2The time, realize Best Coupling (sound wave resonance), must observe the following law of conservation of energy and the law of conservation of momentum
In the formula, for composing bright gram constant, k 1, k 2, k 3Be the sound wave wave number.When satisfying formula (1), can only consider ω 1, ω 2, ω 3The coupling of this three row frequency wave, and can not consider any coupling of this three row frequency wave and all other frequency waves fully.Sound wave is the acoustic energy conservation in the nonlinear interaction process, and when formula (1) was satisfied in acoustic wave energy generation transfer, energy exchange was maximum between the sound wave.
In desirable nonlinear dielectric,
Figure BDA0000099677330000063
weak signal ripple and pump ripple interaction process in nonlinear dielectric can be write as general Burgers equation form when sound Reynolds number
Figure BDA0000099677330000062
In the formula, v is a underwater acoustic wave vibration velocity, B/A is the nonlinear parameter of medium, is the ratio of quadratic term coefficient and linear coefficient in the state equation taylor series expansion, and it is the basic parameter of nonlinear acoustics, c 0Be the static velocity of sound, ρ 0Be the density of aqueous medium,
Figure BDA0000099677330000066
Be time delay, x is a propagation distance.B is the medium coefficient of viscosity,
Figure BDA0000099677330000067
ζ is for cutting the coefficient of viscosity, and η is the body coefficient of viscosity,
Figure BDA0000099677330000068
Be temperature conductivity coefficient, c v, c pBe electric capacity specific heat and voltage specific heat.
Suppose that sound source locates to launch sinusoidal signal at , provide the form that displacement x place sound wave separates according to the Burgers equation and do
Figure BDA00000996773300000610
In the formula, A 1(x), A 2(x), A 3(x) be the complex amplitude of sound wave, c is a constant.
Consider the dissipation effect that high-power sound wave is propagated in nonlinear dielectric, this moment, acoustic energy was non-vanishing to the derivative of displacement x, and the wave amplitude equation does after the interaction that obtains three row sound waves in the formula (3) of deriving
Figure BDA00000996773300000611
In the formula,
Figure BDA00000996773300000612
is the phase mismatch factor.If
Figure BDA00000996773300000613
then three ripples is phase matched, be equivalent to three phonon conservations of momentum. is dissipation factor.
Formula (4) is the Non-Self-Governing form; For make find the solution on this equation mathematics more convenient; Numerical computations also is convenient to handle, and does simple conversion and turns to autonomous form
Figure BDA0000099677330000071
Formula (5) is a form common in the three-wave interaction problem.
Conservation of momentum during three ripple Best Coupling; Get
Figure BDA0000099677330000072
and multiply by
Figure BDA0000099677330000074
successively when x gets arbitrary value when three formulas of the dissipation effect of ignoring sound wave
Figure BDA0000099677330000073
formula (4), formula (4) becomes through after the integral operation
Figure BDA0000099677330000075
Per two equations in the formula (6) subtract each other all can obtain another equation, is not independently each other.Second with formula (6) is multiplied by ω 1, the 3rd is multiplied by ω 2, and utilize After conversion, can obtain law of conservation of energy
Figure BDA0000099677330000077
Formula (7) shows; If total acoustic energy at coordinate x place be
Figure BDA0000099677330000078
following formula (7) just
Figure BDA0000099677330000079
explain that its gross energy is constant in the interactional process of three row frequency sound waves.In other words, the energy value of sound wave just exchanges between each frequency sound wave, and medium is not participated in, and only plays instrumentality, and this is interactional characteristics of all parameters just also.
Formula (6) can further be written as
Figure BDA00000996773300000710
In the formula, E 1(x), E 2(x), E 3(x) be that sound wave is at the acoustic energy value apart from sound source x place.Between formula (8) and the nonlinear optics derivation light wave interactional Manley-Rowe form class seemingly, the nonlinear interaction form of underwater acoustic wave and meaning also can be used for reference in the optics the interactional meaning of light wave and analyze and handle thus.According to the Manley-Rowe theorem, formula (6) differential is got
Figure BDA0000099677330000081
Utilize the physical definition of sound energy flux density, the phonon average flux expression formula for then formula (9) be rewritten as
Figure BDA0000099677330000083
According to formula (10), the acoustic energy transfer relationship of nonlinear interaction between the sound wave is carried out qualitative analysis.Formula (10) shows in the sound wave interaction process that frequency is ω 3Phonon of the every minimizing of sound wave, then frequency is ω 1And ω 2Sound wave all to increase a phonon, this that is to say powerful sound wave ω in accordance with the acoustic energy law of conservation 3Incide in the nonlinear dielectric, making the incident frequency is ω 1Weak sound wave convert frequency to through difference frequency and do
Figure BDA0000099677330000084
(perhaps ) sound wave, this process is the process of sound wave decay, i.e. dN 1, dN 2Increment for just, and dN 3Increment for negative.Can estimate the action effect of nonlinear interaction thus according to formula (9) and formula (10).
The complex amplitude form of expression of three row sound waves does
Figure BDA0000099677330000086
In the formula, B 1(x), B 2(x), B 3(x) and
Figure BDA0000099677330000087
Be frequencies omega 1, ω 2, ω 3The real number amplitude and the phase constant of sound wave.
Below divide two kinds of situation to draw nonlinear interaction and produce the resonance wave amplitude
1 and sound wave frequently
Beginning not have frequency in the sound field is ω 2Sound wave, this component is to do by frequency
Figure BDA0000099677330000088
Sound wave coupling form.If the pump intensity of wave can not have bigger change because of generating sound wave; Then can the pump wave field be approximately a constant field, promptly
Figure BDA0000099677330000089
has the formula (4) of formula (11) substitution three ripple couplings
Figure BDA0000099677330000091
In the formula;
Figure BDA0000099677330000092
is approximately at 0 o'clock in dissipation, and preceding two equations of formula (12) can also be expressed as
Figure BDA0000099677330000093
Wushu (13) generation is to the 3rd equation of formula (12), and through after the integral operation does
In formula (14); When initial condition ; Formula (14) both members is equated,
Figure BDA0000099677330000096
arranged discuss below
Figure BDA0000099677330000097
at this moment and get
Figure BDA0000099677330000098
so formula (13) can further be written as
Because formula (15) is linear, separate the form of being write as following formula to it so long
Figure BDA00000996773300000910
For asking the particular solution of formula (15), with formula (16) substitution formula (15), make that the determinant of equation is zero after the substitution, promptly
Figure BDA00000996773300000911
Formula (17) can be confirmed two value l of l 1, l 2, formula (15) separating and can be rewritten into thus
Figure BDA00000996773300000912
Wherein, C I, jBe undetermined constant, by
Figure BDA00000996773300000913
The boundary condition at place
Confirm.
Because formula (15) only can be found the solution with mathematical method or numerical method,, can suppose when pumping source is strong for trying to achieve its analysis result
Figure BDA0000099677330000102
The back frequency of sound wave that obtains thus interacting is ω 1, ω 2The amplitude analytic solutions
Figure BDA0000099677330000103
Find out B by formula (20) 2(x) amplitude is by weak signal strength B in the sound field 1(0) confirms.
2, difference frequency sound wave
Beginning not have frequency in the sound field is ω 2Sound wave, this component is to do by frequency
Figure BDA0000099677330000104
Sound wave coupling form.When considering pump ripple ω 3During energy variation in the nonlinear interaction process, obtain following ACOUSTIC WAVE EQUATION according to three ripple coupled wave equation formulas (4)
Figure BDA0000099677330000105
Ordinary circumstance following formula (21) only can be found the solution with mathematical method or numerical method; For trying to achieve its analysis result; When the dissipation factor
Figure BDA0000099677330000106
of three ripples; Utilizing
Figure BDA0000099677330000107
locates; Initial condition formula (19) and
Figure BDA0000099677330000108
are similar with the solution of formula (12); Per two equations combination of formula (21) obtains after integral and calculating
Figure BDA0000099677330000109
The form that formula (22) substitution formula (23) obtains variables separation does
Figure BDA0000099677330000111
In the formula;
Figure BDA0000099677330000112
obtains after introducing
Figure BDA0000099677330000113
conversion for asking separating of formula (23)
Figure BDA0000099677330000114
The left-hand component of formula (24) contains the single order ellptic integral, utilizes ellptic integral
Figure BDA0000099677330000115
After variable
Figure BDA0000099677330000116
replacement, obtain
Figure BDA0000099677330000117
In the formula;
Figure BDA0000099677330000118
be κ (y thus; K) also can be write as κ (ψ; K) form; The inverse function
Figure BDA00000996773300001110
of claiming
Figure BDA0000099677330000119
is the Jacobi elliptic function, and two kinds of forms of other of elliptic function can be write
Figure BDA00000996773300001111
Figure BDA00000996773300001112
During as ; The result of integration type (25) is a complete elliptic integral of the first kind; Its value is designated as ∝, and the expression-form of ∝ is:
Figure BDA00000996773300001114
According to the ellptic integral of the Jacobian transform (24) of the complete elliptic function of the first kind, the ripple of the back three row sound waves that interact is separated and can be write as the following formula form
Figure BDA00000996773300001115
Formula (27) provides behind the three-wave interaction sound wave amplitude form of expression that the elliptic function form of substantial connection is arranged with variable
Figure BDA00000996773300001116
, has the obvious periodic feature of elliptic function.
In the formula (20) in
Figure BDA00000996773300001117
and the formula (27)
Figure BDA00000996773300001118
all represent the sound wave phase place; And phase place and propagation distance x are proportional, and the sound wave amplitude is the function of measuring distance.Find out B by formula (27) 2(x) maximum is by stronger in a sound field pumping source intensity B 3(0) confirms.
Below utilize the emulation of Runger-Kutta method that the amplitude or the energy of each sound field of obtaining carried out emulation, obtain amplitude or the energy of each sound field change curve, wherein ω with propagation distance 1Weak signal ripple initial magnitude is 1, and frequency is 1kHz, and the non-linear variable element of aqueous medium is 3.6.
(1) do not considering to have carried out following simulation analysis according to separating (20) under the situation that the pumping source acoustic energy changes.A: measuring distance x is 178.5 meters, and the pumping source initial magnitude is 10, and frequency is 80000rad/s; Amplitude is as shown in Figure 1 with the variable in distance curve, and energy is as shown in Figure 2 with the variable in distance curve.B: measuring distance is 78.5 meters, and the pumping source initial magnitude is taken as 10, and frequency is 1000rad/s~80000rad/s, and energy is as shown in Figure 3 with the pumping source frequency variation curve; When the pumping source amplitude is 1100, the weak signal wave amplitude is as shown in Figure 4 with pumping source amplitude and frequency variation curve; ω 1The energy difference is as shown in Figure 5 with pumping source amplitude and frequency variation curve before and after the ripple nonlinear interaction.
When (2) considering the pumping source energy variation, utilize the Runger-Kutta method that sound wave interaction back acoustic signature is carried out numerical analysis.Parameter: the coefficient of viscosity of aqueous medium does when getting 10 ℃
Figure BDA0000099677330000121
The dissipation parameter
Figure BDA0000099677330000122
A: the pumping source initial magnitude is 10, and frequency is 2kHz, 5kHz, and 10kHz, it is as shown in Figure 6 with the measuring distance change curve to obtain sound wave interaction back amplitude, generates ω 2Wave energy is as shown in Figure 7 with the variable in distance curve.B: the pumping source frequency is taken as 2kHz, and amplitude is taken as 5,10,20, generates ω 1Wave energy is as shown in Figure 8 with the variable in distance curve.
B among each figure 1, B 2, B 3Frequency is ω after representing the sound wave variable element to interact respectively 1, ω 2And ω 3The amplitude of sound wave.DE 1Represent ω 1Energy difference before the wave interaction after energy and the interaction, dE 2After interacting, representative generates ω 2The energy value of ripple.
Weak signal ripple ω among Fig. 1 1With pump ripple ω 3Through and frequently
Figure BDA0000099677330000123
Generate ω 2Ripple, sound wave is the acoustic energy conservation in the process that nonlinear interaction takes place, and the acoustic wave energy that generates mainly comes from weak signal ripple ω 1So the weak signal ripple shifts energy to newly-generated sound wave ω again 2Process in the energy of weak signal ripple will reduce.See newly-generated sound wave ω from expression formula (20) 2Amplitude presents periodic variation, generates thus and frequency ω 2Sound field can not unrestrictedly increase, the power entrained when weak signal is evacuated, ω 2Power increase process gradually in, saturated phenomenon can appear, ω 2Ripple can be passed to ω to energy back again 1, as can be seen from Figure 1 waveform presents the variation of pulsation trend.Fig. 2~5 simulation results see that the energy of weak signal ripple all can reduce as long as between the sound wave nonlinear interaction has taken place, and greater than zero, maximum decrease can reach thirties dB to difference all the time, and magnitude and energy variation all present the cycle and pulse.By (a) (b) (c) figure contrast discovery of Fig. 6, initial sound wave ω 2Acoustic pressure be zero, weak signal sound wave ω 1With pump ripple ω 3Nonlinear interaction takes place
Figure BDA0000099677330000124
Shift and give sound wave ω 2, sound wave ω thus 2Amplitude start from scratch gradually and to increase, want to make weak signal ω 1Magnitudes of acoustic waves is farthest amplified in certain segment distance to be needed to satisfy
Figure BDA0000099677330000131
When
Figure BDA0000099677330000132
The time, added damping parameter, find out that from figure damping parameter can change the range attenuation rate of sound wave, to the not influence of propagation law of sound wave, frequency is high more, and it is big more to decay.Fig. 7 shows the increase with the pumping source frequency, ω 2Acoustic wave energy reduces; When Fig. 8 was illustrated in the pumping source frequency and is the twice of signal frequency, along with the increase of pumping source frequency, the amplification quantity of weak signal increased.Increase or reduce trend and still present the variation of pulsation rule.
Can find out through simulation result:
(1) between the sound wave the interactional process of non-linear variable element takes place in, generate new frequency content and shift with energy.
(2) when generating with sound wave frequently in the non-linear variable element interaction process of sound wave and frequently the energy of ripple mainly comes from the weak signal ripple because the conservation of energy in the interaction process, institute so that the weak signal ripple reduce certain section propagation distance self-energy.
(3) when the non-linear variable element interaction process of sound wave generates the difference frequency sound wave; The energy of difference frequency ripple mainly comes from pumping source; When satisfying between applicator Best Coupling and being two frequencys multiplication, the weak signal wave energy will farthest be amplified in certain section propagation distance.
(4) between the sound wave in when, non-linear variable element taking place interacting, energy shifts to have periodically, thereby causes the energy value periodically-varied of weak signal ripple, has long-term modularity.
Therefore, after non-linear variable element takes place between the sound wave interacting, can realize that the trend that the weak signal wave energy presents pulsation changes, can realize the output energy adjustment of three train waves making the output energy reduce or amplify through the parameter of selecting to be fit to.
The present invention also can have other various embodiments; Under the situation that does not deviate from spirit of the present invention and essence thereof; Those of ordinary skill in the art are when making various corresponding changes and distortion according to the present invention, and these change and be out of shape the protection range that all should belong to the appended claim of the present invention accordingly.

Claims (8)

1. the interact output adjusting method of back acoustic energy of sound wave in the nonlinear dielectric is characterized in that may further comprise the steps:
(a) frequency is ω 3Pump ripple and frequency be ω 1Weak signal ripple generation nonlinear interaction, the generation frequency is ω 2Resonance wave;
(b) according to said pump ripple, weak signal ripple resonant wave frequency ω 3, ω 1And ω 2, the amplitude B at three train wave displacement x places after calculating interacts respectively 1(x), B 2(x) and B 3(x);
(c) according to the pump wave amplitude B that obtains 3(x), weak signal wave amplitude B 1(x) and resonance wave amplitude B 2(x) variation characteristics realize the output energy adjustment to three train waves.
2. the output adjusting method of acoustic energy after sound wave interacts in the nonlinear dielectric according to claim 1 is characterized in that: said resonance wave is and frequency resonance wave or difference frequency resonance wave, i.e. ω 21+ ω 3Perhaps ω 231
3. the output adjusting method of acoustic energy is characterized in that: the amplitude B at said calculating three train wave displacement x places after sound wave interacted in the nonlinear dielectric according to claim 1 1(x), B 2(x) and B 3(x) method is:
Step (b1) is calculated and is obtained at the displacement x place sound wave vibration velocity v (x) according to pump ripple and weak signal ripple interactional Burgers equation in nonlinear dielectric;
The sound wave vibration velocity v (x) that step (b2) obtains according to step (b1) calculates the wave amplitude equation that obtains after three row sound waves interact;
Three train waves after step (b3) is obtained to interact by the wave amplitude Equation for Calculating are at the amplitude B at displacement x place 1(x), B 2(x) and B 3(x).
4. the output adjusting method of acoustic energy after sound wave interacts in the nonlinear dielectric according to claim 3, it is characterized in that: said Burgers equation is:
∂ v ∂ x - β c 0 2 v ∂ v ∂ τ = b 2 c 0 2 ρ 0 ∂ 2 v ∂ τ 2
Wherein, v is a sound wave vibration velocity, β=1+B/2A, and B/A is the nonlinear parameter of medium, is the ratio of quadratic term coefficient and linear coefficient in the state equation taylor series expansion, it is the basic parameter of nonlinear acoustics, c 0Be the static velocity of sound, ρ 0Be the density of nonlinear dielectric, τ=t-x/c 0Be time delay, x is a measuring distance, and b is the medium coefficient of viscosity,
Figure FDA0000099677320000012
ζ is for cutting the coefficient of viscosity, and η is the body coefficient of viscosity,
Figure FDA0000099677320000013
Be temperature conductivity coefficient, c v, c pBe electric capacity specific heat and voltage specific heat;
The sine wave that sound source is sent at the x=0 place, sound wave vibration velocity v (x) is expressed as at the displacement x place to calculate acquisition:
v ( x , t ) = Σ n = 1,2,3 A n ( x ) exp ( i ω n ( t - x c 0 ) ) + c
Wherein, A 1(x), A 2(x), A 3(x) be the complex amplitude of three row sound waves, c is a constant.
5. the output adjusting method of acoustic energy after sound wave interacts in the nonlinear dielectric according to claim 1 is characterized in that: the wave amplitude equation after three row sound waves interact:
d A 1 ( x ) dx + δ 1 A 1 ( x ) = j A 2 * ( x ) A 3 ( x ) e jΔkx βω 1 c 0 2 d A 2 ( x ) dx + δ 2 A 2 ( x ) = j A 1 * ( x ) A 3 ( x ) e jΔkx βω 2 c 0 2 d A 3 ( x ) dx + δ 3 A 3 ( x ) = j A 1 * ( x ) A 2 ( x ) e - jΔkx βω 3 c 0 2
Wherein, k 11/ c 0, k 22/ c 0, k 33/ c 0, Δ k=k 3-k 1-k 2Be the phase mismatch factor.If Δ k=0; Then three ripples are phase matched; Be equivalent to three phonon conservations of momentum,
Figure FDA0000099677320000023
is dissipation factor.
6. the sound wave output adjusting method of back acoustic energy that interacts in the nonlinear dielectric according to claim 2 is characterized in that: when the resonance wave that produces for resonance wave frequently the time, i.e. ω 21+ ω 3The time, the complex amplitude of three train waves is expressed as the form of real amplitude and phase place:
Figure FDA0000099677320000024
Wherein, B 1(x), B 2(x), B 3(x) and
Figure FDA0000099677320000025
Be frequencies omega 1, ω 2, ω 3The real number amplitude and the phase constant of sound wave;
Dissipation factor δ i=0 (i=1,2), pump intensity of wave do not change initial condition B because generate resonance wave 2(0)=0,
Figure FDA0000099677320000026
The time, trying to achieve the back frequency of sound wave that interacts is ω 1, ω 2Amplitude do
B 1 ( x ) = B 1 ( 0 ) cos ( xβ B 3 ( 0 ) ω 1 ω 2 c 0 2 ) B 2 ( x ) = ω 2 ω 1 B 1 ( 0 ) sin ( xβ B 3 ( 0 ) ω 1 ω 2 c 0 2 ) .
7. the sound wave output adjusting method of back acoustic energy that interacts in the nonlinear dielectric according to claim 2 is characterized in that: when the resonance wave that produces is the difference frequency resonance wave, i.e. and ω 231The time, the complex amplitude of three train waves is expressed as the form of real amplitude and phase place:
Figure FDA0000099677320000032
Wherein, B 1(x), B 2(x), B 3(x) and
Figure FDA0000099677320000033
Be frequencies omega 1, ω 2, ω 3The real number amplitude and the phase constant of sound wave;
Dissipation factor δ i=0, i=1,2,3, initial condition B 2(0)=0,
Figure FDA0000099677320000034
The time, the amplitude of trying to achieve the back three row sound waves that interact is:
B 1 ( x ) = B 3 ( 0 ) ω 1 / ω 3 k dn ( κ + β B 3 ( 0 ) x ω 1 ω 3 kc 0 2 ) B 2 ( x ) = ω 2 / ω 3 B 3 ( 0 ) cn ( κ + β B 3 ( 0 ) x ω 1 ω 3 kc 0 2 ) B 3 ( x ) = B 3 ( 0 ) sn ( κ + β B 3 ( 0 ) x ω 1 ω 3 kc 0 2 )
(u k) is the Jacobi elliptic function to y=sn, cn (u)=(1-sn 2) 1/2, dn (u)=(1-k 2Sn 2) 1/2
Wherein, κ = ∫ 0 1 ( 1 - y 2 ) ( 1 - k 2 y 2 ) Dy , k ∈ ( 0,1 ) .
8. the output adjusting method of acoustic energy after sound wave interacts in the nonlinear dielectric according to claim 1 is characterized in that said real root is according to the pump wave amplitude B that obtains 3(x), weak signal wave amplitude B 1(x) and resonance wave amplitude B 2(x) variation characteristics realize that the output energy adjustment to three train waves is one of following:
(1) according to the pump wave amplitude B that obtains 3(x), weak signal wave amplitude B 1(x) resonant wave amplitude B 2(x) variation characteristics realize the output energy adjustment to three train waves, regulate the output energy of weak signal ripple or resonance wave;
(2) according to the back weak signal wave amplitude B that interacts 1(x) resonant wave amplitude B 2(x) with pump wave frequency ω 3The variation characteristics, regulate the output energy of weak signal ripple or resonance wave;
(3) according to the back weak signal wave amplitude B that interacts 1(x)) resonant wave amplitude B 2(x) with weak signal wave frequency ω 1The variation characteristics, regulate the output energy of weak signal ripple or resonance wave;
(4) according to the back weak signal wave amplitude B that interacts 1(x) resonant wave amplitude B 2(x) with pump wave amplitude B 3(0) variation characteristics, the output energy of adjusting weak signal ripple or resonance wave.
CN 201110316125 2011-10-18 2011-10-18 Output adjustment method for sound energy obtained by sound wave interaction in nonlinear medium Expired - Fee Related CN102510548B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110316125 CN102510548B (en) 2011-10-18 2011-10-18 Output adjustment method for sound energy obtained by sound wave interaction in nonlinear medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110316125 CN102510548B (en) 2011-10-18 2011-10-18 Output adjustment method for sound energy obtained by sound wave interaction in nonlinear medium

Publications (2)

Publication Number Publication Date
CN102510548A true CN102510548A (en) 2012-06-20
CN102510548B CN102510548B (en) 2013-12-18

Family

ID=46222594

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110316125 Expired - Fee Related CN102510548B (en) 2011-10-18 2011-10-18 Output adjustment method for sound energy obtained by sound wave interaction in nonlinear medium

Country Status (1)

Country Link
CN (1) CN102510548B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103576574A (en) * 2013-10-10 2014-02-12 哈尔滨工程大学 Phase control method of acoustic energy conversion under nonlinear mutual action of three acoustic wavetrains
CN105301582A (en) * 2015-10-14 2016-02-03 哈尔滨工程大学 Underwater acoustic field weak signal enhancement method
CN106841382A (en) * 2017-01-23 2017-06-13 哈尔滨工程大学 Non-uniform mixing medium nonlinear factor measuring method based on three couple waves interaction principle
CN112954562A (en) * 2021-01-27 2021-06-11 哈尔滨工程大学 Acoustic signal enhancer based on parameter excitation

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1638533A (en) * 1996-07-17 2005-07-13 美国技术公司 Sound heterodyne apparatus and method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1638533A (en) * 1996-07-17 2005-07-13 美国技术公司 Sound heterodyne apparatus and method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨德森等: "《水中声波之间非线性作用时能量变化的分析》", 《哈尔滨工程大学学报》 *
董雷等: "《大振幅波非线性相互作用的谱分解法研究》", 《哈尔滨工程大学学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103576574A (en) * 2013-10-10 2014-02-12 哈尔滨工程大学 Phase control method of acoustic energy conversion under nonlinear mutual action of three acoustic wavetrains
CN103576574B (en) * 2013-10-10 2015-12-23 哈尔滨工程大学 The phase control method of acoustic energy conversion under a kind of three row sound wave nonlinear interactions
CN105301582A (en) * 2015-10-14 2016-02-03 哈尔滨工程大学 Underwater acoustic field weak signal enhancement method
CN105301582B (en) * 2015-10-14 2017-10-31 哈尔滨工程大学 A kind of weak signal enhancement method in sound field under water
CN106841382A (en) * 2017-01-23 2017-06-13 哈尔滨工程大学 Non-uniform mixing medium nonlinear factor measuring method based on three couple waves interaction principle
CN106841382B (en) * 2017-01-23 2019-06-14 哈尔滨工程大学 Based on three couple waves interaction non-uniform mixing medium nonlinear factor measurement method
CN112954562A (en) * 2021-01-27 2021-06-11 哈尔滨工程大学 Acoustic signal enhancer based on parameter excitation

Also Published As

Publication number Publication date
CN102510548B (en) 2013-12-18

Similar Documents

Publication Publication Date Title
Agnon et al. A new approach to high-order Boussinesq models
Shen et al. Low-complexity sinusoidal-assisted EMD (SAEMD) algorithms for solving mode-mixing problems in HHT
Hendy et al. A numerically efficient and conservative model for a Riesz space-fractional Klein–Gordon–Zakharov system
Colosi et al. Statistics of low-frequency normal-mode amplitudes in an ocean with random sound-speed perturbations: Shallow-water environments
CN102510548B (en) Output adjustment method for sound energy obtained by sound wave interaction in nonlinear medium
Colosi et al. Observations of sound-speed fluctuations in the western Philippine Sea in the spring of 2009
Helsing et al. Variants of an explicit kernel-split panel-based Nyström discretization scheme for Helmholtz boundary value problems
Tariq et al. Computational soliton solutions to (3+ 1)-dimensional generalised Kadomtsev–Petviashvili and (2+ 1)-dimensional Gardner–Kadomtsev–Petviashvili models and their applications
Niu et al. Theoretical analysis of warping operators for non-ideal shallow water waveguides
Raftari et al. Analytical solution of second-order hyperbolic telegraph equation by variational iteration and homotopy perturbation methods
CN105353408A (en) Wigner higher-order spectrum seismic signal spectral decomposition method based on matching pursuit
Bhaumik et al. Impulse response and spatio-temporal wave-packets: The common feature of rogue waves, tsunami, and transition to turbulence
Colosi et al. Coupled mode transport theory for sound transmission through an ocean with random sound speed perturbations: Coherence in deep water environments
Whitfield et al. Rotation-induced nonlinear wavepackets in internal waves
Xu et al. Lie symmetry analysis, conservation laws, solitary and periodic waves for a coupled Burger equation
Kim et al. Energy flow model considering near field energy for predictions of acoustic energy in low damping medium
CN102693342A (en) Parameter selecting method for restraining sound wave energy in strong nonlinear medium
Akbari-Moghanjoughi Electrostatic rogue-waves in relativistically degenerate plasmas
CN103200491B (en) A kind of based on bubble motion chaos anti-control method in the sound field of parameter perturbation
Gandarias et al. Conservation laws and symmetries of a generalized Kawahara equation
BEREKETOĞLU et al. Behavior of the solutions of a partial differential equation with a piecewise constant argument
Zhang et al. General third-order solutions for irregular waves in deep water
Gao Numerical oscillation and non-oscillation for differential equation with piecewise continuous arguments of mixed type
Diaz et al. Stabilized perfectly matched layer for advective acoustics
Zheng Singularity formation for the compressible Euler equations with general pressure law

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131218

Termination date: 20191018

CF01 Termination of patent right due to non-payment of annual fee