Acoustic-Structure Simulation of Exhaust Systems by Coupled FEM and Fast BEM
|
|
- Nelson Knight
- 6 years ago
- Views:
Transcription
1 Acoustic-Structure Simuation of Exhaust Systems by Couped FEM and Fast BEM Lothar Gau, Michae Junge Institute of Appied and Experimenta Mechanics, University of Stuttgart Pfaffenwadring 9, 755 Stuttgart, Germany ABSTRACT The sound properties of modern cars pay an increasingy important roe on the overa impression of a vehice. One of the major sound sources is hereby the exhaust system, which is exposed to arge pressure pusations resuting from the periodicay bown out exhaust gas of the car engine. These arge pressure pusations ead to vibrations of the thin was of the exhaust system and thus create another sound source which is additionay contributing to the sound radiation of the exhaust system besides the sound at the orifice. In this work, the sound radiation of an associated simpified exhaust system is anayzed. The sound transmission of the interior exhaust system is simuated by using conforming couping of FE-FE-modes. For the subsequent soution of the infinite exterior, probem the Fast-Mutipoe BEM is empoyed enabing a very fast and efficient computation. This work focuses on the contribution of the vibrating structure to the overa sound radiation. It is found that cose to some resonance frequencies of the in-vacuo structure there is a significant contribution of the vibrating structure to the overa radiated sound power. Nomencature F s, F f oad vectors u dispacement degree of freedom I a acoustic intensity vector W radiated soundpower I unity matrix v n norma veocity M s, K s mass and stiffness matrix of the structure Z specific acoustic impedance M f, K f mass and stiffness matrix of the fuid Γ cosed surface of radiating object M fs, K fs couping matrices κ wave number n unit norma vector on boundary surface ρ, c density and wave speed of fuid P fundamenta soution (V q)(x) singe ayer potentia p pressure degree of freedom (Kp)(x) doube ayer potentia q acoustic fux (K q)(x) adjoint doube ayer potentia R compex refection coefficient (Dp)(x) hyper-singuar operator 1 Introduction On a car there are severa sound sources each dominating at a different speed. At ow speeds up to 5 km/h the major sound source is the car engine. In most cases, the engine noise is very efficienty reduced by an exhaust system. However, the arge pressure pusations on the outet of the piston engine may ead to vibrations of the thin was of the exhaust system. These vibrations in the exhaust structure are additionay contributing to the sound radiation of the exhaust system. Besides, the propagation of the structure-borne sound induced in the was can acousticay short-circuit the siencer instaed in the duct as reported for inet and outet ducts of arge gas turbines by [1]. The excitation of structura vibration modes by the acoustic path is investigated and the sound radiation of an ideaized car exhaust system is anayzed. A specia focus is set on the contribution of the vibrating structure to the overa sound radiation of the exhaust system.
2 2 Simuation mode The simuation mode of the ideaized exhaust system consists of two sub-modes. In the first step, a fuy couped FE-FEmode is used to compute the interior acoustic probem of the exhaust system. The fuid-structure interface accounts for the infuence of the acoustics on the fexibe was. By means of this mode, the excitation of structura vibration modes by the acoustic path is computed. In the next step, the cacuated norma veocities on the boundary of the structura parts of the exhaust system and the pressure on the outet are used as Neumann and Dirichet boundary data respectivey for the BE-mode in order to sove the exterior acoustic probem. 2.1 Ideaized exhaust system In this work, an ideaized exhaust system is used. It consists of one expansion chamber, an ebow, an exhaust pipe and a tai pipe as shown in Fig. 1. The system is fixed at its ends by circuar soft rings. It is excited on the inet with a uniform, harmonic pusating pressure p = ˆpe jωt. This mode does not consider mean fow and vortices. It assumes that the inear acoustics theory is appicabe. The tota pipe ength is 2.4 m, the pipe radius is chosen as r =.34 m and the wa thickness is t w =.1 m. 3 Interior acoustic probem FSI-couped undamped probems are formuated by FE-discretization as presented in [2] [ ][ü ] [ ] [ ] [ Ms Ks K + fs u F s] =, (1) M f p K f p M fs with the noda dispacement degrees of freedom u, and the noda acoustic pressure degrees of freedom p. The mass couping matrix M fs and the stiffness couping matrix K fs are reated by M fs = ρk T fs. At the outet, a radiation impedance condition is appied using the resuts derived in [3]: The specific acoustic impedance Z is expressed in terms of the compex refection coefficient R by F f Z = p v n = ρc 1 + R 1 R. (2) The compex refection coefficient R is characterized by its magnitude R, the end correction factor δ and the wave number k R = R e j(π 2kδ). (3) Figure 1: Structure of the ideaized exhaust system
3 R δ f f Figure 2: Radiation impedance condition: Magnitude and end correction factor of a circuar pipe (radius r =.34 m) according to [3]. The vaues of R and δ for the system at hand are shown in Fig. 2. For ow frequencies the magnitude of the refection coefficient is cose to unity whie its phase ange is amost π, hence an incoming wave is amost fuy refected at the tai pipe s end. In the FE-program ANSYS, the compex impedance condition (2) is modeed in the harmonic case by adding massess surface eements onto the outet, assigned with viscosity and added mass per unit area. In the simuation the expansion chamber shows a typica transmission oss behavior with a maximum vaue of 15 db as shown in Fig. 3. These resuts are in good agreement with those presented in [4]. Transmission oss [db] f [Hz] Figure 3: Transmission oss of ideaized exhaust system with a maximum vaue of 15 db. The moda anaysis of the in-vacuo system (exhaust system without fuid eements) shows resonance frequencies as indicated by the vertica ines in the ower pot of Fig. 6. The eigenmode of the exhaust system at f = Hz is depicted in Fig. 4. It is evident that the exhaust pipe and the expansion chamber are moving in opposite phase. The ampitude of the taipipe is much ower than the ampitude of the exhaust pipe and the expansion chamber.
4 4 Exterior acoustic probem Figure 4: Eigenmode of the in-vacuo exhaust system at f = Hz Having soved the interior acoustic probem, the computed dispacements on the pipe and the pressure vaues on the outet are used as boundary data for the exterior acoustic probem, yieding vaues for the pressure and acoustic fux respectivey on the surface of the exhaust system. This probem is computed by means of the Fast Mutipoe Mutieve Boundary Eement Method (FMM-BEM) [5] which wi be descriped briefy in the next two subsections. 4.1 The Boundary Eement Method (BEM) The Boundary Eement Method is we suited for acoustic computations, offering high accuracy and a simpe mesh generation, especiay for modeing exterior acoustic probems [6]. Having once computed the acoustic pressure and the acoustic fux on the boundary of the acoustic domain, the acoustic fied pressure and the fied fux are computed at any arbitrary fied point with ow numerica costs. For time-harmonic behavior, the pressure p in an acoustic fied is governed by the Hemhotz equation 2 p(x) + κ 2 p(x) =, (4) where κ = ω/c is the acoustic wavenumber and c is the sound veocity. The acoustic fux q on the boundary is prescribed as q(x) := p(x), x Γ, (5) n which is proportiona to the surface veocity and thus to the dispacement ampitudes of a vibrating structure. For the exterior probem, the Hemhotz equation is soved in an exterior domain Ω e, with the boundary Γ = Ω e. The fundamenta soution P (x,y) for the Hemhotz equation with spatia Dirac source is given by P (x,y) = 1 e iκr 4π r, (6) where r = x y is the Eucidean distance between the fied point x and the oad point y. Appication of Green s second theorem to the weak form of the Hemhotz equation Ω e P (x,y) ( 2 p(x) + κ 2 p(x) ) dω e = (7) yieds the representation formua p(x) = Γ P P (x,y) (x,y)q(y)ds y + p(y)ds y, x Ω e. (8) Γ n y By shifting the fied point onto the smooth boundary, one obtains the boundary integra equation p(x) = 1 2 p(x) P P (x,y) (x,y)q(y)dγ y + p(y)dγ y, x Γ. (9) Γ Γ n }{{} y }{{} (V q)(x) (Kp)(x)
5 y a z a D z b x b Figure 5: The action of sources at y a are bunded in z a, then transated to another center z b, and evauated at the fied points x b. The expressions (V q)(x) and (Kp)(x) are caed the singe and doube ayer potentias, respectivey. The hypersinguar boundary integra equation [5] is obtained by taking the norma derivative and the imit x Γ of Eq. (8) q(x) = 1 2 q(x) P (x,y) q(y)dγ y Γ n } x {{} (K q)(x) + 2 P (x,y) p(y)dγ y, Γ n x n y }{{} x Γ, () (Dp)(x) where (K q)(x) and (Dp)(x) are the adjoint doube ayer potentia and the hyper-singuar operator, respectivey. A more detaied description of the Boundary Eement Method for acoustics can be found in [7, 8, 9]. The standard BEM has one significant drawback. For frequencies that correspond to the eigenfrequencies of the associated interior probem, the soution is not unique. Burton and Mier [] proposed a combination of the standard and the hyper-singuar boundary integra equation (Eqs. (9) and ()) to avoid this probem at the so caed critica frequencies. Using an isoparametric Gaerkin formuation with piecewise inear shape functions for the pressure and constant shape functions for the acoustic fux, the Neumann probem considered in this paper is written as Burton-Mier BEM formuation by superposition and i κ -weighting ( 1 2 I + K + i ) κ D p = ( V i 2κ I i ) κ K q. (11) 4.2 Fast Mutipoe-BEM (FBEM) A major imitation of cassica BEM formuations is the numerica compexity of order O(N 2 ) for setting up and storing the fuy popuated system matrices. Thus the appication of cassica BEM methods is imited to the ow frequency regime with moderate numbers of unknowns. Since modern iterative sovers ike GMRES [11] ony require the evauation of matrix-vector products, the fuy popuated system matrices do not need to be set up. The key idea of the fast mutipoe method is the efficient evauation of the matrixvector products by using a series expansion of the fundamenta soution, which combines the effect of sources far away from a fiedpoint in a farfied representation. Foowing Rokhin [12], the mutipoe expansion reads with d = x b z b + z a y a and D as depicted in Figure 5 e iκ x y x y = eiκ D+d D + d = iκ (2 + 1)( 1) j (κ d )h (1) (κ D )P ( ˆD ˆd), D > d, (12) = where j ( ) denotes spherica Besse functions, h ( ) Hanke function, P ( ) Legendre poynomias and ˆ ( ) indicates normaized vectors. A diagona form of the mutipoe expansion is obtained by using the orthonormaity of the Legendre poynomias on the unit sphere and the expansion of spherica waves e iκ D+d D + d = iκ 4π = (2 + 1)i h (1) (κ D ) e iκd s P (s ˆD)ds. (13) S 2 The summation must be truncated to the expansion ength L, since the Hanke functions diverge for arge. Therefor,
6 transation operators are defined M L (s,d) = L = (2 + 1)i h (1) (κ D )P (s ˆD), (14) which ony depend on the distance vector D. Thus, a hierarchica custering strategy is appied to achieve an efficient mutieve scheme. A custers which are within some distance of a custer form the nearfied. Custers, whose father custers fufi the nearfied condition, but are not in each others nearfied themseves, form the interaction ist. In the foowing, the evauation of the matrix-vector product for the singe ayer potentia is discussed. For the singe ayer potentia with constant shape functions, the v-th component of the matrix-vector product is now approximated by G ν v ν (V nearfied u) ν + w ν,j ν j=1 µ farfied w µ,i µ P (x ν,j,y µ,i ), (15) u µ G µ i=1 } {{ } evauation with mutipoe method where µ, w µ,i and y µ,i are the Jacobi determinant, the Gauss weights, and the integration points for eement µ, respectivey. The terms are defined anaogous for the eement ν. The expression on the right side are evauated efficienty using the mutipoe agorithm, which can be sketched as foows: First, a farfied signature F γ max (s) is computed for a custers γ on the owest eve. The farfied signature is transated to the interaction ist using the transation operators (14) N γ (s) = interactionist M L (s,d)f γ (s). (16) At the next step, F γ is shifted to the center of the father custer. The ast two steps must be repeated unti the interaction ist is empty. In the downward pass, the nearfied signatures N γ (s) are shifted to the chid custers. Finay, the soution in the integration points is recovered on the owest eve. For the compete mutipoe agorithm, the computing cost is of order O(N og 2 N) [5]. Thus the FMM agorithm enabes a fast and efficient computation of arge scae BEM systems up to more than degrees of freedom using standard desktop PC. 4.3 Sound radiation The sound radiation of the exhaust system is evauated in a post-processing step by separatey computing the radiated sound power through the cross-section of the outet W ΓD and through the was of the exhaust system W ΓN W Γ = W ΓD + W ΓN = I a n dγ + Γ D I a n dγ Γ N, (17) where I = pv is the oca intensity vector on the surface and n the unit norma vector pointing outwards from the cosed surface Γ = Γ D + Γ N. The soid ine in the upper pot of Fig. 6 shows the tota radiated sound power W Γ scaed to its maximum vaue. The oca maxima and minima every 4-6 Hz resut from the constructive and destructive interference of the refected wave on the outet. The dash-dotted ine indicates the magnitude of the sound power radiated by the vibrating structure W ΓN, whie in the ower pot of Fig. 6 the ratio of the magnitude of W ΓN and W Γ is depicted. Both pots ceary show that W ΓN does not contribute significanty to the tota radiated sound power for the vast majority of frequencies. Yet, for some frequencies marked with an asterisk, energy is absorbed by the was, though the vaues are very sma. In contrast to this, it is found that for some frequencies cose to resonance frequencies of the in-vacuo exhaust system the ratio of the magnitude of W ΓN and W Γ is higher than %. Thus W ΓN significanty contributes to the tota radiated sound power. Some of the acoustic energy inserted into the system by the appied uniform, harmonic pusating acoustic pressure is transferred via the fuid-structure interface onto the was of the exhaust system, generating structura vibrations. So far, the radiated sound power of the exhaust system has been considered. In this paragraph the sound fied of the exhaust system at a certain frequency wi be investigated. Figure 7 shows the computed sound fied at a frequency f = 288 Hz. The vertica pane cose to the outet of the exhaust system shows the circuar radiation of the sound at the outet. In this area the tai pipe sound is the dominating source. The vertica pane in the rear is chosen at a distance of 2 cm from the ongitudina axis of the exhaust system, giving a good impression how the sound pressure fied might ook ike at the underfoor of a car. Two dark spots (one bue spot and one red spot) are observed representing areas of amost the same magnitude but opposite sign resuting from the motion of the exhaust system simiar to the eigenmode shown earier in Fig. 4.
7 normaized sound power -5 WΓ WΓN f [Hz] WΓN /WΓ -1 power ratio f [Hz] Figure 6: Top: Normaized radiated sound power of the exhaust system W Γ and the contribution of the was of exhaust structure W ΓN. Bottom: Ratio of W ΓN to W Γ. The vertica ines represent the resonance frequencies of the in-vacuo exhaust system. W ΓN does not contribute significanty to the overa radiated sound power W Γ, except for the case of pressure induced resonance of the exhaust structure. Figure 7: Sound fied of vibrating exhaust system at f = 288 Hz.
8 5 Concusion In this work a method to efficienty compute the sound radiation of an exhaust system is presented. It consists of a consecutive computation of the acoustic interior and exterior probem. The simuation resuts of the acoustic interior probem soved with a FE-FE-formuation are used as boundary data for a BE-formuation of the acoustic exterior probem. For the computation of the acoustic exterior probem, the Fast Mutipoe Boundary Eement Method is appied enabing a fast and efficient computation of arge scae systems. For the presented ideaized exhaust system the simuation resuts suggest a minor contribution of the sound generated by the vibrating structure to the tota radiated sound power. However, cose to resonance frequencies of the in-vacuo system an energy exchange on the FE-FE interface is observed and structura vibration of the structure of the exhaust system is induced. In the case of a more reaistic exhaust system, especiay for the case of a more effective muffer with a much higher transmission oss, the infuence of the sound radiated by the vibrating structure wi pay a more important roe. In this work ony the excitation of structura vibration by the acoustic path is considered. Vibrations induced by the movement of the attached engine are not taken into account. The presented method proves viabe for the computation of the sound radiation of exhaust systems and can easiy be extended to account for the additiona effects mentioned. 6 Acknowedgement Funding of this project by the Friedrich-und-Eisabeth-Boysen-Stiftung is gratefuy acknowedged. References [1] A. Facian, A. Nisson, L. Feng, and E. Nisson. Propagation of structure-borne sound in siencers used in power pants. In Proceedings of the Eeventh Internationa Congress on Sound and Vibration, pages 69 76, 24. [2] O.C Zienkiewicz and R.L. Tayor. The Finite Eement Method, voume 1. Butterworth-Heinemann, Oxford, 5th edition, 2. [3] H. Levine and J. Schwinger. On the radiation of sound from an unfanged circuar pipe. Physica Review, 73(4):383 46, [4] M.L. Munja. Acoustics of Ducts and Muffers. John Wiey & Sons, New York, [5] M. Fischer. The Fast Mutipoe Boundary Eement Method and its Appications to Structure-Acoustic Fied Interaction. PhD thesis, University of Stuttgart, 24. [6] L. Gau and Wagner M. Kög M. Boundary Eement Methods for Engineers and Scientists. Springer, Berin, 23. [7] O. v. Estorff, editor. Boundary Eements in Acoustics: Advances and Appications. WIT Press, Southampton, UK, 2. [8] T.W. Wu. Boundary Eement Acoustics: Fundamentas and Computer Codes. WIT Press, Southampton, UK, 2. [9] L. Gau and M. Fischer. Large Scae Simuations of Acoustic Structure Interaction Using the Fast Mutipoe BEM. In R. Hemig, A. Mieke, and B. Womuth, editors, Mutified Probems in Soid and Fuid Mechanics, pages , 26. [] A.J. Burton and G.F. Mier. The appication of the integra equation method to the numerica soution of some exterior boundary-vaue probems. Proc. Roya Soc. Lond, A, 323:21 2, [11] Y. Saad. Iterative Methods for Sparse Linear Systems. SIAM, Society for Industria and Appied Mathematics, Phiadephia, 23. [12] V. Rokhin. Diagona forms of transation operators for the hemhotz equation in three dimensions. Appied and Computationa Harmonic Anaysis, 1:82 93, 1993.
Proceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More informationIn-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017
In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationHILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT
Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is
More informationSEA Subsystem Attribute Correction Based on Stiffness Multipliers
Internationa Conference on Eectrica, Mechanica and Industria Engineering (ICEMIE 016) SEA Subsystem Attribute Correction Based on Stiffness Mutipiers Jie Gao1,* and Yong Yang 1 Coege of hysics and Information
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More informationLegendre Polynomials - Lecture 8
Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha
More informationFormulas for Angular-Momentum Barrier Factors Version II
BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A
More informationSECTION A. Question 1
SECTION A Question 1 (a) In the usua notation derive the governing differentia equation of motion in free vibration for the singe degree of freedom system shown in Figure Q1(a) by using Newton's second
More informationDYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE
3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses
More informationTM Electromagnetic Scattering from 2D Multilayered Dielectric Bodies Numerical Solution
TM Eectromagnetic Scattering from D Mutiayered Dieectric Bodies Numerica Soution F. Seydou,, R. Duraiswami, N.A. Gumerov & T. Seppänen. Department of Eectrica and Information Engineering University of
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More informationLECTURE 10. The world of pendula
LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive
More informationModal analysis of a multi-blade system undergoing rotational motion
Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa
More informationPhysics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27
Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that
More informationApproximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation
Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische
More informationModule 22: Simple Harmonic Oscillation and Torque
Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque
More informationarxiv: v1 [hep-th] 10 Dec 2018
Casimir energy of an open string with ange-dependent boundary condition A. Jahan 1 and I. Brevik 2 1 Research Institute for Astronomy and Astrophysics of Maragha (RIAAM, Maragha, Iran 2 Department of Energy
More informationEffect of Oxygen Injection into Argon Induction Plasmas on Chemically Non-Equilibrium Conditions
Proceedings of 17th Internationa Symposium on Pasma Chemistry, Toronto, Canada, August 7-12, 25 Effect of Oxygen Injection into Argon Induction Pasmas on Chemicay Non-Equiibrium Conditions Nobuhiko Atsuchi
More informationIntroduction to LMTO method
1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations
More informationRadiation Fields. Lecture 12
Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,
More informationUNCOMPLICATED TORSION AND BENDING THEORIES FOR MICROPOLAR ELASTIC BEAMS
11th Word Congress on Computationa Mechanics WCCM XI 5th European Conference on Computationa Mechanics ECCM V 6th European Conference on Computationa Fuid Dynamics ECFD VI E. Oñate J. Oiver and. Huerta
More informationNonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationHigh-order approximations to the Mie series for electromagnetic scattering in three dimensions
Proceedings of the 9th WSEAS Internationa Conference on Appied Mathematics Istanbu Turkey May 27-29 2006 (pp199-204) High-order approximations to the Mie series for eectromagnetic scattering in three dimensions
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationRelated Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage
Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is
More informationSEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More informationSupporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers
Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationMore Scattering: the Partial Wave Expansion
More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationSarabjeet Singh (1), Colin H. Hansen (2) and Carl Q. Howard (2)
Acoustics 8 Geeong, Victoria, Austraia 4 to 6 November 8 Acoustics and ustainabiity: How shoud acoustics adapt to meet future demands? A detaied tutoria for evauating in-duct net acoustic power transmission
More informationApplied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation
22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements
More informationLobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z
Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant
More informationVTU-NPTEL-NMEICT Project
MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid
More informationThe influence of temperature of photovoltaic modules on performance of solar power plant
IOSR Journa of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vo. 05, Issue 04 (Apri. 2015), V1 PP 09-15 www.iosrjen.org The infuence of temperature of photovotaic modues on performance
More informationLecture 8 February 18, 2010
Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some
More informationPhysics 506 Winter 2006 Homework Assignment #6 Solutions
Physics 506 Winter 006 Homework Assignment #6 Soutions Textbook probems: Ch. 10: 10., 10.3, 10.7, 10.10 10. Eectromagnetic radiation with eiptic poarization, described (in the notation of Section 7. by
More informationarxiv:nlin/ v2 [nlin.cd] 30 Jan 2006
expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische
More information12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes
Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.
More informationTHE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE
THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE KATIE L. MAY AND MELISSA A. MITCHELL Abstract. We show how to identify the minima path network connecting three fixed points on
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),
More informationMethods for Ordinary Differential Equations. Jacob White
Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa
More informationQuantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18
Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationCandidate Number. General Certificate of Education Advanced Level Examination January 2012
entre Number andidate Number Surname Other Names andidate Signature Genera ertificate of Education dvanced Leve Examination January 212 Physics PHY4/1 Unit 4 Fieds and Further Mechanics Section Tuesday
More informationAn approximate method for solving the inverse scattering problem with fixed-energy data
J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999
More informationTracking Control of Multiple Mobile Robots
Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart
More informationSmoothers for ecient multigrid methods in IGA
Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third
More informationshould the warm BPMs in LHC be coated with a 100 micron copper layer? (question by Gerhard Schneider)
shoud the warm BPMs in LHC be coated with a micron copper ayer? (question by Gerhard Schneider) 46 BPMs per beam (6 BPMSW, 8 BPMW, 4 BPMWA, 8 BPMWB) Average beta Injection Top Horizonta beta Vertica beta
More information18. Atmospheric scattering details
8. Atmospheric scattering detais See Chandrasekhar for copious detais and aso Goody & Yung Chapters 7 (Mie scattering) and 8. Legendre poynomias are often convenient in scattering probems to expand the
More informationUNIVERSITY OF CINCINNATI
UNIVERSITY OF CINCINNATI 5 August 006 Date: I,, hereby submit this work as part of the requirements for the degree of: in: Miind Shrikant Bapat Master of Science Department of Mechanica Engineering It
More informationA nodal collocation approximation for the multidimensional P L equations. 3D applications.
XXI Congreso de Ecuaciones Diferenciaes y Apicaciones XI Congreso de Matemática Apicada Ciudad Rea, 1-5 septiembre 9 (pp. 1 8) A noda coocation approximation for the mutidimensiona P L equations. 3D appications.
More informationBohr s atomic model. 1 Ze 2 = mv2. n 2 Z
Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationA Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)
A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,
More informationElements of Kinetic Theory
Eements of Kinetic Theory Statistica mechanics Genera description computation of macroscopic quantities Equiibrium: Detaied Baance/ Equipartition Fuctuations Diffusion Mean free path Brownian motion Diffusion
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More information1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2
MATEC Web of Conferences 159, 01005 (018) IJCAET & ISAMPE 017 https://doi.org/10.1051/matecconf/01815901005 Dynamic Response of RC Cantiever Beam by Equivaent Singe Degree of Freedom Method on Eastic Anaysis
More informationNumerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process
ISIJ Internationa Vo 54 (2014) No 2 pp 254 258 Numerica Simuation for Optimizing Temperature Gradients during Singe Crysta Casting Process Aeksandr Aeksandrovich INOZEMTSEV 1) Aeksandra Sergeevna DUBROVSKAYA
More informationD. Prémel, J.M. Decitre and G. Pichenot. CEA, LIST, F Gif-sur-Yvette, France
SIMULATION OF EDDY CURRENT INSPECTION INCLUDING MAGNETIC FIELD SENSOR SUCH AS A GIANT MAGNETO-RESISTANCE OVER PLANAR STRATIFIED MEDIA COMPONENTS WITH EMBEDDED FLAWS D. Préme, J.M. Decitre and G. Pichenot
More informationPost-buckling behaviour of a slender beam in a circular tube, under axial load
Computationa Metho and Experimenta Measurements XIII 547 Post-bucking behaviour of a sender beam in a circuar tube, under axia oad M. Gh. Munteanu & A. Barraco Transivania University of Brasov, Romania
More informationResearch on liquid sloshing performance in vane type tank under microgravity
IOP Conference Series: Materias Science and Engineering PAPER OPEN ACCESS Research on iquid soshing performance in vane type tan under microgravity Reated content - Numerica simuation of fuid fow in the
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More informationCoupling of LWR and phase transition models at boundary
Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of
More informationTIME DEPENDENT TEMPERATURE DISTRIBUTION MODEL IN LAYERED HUMAN DERMAL PART
VOL. 8, No. II, DECEMBER, 0, pp 66-76 TIME DEPENDENT TEMPERATURE DISTRIBUTION MODEL IN LAYERED HUMAN DERMAL PART Saraswati Acharya*, D. B. Gurung, V. P. Saxena Department of Natura Sciences (Mathematics),
More informationDIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM
DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing
More informationTHE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationThe Nottingham eprints service makes this work by researchers of the University of Nottingham available open access under the following conditions.
Uku, H. Arda and Bogaert, Ignace and Coos, Kristof and Andriui, Francesco P. and Bagci, Hakan (2017) of the time-domain MFIE at ow frequencies. IEEE Antennas and Wireess Propagation Letters, 16. pp. 1565-1568.
More information14-6 The Equation of Continuity
14-6 The Equation of Continuity 14-6 The Equation of Continuity Motion of rea fuids is compicated and poory understood (e.g., turbuence) We discuss motion of an idea fuid 1. Steady fow: Laminar fow, the
More informationPhysics 505 Fall Homework Assignment #4 Solutions
Physics 505 Fa 2005 Homework Assignment #4 Soutions Textbook probems: Ch. 3: 3.4, 3.6, 3.9, 3.0 3.4 The surface of a hoow conducting sphere of inner radius a is divided into an even number of equa segments
More informationarxiv: v1 [physics.flu-dyn] 2 Nov 2007
A theoretica anaysis of the resoution due to diffusion and size-dispersion of partices in deterministic atera dispacement devices arxiv:7.347v [physics.fu-dyn] 2 Nov 27 Martin Heer and Henrik Bruus MIC
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationHYDROGEN ATOM SELECTION RULES TRANSITION RATES
DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS
More informationUltrasonic Measurements of Kinematic Viscosity for Analize of Engine Oil Parameters
th European Conference on Non-Destructive Testing (ECNDT 04), October 6-0, 04, Prague, Czech Repubic More Info at Open Access Database www.ndt.net/?id=6344 Utrasonic Measurements of Kinematic Viscosity
More informationTwo Kinds of Parabolic Equation algorithms in the Computational Electromagnetics
Avaiabe onine at www.sciencedirect.com Procedia Engineering 9 (0) 45 49 0 Internationa Workshop on Information and Eectronics Engineering (IWIEE) Two Kinds of Paraboic Equation agorithms in the Computationa
More informationAerodynamic Design of Turbine Blades Using an Adjoint Equation Method
43rd AIAA Aerospace Sciences Meeting and Exhibit, January 10-13, 2005, Reno, Nevada Aerodynamic Design of Turbine ades Using an Adjoint Equation Method Hsiao-Yuan Wu, and Feng Liu Department of Mechanica
More informationNumerical simulation of javelin best throwing angle based on biomechanical model
ISSN : 0974-7435 Voume 8 Issue 8 Numerica simuation of javein best throwing ange based on biomechanica mode Xia Zeng*, Xiongwei Zuo Department of Physica Education, Changsha Medica University, Changsha
More informationO9e Fringes of Equal Thickness
Fakutät für Physik und Geowissenschaften Physikaisches Grundpraktikum O9e Fringes of Equa Thickness Tasks 1 Determine the radius of a convex ens y measuring Newton s rings using ight of a given waveength.
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationAND RAMANI DURAISWAMI
MULTIPLE SCATTERING FROM N SPHERES USING MULTIPOLE REEXPANSION NAIL A. GUMEROV AND RAMANI DURAISWAMI Abstract. A semi-anaytica technique for the soution of probems of wave scattering from mutipe spheres
More information2.1. Cantilever The Hooke's law
.1. Cantiever.1.1 The Hooke's aw The cantiever is the most common sensor of the force interaction in atomic force microscopy. The atomic force microscope acquires any information about a surface because
More informationSubstructuring Preconditioners for the Bidomain Extracellular Potential Problem
Substructuring Preconditioners for the Bidomain Extraceuar Potentia Probem Mico Pennacchio 1 and Vaeria Simoncini 2,1 1 IMATI - CNR, via Ferrata, 1, 27100 Pavia, Itay mico@imaticnrit 2 Dipartimento di
More informationThermophoretic interaction of heat releasing particles
JOURNAL OF APPLIED PHYSICS VOLUME 9, NUMBER 7 1 APRIL 200 Thermophoretic interaction of heat reeasing partices Yu Doinsky a) and T Eperin b) Department of Mechanica Engineering, The Pearstone Center for
More informationChapter 2 Theoretical Preliminaries of Acoustics
Chapter 2 Theoretica Preiminaries of Acoustics In this chapter, we review some of the fundamentas of acoustics and introduce the spherica harmonic expansion of a sound fied, which is the basis for the
More informationExperimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure
Experimenta Investigation and Numerica Anaysis of New Muti-Ribbed Sab Structure Jie TIAN Xi an University of Technoogy, China Wei HUANG Xi an University of Architecture & Technoogy, China Junong LU Xi
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationTorsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg
Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................
More informationInternational Journal of Mass Spectrometry
Internationa Journa of Mass Spectrometry 280 (2009) 179 183 Contents ists avaiabe at ScienceDirect Internationa Journa of Mass Spectrometry journa homepage: www.esevier.com/ocate/ijms Stark mixing by ion-rydberg
More informationTarget Location Estimation in Wireless Sensor Networks Using Binary Data
Target Location stimation in Wireess Sensor Networks Using Binary Data Ruixin Niu and Pramod K. Varshney Department of ectrica ngineering and Computer Science Link Ha Syracuse University Syracuse, NY 344
More informationSimple models for rope substructure mechanics: application to electro-mechanical lifts
Journa of Physics: Conference Series PAPER OPEN ACCESS Simpe modes for rope substructure mechanics: appication to eectro-mechanica ifts To cite this artice: I Herrera and S Kaczmarczyk 016 J. Phys.: Conf.
More informationChapter 11: Two-Phase Flow and Heat Transfer Forced Convective Boiling in Tubes
11.5 Forced Convective 11.5.1 Regimes in Horizonta and Vertica Tubes The typica sequence of fow regimes for upward fow forced convective boiing in a uniformy-heated vertica tube (q =const) is shown in
More information