Sensors & Transducers 2015 by IFSA Publishing, S. L.
|
|
- Stewart Berry
- 6 years ago
- Views:
Transcription
1 Sensors & Transucers, Vol. 184, Issue 1, January 15, pp Sensors & Transucers 15 by IFSA Publishing, S. L. Non-invasive an Locally Resolve Measurement of Soun Velocity by Ultrasoun Mario WOLF, Elfgar KÜHNICKE Soli State electronics Laboratory, Technical University of Dresen, 16, Germany Tel.: , fax: Receive: 14 November 14 /Accepte: 15 December 14 /Publishe: 31 January 15 Abstract: A new metho to measure soun velocity an istance simultaneously an locally resolve by an ultrasoun annular array is evelope in meia with constant soun velocity an then applie in meia with continuously changing properties. Instea of using reflectors at known positions the echoes of moving scattering particles are analyze to etermine the focus position. Although the metho reaches a very high accuracy for constant soun velocity, there are systematic eviation between measurements an soun fiel simulations using Fermat s principle. The soun propagation has to be escribe with a moifie wave equation. A new approach etermining GREEN s functions for a half space with continuously changing properties is presente. It combines the high frequency approximation with an integral transform metho. Copyright 15 IFSA Publishing, S. L. Keywors: Ultrasoun annular arrays, Measurement of soun velocity, Locally resolve, High frequency approximation, Integral transform metho. 1. Introuction A locally resolve monitoring of soun velocity allows estimating locally physical quantities like concentration or temperature or material properties like ensity or elasticity. This facilitates investigating an optimizing many inustrial processes, like mixing or chemical reactions, as well as meical therapy like hyperthermia for cancer treatment. In this contribution, a measurement technique is applie, which allows measuring soun velocity locally resolve using an ultrasoun annular array with concentric rings. In contrast to conventional tomographic techniques, it works without any reflectors or aitional receivers at known positions. Instea of evaluating various propagation paths, the focusing of the array is varie an the focus position an the soun velocity to the focus point are etermine simultaneously by analyzing the echoes of moving scattering particles. This is possible because the focus position epens on the soun velocity an the parameters of the use transucer. Therefore, the time of flight to the focus is use with calibration curves for the simultaneous etermination of soun velocity an focus position. The time of flight to the focus point is etermine from the average amplitue of the echo signals: The emitte wave is reflecte at each particle while the amplitue of the reflecte signal is proportional to the amplitue of the incient wave (Fig. 1a). Therefore the echoes from particles within the focus area are strongest. Although a single echo nearly appears like noise averaging over a sufficient number of signals generates a clear maximum (see Fig. 1b). As particles are in movement it is possible to consier a uniform istribution of particles in average time. So the 53
2 Sensors & Transucers, Vol. 184, Issue 1, January 15, pp average echo signal amplitue is a measure for the pressure along the acoustic axis of the transucer.. Metho to Measure soun Velocity an Distance Simultaneously.1. Soun Fiel Calculation by GREEN S Functions GREEN s functions are a beneficial tool for moelling soun propagation in homogeneous meia, which escribe the impulse response of a point source in a space with a specific geometry. Exact solutions ha been foun for half spaces, plates an layere meia by means of integral transform [3], [4]. Also approximate harmonic GREEN s functions are calculate giving the transfer function in the Fourier transforme omain [5]. GREEN s functions for a half space are also use in this contribution 1 jkr G ω ( R, θ ) = S( θ ) e, (1) R Fig. 1. Normalize average echo signal amplitue for ifferent numbers of averages in water of 3 C. The basic set-up for the measurements is quiet simple. Only a focusing single probe an a evice for pulsing an ata acquisition are require. The probe can be locate in the examine meium without any ajustment. The meium has to contain scattering particles much smaller than wavelength an in a low concentration, so that the properties of the propagation meium are not influence. The metho ha been introuce by Lenz [1] for meia with constant soun velocity using a focusing single probe. This contribution enhances the metho to annular arrays, which allow moving the focus along the acoustic axis an so measuring in ifferent epths. At first this is qualifie for meia with constant soun velocity an than employe for meia with nonconstant soun velocity. A new metho for moelling soun propagation in meia with continuously changing material parameters is introuce for the evaluation. The moifie wave equation cause by non-constant material parameters is solve for a point source with an integral transform metho an a high frequency approximation []. This paper is ivie in five sections. Section explains how to measure soun velocity an focus position simultaneously. It motivates the moelling of soun propagation in meia with non-constant material properties which is iscusse in Section 3. Section 4 shows some measurement results for meia with constant soun velocity an a linear temperature graient in water. Section 5 gives a summary an perspectives. with the wave number k, ie istance R between source an observation point an the angle θ between the vector between source an observation point an the normal vector of the surface where the source is locate. The irectivity pattern S(θ) concerns the bounary conitions at the interface. An example for an interface of water an PZT (lea- zirconatetitanate) is shown in Fig.. Fig.. Directivity pattern of a point source on a water-pzt-interface... Annular Arrays for Focusing Annular arrays nee far less elements for a well focusing as the wiely sprea linear array. The array use in this contribution consists of 6 elements (see Fig. 3 for structure an imensions) an reaches an extension of the soun beam in the range of the wavelength. Focusing works by a superposition of measure signals where each is elaye with a calculate time: The soun paths from each element to the esigne focus point are etermine (Fig. 4). With a known soun velocity of the propagation meium these paths allow to calculate the ifferences in time of flight, which are use as elay times for each element so that all waves arrive at the focus point at the same time. 54
3 Sensors & Transucers, Vol. 184, Issue 1, January 15, pp Element ri [mm] ra [mm] c = 143 m/s c = 156 m/s Fig. 3. Structure an imensions of the curve annular array, centre frequency: fm = 9 MHz, Curvature: R = 5 mm; ri an ra are the insie an outsie raii of the rings. As there is positive interference of course the pressure at this point becomes maximal. Fig. 5. Soun fiels for ifferent soun velocities. Fig. 4. Focusing soun paths to focus points for ifferent soun velocities. Of course if e.g. the temperature in the meium changes this causes a change of soun velocity. The focus moves because the positive interference arises at another point (Fig. 4, re an blue lines). As calculate soun fiels are in the Fourier transforme omain also the elay times have to be transforme into phase shifts. A superposition of the phase shifte soun fiels generates the focus at the esigne point. Fig. 5 shows the soun fiels of the use array, riven with the same elay times for ifferent soun velocities. They were chosen so that the focus woul arise at a istance of 4 mm in a meium with a soun velocity of c = 15 m/s..3. Simultaneous Determination of Soun Velocity an Focus Position It has been shown, that a change of soun velocity isplaces the focus position resulting in a ifferent time of flight to the focus. If the parameters of the transucer are known the time of flight can be calculate as a function of soun velocity. Fig. 6 shows functions for several sets of elay times. The labels mean that the use elay times woul cause a focus in a calibration meium with c = 15 m/s at the corresponing epth. Fig. 6. Calibration curves. The measure time of flight to the focus point can be use to etermine the soun velocity by these curves. Of course then the real focus position emerges. If the elay times are varie the focus is move along the acoustic axis of the array. So the soun velocity can be measure at each point, which allows measuring the istribution of soun velocity. 3. Soun Propagation in Meia with Non-constant Properties As mentione in the previous section the parameters of the transucer have to be known to calculate the soun fiel of the transucer. The focus position also has to be preicte for meia with continuously changing properties where aitional consierations are require Simultaneous Determination of Soun Velocity an Focus Position The most obvious way to moel the soun propagation in meia with continuously changing properties is applying Fermat s principle. It says that a wave propagates along that path which results in a minimal propagation time. 55
4 Sensors & Transucers, Vol. 184, Issue 1, January 15, pp It can be use to erive Snell s law but if soun velocity changes continuously the propagation path becomes curve an has to be etermine by a calculus of variation [6]. The time of flight T to a point P(x,y) an the path length S to this point are: with a = 1 c( x) 1+ 1 a c T = ln[ ], () c 1+ 1 a c ( x) 1 S = [arcsin( c ( x ) a ) arcsin( c a )], (3) a ( y c x x ) y + 4 y ( c + x) for a linear soun velocity graient (c(x) = c+ x,) an a point source in the origin of coorinates. This ha been use to erive fiels of point sources for linear soun velocity graients. The impulse response is assume by 1 Gδ ( x, y, t) = δ ( t T ( x, y)) (4) S( x, y) with the Dirac function δ. In the Fourier transforme omain this is 1 jωt ( x, y) Gω ( x, y, ω ) = e (5) S( x, y) Fig. 7 shows the phase plots of two examples where soun velocity increases (left) an ecreases (right) in positive irection x. p = v, χ p= v, (6) with pressure p, particle velocity v, an the material parameters mass ensity an elasticity χ. If the material parameters are constant, ifferentiating these equations with respect to location an time leas to the well-known wave equation. If these parameters are functions of location an aitional term appears in the graient of (6). This aitional term also appears in the wave equation: 1 Δp + p χp = (7) Also the efinition of potential Φ is moifie: 1 v = graφ (8) The following examination shall be one for a one-imensional epenence of the material properties in irection z One-imensional Solution Consiering a plane wave propagating in the irection of the graient of the material properties the wave equation for the potential is obtaine Φ+ z z 1 z Φ χ Φ = t (9) If the mass ensity epens linearly on z this equation has a solution in the form of a generalize power series [7]: n+ K ( z z ) (1) Φ = cn n= 1 Fig. 7. Phase plot for positive an negative soun velocity graient. However, there are alreay numerical problems in evaluating the coefficients an the solution with respect to the convergence for this one-imensional case. So it seems not to be feasible fining an exact solution for a two- or three-imensional problem. Although Fermat s principle gives the correct time of flight of the wave there is no information about the amplitue of the wave. 3.. Derivation of the Moifie Wave Equation The soun propagation in liquis is base on two funamental equations, the equations of motion an elasticity: 3.4. High Frequency Approximation an Integral Transform High frequency approximation ha been evelope in geophysics an is applie in techniques like ray tracing [8]. In this contribution, harmonic GREEN's functions shall be erive with this approach in combination with an integral transform metho. This allows calculating a transfer function for a point source for a specific geometry. Just the 56
5 Sensors & Transucers, Vol. 184, Issue 1, January 15, pp axially symmetric problem is solve because it is much easier manageable than the general problem. The wave equation in cylinrical coorinates (r, θ, z) is use: 1 = r + Φ +... r r r z 1 Φ χ Φ z z t (11) The single erivative to r in the secon term vanishes ue to its scalar multiplication with the graient of the mass ensity, having only a component in z irection. Applying a Hankel transform as escribe in [3] with respect to r Φ H ( ξ,z,t) = Φ( r,z,t) J ( jωξr)rr (1) leas to a one-imensional wave equation in the transforme omain, which is enote by the inex H: = z ξ Φ H H ω Φ 1 + Φ z z H χ Φ t H... (13) The high frequency approximation assumes that equation (8) can be solve by the following ansatz Φ = A jω t T ( x,z ) ( x,z) e ( ) (14) Now, Φ is replace with this ansatz an the terms are arrange accoring to its powers of ω. = Aω T 1 A T 1 T jω z A z z z (15) z... + T χ z z ξ... 1 A A + z z A ( ξ,z) = A ( ( z ) χ( z ) ξ ) ( z) ( z) χ( z) ξ z ( ) ( ), (17) with the solely free parameter A. Note that this leas to the solution for a homogeneous meium if an χ are constants. All methos of generalize ray theory explaine in [3] like the erivation of source functions for point source acting on an interface consiering the bounary conitions can be applie to this solution. Finally, the inverse transformation has to be one: jω t T ( ξ,z ) ( ) = ω A( ξ,z) e Φ r,z ( ) J ( jωξr )ξξ (18) Current work is on an evaluation of this integral with a steepest escent approximation as it is escribe in [9]. The approximation facilitates the integral into a solvable form. This causes a neglect of surface waves being not of interest for the presente application. However, the metho is complicate because of the integral expression of T in the exponent. Though, the integral can be evaluate by a finite series expansion resulting in aitional terms containing higher powers of ξ. 4. Measurement Results 4.1. Measurements for Constant Soun Velocity Primarily, soun velocity was measure in ifferent meia with constant temperature, as escribe in Section II, to compare measurements an simulations. It is striking that the measure curve for Θ = 6 C fits very well to the simulate one, whereas the measure curve for Θ = 6 C eviates (Fig. 8). Therefore, aitional measurements with water-ethanol solutions were one. Different soun velocities can be ajuste with ifferent concentrations of ethanol, without a change in temperature [1]. Assuming large ω the first two lines are equate to zero inepenently an the frequency-inepenent thir line is neglecte. So T can be etermine irectly from the first line an with this solution A is etermine from the secon line T z ( ξ,z) = ± z) χ( z + z ( ξ z, (16) ) Fig. 8. Comparison of measure an calculate times of flight as a function of the use set of elay times for ifferent meia. 57
6 Sensors & Transucers, Vol. 184, Issue 1, January 15, pp It can be recognize that the measure curve of a water-ethanol solution with a soun velocity of 1565 m/s in Fig. 8 agrees to the simulate curve. This implicates that the working temperature of the probe has to be taken into account calculating the calibration curves. To analyze the accuracy measurements with ifferent numbers of averages have been one. Table 1 shows that the relative error of soun velocity is less than 1 %. Principle. Fig. 1 shows the comparison of calculate an measure times of flight to the focus as a function of the use set of elay times, corresponing with a focus (Fok) in the calibration meium water of C (soun velocity graient in water from 4 C at the transucer to 6 C at a istance of 5 mm). Table 1. Measurement capability. Number of average signals Stanar eviation ± 3σ of time of flight Absolute error of soun velocity 5 5 ns 13 m/s 5 ns 6 m/s 1 15 ns 4 m/s 5 1 ns.5 m/s 4.. Measurements for a Soun Velocity Graient An evient way to achieve a soun velocity graient is to generate a temperature graient with water, because the soun velocity as a function of temperature is well known for water [11] an it can be generate in a stable state. Fig. 9 shows the general set-up. Water in a basin locate at the bottom is kept at a temperature of 6 C. A secon smaller basin is place above. It contains a metal plate at its bottom for goo thermal conuction an a heat source at its top. This generates a vertical layere arrangement of warm water above cooler water, whereby a flui flow is avoie. The temperature is measure by an array of temperature sensors to etermine the soun velocity profile in the experimental set-up. Fig. 1. Comparison of measure an calculate times of flight as a function of the use set of elay times for a temperature graient. Although the notable ifference is just in the range of one microsecon this woul cause an error of more than 1 m/s in etermining the soun velocity. For aitional examinations, the soun velocity was measure conventionally via measuring the time of flight to a reflector at a known position. Moving the reflector along the acoustic axis allows a stepwise reconstruction of the soun velocity profile. Aitionally the soun velocity was etermine from temperature sensors again. Fig. 11 shows a comparison of the two conventional methos to etermine the mean soun velocity between the probe an a reflector at istance z. First, the time of flight for various reflector istances is measure (blue line). Secon, the temperatures are measure at various locations an converte to a soun velocity accoring to [11] an average over the propagation path (green line). The systematic eviation can be seen here, too Measurement of mean soun velocity c via time of flight c via temperatur c mean [m/s] Fig. 9. Experimental set-up for a soun velocity graient. The soun fiel for this graient an the time of flight to the focus were calculate applying Fermat's z [mm] Fig. 11. Comparison of two conventional methos to etermine the mean soun velocity between the probe an a reflector at istance z. 58
7 Sensors & Transucers, Vol. 184, Issue 1, January 15, pp Both eviations (see Fig. 1 an Fig. 11) result from the eficient assumption that wave propagation can be escribe with the wave equation for homogeneous meia. 5. Summary an Perspective This contribution iscusses a non-invasive metho to measure the soun velocity locally resolve along the acoustic axis of the use annular array. The three-imensional istribution of soun velocity can be obtaine by scanning. The metho has been qualifie for meia with constant soun velocity an then extene to meia with continuously changing properties. It has been shown that the continuous change of material properties has to be taken into account for the moelling of soun propagation. Fermat s principle gives the correct time of flight, but to obtain information about the amplitue, the moifie wave equation has to be solve. The potential of a point source has been calculate in the Hankel transforme omain. The inverse transform is actually realize an will allow calculating GREEN s functions for meia with continuously changing properties. Due to the assume change of material properties in only one imension a change of these properties in other imensions woul cause a lateral eviation of the focus position. This effect has to be examine in further works. Acknowlegements The authors woul like to thank Deutsche Forschungsgemeinschaft (DFG) for the financial support of the ongoing research project KU175/17-1. References [1]. M. Lenz, M. Bock, E. Kühnicke, J. Pal, A. Cramer, Measurement of the soun velocity in fluis using the echo signals from scattering particles, Ultrasonics, Vol. 5, Issue 1, January 1, pp []. M. Wolf, E. Kühnicke, Monitoring of Temperature Distribution in Liquis with Ultrasoun by Locally Resolve Measuring of Soun Velocity, in Proceeings of the 5 th International Conference on Sensor Device Technologies an Applications (SENSORDEVICES' 14) Lisbon, Portugal, 16- November 14, pp [3]. Y. H. Pao, R. R. Gajewski, The generalize ray theory an transient responses of layere elastic solis, in Physical Acoustics (W. P. Mason, e.), Acaemic Press, New York, Vol. 13, [4]. A. N. Ceranoglu, Y. H. Pao, Propagation of Elastic Pulses an Acoustic Emission in a Plate Part 1, Journal of Applie Mechanics, Vol. 48, Issue 1, 1981, pp [5]. E. Kuhnicke, Calculation of three-imensional harmonic waves in layere meia, in Proceeings of the IEEE International Ultrasonics Symposium (IUS), 7-1 Nov. 1995, pp [6]. M. Wolf, E. Kühnicke, M. Lenz, Moelling of soun propagation in meia with continuously changing properties towars a locally resolve measurement of soun velocity, in Proceeings of the IEEE International Ultrasonics Symposium (IUS), 1-5 July 13, pp [7]. L. Fuchs, Zur Theorie er linearen Differentialgleichungen mit veränerlichen Coefficienten, Journal für ie Reine un Angewante Mathematik, Vol. 1866, No. 66, 1995, pp [8]. V. Cerveny, Seismic Ray Theory, Cambrige University Press, 5. [9]. E. Kühnicke, Plane arrays Funamental investigations for correct steering by means of soun fiel calculations, Wave Motion, Vol. 44, Issue 4, March 7, pp [1]. C Burton, A Stuy of Ultrasonic Velocity an Absorption in liqui mixtures, Journal of the Acoustical Society of America, Vol., Issue, [11]. W. Marczak, Water as a stanar in the measurements of spee of soun in liquis, Journal of the Acoustical Society of America, 1997, pp Copyright, International Frequency Sensor Association (IFSA) Publishing, S. L. All rights reserve. ( 59
Table of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationCharacterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy
JOURNAL OF APPLIED PHYSICS VOLUME 85, NUMBER 1 15 JUNE 1999 Characterization of lea zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy Haifeng Wang, Wenhua Jiang, a) an Wenwu
More informationSimulation of Angle Beam Ultrasonic Testing with a Personal Computer
Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationP. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
Time-harmonic torsional waves in a composite cyliner with an imperfect interface J. R. Berger a) Division of Engineering, Colorao School of Mines, Golen, Colorao 80401 P. A. Martin b) Department of Mathematics,
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationModelling of Transit-Time Ultrasonic Flow Meters Under Multi-phase Flow Conditions
Moelling of Transit-Time Ultrasonic Flo Meters Uner Multi-phase Flo Conitions Matej Simura an Lars Duggen Mas Clausen Institute University of Southern Denmark Sønerborg, 64, Denmark Email: simura@mci.su.k
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationAn Approach for Design of Multi-element USBL Systems
An Approach for Design of Multi-element USBL Systems MIKHAIL ARKHIPOV Department of Postgrauate Stuies Technological University of the Mixteca Carretera a Acatlima Km. 2.5 Huajuapan e Leon Oaxaca 69000
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More information5-4 Electrostatic Boundary Value Problems
11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions
More informationStudy on aero-acoustic structural interactions in fan-ducted system
Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,
More informationEvaporating droplets tracking by holographic high speed video in turbulent flow
Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationApproaches for Predicting Collection Efficiency of Fibrous Filters
Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationCHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić
Sensors an Actuators A: Physical 168 (211) 149-154 oi: 1.116/j.sna.211..7 211 Elsevier B.V. CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS Jože Kutin *, Ivan Bajsić Laboratory of
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationRelation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function
Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives
More informationA Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion
Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationGet Solution of These Packages & Learn by Video Tutorials on WAVE OPTICS
FREE Downloa Stuy Package from website: wwwtekoclassescom & wwwmathsbysuhagcom Get Solution of These Packages & Learn by Vieo Tutorials on wwwmathsbysuhagcom PRINCIPLE OF SUPERPOSITION When two or more
More informationElectromagnet Gripping in Iron Foundry Automation Part II: Simulation
www.ijcsi.org 238 Electromagnet Gripping in Iron Founry Automation Part II: Simulation Rhythm-Suren Wahwa Department of Prouction an Quality Engineering, NTNU Tronheim, 7051, Norway Abstract This paper
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationRETROGRADE WAVES IN THE COCHLEA
August 7, 28 18:2 WSPC - Proceeings Trim Size: 9.75in x 6.5in retro wave 1 RETROGRADE WAVES IN THE COCHLEA S. T. NEELY Boys Town National Research Hospital, Omaha, Nebraska 68131, USA E-mail: neely@boystown.org
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationSituation awareness of power system based on static voltage security region
The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran
More informationHEAT TRANSFER ENHANCED PARABOLIC TROUGH RECEIVER FOR DSG WITH CAPILLAR SYSTEMS
HEAT TRANSFER ENHANCED PARABOLIC TROUGH RECEIVER FOR DSG WITH CAPILLAR SYSTEMS Rojas, M.E. Parabolic Trough Technology Group, PSA / CIEMAT, Av Complutense,, 8040 Mari, Spain Phone: 34 91 346 6049; Fax:
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationELECTRON DIFFRACTION
ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the
More informationAsymptotics of a Small Liquid Drop on a Cone and Plate Rheometer
Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity
More informationIPMSM Inductances Calculation Using FEA
X International Symposium on Inustrial Electronics INDEL 24, Banja Luka, November 68, 24 IPMSM Inuctances Calculation Using FEA Dejan G. Jerkan, Marko A. Gecić an Darko P. Marčetić Department for Power,
More informationDistribution Theory for Discontinuous Test Functions and Differential Operators with Generalized Coefficients
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 1, 9733 1996 ARTICLE NO 56 Distribution Theory for Discontinuous Test Functions an Differential Operators with Generalize Coefficients P Kurasov* Department
More informationThe Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis
The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationThe Principle of Least Action and Designing Fiber Optics
University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic
More informationA General Analytical Approximation to Impulse Response of 3-D Microfluidic Channels in Molecular Communication
A General Analytical Approximation to Impulse Response of 3- Microfluiic Channels in Molecular Communication Fatih inç, Stuent Member, IEEE, Bayram Cevet Akeniz, Stuent Member, IEEE, Ali Emre Pusane, Member,
More informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More informationIn the usual geometric derivation of Bragg s Law one assumes that crystalline
Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this
More information(3-3) = (Gauss s law) (3-6)
tatic Electric Fiels Electrostatics is the stuy of the effects of electric charges at rest, an the static electric fiels, which are cause by stationary electric charges. In the euctive approach, few funamental
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationCrack-tip stress evaluation of multi-scale Griffith crack subjected to
Crack-tip stress evaluation of multi-scale Griffith crack subjecte to tensile loaing by using periynamics Xiao-Wei Jiang, Hai Wang* School of Aeronautics an Astronautics, Shanghai Jiao Tong University,
More informationOptimized Signal De-noising Algorithm for Acoustic Emission Leakage
1009 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 46, 2015 Guest Eitors: Peiyu Ren, Yancang Li, Huiping Song Copyright 2015, AIDIC Servizi S.r.l., ISBN 978-88-95608-37-2; ISSN 2283-9216 The
More informationComparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope
Proceeings of the Worl Congress on Engineering Vol WCE, July 6-8,, Lonon, U.K. Comparative Approaches of Calculation of the Back Water Curves in a Trapezoial Channel with Weak Slope Fourar Ali, Chiremsel
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationThe Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationThe influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section
IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A
More information1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity
AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationSYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d
SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,
More information(2012) , ISBN
Cregan, V. an O'Brien, Stephen B.G. an McKee, Sean (2012) Asymptotics of a small liqui rop on a cone an plate rheometer. In: Progress in Inustrial Mathematics at ECMI 2010. Mathematics in Inustry: Progress
More informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationImproving the Lorentz Force Amplitude of an EMAT Using Stacked Coil Configuration
Sensors & Transucers, Vol. 155, Issue 8, August 213, pp. 262-27 Sensors & Transucers 213 by IFSA http://www.sensorsportal.com Improving the Lorentz Force Amplitue of an MAT Using Stacke Coil Configuration
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationTrack Initialization from Incomplete Measurements
Track Initialiation from Incomplete Measurements Christian R. Berger, Martina Daun an Wolfgang Koch Department of Electrical an Computer Engineering, University of Connecticut, Storrs, Connecticut 6269,
More informationWavelet-based Local Tomography
Inian Society for Non-Destructive Testing Hyeraba Chapter Proc. National Seminar on Non-Destructive Evaluation Dec. 7-9, 006, Hyeraba Wavelet-base Local Tomography M.V. Gopala Rao, E.M.L. Tanua an S. Vathsal
More informationPAijpam.eu RELATIVE HEAT LOSS REDUCTION FORMULA FOR WINDOWS WITH MULTIPLE PANES Cassandra Reed 1, Jean Michelet Jean-Michel 2
International Journal of Pure an Applie Mathematics Volume 97 No. 4 2014 543-549 ISSN: 1311-8080 (printe version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu oi: http://x.oi.org/10.12732/ijpam.v97i4.13
More informationTwo Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices
Journal of Physics: Conference Series PAPER OPEN ACCESS Two Dimensional Numerical Simulator for Moeling NDC Region in SNDC Devices To cite this article: Dheeraj Kumar Sinha et al 2016 J. Phys.: Conf. Ser.
More informationSingle Arm, Centrifugal, Water Turbine for Low Head and Low Flow Application: Part 1- Theory and Design
Energy an Power 2018, 8(2): 51-55 DOI: 10.5923/j.ep.20180802.03 Single Arm, Centrifugal, Water Turbine for Low ea an Low Flow Application: Part 1- Theory an Design Kiplangat C. Kononen 1, Augustine B.
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationGaussian processes with monotonicity information
Gaussian processes with monotonicity information Anonymous Author Anonymous Author Unknown Institution Unknown Institution Abstract A metho for using monotonicity information in multivariate Gaussian process
More informationarxiv: v2 [cond-mat.stat-mech] 11 Nov 2016
Noname manuscript No. (will be inserte by the eitor) Scaling properties of the number of ranom sequential asorption iterations neee to generate saturate ranom packing arxiv:607.06668v2 [con-mat.stat-mech]
More information24th European Photovoltaic Solar Energy Conference, September 2009, Hamburg, Germany
4th European hotovoltaic Solar Energy Conference, 1-5 September 9, Hamburg, Germany LOCK-IN THERMOGRAHY ON CRYSTALLINE SILICON ON GLASS (CSG) THIN FILM MODULES: INFLUENCE OF ELTIER CONTRIBUTIONS H. Straube,
More informationResearch Article Numerical Analysis of Inhomogeneous Dielectric Waveguide Using Periodic Fourier Transform
Microwave Science an Technology Volume 2007, Article ID 85181, 5 pages oi:10.1155/2007/85181 Research Article Numerical Analysis of Inhomogeneous Dielectric Waveguie Using Perioic Fourier Transform M.
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.
More informationGravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann
Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)
More informationBalance laws on domains with moving interfaces. The enthalpy method for the ice melting problem.
Balance laws on omains with moving interfaces. The enthalpy metho for the ice melting problem. Gowin Kakuba 12th March 2008 Outline 1 Balance equations on omains with moving bounaries Introuction Rankine-Hugoniot
More information