DE4030753A1 - Calculating sound dispersion in flowing media - separating surface into partial areas small compared with wavelengths produced by sound pressure of vibrating surface - Google Patents
Calculating sound dispersion in flowing media - separating surface into partial areas small compared with wavelengths produced by sound pressure of vibrating surfaceInfo
- Publication number
- DE4030753A1 DE4030753A1 DE19904030753 DE4030753A DE4030753A1 DE 4030753 A1 DE4030753 A1 DE 4030753A1 DE 19904030753 DE19904030753 DE 19904030753 DE 4030753 A DE4030753 A DE 4030753A DE 4030753 A1 DE4030753 A1 DE 4030753A1
- Authority
- DE
- Germany
- Prior art keywords
- sound
- calculated
- point
- waves
- spherical waves
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/66—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/02—Analysing fluids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02836—Flow rate, liquid level
Abstract
Description
Die vorliegende Erfindung betrifft ein Verfahren zur Berech nung der Schallausbreitung in strömenden Medien.The present invention relates to a method for calculation Sound propagation in flowing media.
Die Kenntnis über die Schallausbreitung in strömenden Medien ist beispielsweise wesentlich für die Entwicklung und Opti mierung von Ultraschall-Durchflußmessern.Knowledge of sound propagation in flowing media is, for example, essential for development and opti mation of ultrasonic flow meters.
Bisher sind Berechnungsverfahren der folgenden Arten bekannt geworden:So far, calculation methods of the following types are known become:
- a) Messung der Schallausbreitung mit Sondenmikrofonen in maßstäblich vergrößerten Modellen und Umsetzung der betref fenden Meßergebnisse in entsprechende Konstruktionsvorgaben,a) Measurement of sound propagation with probe microphones in models enlarged to scale and implementation of the subject Find measurement results in corresponding design specifications,
- b) Berechnungsverfahren, z. B. nach der sog. Finite-Elemen te-Methode (hierbei ist eine sehr große Anzahl von Elementen erforderlich, da die Elemente klein gegen die Wellenlänge sein müssen; die Strömungseinflüsse sind in auf dem Markt angebotenen Rechnerprogrammen nicht implementiert),b) calculation methods, e.g. B. according to the so-called finite elements te method (here is a very large number of elements required because the elements are small against the wavelength have to be; the flow influences are in on the market offered computer programs not implemented),
- c) Berechnung der Schallausbreitung nach dem Huygenschen Prinzip ohne Berücksichtigung der Strömungseinflüsse.c) Calculation of sound propagation according to Huygens Principle without taking flow influences into account.
Der vorliegenden Erfindung liegt die Aufgabe zugrunde, ein Verfahren der eingangs genannten Art zu schaffen, das ge genüber dem Stand der Technik eine zeitsparende und damit kostengünstige Berechnung der Schallausbreitung in strömen den Medien und damit eine Optimierung von Ultraschall- Durchflußmessern gestattet. The present invention is based on the object To create procedures of the type mentioned, the ge Compared to the state of the art a time-saving and therefore cost-effective calculation of sound propagation in flows the media and thus an optimization of ultrasound Flow meters allowed.
Zur Lösung dieser Aufgabe wird ein Verfahren nach dem Ober begriff des Patentanspruchs 1 vorgeschlagen, das erfin dungsgemäß dadurch gekennzeichnet ist, daß die durch das strömende Medium hervorgerufene Verwehung jeder dieser Ku gelwellen berechnet wird und daß die Überlagerung der er zeugten Kugelwellen erst am Empfangsort durchgeführt wird.To solve this problem, a procedure according to the Ober Concept of claim 1 proposed that invented is accordingly characterized in that the flowing medium-induced drift of each of these Ku gel waves is calculated and that the superposition of the he generated spherical waves is only carried out at the receiving location.
Vorteilhafte Weiterbildungen der Erfindung sind durch die in den Unteransprüchen angegebeben Merkmale gekennzeichnet.Advantageous developments of the invention are characterized by the in Characteristics indicated in the subclaims.
Im folgenden wird die Erfindung anhand meherer Figuren im einzelnen beschrieben.In the following the invention with reference to several figures in described.
Fig. 1 zeigt schematisch die Ausbreitungscharakteristik einer Kugelwelle. Fig. 1 shows schematically the propagation characteristic of a spherical wave.
Fig. 2 zeigt schematisch die Verwehung der Kugelwelle bei ortsunabhängiger Strömungsgeschwindigkeit. Fig. 2 shows schematically the drift of the spherical wave at a location-independent flow rate.
Fig. 3 zeigt schematisch die Verwehung der Kugelwelle bei ortsabhängiger Strömungsgeschwindigkeit. Fig. 3 shows schematically the drift of the spherical wave at location-dependent flow speed.
Fig. 4 zeigt ein bevorzugtes Ausführungsbeispiel eines Flußdiagramms, das die einzelnen Funktionsblöcke bzw. Funktionsschritte zur Durchführung des erfindungsge mäßen Verfahrens angibt. Fig. 4 shows a preferred embodiment of a flow chart indicating the individual function blocks or functional steps for performing the method according to the invention.
Für das Verfahren wird zur Berechnung der Schallausbreitung in strömenden Medien vom Huygenschen Prinzip ausgegangen, gemäß dem der von einer schwingenden Oberfläche erzeugte Schalldruck durch Zerlegung der Oberfläche in Teilflächen, die klein gegen die Wellenlänge sind, und durch phasenrich tige Überlagerung der von diesen Teilflächen abgestrahlten Kugelwellen am Empfangsort berechenbar ist: Erfindungsgemäß ist vorgesehen, daß die durch das strömende Medium hervor gerufene Verwehung jeder dieser Kugelwellen berechnet wird und daß die Überlagerung der erzeugten Kugelwellen erst am Empfangsort durchgeführt wird. Es wird eine endliche Anzahl von beliebig geformten Reflektoren in dem Schallausbrei tungsweg angenommen, deren Oberflächen wiederum als Aus gangspunkt von Kugelwellen zu betrachten sind.For the procedure is used to calculate the sound propagation based on Huygens principle in flowing media, according to that generated by a vibrating surface Sound pressure by dividing the surface into partial areas, which are small compared to the wavelength, and by phase term superposition of the radiated from these partial areas Spherical waves can be calculated at the receiving location: According to the invention it is intended that the flowing through the medium called drift of each of these spherical waves is calculated and that the superposition of the spherical waves generated only on Receiving location is carried out. It will be a finite number of arbitrarily shaped reflectors in the sound spread path assumed, their surfaces in turn as Aus point of view of spherical waves are to be considered.
Die betreffenden Oberflächen werden in Teilflächen zerlegt, die klein gegen die Wellenlänge des Schalls sind, wobei in Ansatz gestellt wird, daß jede dieser Teilflächen eine Ku gelwelle abstrahlt, deren Intensität proportional zu ihrer Ausbreitungsgeschwindigkeit und zu ihrer Fläche ist. Für jeweils zwei Oberflächen, die über einen Schallweg mitein ander verbunden sind, wird eine Laufzeitmatrix aufgestellt, wozu für jeden Punkt des "Schallsenders" und für jeden Punkt des "Schallempfängers" die Laufzeit des Schallsignals be rechnet wird, wobei ausgehend von der direkten Verbindung der Punkte die Verwehung während der Laufzeit iterativ kom pensiert wird. Auf der Grundlage der berechneten Laufzeit matrizen werden sowohl stationäre als auch zeittransiente Vorgänge durch Multiplikation der Eingangssignale mit den Matrizen berechnet.The surfaces concerned are broken down into sub-areas, which are small compared to the wavelength of the sound, where in Approach is made that each of these sub-areas a Ku emits a wave whose intensity is proportional to its Speed of propagation and its area is. For two surfaces each, which coexist over a sound path are connected, a runtime matrix is drawn up, why for every point of the "sound transmitter" and for every point the "sound receiver" be the duration of the sound signal is calculated, starting from the direct connection the points the drift iteratively com during runtime is pensed. Based on the calculated transit time Matrices become both stationary and time-transient Operations by multiplying the input signals by Matrices calculated.
Die Berechnungsschritte werden parametrisiert und als Grundlage für ein geometrieunabhängiges Simulationsprogramm benutzt, wonach die betreffende Geometrie durch den Pro grammbenutzer festlegbar ist.The calculation steps are parameterized and as Basis for a geometry-independent simulation program used, according to which the geometry in question by the Pro gram user can be determined.
Die betreffenden Geometriedaten der schallerzeugenden Ober flächen, sämtlicher Reflektoren und der schallempfangenden Oberflächen werden in einen mit dem Simulationsprogramm geladenen Rechner eingegeben, es werden die Schallwege (z. b. Schallsender - Reflektor 1 - Reflektor 2 - Schallempfän ger), und es werden die Strömungsprofile abschnittsweise entweder als Ortsfunktionen oder als Rasterdaten eingegeben.The relevant geometry data of the sound-generating surfaces, all reflectors and the sound-receiving surfaces are entered into a computer loaded with the simulation program, the sound paths (e.g. sound transmitter - reflector 1 - reflector 2 - sound receiver), and the flow profiles are either sectionally as Location functions or entered as raster data.
Schließlich kann vorteilhafterweise vorgesehen sein, daß eine graphische Ausgabe der berechneten Schalldruckverläufe sowohl als Zeitsignal für jeden Ort als auch als Ortsfunk tion für jeden Zeitpunkt bewirkt wird.Finally, it can advantageously be provided that a graphical output of the calculated sound pressure curves both as a time signal for each location and as a local radio tion is effected for every point in time.
Fig. 1 zeigt, wie bereits angegeben, schematisch die Aus breitungscharakteristik einer Kugelwelle. Fig. 2 u. Fig. 3 zeigen jeweils schematisch die Verwehung der Kugelwelle bei ortsunabhängiger Strömungsgeschwindigkeit bzw. bei ortsab hängiger Strömungsgeschwindigkeit. Fig. 1 shows, as already indicated, schematically from the propagation characteristics of a spherical wave. Fig. 2 u. Fig. 3 each schematically show the drift of the spherical shaft at a location-independent flow rate or at a location-dependent flow rate.
Die von einem Punktstrahler ausgesandte Kugelwelle wird unter dem Einfluß der Strömung dergestalt verformt, daß die Ausbreitungsrichtung erhalten bleibt, die Punkte an der Oberfläche jedoch verschoben werden (Mitnahmeeffekt).The spherical wave emitted by a spotlight is deformed under the influence of the flow so that the Direction of propagation remains, the points at the However, the surface can be moved (entrainment effect).
Der Oberflächenpunkt der Kugelwelle sei im Winkel α zur x- Achse angenommen. Ohne Strömung können die Koordianten des Punkts wie folgt beschrieben werden:The surface point of the spherical wave is at an angle α to the x- Axis accepted. Without current, the coordinators of the Can be described as follows:
x(t) = c * cos (α) * t
y(t) = c * sin (α) * t (1)x (t) = c * cos (α) * t
y (t) = c * sin (α) * t (1)
Mit Strömung können die Koordinaten des Punkts wie folgt beschrieben werden:With flow, the coordinates of the point can be as follows to be discribed:
x(t) = c * cos (α) * t + vx (α, t) dt
y(t) = c * sin (α) * t + vy (α, t) dt (2)x (t) = c * cos (α) * t + v x (α, t) dt
y (t) = c * sin (α) * t + v y (α, t) dt (2)
Die Integrale bedeuten hierbei ein Zeitintegral über einen vorgegebenen zeitabhängigen Weg.The integrals here mean a time integral over one given time-dependent path.
Aus Fig. 3 geht die Verwehung der Kugelwelle bei ortabhän giger Strömungsgeschwindigkeit hervor.From Fig. 3 shows the drift of the spherical shaft at location-dependent flow velocity.
Da nur der Endpunkt des Strahlweges bekannt ist, nicht je doch der Winkel, unter dem abgestrahlt wird, um den Endpunkt zu erreichen, wird in einer ersten Iteration die mittlere Strömungsgeschwindigkeit entlang der direkten Verbindung zwischen "Sender" und "Empfänger" ermittelt. Der Signalweg beträgt somitSince only the end point of the beam path is known, not ever but the angle at which radiation is emitted around the end point to achieve, the middle will be achieved in a first iteration Flow speed along the direct connection determined between "sender" and "receiver". The signal path is thus
x(t) = c * cos (α) + vx * t
y(t) = c * sin (α) + vy * t (3)x (t) = c * cos (α) + v x * t
y (t) = c * sin (α) + v y * t (3)
Hieraus kann auf den Abstrahlwinkel α und die tatsächliche Laufzeit t zurückgerechnet werden. Diese erste Iteration liefert bei ortsunabhängiger Strömung, die nur Funktion einer Ortskoordinate ist (z. B. laminare Strömung in einem geraden Rohr) bereits konkrete Ergebnisse.From this, the radiation angle α and the actual Term t can be calculated back. This first iteration delivers with location-independent flow, the only function is a location coordinate (e.g. laminar flow in one straight tube) already have concrete results.
Ausgehend vom Abstrahlwinkel α der ersten Iteration wird in weiteren Iterationen der Signalweg verfolgt und so modifi ziert, daß der Endpunkt getroffen wird. Beim Einsatz des Newton′schen Iterationsverfahren setzt dieses Vorgehen selbstverständlich eindeutige Signalwege voraus (z. B. keine strömungsbedingten akustischen Linsen).Starting from the radiation angle α of the first iteration, in followed further iterations of the signal path and so modifi adorns that the end point is hit. When using the Newton's iteration method implements this procedure of course, clear signal paths ahead (e.g. none flow-related acoustic lenses).
Das erfindungsgemäße Verfahren bietet gegenüber bekannten betreffenden Verfahren den Vorteil einer wesentlichen Zeit- und damit Kostenersparnis.The method according to the invention offers compared to known relevant process the advantage of a substantial time and thus cost savings.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE19904030753 DE4030753A1 (en) | 1990-09-28 | 1990-09-28 | Calculating sound dispersion in flowing media - separating surface into partial areas small compared with wavelengths produced by sound pressure of vibrating surface |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE19904030753 DE4030753A1 (en) | 1990-09-28 | 1990-09-28 | Calculating sound dispersion in flowing media - separating surface into partial areas small compared with wavelengths produced by sound pressure of vibrating surface |
Publications (1)
Publication Number | Publication Date |
---|---|
DE4030753A1 true DE4030753A1 (en) | 1992-04-02 |
Family
ID=6415197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE19904030753 Withdrawn DE4030753A1 (en) | 1990-09-28 | 1990-09-28 | Calculating sound dispersion in flowing media - separating surface into partial areas small compared with wavelengths produced by sound pressure of vibrating surface |
Country Status (1)
Country | Link |
---|---|
DE (1) | DE4030753A1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1914008A3 (en) * | 2006-10-20 | 2010-06-02 | Kawasaki Jukogyo Kabushiki Kaisha | Jet noise modeling method, jet noise analyzing method,and aircraft designing method |
EP2216316A1 (en) | 2004-11-04 | 2010-08-11 | Bayer CropScience Aktiengesellschaft | Phenyl acetic acid derivatives |
CN117310671A (en) * | 2023-11-29 | 2023-12-29 | 中国海洋大学 | Shallow sea sound source distance environment self-adaptive estimation method applying frequency dispersion elimination transformation |
-
1990
- 1990-09-28 DE DE19904030753 patent/DE4030753A1/en not_active Withdrawn
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2216316A1 (en) | 2004-11-04 | 2010-08-11 | Bayer CropScience Aktiengesellschaft | Phenyl acetic acid derivatives |
EP2216317A1 (en) | 2004-11-04 | 2010-08-11 | Bayer CropScience AG | Phenyl acetic acid halogenides |
EP2253207A1 (en) | 2004-11-04 | 2010-11-24 | Bayer CropScience AG | Substituted phenyl acetic acid esters |
EP1914008A3 (en) * | 2006-10-20 | 2010-06-02 | Kawasaki Jukogyo Kabushiki Kaisha | Jet noise modeling method, jet noise analyzing method,and aircraft designing method |
US8219363B2 (en) | 2006-10-20 | 2012-07-10 | Kawasaki Jukogyo Kabushiki Kaisha | Jet noise modeling method, jet noise analyzing method, and aircraft designing method |
CN117310671A (en) * | 2023-11-29 | 2023-12-29 | 中国海洋大学 | Shallow sea sound source distance environment self-adaptive estimation method applying frequency dispersion elimination transformation |
CN117310671B (en) * | 2023-11-29 | 2024-03-01 | 中国海洋大学 | Shallow sea sound source distance environment self-adaptive estimation method applying frequency dispersion elimination transformation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11763474B2 (en) | Method for generating simulated point cloud data, device, and storage medium | |
EP1878308B1 (en) | Device and method for generation and processing of sound effects in spatial audio reproduction systems using a graphical user interface | |
AT521120B1 (en) | Method and device for determining a radar cross section, method for training an interaction model and radar target emulator and test bench | |
DE2633003A1 (en) | FLOW METER | |
DE3007157A1 (en) | DOPPLER METHOD FOR FLOW MEASUREMENT | |
DE19548433C2 (en) | Multi-channel ultrasonic flow measuring device and method for measuring the flow rate of a flow or an open channel | |
DE4030753A1 (en) | Calculating sound dispersion in flowing media - separating surface into partial areas small compared with wavelengths produced by sound pressure of vibrating surface | |
DE102014118622A1 (en) | Method for simulatively determining an interaction between a sensor of a motor vehicle and a virtual object in a virtual surrounding area of the motor vehicle and computing device | |
JP2022110069A (en) | Generation method of simulation scene, device, electronic device and storage medium | |
EP3729044A1 (en) | Method for analysing an automation system of an installation, emulator for at least partial virtual operation of an automation system of an installation, and system for analysing an automation system of an installation | |
JP3309067B2 (en) | In-plane flow display method | |
CN107462883A (en) | Method and laser range finder of the measurement point to any dimension plan range in space | |
KR102418929B1 (en) | Method for computer-implemented simulation of radar raw data | |
Kagawa et al. | Finite element simulation of non-linear sound wave propagation | |
JP3492347B2 (en) | Apparatus, method and program for simulating underwater sound propagation characteristics | |
CN116996134A (en) | Vehicle antenna performance simulation determination method and device, electronic equipment and storage medium | |
Lundén et al. | Psychoacoustic evaluation as a tool for optimization in the development of an urban soundscape simulator | |
DE60224864T2 (en) | A method of determining the values of an electromagnetic field generated by a radio base station in an urban environment | |
Kagawa et al. | DISCRETE HUYGENS'MODELLING SIMULATION OF SOUND WAVE PROPAGATION IN VELOCITY VARYING ENVIRONMENTS | |
EP0917645A1 (en) | Method and device for ultrasonic flow measurement | |
Fischer et al. | Time-domain imaging techniques for aeroacoustic sources | |
CN109143333A (en) | Forward modeling method and computer readable storage medium based on triangulation model | |
Faure et al. | On the pressure field of a transducer in the form of a curved strip | |
Schmitz et al. | Merging software for sound reinforcement systems and for room acoustics | |
Ende et al. | Auralization and Visualization for Infrastructure Planning in the Joint Project |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
8139 | Disposal/non-payment of the annual fee |