CN117375577B - Numerical filtering method and device for sound propagation problem, electronic equipment and storage medium - Google Patents
Numerical filtering method and device for sound propagation problem, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN117375577B CN117375577B CN202311659770.2A CN202311659770A CN117375577B CN 117375577 B CN117375577 B CN 117375577B CN 202311659770 A CN202311659770 A CN 202311659770A CN 117375577 B CN117375577 B CN 117375577B
- Authority
- CN
- China
- Prior art keywords
- filtering
- wave number
- sound field
- acoustic
- local
- 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.)
- Active
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 178
- 238000000034 method Methods 0.000 title claims abstract description 148
- 230000008569 process Effects 0.000 claims abstract description 38
- 238000012545 processing Methods 0.000 claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 238000004590 computer program Methods 0.000 claims description 9
- 230000036961 partial effect Effects 0.000 claims description 5
- 238000012546 transfer Methods 0.000 abstract description 11
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 12
- 230000002829 reductive effect Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000010355 oscillation Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000013473 artificial intelligence Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0248—Filters characterised by a particular frequency response or filtering method
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0223—Computation saving measures; Accelerating measures
- H03H17/0238—Measures concerning the arithmetic used
Abstract
The application discloses a numerical filtering method, a numerical filtering device, electronic equipment and a storage medium for an acoustic propagation problem, which are applied to the technical field of acoustic propagation. Generating grids of the acoustic wave field to be processed according to the initial value information of the acoustic field of the acoustic wave field to be processed; determining local wave number information of the current grid point according to the sound field quantity on the adjacent grid point to serve as a nonlinear multiplier of a dissipation term in the filtering process; and performing filtering processing on the acoustic wave field to be processed once based on the local wave number information according to a preset time format every time, until reaching a preset calculation time, and obtaining sound field information changing with time. The method and the device can solve the defect that the filtering of the related technology uses a wide template, can improve the credibility of a filtering method adopted in the acoustic propagation process, and reduce the difficulty of boundary processing and the data transfer quantity among parallel computing blocks.
Description
Technical Field
The present invention relates to the field of sound propagation technologies, and in particular, to a method and apparatus for filtering a numerical value of a sound propagation problem, an electronic device, and a readable storage medium.
Background
Acoustic propagation problems are typically calculated in a numerical simulation process using a center differential format without dissipation. The format has no dissipation, but problems such as numerical dispersion, grid unevenness, boundary processing and the like in numerical simulation can cause that numerical oscillation can not be effectively inhibited or eliminated, and the accuracy of a numerical simulation result is seriously affected. In order to improve the accuracy of the numerical simulation result, the related technology adds filtering processing to the acoustic propagation numerical simulation to eliminate the numerical oscillation, so that the final numerical solution is more accurate. When the central differential format is adopted to carry out the numerical simulation of the acoustic propagation problem, in the numerical simulation step, a Shapiro explicit filtering method is generally adopted, the current sound field needs to be filtered every time step or every time advance of a plurality of time steps, and dissipation items of the method have larger dissipation in a high wave number region and small dissipation in a low wave number region, so that high-frequency oscillation can be well restrained, and a better filtering effect is achieved.
However, the filtering template adopted by the numerical simulation explicit filtering method in the related art is wide, the use of the wide template can reduce the reliability of the filtering result, and the increase of the difficulty of boundary processing and the increase of the data transfer quantity between blocks in parallel computing are caused.
In view of this, reducing the width of the filtering method template is a technical problem that needs to be solved by those skilled in the art.
Disclosure of Invention
The application provides a numerical filtering method, a device, electronic equipment and a readable storage medium for an acoustic propagation problem, which can reduce the width of a filtering method template, improve the credibility of a filtering method adopted in an acoustic propagation process, and reduce the difficulty of boundary processing and the data transfer quantity among parallel computing blocks.
In order to solve the technical problems, the application provides the following technical scheme:
in one aspect, the present application provides a method for filtering a numerical value of an acoustic propagation problem, including:
generating grids of the acoustic wave field to be processed according to the initial value information of the acoustic field to be processed;
determining local wave number information of the current grid point according to the sound field quantity on the adjacent grid point to serve as a nonlinear multiplier of a dissipation term in the filtering process;
and performing filtering processing on the acoustic wave field to be processed once based on the local wave number information according to a preset time format every time, until reaching a preset calculation time, and obtaining sound field information changing with time.
Illustratively, the determining the local wave number information of the current lattice point according to the sound field quantity on the adjacent lattice point includes:
calculating a local wave number indicator according to sound field quantity on the adjacent network point;
and taking the smaller value of the local wave number indicator and the preset value as local wave number information.
Illustratively, the selecting the smaller value of the local wave number indicator and the preset value as the local wave number information includes:
and selecting a smaller value between the local wave number indicator and the numerical value 1 as the local wave number information.
Illustratively, the filtering the acoustic wave field to be processed once based on the local wave number information includes:
and sequentially carrying out one-time filtering processing on the acoustic wave field to be processed based on the local wave number information for each dimension direction.
Illustratively, the determining the local wave number information of the current lattice point according to the sound field quantity on the adjacent lattice point includes:
invoking a local wave number indicator calculation relation, and calculating a local wave number indicator of each grid point; the partial wave number indicator calculation relation is as follows:
；
wherein,α j is the firstjThe local wavenumber indicators of the individual grid points,δin order to be an operator,，/>is the firstjGrid points, thjSound field variable values at midpoints of +1 lattice points, +.>Is the firstjGrid points, thj-sound field variable values at the midpoints of 1 grid points, Δ being an operator, Δf j =f j+1 -f j-1 ，f j Is the firstjThe amount of sound field at each grid point,f j+1 is the firstjSound field quantity of +1 grid points,f j-1 is the firstjThe sound field quantity of 1 grid point,εis constant.
Illustratively, the filtering the acoustic wave field to be processed once based on the local wave number information includes:
invoking a first filtering relation, and performing primary filtering treatment on the acoustic wave field to be treated; the first filtering relation is:
;
in the method, in the process of the invention,is the firstjThe amount of sound field of the individual grid points after filtering,α j is the firstjThe local wavenumber indicators of the individual grid points,f j is the firstjSound field quantity of each grid point.
Illustratively, the filtering the acoustic wave field to be processed once based on the local wave number information includes:
invoking a second filtering relation, and performing primary filtering treatment on the acoustic wave field to be treated; the second filtering relation is:
；
in the method, in the process of the invention,is the firstjThe amount of sound field of the individual grid points after filtering,α j is the firstjThe local wavenumber indicators of the individual grid points,f j is the firstjSound field quantity of each grid point.
Another aspect of the present application provides a numerical filtering apparatus for an acoustic propagation problem, comprising:
the grid generation module is used for generating grids of the acoustic wave field to be processed according to the initial value information of the acoustic field of the acoustic wave field to be processed;
the local wave number information determining module is used for determining the local wave number information of the current grid point according to the sound field quantity on the adjacent grid point to serve as a nonlinear multiplier of a dissipation term in the filtering process;
and the filtering module is used for carrying out filtering processing on the acoustic wave field to be processed once based on the local wave number information according to a preset time format every time, and obtaining sound field information changing with time until reaching a preset calculation time.
The present application also provides an electronic device comprising a processor for implementing the steps of the numerical filtering method of the acoustic propagation problem of any one of the preceding claims when executing a computer program stored in a memory.
The present application finally provides a readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the method for numerical filtering of acoustic propagation problems as described in any of the preceding claims.
The technical scheme provided by the application has the advantages that the local wave number information of the current grid point is determined according to the sound field quantity on the adjacent grid point and is used as a nonlinear multiplier of a dissipation term in the filtering process, so that the convergence precision of the filtering method can be effectively improved by adding a nonlinear coefficient to the dissipation term of the original filtering method, the number of filtering template points is reduced under the requirement of the same precision, the requirement on the width of the template is lower, namely, compared with the prior art, the template with smaller width can be used for meeting the requirement on the same filtering precision, the purposes of small dissipation in a low wave number area to keep sound field effective information and larger dissipation in a high wave number area are achieved, the high-frequency numerical oscillation is effectively restrained, the reliability of the filtering method is improved, and the boundary processing difficulty and the data transfer quantity among parallel computing blocks are effectively reduced.
In addition, the application provides a corresponding implementation device, electronic equipment and a readable storage medium for the numerical filtering method of the sound propagation problem, so that the method is more practical, and the device, the electronic equipment and the readable storage medium have corresponding advantages.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the application.
Drawings
For a clearer description of the technical solutions of the present application or of the related art, the drawings that are required to be used in the description of the embodiments or of the related art will be briefly described, and it is obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive effort for a person of ordinary skill in the art.
Fig. 1 is a flow chart of a numerical filtering method for acoustic propagation problem provided in the present application;
FIG. 2 is a graph of sound pressure distribution of vortex sound scattering problem numerical simulation results for an exemplary application scenario provided herein;
FIG. 3 is a block diagram of one embodiment of a numerical filter device for acoustic propagation problems provided herein;
fig. 4 is a block diagram of an embodiment of an electronic device provided in the present application.
Detailed Description
In order to provide a better understanding of the present application, those skilled in the art will now make further details of the present application with reference to the drawings and detailed description. It will be apparent that the described embodiments are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
The terms first, second and the like in the description and in the claims of the present application and in the above-described figures, are used for distinguishing between different objects and not necessarily for describing a particular sequential or chronological order. Furthermore, the terms "comprise" and "have," as well as any variations of the two, are intended to cover a non-exclusive inclusion. The term "exemplary" means "serving as an example, embodiment, or illustration. Any embodiment described herein as "exemplary" is not necessarily to be construed as preferred or advantageous over other embodiments. Various non-limiting embodiments of the present application are described in detail below. Numerous specific details are set forth in the following description in order to provide a better understanding of the invention. It will be understood by those skilled in the art that the present application may be practiced without these specific details. In other instances, well-known methods, procedures, components, and circuits have not been described in detail so as not to obscure the present application.
Referring first to fig. 1, fig. 1 is a flow chart of a numerical filtering method for an acoustic propagation problem provided in the present application, where the present application may include the following:
s101: and generating grids of the acoustic wave field to be processed according to the initial value information of the acoustic field to be processed.
In this step, the acoustic wave field to be processed is the acoustic wave field of the current research acoustic propagation problem, for example, in the application scenario of the acoustic wave field when the acoustic wave passes through the vortex, the acoustic wave field to be processed is the acoustic wave field when the acoustic wave passes through the vortex. It will be appreciated that in the sound propagation problem, sound field initial value information is typically given, and then time-varying sound field information needs to be acquired. The problems are unsteady problems, and a numerical simulation process of the problems firstly sets a calculation area to generate grids so as to discretize an initial value of a sound field. The sound field initial value information in this step includes at least grid parameters including, but not limited to, grid spacing and time step, and a calculation region. The user issues the self-defined initial value information of the sound field, a calculation area is determined according to the received initial value information of the sound field, and a uniform grid is generated according to grid parameters, wherein the calculation area can be [ -20, 20] × [ -20, 20], the grid interval and the time step can be respectively 0.1 and 0.05, and the uniform grid of the acoustic wave field to be processed can be respectively generated according to the grid interval and the time step.
S102: and determining the local wave number information of the current grid point according to the sound field quantity on the adjacent grid point to serve as a nonlinear multiplier of a dissipation term in the filtering process.
In this step, the sound field amount on the grid point includes, but is not limited to, density pulsation and pressure pulsation, and the adjacent grid point of the current grid point may include a plurality of grid points above, below, left and right, which can be flexibly determined by those skilled in the art according to actual needs, and the present application is not limited in any way. The local wave number information is used as a nonlinear multiplier of a dissipation term in the filtering process, and the value is related to the value of the sound field quantity of the nearby grid points, so that the method belongs to a nonlinear method, the convergence precision of the filtering method can be improved, and the number of points required by the filtering method suitable for the numerical simulation method with the same precision can be smaller through the local wave number information, namely, a large filtering template is not needed. It is generally believed that the filtering method requires a higher two-order precision than the spatial value dispersion, such as when using a center difference method with 7-point 6-order precision for simulation, the filtering method is required to have 8-order precision, while the linear 8-order explicit filtering format requires that the template contains at least 9 points. An excessively wide template reduces the reliability of the filtering result and leads to an increase in the difficulty of boundary processing and an increase in the amount of data transfer between blocks in parallel computation. The nonlinear multiplier of the local wave number information is added to the dissipation term in the filtering process, the 6-order explicit filtering method can be adapted to a 7-point 6-order precision center difference method, obviously, compared with an 8-order explicit filtering method, the 6-order explicit filtering method can reduce templates, further improve the reliability of the filtering method, and reduce the boundary processing difficulty and the data transfer quantity among parallel computing blocks.
S103: and performing filtering processing on the acoustic wave field to be processed once based on the local wave number information according to the preset time format every time, until the preset calculation time is reached, and obtaining the sound field information changing with time.
At step S101, discretizing an initial sound field value of the acoustic wave field to be processed, determining a sound field equation and a boundary condition of the acoustic wave field to be processed according to an application scene of the acoustic wave field to be processed after determining a dissipation term adopted in a filtering process by step S102, performing discretization processing on the sound field equation and the boundary condition, and performing time propulsion by combining a proper time format, namely a preset time format, wherein the preset time format can be one step or a plurality of steps, and accordingly, performing one filtering processing on the sound field every time one step or a plurality of steps is performed. When the filtering processing is performed, any filtering method may be used, for example, a shape filtering method may be used in the embodiment, and the local wave number information selected in S102 needs to be combined in the filtering process, and when the preset computing time is reached, the computing is terminated, so that the sound field information of the sound wave field to be processed changing along with the time can be obtained.
In the technical scheme provided by the application, the local wave number information of the current grid point is determined according to the sound field quantity on the adjacent grid point and is used as a nonlinear multiplier of a dissipation term in the filtering process, so that the convergence precision of the filtering method can be effectively improved by adding a nonlinear coefficient to the dissipation term of the original filtering method, the number of filtering template points is reduced under the requirement of the same precision, the requirement on the template width is lower, namely, compared with the prior art, the template with smaller width can be used for meeting the requirement on the same filtering precision, the purposes of not only realizing small dissipation in a low wave number region to keep sound field effective information and larger dissipation in a high wave number region are achieved, thereby effectively inhibiting high-frequency numerical oscillation, but also improving the reliability of the filtering method and effectively reducing the difficulty of boundary processing and the data transfer quantity among parallel computing blocks.
It should be noted that, in the present application, the steps may be executed simultaneously or in a certain preset order as long as the steps conform to the logic order, and fig. 1 is only a schematic manner and does not represent only such an execution order.
In the above embodiment, how to determine the local wave number information is not limited, and this embodiment provides an efficient and simple calculation method of the local wave number information, which may include the following:
invoking a local wave number indicator calculation relation, and calculating a local wave number indicator of each grid point; the local wavenumber indicator calculation relationship can be expressed as:
；
wherein,α j is the firstjThe local wavenumber indicators of the individual grid points,δin order to be an operator,，/>is the firstjGrid points, thjSound field variable values at midpoints of +1 lattice points, +.>Is the firstjGrid points, thj-sound field variable values at the midpoints of 1 grid points, Δ being an operator, Δf j =f j+1 -f j-1 ，f j Is the firstjThe amount of sound field at each grid point,f j+1 is the firstjSound field quantity of +1 grid points,f j-1 is the firstjThe sound field quantity of 1 grid point,εis a constant, which is a small positive number, and the denominator for placing the above-mentioned local wave number indicator is 0, for example, 10 -40 Accordingly, the above-described partial wave number indicator calculation relation may be expressed as:
。
in this embodiment, accuracy analysis can be performed based on the above-described partial wave number indicator calculation relation, if the grid pitch ishThenV 6 =O(h 6 ). Wherein,O(h m ) Indicating whenhToward 0, the amount is less thanh m Is a constant multiple of (a). At a point where the first order derivative is not 0,α j =O(h 2 ) The accuracy of the filtering method is 8 th order at this time; whereas at the level one lead of 0,α j =O(h) The accuracy of the filtering method is 7 th order. Since the first order 0 points are only individual points, the filtering method is finallyL 1 The error accuracy is 8 th order. Wherein,L 1 the error is a value obtained by multiplying the absolute value of each grid point error by the grid pitch and summing the multiplied absolute value. Therefore, the filtering method can be used in combination with a 7-point 6-order precision center difference method without a larger template, so that the width of the adopted template can be reduced, the reliability of the filtering method is improved, and the boundary processing difficulty and the data transfer quantity among parallel computing blocks are effectively reduced.
It can be appreciated that the value of the local wave number indicator may be relatively large, which results in that the dissipation term in the filtering process is larger than the original dissipation term, so that the width of the finally used filtering template cannot be reduced, based on this, in order to effectively and accurately implement the reduction of the width of the filtering template to improve the reliability of the filtering method, based on the foregoing embodiment, the present application further provides an efficient and simple local wave number information determining manner, which may include the following:
calculating a local wave number indicator according to sound field quantity on the adjacent network point; and taking the smaller value of the local wave number indicator and the preset value as the local wave number information.
The preset value is used to ensure that the dissipation term adopted in the filtering process of the application is not greater than the original dissipation term, and accordingly, after the local wave number indicator is obtained by calculation in any embodiment, a further selection is needed, that is, a smaller value between the local wave number indicator and the preset value is selected. For example, as a simple embodiment, the preset value may be 1, and correspondingly, a smaller value may be selected between the local wave number indicator and the numerical value 1 as the local wave number information. Of course, the preset value may be 1, which fluctuates within any small fluctuation range, such as 0.9, which does not affect the implementation of the present application.
Taking the filtering method as an example, shapiro filtering is adopted, and the wave number iskOf (2), i.e. the frequency in the spatial direction iskThe grid spacing ishCalculated to obtainα j =4sin 4 (khAnd 2), the formula is smaller in the low wave number area, and larger in the high wave number area, so that the filtering requirement is well met. By means of oppositeα j And the filter effect of the filter method on the high-frequency wave is equivalent to that of the original shape filter method by taking a smaller value with a preset value 1.
It is to be understood that the filtering in the above embodiment may be one-dimensional filtering or multidimensional filtering, and for the one-dimensional filtering, the filtering processing may be directly performed according to the manner described in the above embodiment, and for the multidimensional filtering, the filtering may be sequentially performed on each dimension direction according to the above method, that is, for each dimension direction sequentially, the filtering processing may be performed on the acoustic wave field to be processed once based on the local wave number information.
The above embodiment does not limit the filtering method, and the present embodiment also provides a plurality of nonlinear filtering methods, which may include the following:
the nonlinear filtering method adopted in the embodiment can be obtained by adding a nonlinear multiplier to a dissipation term of a 6-order shape filtering method, and can be used for carrying out primary filtering processing on an acoustic wave field to be processed by calling a first filtering relation; the first filter relation may be expressed as:
；
in the method, in the process of the invention,is the firstjThe amount of sound field of the individual grid points after filtering,α j is the firstjThe local wavenumber indicators of the individual grid points,f j is the firstjSound field quantity of each grid point.
The nonlinear filtering method adopted in the embodiment can be obtained by adding a nonlinear multiplier to a dissipation term of an 8-order shape filtering method, and a second filtering relation can be called to perform primary filtering treatment on an acoustic wave field to be treated; the second filter relation may be expressed as:
；
in the method, in the process of the invention,is the firstjThe amount of sound field of the individual grid points after filtering,α j is the firstjThe local wavenumber indicators of the individual grid points,f j is the firstjSound field quantity of each grid point.
The embodiment can adopt a 7-point filtering method, and can also be popularized to the conditions of 9 points or even higher odd points, which do not influence the realization of the application, and the technical personnel in the field can flexibly select according to the actual conditions.
To make the implementation of the present application more clear to those skilled in the art, the present application further provides an exemplary implementation manner, where the acoustic wave field to be processed in this embodiment is an acoustic wave field when an acoustic wave passes through a vortex, and may include the following:
when the sound wave passes through the vortex, strong nonlinear scattering phenomenon occurs, and the amplitude and the phase of the sound wave are changed obviously. The present embodiment uses an enhanced optimized differential format in combination with the first filtering relation to calculate the scattering problem of the plane acoustic wave passing through the isentropic vortex. The problem is dimensionless at infinity speed of sound and density. The vortex center of the two-dimensional isentropic Taylor vortex in the background flow is positioned at the center (0, 0) of the calculation region, and the flow parameters can be as follows:
；
wherein,ρin order to achieve a density of the particles,u θ for the background flow tangential velocity,u r in order to be a radial velocity,pis the pressure. Maximum Mach numberM v =0.125，γFor the gas insulation index, the air content is generally 1.4. The left boundary incident plane wave form is:
；
wherein,ρ＇、u＇、ν＇、p' is the density of the incident sound wave, the velocity in the x direction, the velocity in the y direction, and the pressure, respectively. Sound pressure amplitude ε' =10 -5 The wavelength lambda takes 1.
In this embodiment, the initial value information of the sound field is calculated by the calculation area of [ -20, 20] × [ -20, 20], the grid interval and the time step are respectively 0.1 and 0.05, in fig. 2, the sound pressure distribution of the calculation result of the preset calculation time is shown by the x-axis, the speed of the incident sound wave in the x-direction is shown by the y-axis, the speed of the incident sound wave in the y-direction is shown by the y-axis, as can be seen from fig. 2, the numerical result better reflects the scattering condition of the vortex on the sound wave, and the sound wave can be obviously distorted after passing through the vortex. There is a significant phase shift just behind the vortex, the sound wave amplitude decays rapidly, several noise interference bands are formed behind the vortex, and the up and down asymmetry is achieved.
From the above, in this embodiment, by adding a nonlinear coefficient to the dissipation term of the shape filtering method, the convergence accuracy of the filtering method is improved, so that the requirement on the width of the template is reduced, and the adverse effect caused by a large template is solved. The local wave number indicator and 1 are smaller, so that the filtering effect of the filtering method on the high-frequency wave is equivalent to that of the original shape filtering method, the credibility of the filtering method adopted in the process that the simulated sound wave passes through the vortex can be improved, and the boundary processing difficulty and the data transfer quantity among parallel computing blocks can be reduced.
The application also provides a corresponding device for the numerical filtering method of the sound propagation problem, so that the method is more practical. Wherein the device may be described separately from the functional module and the hardware. In the following description, a numerical filtering apparatus for an acoustic propagation problem provided in the present application is described, which is used to implement the numerical filtering method for an acoustic propagation problem provided in the present application, and in this embodiment, the numerical filtering apparatus for an acoustic propagation problem may include or be divided into one or more program modules, where the one or more program modules are stored in a storage medium and executed by one or more processors, to implement the numerical filtering method for an acoustic propagation problem disclosed in the first embodiment. Program modules in the present application refer to a series of computer program instruction segments capable of performing a specified function, more suitable than the program itself for executing a numerical filtering device describing sound propagation problems in a storage medium. The following description will specifically describe functions of the program modules of the present embodiment, and a numerical filtering apparatus for an acoustic propagation problem described below and a numerical filtering method for an acoustic propagation problem described above may be referred to correspondingly to each other.
Based on the angle of the functional modules, referring to fig. 3, fig. 3 is a block diagram of a numerical filtering device for acoustic propagation problem provided in the present application under a specific embodiment, where the device may include:
the grid generating module 301 is configured to generate a grid of the acoustic wave field to be processed according to the acoustic field initial value information of the acoustic wave field to be processed.
The local wave number information determining module 302 is configured to determine local wave number information of a current grid point according to the sound field amounts on adjacent grid points, so as to be a nonlinear multiplier of a dissipation term in the filtering process.
The filtering module 303 is configured to perform a filtering process on the acoustic wave field to be processed based on the local wave number information once per time of time advancing according to a preset time format until reaching a preset calculation time, and obtain sound field information that varies with time.
Optionally, in some implementations of this embodiment, the local wave number information determining module 302 may be further configured to: calculating a local wave number indicator according to sound field quantity on the adjacent network point; and taking the smaller value of the local wave number indicator and the preset value as the local wave number information.
As an exemplary implementation of the foregoing embodiment, the foregoing local wave number information determining module 302 may be further configured to: a smaller value is selected between the local wave number indicator and the value 1 as the local wave number information.
Optionally, in other implementations of this embodiment, the filtering module 303 may be further configured to: and sequentially carrying out primary filtering processing on the acoustic wave field to be processed based on the local wave number information in each dimension direction.
Illustratively, in some implementations of this embodiment, the local wavenumber information determining module 302 may be further configured to:
invoking a local wave number indicator calculation relation, and calculating a local wave number indicator of each grid point; the local wavenumber indicator calculation relationship is:
；
wherein,α j is the firstjThe local wavenumber indicators of the individual grid points,δin order to be an operator,，/>is the firstjGrid points, thjSound field variable values at midpoints of +1 lattice points, +.>Is the firstjGrid points, thj-sound field variable values at the midpoints of 1 grid points, Δ being an operator, Δf j =f j+1 -f j-1 ，f j Is the firstjThe amount of sound field at each grid point,f j+1 is the firstjSound field quantity of +1 grid points,f j-1 is the firstjThe sound field quantity of 1 grid point,εis constant.
Illustratively, in other implementations of the present embodiment, the filtering module 303 may further be configured to:
invoking a first filtering relation, and performing primary filtering treatment on the acoustic wave field to be treated; the first filtering relation is:
;
in the method, in the process of the invention,is the firstjThe amount of sound field of the individual grid points after filtering,α j is the firstjThe local wavenumber indicators of the individual grid points,f j is the firstjSound field quantity of each grid point.
Illustratively, in still other implementations of this embodiment, the filtering module 303 may be further configured to: invoking a second filtering relation, and performing primary filtering treatment on the acoustic wave field to be treated; the second filtering relation is:
；
in the method, in the process of the invention,is the firstjThe amount of sound field of the individual grid points after filtering,α j is the firstjThe local wavenumber indicators of the individual grid points,f j is the firstjSound field quantity of each grid point.
The functions of each functional module of the numerical filter device for the acoustic propagation problem can be specifically implemented according to the method in the above method embodiment, and the specific implementation process of the functional module can refer to the related description of the above method embodiment, which is not repeated herein.
As can be seen from the above, the present embodiment can reduce the width of the filtering method template, not only can improve the reliability of the filtering method adopted in the acoustic propagation process, but also can reduce the difficulty of boundary processing and the data transfer amount between parallel computing blocks.
The above-mentioned numerical filtering device for the acoustic propagation problem is described from the viewpoint of a functional module, and further, the present application also provides an electronic device, which is described from the viewpoint of hardware. Fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present application in an implementation manner. As shown in fig. 4, the electronic device comprises a memory 40 for storing a computer program; a processor 41 for implementing the steps of the numerical filtering method of the acoustic propagation problem as mentioned in any of the embodiments above when executing a computer program.
Processor 41 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and processor 41 may also be a controller, microcontroller, microprocessor, or other data processing chip, among others. The processor 41 may be implemented in at least one hardware form of DSP (Digital Signal Processing ), FPGA (Field-Programmable Gate Array, field programmable gate array), PLA (Programmable Logic Array ). The processor 41 may also comprise a main processor, which is a processor for processing data in an awake state, also called CPU (Central Processing Unit ); a coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 41 may be integrated with a GPU (Graphics Processing Unit, image processor) for taking care of rendering and drawing of the content that the display screen is required to display. In some embodiments, the processor 41 may also include an AI (Artificial Intelligence ) processor for processing computing operations related to machine learning.
Memory 40 may include one or more computer-readable storage media, which may be non-transitory. Memory 40 may also include high-speed random access memory as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. Memory 40 may be an internal storage unit of the electronic device, such as a hard disk of a server, in some embodiments. The memory 40 may in other embodiments also be an external storage device of the electronic device, such as a plug-in hard disk provided on a server, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash Card (Flash Card) or the like. Further, the memory 40 may also include both internal storage units and external storage devices of the electronic device. The memory 40 may be used to store not only application software installed in the electronic device, but also various types of data, such as: code of a program or the like during the numerical filtering method that performs the acoustic propagation problem can also be used to temporarily store data that has been output or is to be output.
In this embodiment, the memory 40 is at least used for storing a computer program 401, which, when loaded and executed by the processor 41, is capable of implementing the relevant steps of the numerical filtering method for the acoustic propagation problem disclosed in any of the foregoing embodiments. In addition, the resources stored in the memory 40 may further include an operating system 402, data 403, and the like, where the storage manner may be transient storage or permanent storage. Operating system 402 may include, among other things, windows, unix, linux. The data 403 may include, but is not limited to, data corresponding to the numerical filtering results of the acoustic propagation problem, and the like.
In some embodiments, the electronic device may further include a display 42, an input/output interface 43, a communication interface 44, or referred to as a network interface, a power supply 45, and a communication bus 46. Among other things, the display 42, input-output interface 43 such as a Keyboard (Keyboard) belong to a user interface, which may alternatively include a standard wired interface, a wireless interface, etc. Alternatively, in some embodiments, the display may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, an OLED (Organic Light-Emitting Diode) touch, or the like. The display may also be referred to as a display screen or display unit, as appropriate, for displaying information processed in the electronic device and for displaying a visual user interface. The communication interface 44 may optionally include a wired interface and/or a wireless interface, such as a WI-FI interface, a bluetooth interface, etc., typically used to establish a communication connection between the electronic device and other electronic devices. The communication bus 46 may be a peripheral component interconnect standard (peripheral component interconnect, PCI) bus, or an extended industry standard architecture (extended industry standard architecture, EISA) bus, among others. The bus may be classified as an address bus, a data bus, a control bus, etc. For ease of illustration, only one thick line is shown in fig. 4, but not only one bus or one type of bus.
Those skilled in the art will appreciate that the configuration shown in fig. 4 is not limiting of the electronic device and may include more or fewer components than shown, for example, may also include sensors 47 to perform various functions.
The functions of each functional module of the electronic device may be specifically implemented according to the method in the above method embodiment, and the specific implementation process may refer to the related description of the above method embodiment, which is not repeated herein.
As can be seen from the above, the present embodiment can reduce the width of the filtering method template, not only can improve the reliability of the filtering method adopted in the acoustic propagation process, but also can reduce the difficulty of boundary processing and the data transfer amount between parallel computing blocks.
It will be appreciated that the numerical filtering method of the acoustic propagation problem in the above embodiments may be stored in a computer readable storage medium if implemented in the form of a software functional unit and sold or used as a stand alone product. Based on such understanding, the technical solution of the present application, or a part contributing to the related art, or all or part of the technical solution, may be embodied in the form of a software product stored in a storage medium, performing all or part of the steps of the methods of the embodiments of the present application. And the aforementioned storage medium includes: a U disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), an electrically erasable programmable ROM, registers, a hard disk, a multimedia card, a card-type Memory (e.g., SD or DX Memory, etc.), a magnetic Memory, a removable disk, a CD-ROM, a magnetic disk, or an optical disk, etc., that can store program code.
Based on this, the present application further provides a readable storage medium storing a computer program which, when executed by a processor, performs the steps of the numerical filtering method of the acoustic propagation problem of any of the above embodiments.
The functions of each functional module of the readable storage medium may be specifically implemented according to the method in the above method embodiment, and the specific implementation process may refer to the relevant description of the above method embodiment, which is not repeated herein.
In this specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, so that the same or similar parts between the embodiments are referred to each other. For the hardware including the device and the electronic equipment disclosed in the embodiments, the description is relatively simple because the hardware includes the device and the electronic equipment corresponding to the method disclosed in the embodiments, and relevant places refer to the description of the method.
Those of skill would further appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative elements and steps are described above generally in terms of functionality in order to clearly illustrate the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
The above describes in detail a method, an apparatus, an electronic device and a readable storage medium for filtering a sound propagation problem provided in the present application. Specific examples are set forth herein to illustrate the principles and embodiments of the present application, and the description of the examples above is only intended to assist in understanding the methods of the present application and their core ideas. It should be noted that it would be obvious to those skilled in the art that various improvements and modifications can be made to the present application without departing from the principles of the present application, and such improvements and modifications fall within the scope of the claims of the present application.
Claims (9)
1. A method of numerical filtering of an acoustic propagation problem, comprising:
generating grids of the acoustic wave field to be processed according to the initial value information of the acoustic field to be processed;
determining local wave number information of the current grid point according to the sound field quantity on the adjacent grid point to serve as a nonlinear multiplier of a dissipation term in the filtering process;
performing a filtering process on the acoustic wave field to be processed based on the local wave number information every time according to a preset time format until reaching a preset calculation time to obtain sound field information changing with time;
wherein the determining the local wave number information of the current grid point according to the sound field quantity on the adjacent grid point comprises the following steps:
invoking a local wave number indicator calculation relation, and calculating a local wave number indicator of each grid point; the partial wave number indicator calculation relation is as follows:
；
wherein,α j is the firstjThe local wavenumber indicators of the individual grid points,δin order to be an operator,，/>is the firstjGrid points, thjSound field variable values at midpoints of +1 lattice points, +.>Is the firstjGrid points, thj-sound field variable values at the midpoints of 1 grid points, Δ being an operator, Δf j = f j+1 - f j-1 ，f j Is the firstjThe amount of sound field at each grid point,f j+1 is the firstjSound field quantity of +1 grid points,f j-1 is the firstjThe sound field quantity of 1 grid point,εis constant.
2. The method of numerical filtering of an acoustic propagation problem according to claim 1, wherein the determining the local wave number information of the current lattice point from the sound field amounts on the neighboring lattice points includes:
calculating a local wave number indicator according to sound field quantity on the adjacent network point;
and taking the smaller value of the local wave number indicator and the preset value as local wave number information.
3. The method of numerical filtering of an acoustic propagation problem according to claim 2, wherein said using the smaller of the local wave number indicator and a preset value as the local wave number information includes:
and selecting a smaller value between the local wave number indicator and the numerical value 1 as the local wave number information.
4. The numerical filtering method of the acoustic propagation problem according to claim 1, wherein the performing a filtering process on the acoustic wave field to be processed based on the local wave number information includes:
and sequentially carrying out one-time filtering processing on the acoustic wave field to be processed based on the local wave number information for each dimension direction.
5. The numerical filtering method of the acoustic propagation problem according to any one of claims 1 to 4, wherein the performing a filtering process on the acoustic wave field to be processed based on the local wave number information includes:
invoking a first filtering relation, and performing primary filtering treatment on the acoustic wave field to be treated; the first filtering relation is:
;
in the method, in the process of the invention,is the firstjThe amount of sound field of the individual grid points after filtering,α j is the firstjThe local wavenumber indicators of the individual grid points,f j is the firstjSound field quantity of each grid point.
6. The numerical filtering method of the acoustic propagation problem according to any one of claims 1 to 4, wherein the performing a filtering process on the acoustic wave field to be processed based on the local wave number information includes:
invoking a second filtering relation, and performing primary filtering treatment on the acoustic wave field to be treated; the second filtering relation is:
；
in the method, in the process of the invention,is the firstjThe amount of sound field of the individual grid points after filtering,α j is the firstjThe local wavenumber indicators of the individual grid points,f j is the firstjSound field quantity of each grid point.
7. A numerical filter apparatus for acoustic propagation problems, comprising:
the grid generation module is used for generating grids of the acoustic wave field to be processed according to the initial value information of the acoustic field of the acoustic wave field to be processed;
the local wave number information determining module is used for determining the local wave number information of the current grid point according to the sound field quantity on the adjacent grid point to serve as a nonlinear multiplier of a dissipation term in the filtering process;
the filtering module is used for carrying out filtering processing on the acoustic wave field to be processed once based on the local wave number information according to a preset time format every time, and obtaining sound field information changing with time until reaching a preset calculation time;
wherein the local wave number information determining module is further configured to:
invoking a local wave number indicator calculation relation, and calculating a local wave number indicator of each grid point; the partial wave number indicator calculation relation is as follows:
；
wherein,α j is the firstjThe local wavenumber indicators of the individual grid points,δin order to be an operator,，/>is the firstjGrid points, thjSound field variable values at midpoints of +1 lattice points, +.>Is the firstjGrid points, thj-sound field variable values at the midpoints of 1 grid points, Δ being an operator, Δf j = f j+1 - f j-1 ，f j Is the firstjThe amount of sound field at each grid point,f j+1 is the firstjSound field quantity of +1 grid points,f j-1 is the firstjThe sound field quantity of 1 grid point,εis constant.
8. An electronic device comprising a processor and a memory, the processor being adapted to implement the steps of the method for numerical filtering of acoustic propagation problems according to any one of claims 1 to 6 when executing a computer program stored in the memory.
9. A readable storage medium, characterized in that the readable storage medium has stored thereon a computer program which, when executed by a processor, implements the steps of the method for numerical filtering of acoustic propagation problems according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311659770.2A CN117375577B (en) | 2023-12-06 | 2023-12-06 | Numerical filtering method and device for sound propagation problem, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311659770.2A CN117375577B (en) | 2023-12-06 | 2023-12-06 | Numerical filtering method and device for sound propagation problem, electronic equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117375577A CN117375577A (en) | 2024-01-09 |
| CN117375577B true CN117375577B (en) | 2024-03-12 |
Family
ID=89389627
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311659770.2A Active CN117375577B (en) | 2023-12-06 | 2023-12-06 | Numerical filtering method and device for sound propagation problem, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117375577B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102008014575A1 (en) * | 2008-03-13 | 2009-09-17 | Volkswagen Ag | Acoustic sources locating method for use in automobile, involves performing local area-wave number range-transform based on acoustic pressure signals, and deriving directions of arrival of sound using wave number spectrum |
| CN109155885A (en) * | 2016-05-30 | 2019-01-04 | 索尼公司 | Local sound field forms device, local sound field forming method and program |
| KR101919241B1 (en) * | 2017-06-20 | 2019-02-08 | 주식회사 한화 | Apparatus and Method for optimum smoothing filtering to eliminate partial continuous filtering |
| CN113239573A (en) * | 2021-06-05 | 2021-08-10 | 西北工业大学 | Closed space sound field reconstruction method based on grid fluctuation-free modeling |
| CN113758559A (en) * | 2021-09-15 | 2021-12-07 | 哈尔滨工程大学 | Method for measuring shallow sea channel based on scattering sound field separation algorithm |
| CN113987691A (en) * | 2021-11-25 | 2022-01-28 | 中国人民解放军国防科技大学 | High-precision hybrid calculation method, device, equipment and storage medium for shock wave instability |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101702330B1 (en) * | 2010-07-13 | 2017-02-03 | 삼성전자주식회사 | Method and apparatus for simultaneous controlling near and far sound field |
| GB2547940A (en) * | 2016-03-04 | 2017-09-06 | Robertsson Johan | Source separation method |
| US11050591B2 (en) * | 2019-04-15 | 2021-06-29 | Alexei V. Nikitin | Method and apparatus for nonlinear filtering and for secure communications |
| CN114151066B (en) * | 2021-10-09 | 2023-04-25 | 电子科技大学 | Reverse time migration imaging method for acoustic interface of ultrasonic Lamb wave logging well wall |
| CN114895351A (en) * | 2022-04-29 | 2022-08-12 | 中国科学院地质与地球物理研究所 | Method and device for medium modeling for simulating seismic wave propagation at any discontinuous interface |
-
2023
- 2023-12-06 CN CN202311659770.2A patent/CN117375577B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102008014575A1 (en) * | 2008-03-13 | 2009-09-17 | Volkswagen Ag | Acoustic sources locating method for use in automobile, involves performing local area-wave number range-transform based on acoustic pressure signals, and deriving directions of arrival of sound using wave number spectrum |
| CN109155885A (en) * | 2016-05-30 | 2019-01-04 | 索尼公司 | Local sound field forms device, local sound field forming method and program |
| KR101919241B1 (en) * | 2017-06-20 | 2019-02-08 | 주식회사 한화 | Apparatus and Method for optimum smoothing filtering to eliminate partial continuous filtering |
| CN113239573A (en) * | 2021-06-05 | 2021-08-10 | 西北工业大学 | Closed space sound field reconstruction method based on grid fluctuation-free modeling |
| CN113758559A (en) * | 2021-09-15 | 2021-12-07 | 哈尔滨工程大学 | Method for measuring shallow sea channel based on scattering sound field separation algorithm |
| CN113987691A (en) * | 2021-11-25 | 2022-01-28 | 中国人民解放军国防科技大学 | High-precision hybrid calculation method, device, equipment and storage medium for shock wave instability |
Non-Patent Citations (3)
| Title |
|---|
| Direct Computation of Sound Generation by the Symmetric Vortex Merger;Xuliang Liu等;《2021 4th International Conference on Information Communication and Signal Processing (ICICSP)》;213-218 * |
| 基于声场空间分布特征的机械故障诊断方法及其应用研究;鲁文波等;《中国博士学位论文全文数据库信息科技辑》(第10（2012年）期);I136-9 * |
| 旋涡声散射的空间尺度特性数值研究;王益民等;《物理学报》;第70卷(第19期);105-117 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117375577A (en) | 2024-01-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xu et al. | An information-theoretic framework for flow visualization | |
| TWI591490B (en) | Vector computation unit in a neural network processor | |
| CN103390262B (en) | The acquisition methods of weight coefficient of digital filter and device | |
| CN1987928A (en) | Device for real time quick realizing Gauss formwork convolution | |
| JP7435951B2 (en) | Floating point number generation method, apparatus, electronic device, storage medium and computer program for integrated circuit chip verification | |
| CN107038302A (en) | A kind of semi-analytical calculation method in the Eddy Current Nondestructive Testing magnetic field containing columnar defects | |
| CN117375577B (en) | Numerical filtering method and device for sound propagation problem, electronic equipment and storage medium | |
| CN109960841B (en) | Fluid surface tension simulation method, terminal equipment and storage medium | |
| CN108744520B (en) | Method and device for determining placement position of game model and electronic equipment | |
| CN113688552A (en) | Camera module falling simulation method and device, electronic equipment and medium | |
| Soares et al. | An energy-efficient and approximate accelerator design for real-time Canny edge detection | |
| CN114880890B (en) | Far-field boundary condition multi-window processing method, device and medium | |
| EP4357924A1 (en) | Application performance testing method, method and apparatus for establishing performance testing model | |
| US10331837B1 (en) | Device graphics rendering for electronic designs | |
| CN109408566A (en) | A kind of intelligence chart recommended method and device | |
| CN115659766B (en) | Method and device for determining current density vector in multi-turn coil | |
| CN113112488A (en) | Road center line extraction method and device, storage medium and terminal equipment | |
| CN110990963A (en) | Fan interval optimization method and device and computer readable storage medium | |
| CN106874925A (en) | object grouping method, model training method and device | |
| TWI830940B (en) | Information processing circuits and design methods of information processing circuits | |
| CN116881655B (en) | Wind power data processing method, system and storage medium based on scatter diagram linkage | |
| CN116227308B (en) | Numerical simulation method and system for shallow logging natural electric field | |
| CN108229668A (en) | Operation implementation method, device and electronic equipment based on deep learning | |
| JP4855853B2 (en) | Analysis device, computer control method, and model creation program | |
| Singh et al. | Hessian-based robust ray-tracing of implicit surfaces on gpu |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |