Sensors & Transducers 2015 by IFSA Publishing, S. L.

Size: px
Start display at page:

Download "Sensors & Transducers 2015 by IFSA Publishing, S. L."

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

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 information

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy

Characterization 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 information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation 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 information

Homework 7 Due 18 November at 6:00 pm

Homework 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 information

05 The Continuum Limit and the Wave Equation

05 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 information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 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 information

Separation of Variables

Separation 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 information

12.11 Laplace s Equation in Cylindrical and

12.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 information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE 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 information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal 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 information

P. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom

P. 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 information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 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 information

Modelling of Transit-Time Ultrasonic Flow Meters Under Multi-phase Flow Conditions

Modelling 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 information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-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 information

An Approach for Design of Multi-element USBL Systems

An 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 information

The total derivative. Chapter Lagrangian and Eulerian approaches

The 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 information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture 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 information

5-4 Electrostatic Boundary Value Problems

5-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 information

Study on aero-acoustic structural interactions in fan-ducted system

Study 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 information

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating 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 information

6. Friction and viscosity in gasses

6. 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 information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The 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 information

Approaches for Predicting Collection Efficiency of Fibrous Filters

Approaches 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

'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 information

Convective heat transfer

Convective 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 information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: 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 information

Chapter 2 Lagrangian Modeling

Chapter 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 information

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić

CHARACTERISTICS 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 information

PDE Notes, Lecture #11

PDE 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 information

Euler equations for multiple integrals

Euler 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 information

Schrödinger s equation.

Schrö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 information

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function

Relation 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 information

A Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion

A 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 information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math 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 information

Introduction to the Vlasov-Poisson system

Introduction 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 information

Quantum Mechanics in Three Dimensions

Quantum 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 information

Physics 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 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 information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 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 information

Get Solution of These Packages & Learn by Video Tutorials on WAVE OPTICS

Get 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 information

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation

Electromagnet 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 information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon 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 information

RETROGRADE WAVES IN THE COCHLEA

RETROGRADE 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 information

Generalization of the persistent random walk to dimensions greater than 1

Generalization 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 information

The effect of nonvertical shear on turbulence in a stably stratified medium

The 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 information

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA

TOWARDS 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 information

Situation awareness of power system based on static voltage security region

Situation 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 information

HEAT TRANSFER ENHANCED PARABOLIC TROUGH RECEIVER FOR DSG WITH CAPILLAR SYSTEMS

HEAT 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 information

Applications of First Order Equations

Applications 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 information

Average value of position for the anharmonic oscillator: Classical versus quantum results

Average 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 information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization 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 information

ELECTRON DIFFRACTION

ELECTRON 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 information

Asymptotics of a Small Liquid Drop on a Cone and Plate Rheometer

Asymptotics 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 information

IPMSM Inductances Calculation Using FEA

IPMSM 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 information

Distribution Theory for Discontinuous Test Functions and Differential Operators with Generalized Coefficients

Distribution 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 information

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

The 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 information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE 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 information

The Principle of Least Action and Designing Fiber Optics

The 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 information

A General Analytical Approximation to Impulse Response of 3-D Microfluidic Channels in Molecular Communication

A 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 information

arxiv:nlin/ v1 [nlin.cd] 21 Mar 2002

arxiv: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 information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In 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)

(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 information

The Principle of Least Action

The 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 information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic 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 information

Predictive Control of a Laboratory Time Delay Process Experiment

Predictive 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 information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic 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 information

Crack-tip stress evaluation of multi-scale Griffith crack subjected to

Crack-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 information

Optimized Signal De-noising Algorithm for Acoustic Emission Leakage

Optimized 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 information

Comparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope

Comparative 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 information

Vectors in two dimensions

Vectors 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 information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-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 information

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem

Topological 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

θ 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 information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic 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 information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing 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 information

The Press-Schechter mass function

The 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 information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching 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 information

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section

The 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 information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. 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 information

ELEC3114 Control Systems 1

ELEC3114 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 information

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d

SYNTHESIS 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

(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 information

Transmission 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 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 information

Conservation laws a simple application to the telegraph equation

Conservation 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 information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized 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 information

Improving the Lorentz Force Amplitude of an EMAT Using Stacked Coil Configuration

Improving 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 information

Dusty Plasma Void Dynamics in Unmoving and Moving Flows

Dusty 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 information

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS

CONSERVATION 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 information

Track Initialization from Incomplete Measurements

Track 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 information

Wavelet-based Local Tomography

Wavelet-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 information

PAijpam.eu RELATIVE HEAT LOSS REDUCTION FORMULA FOR WINDOWS WITH MULTIPLE PANES Cassandra Reed 1, Jean Michelet Jean-Michel 2

PAijpam.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 information

Two Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices

Two 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 information

Single Arm, Centrifugal, Water Turbine for Low Head and Low Flow Application: Part 1- Theory and Design

Single 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 information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust 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 information

Gaussian processes with monotonicity information

Gaussian 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 information

arxiv: v2 [cond-mat.stat-mech] 11 Nov 2016

arxiv: 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 information

24th European Photovoltaic Solar Energy Conference, September 2009, Hamburg, Germany

24th 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 information

Research Article Numerical Analysis of Inhomogeneous Dielectric Waveguide Using Periodic Fourier Transform

Research 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 information

Centrum voor Wiskunde en Informatica

Centrum 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 information

Gravitation 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 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 information

Balance laws on domains with moving interfaces. The enthalpy method for the ice melting problem.

Balance 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