INTERNAL RESISTANCE OPTIMIZATION OF A HELMHOLTZ RESONATOR IN NOISE CONTROL OF SMALL ENCLOSURES. Ganghua Yu, Deyu Li and Li Cheng 1 I.
|
|
- Lesley Bond
- 5 years ago
- Views:
Transcription
1 ICV4 Cairns Australia 9- July, 7 ITAL ITAC OPTIMIZATIO OF A HLMHOLTZ OATO I OI COTOL OF MALL CLOU Gangua Yu, Deyu Li and Li Ceng Department of Mecanical ngineering, Te Hong Kong Polytecnic University Hung Home, Kowloon, Hong Kong, A of Cina mail address: mmlceng@polyu.edu. Abstract Tis paper examines te influence of te internal resistance of a Helmoltz resonator on te noise control in a small enclosure. Te absorptive process mainly occurring witin te neck of a Helmoltz resonator provides te resonator wit a damping (internal resistance) property, wic directly dissipates te input energy in te resonator. Te remaining non-dissipated energy is re-radiated back to te enclosure, suc forming an effective secondary sound source, and resulting in acoustic interaction wit te primary source. If te internal resistance of te resonator is low, te acoustic interaction between te enclosure and te resonator sarply splits te targeted resonance peak of te enclosure into two parts, and te peak response is significantly attenuated witin a very narrow frequency band. By appropriately increasing te internal resistance at te resonance of te resonator, te working bandwidt can be enlarged at te expense of sacrificing te control performance due to te decreased amplitude. However, if te resistance is over-increased, te strengt of volume velocity out of te resonator aperture becomes too low; compromising te effective acoustic interaction wit te enclosure and resulting in insignificant control at te targeted resonance peak. In tis paper, a matematical model describing te acoustic interaction of a resonator and a small enclosure is presented. An analytical solution is obtained on te pressure field inside te enclosure and te radiation of te resonator. Based on te analytical solutions, an energy reduction index describing te strengt of acoustic interaction in te enclosure wit a resonator is defined. eries of numerical simulations are conducted to illustrate te influence of te internal resistance on te energy reduction and on te dissipated and re-radiated energy. Finally, te optimal internal resistance is obtained. xperimental results using one resonator are also carried out and compared wit simulation results. I. ITODUCTIO Low frequency noise inside small enclosures always as always been a concern in many engineering problems. Active noise control is a promising approac in low-frequency noise controls. However, its cost of implementation, robustness of te control system and te
2 ICV4 9- July 7 Cairns Australia requirement of a large space amper its application in small enclosures. Passive tecniques suc as Helmoltz resonators still sow many appealing features in terms of implementation, control quality and cost effectiveness. Te acoustic interaction between an enclosure and wit Helmoltz resonators as been extensively investigated in te past. Fay and cofield [] built an advanced model to explore some underlying pysics of te resonator. Te experimental results on internal resistance [] sowed tat tere exists an optimal damping value wen oter parameters were fixed. Wit a view to broaden te application of resonators, Cummings [] extended Fay and cofield s single resonator and single room mode model to multiple resonator coupled wit multi room-mode. In is multi-mode model, resonators were taken as pseudo point sources. In order to solve te singularity problem caused by point source, te sound pressure of tat pseudo point source at its own location was calculated by te averaged sound pressure at te surface of a small equivalent pulsating spere. Li and Ceng [3] recently proposed a new model to analyze te interaction between acoustic resonator array and a room wit multiple modes. In tat work, te resonator was treated as a point source wit volume velocity influenced by enclosure. Instead of building te linear equation at te aperture of resonators, tey considered te acoustic equilibrium in te volume of room to avoid te singularity problem previously encountered wen te linear equation set was solved. By comparing te teoretical results wit experimental results, te model sows very good approximation to real multi-mode setup. Despite te consistent effort tat was made in Helmoltz resonator design, researcers are still looking into approaces to improve te noise control performance using a Helmoltz resonator array, wic involves eavy experimental measurements on a trial-and-error basis. Tis paper presents a systemic design tool to optimize te internal resistance of a Helmoltz resonator. It is divided into four sections. Te matematical model of acoustic interaction between an enclosure and one single resonator, te teory of energy optimization, and investigation of te internal resistance of te resonator are presented in ection II. ection III gives numerical and experimental results to examine te internal resistance of resonators. Finally, conclusions for optimal design of a resonator are presented in ection IV. II. THOY A general model for describing te acoustic interaction between an enclosure and only one Helmoltz resonator is presented, and ten, te energy optimization teory is given. Formulas of dissipated energy and radiated energy by an acoustic resonator, wic depends on te internal resistance of te resonator, are derived. Trougout te paper, te superscripts and subscripts,, and indicate te variables associated wit nclosure, esonator, and ource, respectively. A. Acoustic interaction between an enclosure and a single Helmoltz resonator Te inomogeneous wave equation governing te pressure field in te enclosure reads:
3 ICV4 9- July 7 Cairns Australia p(,) r t p(,) r t = ρ q (,) r t, () were p(r, t) is te acoustic pressure, q is te volume velocity source strengt density distribution in te volume and surface of te enclosure. Assuming tat a set of armonic sources wit volume velocity source strengt density q q, q,, q are located at te points, r,, r r form te primary sound field in te enclosure, and te single resonator wit volume velocity source strengt density q centered at te point r (center of te resonator aperture) forms te secondary sound field in te enclosure. otice tat te volume velocity directed out of te resonator as te same sign as tat of te primary sound source, i.e., te positive sign is directed out off te source into te enclosure. Te volume velocity source strengt density out of te acoustic resonator is computed from q (t)= p(r,t)/z, in wic Z is defined as te output acoustic impedance at te aperture of te resonator, and p(r,t) is te sound pressure at r provided tat te largest dimension of te resonator aperture is muc smaller tan te sound wavelengt of interest, equation () becomes p ( r, t) δ ( r r ) (,) (,) () ( r n ), () p r t p r t = ρ qn t δ c + r Z n= were δ(r-r ) is a tree dimensional Dirac delta function.. Acoustic pressure p(r,t) can be decomposed on te basis of eigenfunctions of te enclosure: p(r,t)=σψ (t)φ (r),ψ (t) is te t modal response, and φ (r) is te t eigenfunction. ubstituting tis modal expansion into q. () and applying ortogonality properties of te eigenfunctions yield a discrete acoustic equation ( r ) ( r ) cz ϕ cz ϕ ψ () ϕ ( r ) ψ () ( γ ) ψ () q n() t, (=,, 3,, J) (3) n t t + t = V Z Λ V n= Λ were z =ρ c; z is te caracteristic impedance of te fluid, V is te volume of te enclosure, Λ is te mode normalization factor, given by Λ = ϕ () dv V V r, ϕ ( rn ) is te averaged ϕ ( rn ) over te volume of te nt source (first source), andγ is te t eigenvalue of te enclosures, wic olds on te omogeneous wave equation. Te eigenvalue can be expressed as γ = ω + ic, in wic te real part is te angular frequency and te imaginary part is an equivalent ad oc damping coefficient. Assuming all time dependent variables are armonic, i.e. Ψ (t)=p e iωt and q () t Q e ω n i t = n, ten, in te absence of te resonators, te modal response P can be solved from q. (3) as ( P) iωcz = witout resonator ω γ Λ V n= ( ) ( rn ) ϕ Q, (=,, 3,, J). (4) n 3
4 ICV4 9- July 7 Cairns Australia Wen a resonator is installed, te modal response of P can be solved from q. (3) as ( ) ( ) ( ) ϕ ϕ ϕ n Q czω r r r n n= ( ) Λ Z V ( ) P = iωcz ϕ ( ) n Q γ ω γ ω Λ r n n= ( γ ) ω V cz ( ) ω ϕ Λ i r Contribution of te first sound field Z V Λ ( γ ) ω Contribution of te acoustic resonator, (=,, 3,, J). (5) Te first term on te rigt and side of q. (5) is te contribution coming from te primary sound field, and te second term is te effect of inserting an acoustic resonator into te enclosure. B. nergy optimization Wen inserting a designed resonator to target a resonance peak of te enclosure, two new coupled frequencies were produced on te eiter side of original resonance frequency []. It is ard to carry out parameters optimizations of resonator by selecting ust a single frequency as control target since te presence of tose two new coupled frequencies. It is necessary to involve maor frequency components around te resonance peak [4] for optimizations of resonator. To tis end, energy in a frequency band is determined as te obective function for te optimization of resonator location and damping. Te teory of te energy optimization is derived as below. By using p(r,t)=ψ (t)φ (r) andψ (t)=p e iωt, te energy witin a bandwidt [ω,ω] is calculated based on spatially averaged mean-square pressure : ω T [ ω, ω] = [ e( p) ] dt dv d e( ) V P ω ω ω ωi V = Δ T Λ ( ), (6) i were Δω is te numerical integration step. Te energy reduction is defined as te ratio of energy level [ω,ω ]wit resonator to tat [ω,ω ] witout resonator. [ ω ω ] [ ω, ω] [ ω, ω ], = log = log e ( Λ ) P ( ωi) i ( Λ ) e P ( ωi) i, (7) wereω i varies in bandwidt[ω,ω ]. If local noise control is expected, a small V around tis region can be selected to calculate [ω,ω ]. However, if global noise control is expected, te wole enclosure volume sould be selected for te calculation of [ω,ω ]. C. Damping optimization for Helmoltz resonators As mentioned before, inserting a ligtly damped resonator into a ligtly damped enclosure can split te controlled armonic peak to two peaks. In order to enlarge te working bandwidt of a resonator, damping material is introduced to abase tose armonic peaks between te two coupled frequencies. However, excessive damping is not effective, as 4
5 ICV4 9- July 7 Cairns Australia reported in reference []. Terefore, tere exists an optimal damping to implement optimal control performance of a resonator. Te optimization of internal resistance will be carried out based on te energy reduction described in q. (7). Wen te internal resistance of te inserted resonator varies, te modal response P of enclosure varies, correspondingly. Two important types of energy, dissipated energy and radiated energy by resonator, are te key parameters to evaluate te control performance of a resonator, wic are influenced by resonator resistance. Firstly, te relationsip between dissipated energy and resistance is investigated. Te motion of lumped mass in te resonator neck follows te ewton s second law: ρc ( ) ρ L xt () + z Lxt () + xt () = p( r,) t, (8) V were x(t) is te particle displacement in te resonator neck, wic is assumed positive wen it points to te enclosure, and V te body volume of resonator, te cross sectional area of te neck, L te effective neck lengt, and te internal resistance. Assuming te movement of tis lumped mass is a armonic variable, namely x(t)=xe iωt and considering te decomposition of pressure based on mode sape function sown above, te amplitude of displacement is solved from q. (8) as: J X = ϕ ( ) ( ) P ω ω icω r. (9) + ρl = Because energy is defined in a bandwidt [ω,ω], parameters used to calculate energy must be mean square values in te same band. Tus, based on q. (9), te mean square velocity in tis band can be given as: ( ) J ω ( ) ϕ r P ω ω T = [ e( )] ω, ω ω ω ω v = v dt d = Δ T ( ρl ) ω ω cω ( ) + ( ). () By using te definition of mean square velocity in a band in q. (), te averaged dissipated energy of a resonator in te band [ω,ω] is: cω v P d = i = Δω ω, ω ρl J ϕ ( r ) ( ω ), () = ( ω ) ω + ( cω ) were i =z L is te specific acoustic resistance of a resonator. After discussing te dissipated energy by a resonator, te influence of damping effects on te sound field radiated by a resonator will be investigated as te follows. ince sound pressure in te enclosure is te summation of te primary and secondary sources formed by te inserted resonator, te sound pressure radiated from te resonator can be directly analyzed wen te primary field source is turned off. Assuming tat te sound source velocity strengt of primary source is zero, and te lumped mass in te neck of resonator vibrates wit te displacement x(t)=xe iωt [X can be calculated by q. (9)], te modal response is calculated by: 5
6 ( ) ( γ ) ICV4 9- July 7 Cairns Australia ( r ) czω ϕ P = i ϕ ( ) only _ resonator r P, () ω V Z Λ were P is calculated by q. (5). ubstituting (P ) only_resonator into q. (6), te space averaged energy radiated only by a resonator can be calculated. III. UMICAL IMULATIO AD XPIMT umerical simulations on internal resistance are conducted in te following. Te coupling model between te enclosure and a resonator as been previously validated [3]. One room, wit dimensions l x = m, l y =.7 m and l z =. m, was used as enclosure for experiments and simulations. All pysical parameters used are tabulated in TABL I. umerically 6 enclosure modes are used in modal superposition. atural frequencies wit damping are calculated by formulaγ = ω +, in wic te real part ω is te angular frequency of te ic t enclosure mode witout damping and te image part is te t ad oc damping coefficientc ω Q. Q is set to 56 for te rigid mode () and 46 for oter modes. Te = eigenfunctions φ (r ) of te enclosure for termalviscous boundary conditions were presented in reference []. One square source wit dimensions of mm in te x- and y-directions and zero in z-direction was located at (, 55, ) mm to drive te first sound field. One Brüel & Kær Type 489 ½ micropone located at (.84,.3,.6) m was used to measure PL in te enclosure. TABL I. Pysical parameters Pysical parameter value Ambient temperature, T ( o C) 5 peed of sound, c (m/s) 34.3 Density of air, ρ (kg/m 3 ). pecific eat ratio of air, γ.4 Termal conductivity of air, κ [W/(m K)].63 pecific eat at constant pressure of te air, C p [J/(kg K)]. 3 Coefficient of sear viscosity, µ (Pa. s).85-5 For optimization of internal resistance, mode () was selected as targeted mode. Te calculated resonance frequency for tis mode is 9.8Hz and te experimental value is 4Hz. A Helmoltz resonator was designed to target tis mode and fixed at (.,.3, ) m. Te internal diameter of its neck is mm and te pysical neck lengt is 58mm. Te internal diameter and lengt of its body are 74mm and 67mm, respectively. Te optimization frequency band was from Hz to 4Hz. Te Q-factor of te resonator was varied from to, corresponding to a variation range of te resistance from.5 to 5. mks ayls. Te energy reduction was calculated by q. (7) wit different damping values. Te dissipated energy was calculated by q. (). Te radiated energy was calculated by q. (6). FIG.(a) sows te variation of energy reduction versus te internal resistance. Wen te resistance is equal to 4. mks ayls, te energy reduction in te enclosure reaces te maximum value.34 db. Tus tis resistance is te optimal one based on te energy 6
7 ICV4 9- July 7 Cairns Australia approac. FIG.(b) sows te variation of radiated energy from te resonator and FIG.(c) sows te energy dissipated by resonator. It is obvious tat te variation tendency in FIG.(c) is similar to tat in FIG.(a). Terefore, te influence of resonator resistance on te energy reduction in enclosure is dominated by te dissipation ability of resonator. From FIG. (a, b, c), it is observed tat wen te resistance is muc small, te radiation from te resonator is efficient and te dissipated energy is very small. Terefore, more energy is returned back to te enclosure from te resonator aperture at tis circumstance. Wit te increase of resistance, more sound energy is dissipated by te resonator and te radiated energy is small. However, bot te dissipated energy and radiation energy decreases after te internal resistance of te resonator increases over te optimal value. FIG.. (a) nergy reduction in enclosure measured at (.84,.3,.6) m; (b) adiated energy by resonator; (c) Dissipated energy by resonator. FIG. sows te experimental results in terms of sound pressure level (PL) at te micropone position (.84,.3,.6) m wen resonator resistances are equal to 3.6 and 8. mks ayls. Te resonator resistance 3.6 mks ayls is near te simulated optimal resistance 4. mks ayls, and a large PL reduction up to 6.9dB at 4Hz was observed. A control over te off-target frequency in a relatively broadband is also obtained from FIG. since te resonator is coupled wit all enclosure modes. Te internal resistance 8. mks ayls, muc larger tan te optimal resistance, results in a very low vibrating velocity of te lumped mass inside te resonator neck and subsequently lower energy dissipation inside te enclosure and low radiated energy out of te resonator aperture. Only a.66db PL reduction at 4Hz was 7
8 ICV4 9- July 7 Cairns Australia obtained wen te internal resistance of te resonator is 8. mks ayls. FIG.. PL curves at (..84,.,.6) m. xperiment: witout resonator ; wit a damped resonator r = 3.6 mks ayls ; wit a damped resonator, r = 8. mks ayls.. IV. COCLUIO A teoretical metod for te optimal design of a Helmoltz resonator internal-resistance is presented based on te maximization of energy reduction. umerical results sow tat te dissipated energy dominates te control performance of te resonator (a good cutoff between te peak reduction and te resonator working bandwidt), wic greatly depends on bot te internal resistance in te resonator neck and te vibrating velocity of te lumped mass in te neck. A large vibrating velocity of te lumped mass (wen te internal resistance is small) can provide a large re-radiation volume velocity, wic reacts to te primary sound field to sarply reduce te targeted peak by means of acoustic interaction. Te measured results sow tat te optimally designed internal resistance can provide a better control performance for te resonator. Using a Helmoltz resonator wit an optimal internal resistance, a PL reduction up to 6.9dB was obtained around te targeted resonance wit 5 Hz bandwidt. Acknowledgements Te autors wis to acknowledge a grant from esearc Grants Council of Hong Kong pecial Administrative egion, Cina (Proect o. PolyU 537/6) and support by te Central esearc Grant of Te Hong Kong Polytecnic University troug Grant G-Y67. FC [] F. J. Fay and C. cofield, A note on te interaction between a Helmoltz resonator and an acoustic mode of an enclosure, Journal of ound and Vibration, Vol. 7, (98). [] A. Cummings, Te effects of a resonator array on te sound field in a cavity, Journal of ound and Vibration, Vol. 54, 5-44 (99). [3] Deyu. Li and L. Ceng, Teoretical study of noise control in small enclosures using T-saped acoustic resonators, J. Acoust. oc. Am, to be publised. [4] L. Ceng and. Lapointe, Vibration attenuation of panel structures by optimally saped viscoelastic coating wit added weigt considerations, Tin-walled tru (995). 8
A = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationMicrostrip Antennas- Rectangular Patch
April 4, 7 rect_patc_tl.doc Page of 6 Microstrip Antennas- Rectangular Patc (Capter 4 in Antenna Teory, Analysis and Design (nd Edition) by Balanis) Sown in Figures 4. - 4.3 Easy to analyze using transmission
More informationDiffraction. S.M.Lea. Fall 1998
Diffraction.M.Lea Fall 1998 Diffraction occurs wen EM waves approac an aperture (or an obstacle) wit dimension d > λ. We sall refer to te region containing te source of te waves as region I and te region
More informationThe entransy dissipation minimization principle under given heat duty and heat transfer area conditions
Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area
More informationFabric Evolution and Its Effect on Strain Localization in Sand
Fabric Evolution and Its Effect on Strain Localization in Sand Ziwei Gao and Jidong Zao Abstract Fabric anisotropy affects importantly te overall beaviour of sand including its strengt and deformation
More information6. Non-uniform bending
. Non-uniform bending Introduction Definition A non-uniform bending is te case were te cross-section is not only bent but also seared. It is known from te statics tat in suc a case, te bending moment in
More informationAPPENDIX A. Lukas Novotny. University ofrochester. Rochester, NY Bert Hecht. Institute of Physical Chemistry
OPTICS ON THE NANOMETER SCALE 1 : APPENDIX A Lukas Novotny Te Institute of Optics University ofrocester Rocester, NY 14627 Bert Hect Institute of Pysical Cemistry Swiss Federal Institute of Tecnology,
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationAN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES
Ice in te Environment: Proceedings of te 1t IAHR International Symposium on Ice Dunedin, New Zealand, nd t December International Association of Hydraulic Engineering and Researc AN ANALYSIS OF AMPLITUDE
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationTheoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump
TECHNICAL PAPE Teoretical Analysis of Flow Caracteristics and Bearing Load for Mass-produced External Gear Pump N. YOSHIDA Tis paper presents teoretical equations for calculating pump flow rate and bearing
More informationPECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER
Materials Pysics and Mecanics (5) 5- Received: Marc 7, 5 PECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER Т.I. Belyankova *, V.V. Kalincuk Soutern Scientific Center of Russian
More informationDesalination by vacuum membrane distillation: sensitivity analysis
Separation and Purification Tecnology 33 (2003) 75/87 www.elsevier.com/locate/seppur Desalination by vacuum membrane distillation: sensitivity analysis Fawzi Banat *, Fami Abu Al-Rub, Kalid Bani-Melem
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationModel development for the beveling of quartz crystal blanks
9t International Congress on Modelling and Simulation, Pert, Australia, 6 December 0 ttp://mssanz.org.au/modsim0 Model development for te beveling of quartz crystal blanks C. Dong a a Department of Mecanical
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare ttp://ocw.mit.edu 5.74 Introductory Quantum Mecanics II Spring 9 For information about citing tese materials or our Terms of Use, visit: ttp://ocw.mit.edu/terms. Andrei Tokmakoff, MIT
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationWhy gravity is not an entropic force
Wy gravity is not an entropic force San Gao Unit for History and Pilosopy of Science & Centre for Time, SOPHI, University of Sydney Email: sgao7319@uni.sydney.edu.au Te remarkable connections between gravity
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationEffects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates with Ramped Wall Temperature
Volume 39 No. February 01 Effects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates wit Ramped Wall Temperature S. Das Department of Matematics University of Gour Banga Malda 73
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationGrade: 11 International Physics Olympiad Qualifier Set: 2
Grade: 11 International Pysics Olympiad Qualifier Set: 2 --------------------------------------------------------------------------------------------------------------- Max Marks: 60 Test ID: 12111 Time
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationQuantum Theory of the Atomic Nucleus
G. Gamow, ZP, 51, 204 1928 Quantum Teory of te tomic Nucleus G. Gamow (Received 1928) It as often been suggested tat non Coulomb attractive forces play a very important role inside atomic nuclei. We can
More informationContinuous Stochastic Processes
Continuous Stocastic Processes Te term stocastic is often applied to penomena tat vary in time, wile te word random is reserved for penomena tat vary in space. Apart from tis distinction, te modelling
More informationChapter 4 Optimal Design
4- Capter 4 Optimal Design e optimum design of termoelectric devices (termoelectric generator and cooler) in conjunction wit eat sins was developed using dimensional analysis. ew dimensionless groups were
More informationBallistic electron transport in quantum point contacts
Ballistic electron transport in quantum point contacts 11 11.1 Experimental observation of conductance quantization Wen we discussed te self-consistent calculation of te potential and te modes in an infinite
More informationProblem Solving. Problem Solving Process
Problem Solving One of te primary tasks for engineers is often solving problems. It is wat tey are, or sould be, good at. Solving engineering problems requires more tan just learning new terms, ideas and
More informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
More informationDevelopment of new and validation of existing convection correlations for rooms with displacement ventilation systems
Energy and Buildings 38 (2006) 163 173 www.elsevier.com/locate/enbuild Development of new and validation of existing convection correlations for rooms wit displacement ventilation systems Atila Novoselac
More informationHT TURBULENT NATURAL CONVECTION IN A DIFFERENTIALLY HEATED VERTICAL CHANNEL. Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008
Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008 August 10-14, 2008, Jacksonville, Florida USA Proceedings of HT2008 2008 ASME Summer Heat Transfer Conference August 10-14, 2008, Jacksonville,
More informationExercise 19 - OLD EXAM, FDTD
Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms
More informationJournal of Mechanical Science and Technology 23 (2009) 2072~2084. M. M. Najafizadeh * and M. R. Isvandzibaei
Journal of Mecanical Science and Tecnology (9 7~84 Journal of Mecanical Science and Tecnology www.springerlink.com/content/78-494x DOI.7/s6-9-4- Vibration of functionally graded cylindrical sells based
More informationNon-linear Analysis Method of Ground Response Using Equivalent Single-degree-of-freedom Model
Proceedings of te Tent Pacific Conference on Eartquake Engineering Building an Eartquake-Resilient Pacific 6-8 November 25, Sydney, Australia Non-linear Analysis Metod of Ground Response Using Equivalent
More informationFlapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using the differential transform method
Meccanica 2006) 41:661 670 DOI 10.1007/s11012-006-9012-z Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using te differential transform metod Ozge Ozdemir Ozgumus
More informationA Reconsideration of Matter Waves
A Reconsideration of Matter Waves by Roger Ellman Abstract Matter waves were discovered in te early 20t century from teir wavelengt, predicted by DeBroglie, Planck's constant divided by te particle's momentum,
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationProblem Set 4 Solutions
University of Alabama Department of Pysics and Astronomy PH 253 / LeClair Spring 2010 Problem Set 4 Solutions 1. Group velocity of a wave. For a free relativistic quantum particle moving wit speed v, te
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationNumerical Simulations of the Physical Process for Hailstone Growth
NO.1 FANG Wen, ZHENG Guoguang and HU Zijin 93 Numerical Simulations of te Pysical Process for Hailstone Growt FANG Wen 1,3 ( ), ZHENG Guoguang 2 ( ), and HU Zijin 3 ( ) 1 Nanjing University of Information
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationHARMONIC ALLOCATION TO MV CUSTOMERS IN RURAL DISTRIBUTION SYSTEMS
HARMONIC ALLOCATION TO MV CUSTOMERS IN RURAL DISTRIBUTION SYSTEMS V Gosbell University of Wollongong Department of Electrical, Computer & Telecommunications Engineering, Wollongong, NSW 2522, Australia
More informationComment on Experimental observations of saltwater up-coning
1 Comment on Experimental observations of saltwater up-coning H. Zang 1,, D.A. Barry 2 and G.C. Hocking 3 1 Griffit Scool of Engineering, Griffit University, Gold Coast Campus, QLD 4222, Australia. Tel.:
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationChapter 2 Ising Model for Ferromagnetism
Capter Ising Model for Ferromagnetism Abstract Tis capter presents te Ising model for ferromagnetism, wic is a standard simple model of a pase transition. Using te approximation of mean-field teory, te
More informationCALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY
CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY L.Hu, J.Deng, F.Deng, H.Lin, C.Yan, Y.Li, H.Liu, W.Cao (Cina University of Petroleum) Sale gas formations
More informationFinite Element Analysis of J-Integral for Surface Cracks in Round Bars under Combined Mode I Loading
nternational Journal of ntegrated Engineering, Vol. 9 No. 2 (207) p. -8 Finite Element Analysis of J-ntegral for Surface Cracks in Round Bars under Combined Mode Loading A.E smail, A.K Ariffin 2, S. Abdulla
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
More informationCFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E.
CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E. Publised in: IEA Annex 41 working meeting, Kyoto, Japan Publised:
More informationA Multiaxial Variable Amplitude Fatigue Life Prediction Method Based on a Plane Per Plane Damage Assessment
American Journal of Mecanical and Industrial Engineering 28; 3(4): 47-54 ttp://www.sciencepublisinggroup.com/j/ajmie doi:.648/j.ajmie.2834.2 ISSN: 2575-679 (Print); ISSN: 2575-66 (Online) A Multiaxial
More informationCHAPTER 2 MODELING OF THREE-TANK SYSTEM
7 CHAPTER MODELING OF THREE-TANK SYSTEM. INTRODUCTION Te interacting tree-tank system is a typical example of a nonlinear MIMO system. Heiming and Lunze (999) ave regarded treetank system as a bencmark
More informationLarge eddy simulation of turbulent flow downstream of a backward-facing step
Available online at www.sciencedirect.com Procedia Engineering 31 (01) 16 International Conference on Advances in Computational Modeling and Simulation Large eddy simulation of turbulent flow downstream
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationThe development of contact and noncontact technique to study the heat dissipation in metals under loading
Te development of contact and noncontact tecnique to study te eat dissipation in metals under loading More info about tis article: ttp://www.ndt.net/?id=73 Abstract * ICMM UB RAS, Ac. Koroleva Str., 63
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More informationEFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS
Statistica Sinica 24 2014, 395-414 doi:ttp://dx.doi.org/10.5705/ss.2012.064 EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Jun Sao 1,2 and Seng Wang 3 1 East Cina Normal University,
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More information10.1 VIBRATIONAL RELAXATION *
Andrei Tokmakoff, MIT Department of Cemistry, 3//009 p. 0-0. VIRATIONAL RELAXATION * Here we want to address ow a quantum mecanical vibration undergoes irreversible energy dissipation as a result of interactions
More informationVibro-Acoustics of a Brake Rotor with Focus on Squeal Noise. Abstract
Te International Congress and Exposition on Noise Control Engineering Dearbor MI, USA. August 9-, Vibro-Acoustics of a Brake Rotor wit Focus on Sueal Noise H. Lee and R. Sing Acoustics and Dynamics Laboratory,
More informationInternational Power, Electronics and Materials Engineering Conference (IPEMEC 2015)
International Power, Electronics and Materials Engineering Conference (IPEMEC 2015) Calculation of Stray Losses and Temperature Distribution in Power Transformer Using Coupled Electro-Termal Field Analysis
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationChapters 19 & 20 Heat and the First Law of Thermodynamics
Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,
More informationElmahdy, A.H.; Haddad, K. NRCC-43378
Experimental procedure and uncertainty analysis of a guarded otbox metod to determine te termal transmission coefficient of skyligts and sloped glazing Elmady, A.H.; Haddad, K. NRCC-43378 A version of
More informationStudy of Convective Heat Transfer through Micro Channels with Different Configurations
International Journal of Current Engineering and Tecnology E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rigts Reserved Available at ttp://inpressco.com/category/ijcet Researc Article Study of
More informationACCURATE SYNTHESIS FORMULAS OBTAINED BY USING A DIFFERENTIAL EVOLUTION ALGORITHM FOR CONDUCTOR-BACKED COPLANAR WAVEG- UIDES
Progress In Electromagnetics Researc M, Vol. 10, 71 81, 2009 ACCURATE SYNTHESIS FORMULAS OBTAINED BY USING A DIFFERENTIAL EVOLUTION ALGORITHM FOR CONDUCTOR-BACKED COPLANAR WAVEG- UIDES S. Kaya, K. Guney,
More informationDistribution of reynolds shear stress in steady and unsteady flows
University of Wollongong Researc Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 13 Distribution of reynolds sear stress in steady
More information1 Power is transferred through a machine as shown. power input P I machine. power output P O. power loss P L. What is the efficiency of the machine?
1 1 Power is transferred troug a macine as sown. power input P I macine power output P O power loss P L Wat is te efficiency of te macine? P I P L P P P O + P L I O P L P O P I 2 ir in a bicycle pump is
More informationConsider the element shown in Figure 2.1. The statement of energy conservation applied to this element in a time period t is that:
. Conduction. e General Conduction Equation Conduction occurs in a stationary medium wic is most liely to be a solid, but conduction can also occur in s. Heat is transferred by conduction due to motion
More informationAn optimum design of robotic food handling by using Burger model
DOI 10.1007/s11370-008-0032-5 ORIGINAL RESEARCH PAPER An optimum design of robotic food andling by using Burger model Naoki Sakamoto Mitsuru Higasimori Tosio Tsuji Makoto Kaneko Received: 28 February 2008
More informationNotes: Most of the material in this chapter is taken from Young and Freedman, Chap. 12.
Capter 6. Fluid Mecanics Notes: Most of te material in tis capter is taken from Young and Freedman, Cap. 12. 6.1 Fluid Statics Fluids, i.e., substances tat can flow, are te subjects of tis capter. But
More informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More informationChapter 6: Adiabatic reaction dynamics and Transition State Theory Bimolecular reactions Transition State theory... 72
Winter 204 Capter 6: Adiabatic reaction dynamics and Transition State Teory... 70 imolecular reactions... 7 Transition State teory... 72 Derivation of universal factor T... 74 Kinetic Isotope effect...
More informationAnalysis of Static and Dynamic Load on Hydrostatic Bearing with Variable Viscosity and Pressure
Indian Journal of Science and Tecnology Supplementary Article Analysis of Static and Dynamic Load on Hydrostatic Bearing wit Variable Viscosity and Pressure V. Srinivasan* Professor, Scool of Mecanical
More informationAN IMPROVED WEIGHTED TOTAL HARMONIC DISTORTION INDEX FOR INDUCTION MOTOR DRIVES
AN IMPROVED WEIGHTED TOTA HARMONIC DISTORTION INDEX FOR INDUCTION MOTOR DRIVES Tomas A. IPO University of Wisconsin, 45 Engineering Drive, Madison WI, USA P: -(608)-6-087, Fax: -(608)-6-5559, lipo@engr.wisc.edu
More informationThe Basics of Vacuum Technology
Te Basics of Vacuum Tecnology Grolik Benno, Kopp Joacim January 2, 2003 Basics Many scientific and industrial processes are so sensitive tat is is necessary to omit te disturbing influence of air. For
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationComputational Method of Structural Reliability Based on Integration Algorithms
Sensors & ransducers, Vol. 54, Issue 7, July 03, pp. 5-59 Sensors & ransducers 03 by IFSA ttp://www.sensorsportal.com Computational Metod of Structural Based on Integration Algoritms * Cong Cen, Yi Wan
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationRotational-Vibrational Levels of Diatomic Molecules Represented by the Tietz-Hua Rotating Oscillator
J. Pys. Cem. A 1997, 101, 1595-160 1595 Rotational-Vibrational Levels of Diatomic Molecules Represented by te Tietz-Hua Rotating Oscillator Josep A. Kunc* Departments of Aerospace Engineering Pysics, UniVersity
More informationA finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations
ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationRobotic manipulation project
Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationThe gyrotropic characteristics of hexaferrite ceramics.
Te gyrotropic caracteristics of exaferrite ceramics. D H Martin, Bin Yang and R S Donnan. July 6. 1. Introduction. Hig-performance non-reciprocal devices aving quasi-optical structures and operating at
More informationLast lecture (#4): J vortex. J tr
Last lecture (#4): We completed te discussion of te B-T pase diagram of type- and type- superconductors. n contrast to type-, te type- state as finite resistance unless vortices are pinned by defects.
More informationarxiv: v1 [math.na] 8 Nov 2016
SEOND-HRMONI IMGING IN RNDOM MEDI LILIN BORE, WEI LI, LEXNDER V. MMONOV, ND JOHN. SHOTLND arxiv:6.774v [mat.n] 8 Nov 6 bstract. We consider te problem of optical imaging of small nonlinear scatterers in
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationUNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics
UNIVERSITY O SASKATCHEWAN Department of Pysics and Engineering Pysics Pysics 117.3 MIDTERM EXAM Regular Sitting NAME: (Last) Please Print (Given) Time: 90 minutes STUDENT NO.: LECTURE SECTION (please ceck):
More informationarxiv: v3 [math.na] 15 Dec 2009
THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING M.A. EBRAHIMI, MICHAEL HOLST, AND EVELYN LUNASIN arxiv:91.4548v3 [mat.na] 15 Dec 9 ABSTRACT. In tis paper we investigate te use of te D Navier-Stokes-Voigt
More information