Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods
|
|
- Coral Hampton
- 5 years ago
- Views:
Transcription
1 TEAM 2007, Sept , 2007,Yokohama, Japan Hydoelastic Analysis of a 1900 TEU Containe Ship Using Finite Element and Bounday Element Methods Ahmet Egin 1)*, Levent Kaydıhan 2) and Bahadı Uğulu 3) 1) Pofesso, Faculty of Naval Achitectue and Ocean Engineeing, Istanbul Technical Univesity, Maslak, 34469, Istanbul, Tukey egina@itu.edu.t 2) Depatment Manage, Delta Maine, Okul Sk., Altunizade Sitesi, E/20 Altunizade, 34662, Üsküda, Istanbul, Tukey l.kaydihan@mesh-e.com 3) D., Faculty of Naval Achitectue and Ocean Engineeing, Istanbul Technical Univesity, Maslak, 34469, Istanbul, Tukey bugulu@itu.edu.t Abstact This pape pesents a hydoelastic analysis of a 1900 TEU containe ship using finite element and bounday elements. The method of analysis is sepaated into two pats. In the fist pat, the in vacuo dynamic popeties of the containe ship wee obtained by using a standad finite element method. In the second pat of the analysis, the ship stuctue was assumed vibating in its in vacuo modes when it is in contact with fluid, and the pessue distibution on the wetted suface was calculated sepaately fo each mode. The fluid-stuctue inteaction effects wee calculated in tems of the genealized added mass tems. The wetted suface of the containe ship was idealized by using appopiate bounday elements, efeed to as hydodynamic panels. A highe ode panel method (linea distibution) was adopted fo the calculations. In a futhe analysis, the wet calculations wee epeated using a finite element softwae (ABAQUS). The wet fequencies calculated fom the both analysis wee compaed, and a vey good compaison was obtained between the esults. Keywod: Hydoelasticity, bounday element method, finite element method, containe ship, vibation. 1. INTRODUCTION This pape pesents the dynamic chaacteistics (e.g., natual fequencies and mode shapes) of a 1900 TEU containe ship in full load and ballast conditions. The method pesented is based on a bounday integal equation method togethe with the method of images in ode to impose appopiate bounday condition on the fluid s fee suface. The method poposed in this study has aleady been successfully applied to stuctues patially filled with o patially submeged in quiescent and flowing fluids (see, fo instance, [1] [4]). In this investigation, it is assumed that the fluid is ideal, i.e., inviscid, incompessible and its motion is iotational. It is also assumed that the ship hull vibates in its in vacuo eigenmodes when it is in contact with fluid, and that each in vacuo mode gives ise to a coesponding suface pessue distibution on the wetted suface of the ship stuctue. The in vacuo dynamic analysis entails the vibation of the ship hull in the absence of any extenal foce and stuctual damping, and the coesponding dynamic chaacteistics of the containe ship ae obtained using a standad finite element softwae. At the fluid-stuctue inteface, continuity consideations equie that the nomal velocity of the fluid is equal to that of the stuctue. The nomal velocities on the wetted suface ae expessed in tems of the modal stuctual displacements. By using a bounday integal equation method the fluid pessue is eliminated fom the poblem, and that the fluid-stuctue inteaction foces ae calculated in tems of the genealized hydodynamic added mass coefficients. Duing this analysis, the wetted suface is idealized by using appopiate bounday elements, efeed to as hydodynamic panels. To assess the influence of the fluid on the dynamic esponse behavio of the containe ship, the wet fequencies and associated mode shapes ae calculated.
2 In a futhe analysis, a standad finite element softwae (ABAQUS [5]) was used fo modeling both the stuctual and fluid domains. The fluid domain suounding the ship hull was idealized by using acoustic fluid elements. The cago loads, ballast weight, etc. applied as inetia mass elements ove the cago aea, etc. The wet fequencies and mode shapes wee calculated by solving the fluid-stuctue inteaction poblem. The wet fequency values obtained fom the bounday element method (BEM) compae vey well with those calculated by ABAQUS. 2. MATHEMATICAL MODEL 2.1 Fluid stuctue inteaction poblem The fluid is assumed ideal, i.e., inviscid and incompessible, and its motion is iotational and thee exists a fluid velocity vecto, v, which can be defined as the gadient of the velocity potential function φ as v ( x, y, z, t) = φ ( x, y, z, t). (1) 2 whee φ satisfies the Laplace s equation, φ = 0, thoughout the fluid domain. Befoe descibing the esponses of the flexible stuctue, it is necessay to assign coodinates to the deflections at vaious degees of feedom. One paticula set of genealized coodinates having the advantage of being unambiguous and easily commended is the pincipal coodinates of the dy stuctue (see, fo example, Egin et al. [6]). Fo a stuctue vibating in an ideal fluid, with fequency ω, the pinciple coodinate, descibing the esponse of the stuctue in the th modal vibation, may be expessed by p ( t) = p e 0 iω t (2) The velocity potential function due to the distotion of the stuctue in the th in vacuo vibational mode may be witten as follows (,,, ) i (,, ) e i ω φ x y z t = ωφ x y z p t 0 = 1, 2,, M (3) whee M epesents the numbe of modes of inteest, and p 0 is an unknown amplitude fo the th pincipal coodinate. On the wetted suface of the vibating stuctue the nomal fluid velocity must equal to the nomal velocity on the stuctue and this condition fo the th modal vibation of the elastic stuctue submeged in fluid can be expessed as φ = u n, n whee n is the unit nomal vecto on the wetted suface and points into the egion of inteest. The vecto u denotes the displacement esponse of the stuctue in the th pincipal coodinate and it may be witten as (4) (,,, ) (,, ) e i ω u x y z t = u x y z p t 0 (5) whee u (x, y, z) is the th modal displacement vecto of the median suface of the stuctue, and it is obtained fom the in vacuo analysis. It is assumed that the elastic stuctue vibates at elatively high fequencies so that the effect of suface waves can be neglected. Theefoe, the fee suface condition (infinite fequency limit condition) (see, fo instance, [3]) fo the petubation potential can be appoximated by φ = 0, on the fee suface. (6) The method of images may be used to satisfy this bounday condition. By adding an imaginay bounday egion, the condition given by equation (6) at the hoizontal suface can be omitted; thus the poblem is educed to a classical Neumann case. The details of the method of images may be found in Egin and Temael [4].
3 2.2 Numeical evaluation of petubation potential φ The bounday value poblem fo the petubation potential,φ, may be expessed in the following fom: c s q s s q s ds * * ( ξ) φ ( ξ) = ( φ (, ξ) ( ) φ ( ) (, ξ)) S W (7) * whee ξ and s denote, espectively, the evaluation and field points on the wetted suface. φ is the fundamental solution and expessed as follows, φ ( s, ξ ) = 4π * 1 q = φ / n is the flux, and the distance between the evaluation and field points. The fee tem c (ξ ) is defined as * the faction of φ (ξ ) that lies inside the domain of inteest. Moeove, q ( s, ξ) can be witten as (8) q ( s, ξ ) = ( / n) / 4π * 3 (9) Fo the solution of equation (7) with bounday condition (4), the wetted suface can be idealized by using bounday elements, efeed to as hydodynamic panels, and the distibution of the potential function and its flux ove each hydodynamic panel may be descibed in tems of the shape functions and nodal values as ne ne φ = N φ q = N q e e e, e e e = 1 = 1 (10) Hee, n e epesents the numbe of nodal points assigned to each hydodynamic panel, and N e the shape function adopted fo the distibution of the potential function. e and indicates the numbes of the hydodynamic panels and nodal points, espectively. In the case of a linea distibution adopted in this study, the shape functions fo a quadilateal panel may be expessed as (see Wobel [7]) N e1 = ((1 ς )(1 η)) / 4 N e 2 = ((1 + ς )(1 η)) / 4 (11) N e 3 = ((1 + ς )(1 + η)) / 4 N e 4 = ((1 ς )(1 + η)) / 4 Afte substituting equations (10) and (11) into equation (7) and applying the bounday condition given in equation (4), the unknown potential function values can be detemined fom the following set of algebaic equations m nm m nm * * k φ k + ( φ i = ( i φ ), i= 1 = 1 S i= 1 = 1 i Si c N q ds u n N ds k =1, 2,, m (12) whee m denotes the numbe of nodal points used in the discetization of the stuctue. φ i and u i epesent, espectively, the potential and displacement vecto fo the th nodal point of the ith hydodynamic panel. 2.3 Calculation of genealized fluid-stuctue inteaction foces Using the Benoulli s equation and neglecting the second ode tems, the dynamic fluid pessue on the elastic stuctue due to the th modal vibation becomes
4 P ( x, y, z, t) = ρ φ, t (13) whee ρ is the fluid density. Substituting equation (3) into (13), the following expession fo the pessue is obtained, P x y z t = ρ ω φ (14) 2 iω t (,,, ) ( x, y, z) p0 e. The kth component of the genealized fluid-stuctue inteaction foce due to the th modal in-vacuo vibation of the elastic stuctue can be expessed in tems of the pessue acting on the wetted suface of the stuctue as Z = P ( x, y, z, t) u n ds k k SW = p e ρ ω φ u n ds iω t 2 0 k SW (15) The genealized added mass, A k tem can be defined as A = ρ φ u n ds, k k SW (16) Theefoe, the genealized fluid-stuctue inteaction foce component, Z k, can be witten as Z t = A ω p = A && p t 2 iωt k ( ) k 0 e k ( ) (17) 2.4 Calculation of wet fequencies and mode shapes It should be noted that, in the case when a body oscillates in o nea a fee suface, the hydodynamic coefficients exhibit fequency dependence in the low fequency egion, but show a tendency towads a constant value in the high fequency egion. In this study, it is assumed that the stuctue vibates in a elatively high fequency egion so that the genealized added mass values ae constants and evaluated by use of equation (16). Hence the genealized equation of motion fo the dynamic fluid stuctue inteaction system (see, e.g., Egin et al. [6]), assuming fee vibations with no stuctual damping, is ω 2 ( a + A ) + c p = 0, (18) whee a and c denote the genealized stuctual mass and stiffness matices, espectively. The matix A epesents the infinite fequency genealized added mass coefficients. Solving the eigenvalue poblem, expessed by (18), yields the wet fequencies and associated wet mode shapes of the stuctue in contact with fluid. To each wet fequency ω, thee is a coesponding wet eigenvecto p = p, p, K, p. The coesponding uncoupled wet mode shapes fo the stuctue patially and totally in { } m contact with fluid ae obtained as whee u ( x, y, z) = { u, v, w} = { } = u ( x, y, z) u, v, w u ( x, y, z) p, (19) = 1 M denote the in vacuo mode shapes of the elastic stuctue and M the numbe of mode shapes included in the analysis. It should be noted that the fluid stuctue inteaction foces associated with the inetial effect of the fluid do not have the same spatial distibution as those of the in vacuo modal foms. Consequently, this poduces hydodynamic coupling between the in vacuo modes of the stuctu. This coupling effect is intoduced into equation (18) though the genealized added mass matix A.
5 3. NUMERICAL RESULTS A containe ship is adopted fo the numeical calculations, and its main paticulas ae given in Table 1. Two diffeent loading conditions ae consideed: full load and ballast conditions. The cago loads ae distibuted as inetia mass elements ove the cago aea inne bottom plating. Ballast weight, heavy fuel oil and othe tank weights ae also applied as inetia mass elements. The light weight of the ship (LWT) is 9000 tons. The load goups ae given in Table 2 fo the full load and ballast conditions. The stuctual finite element model of the containe ship is shown in Figue 1. Table 1. The main paticulas of the containe ship Length (L bp ) 171 m Beadth (B) 28 m Depth (D) m Daught (T) m Speed (V) knot The in vacuo dynamic chaacteistics of the containe ship wee obtained using the ABAQUS finite element softwae [5]. This poduces infomation on natual fequencies and the nomal mode shapes of the dy containe ship in vacuum. In these calculations, the ship hull was discetized by using stuctual elements. The in vacuo dynamic chaacteistics of the ship stuctue ae scaled to a genealized stuctual mass of 1 ton m 2. Table 2. Load goups in full load and ballast conditions (tons). Full Load Condition Ballast Condition Cago loads Ballast weight Heavy Fuel Oil Maine Diesel Oil Fesh wate Othe weights Fig.1 The stuctual finite element model of the containe ship.
6 It should be noted that the finite element calculations wee caied out by the Delta Maine of Tukey. The dy natual fequencies of the containe ship ae given in Tables 3 and 4 fo the full load and ballast conditions. In the wet pat of the calculations, the fluid is assumed ideal, i.e., inviscid and incompessible, and its motion is iotational. The wetted suface of the containe ship was idealized by using bounday elements, efeed to as hydodynamic panels and 6739 hydodynamic panels wee distibuted ove the wetted suface of the ship hull, espectively, fo the full load and ballast conditions. A linea distibution of unknown potentials was assumed ove the each hydodynamic panel. These unknown potentials wee calculated using the kinematic bounday conditions imposed at each nodal point. The discitized wetted suface is shown in Figue 2 fo the full load condition. The kinematic bounday conditions wee obtained fom the in vacuo finite element calculations. By using Eq.(18), the eigenvalues and eigenvectos of the fluid-stuctue inteaction system wee calculated, and the wet fequencies ae pesented in Tables 3 and 4 fo the fist five global ship hull vibations togethe with the dy esults, espectively, fo the full load and ballast conditions. A numbe of 12 in vacuo modes was adopted fo the wet calculations. In a futhe analysis, a standad finite element pogam (ABAQUS) was adopted fo the wet calculations. The ship stuctue and fluid suounding the ship hull wee discetized by finite elements. Fo the fluid medium, acoustic finite elements wee employed fo the full load condition. Fig.2 The discetized wetted suface of the containe ship fo the full load condition. Table 3. Calculated fequencies fo the full load condition. In vacuo Wet Wet Diffeence Mode FEM FEM BEM % 1. Tosion Bending Bend.&Tos Tosion Bending Table 4. Calculated fequencies fo the ballast condition. In vacuo Wet Wet Diffeence Mode FEM FEM BEM % 1. Tosion Bending Bend.&Tos Tosion Bending
7 It can be seen fom Tables 3 and 4 that thee ae vey good ageement between the esults obtained fom the finite element analysis (ABAQUS) and bounday element analysis poposed in this study. The lagest diffeences obtained ae 2.5% and 6.2%, espectively, fo the full load and ballast conditions. On the othe hand, it can be obseved fom the tables that the fequencies decease with inceasing aea of contact with the fluid. The ship hull has the lagest aea of contact fo the full load condition, and theefoe, the lowest fequencies wee obtained fo this case. The mode shapes indicated in Tables 3 and 4 ae the global vibation modes of the ship hull, and the fist fou global mode shapes ae shown in Figue 3. (a) (b) (c) (d) Fig.3 The fist fou global modes of the containe ship: (a) 1 st tosion, (b) 1 st bending, tosion, (d) 2 nd bending. (c) hoizontal bending & CONCLUSIONS The dynamic chaacteistics (wet fequencies and associated mode shapes) of a containe ship fo two diffeent load conditions wee calculated by an appoach based on the bounday integal equation method, and the fequencies obtained wee compaed with the esults of the finite element calculations. A vey good compaison was obtained between the pedictions of the bounday element and finite element methods. It can be concluded fom the esults that the bounday integal equation method poposed in this study is suitable fo the dynamic analysis of ship hull stuctues and also fo othe elastic stuctues fully o patially in contact with fluid. ACKNOWLEDGEMENTS The authos acknowledge the suppot of a eseach gant fom the Scientific and Technological Reseach Council of Tukey (TUBITAK); Poect No: 105M041.
8 REFERENCES [1] Uğulu B., Egin, A., A hydoelasticity method fo vibating stuctues containing and/o submeged in flowing fluid, Jounal of Sound and Vibation, 290, (2006), [2] Egin, A., Uğulu, B., Hydoelastic analysis of fluid stoage tanks by using a bounday integal equation method, Jounal of Sound and Vibation, 275, (2004), [3] Egin, A., Uğulu, B., Linea vibation analysis of cantileve plates patially submeged in fluid, Jounal of Fluids and Stuctues, 17, (2003), [4] Egin, A., Temael, P., Fee vibation of a patially liquid-filled and submeged, hoizontal cylindical shell, Jounal of Sound and Vibation, 254, (2002), [5] ABAQUS, Use s Manual (2004). [6] Egin, A., Pice, W.G., Randall, R., Temael, P., Dynamic chaacteistics of a submeged, flexible cylinde vibating in finite wate depths, Jounal of Ship Reseach, 36, (1992), [7] L.C. Wobel, The Bounday Element Method, Applications in Themo-Fluid and Acoustics, Vol.1, John Wiley and Sons, New Yok, ( 2002).
Identification of the degradation of railway ballast under a concrete sleeper
Identification of the degadation of ailway ballast unde a concete sleepe Qin Hu 1) and Heung Fai Lam ) 1), ) Depatment of Civil and Achitectual Engineeing, City Univesity of Hong Kong, Hong Kong SAR, China.
More informationAPPLICATION OF MAC IN THE FREQUENCY DOMAIN
PPLICION OF MC IN HE FREQUENCY DOMIN D. Fotsch and D. J. Ewins Dynamics Section, Mechanical Engineeing Depatment Impeial College of Science, echnology and Medicine London SW7 2B, United Kingdom BSRC he
More informationFREE TRANSVERSE VIBRATIONS OF NON-UNIFORM BEAMS
Please cite this aticle as: Izabela Zamosa Fee tansvese vibations of non-unifom beams Scientific Reseach of the Institute of Mathematics and Compute Science Volume 9 Issue pages 3-9. The website: http://www.amcm.pcz.pl/
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Seies UG Examination 2015 16 FLUID DYNAMICS WITH ADVANCED TOPICS MTH-MD59 Time allowed: 3 Hous Attempt QUESTIONS 1 and 2, and THREE othe questions.
More informationis the instantaneous position vector of any grid point or fluid
Absolute inetial, elative inetial and non-inetial coodinates fo a moving but non-defoming contol volume Tao Xing, Pablo Caica, and Fed Sten bjective Deive and coelate the govening equations of motion in
More informationELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS
THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,
More informationA 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions
EN10: Continuum Mechanics Homewok 7: Fluid Mechanics Solutions School of Engineeing Bown Univesity 1. An ideal fluid with mass density ρ flows with velocity v 0 though a cylindical tube with cosssectional
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More informationA dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids
APCOM & ISCM 11-14 th Decembe, 013, Singapoe A dual-ecipocity bounday element method fo axisymmetic themoelastodynamic defomations in functionally gaded solids *W. T. Ang and B. I. Yun Division of Engineeing
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationLINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE
LINEAR AND NONLINEAR ANALYSES O A WIND-TUNNEL INTRODUCTION BALANCE R. Kakehabadi and R. D. Rhew NASA LaRC, Hampton, VA The NASA Langley Reseach Cente (LaRC) has been designing stain-gauge balances fo utilization
More informationLevitation force analysis of ring and disk shaped permanent magnet-high temperature superconductor
Inn Jounal of Pue & Applied Physics Vol. 55, Apil 017, pp. 61-68 Levitation foce analysis of ing and disk shaped pemanent magnet-high tempeatue supeconducto Sinan Basaan & Selim Sivioglu* Depatment of
More informationA NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM
Poceedings of the ASME 2010 Intenational Design Engineeing Technical Confeences & Computes and Infomation in Engineeing Confeence IDETC/CIE 2010 August 15-18, 2010, Monteal, Quebec, Canada DETC2010-28496
More informationChapter 3 Optical Systems with Annular Pupils
Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The
More informationTHE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN
THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN LIVIU NEAMŢ 1, ALINA NEAMŢ, MIRCEA HORGOŞ 1 Key wods: Magnetostatic shields, Magnetic non-lineaity, Finite element method.
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More informationLiquid gas interface under hydrostatic pressure
Advances in Fluid Mechanics IX 5 Liquid gas inteface unde hydostatic pessue A. Gajewski Bialystok Univesity of Technology, Faculty of Civil Engineeing and Envionmental Engineeing, Depatment of Heat Engineeing,
More informationChapter Introduction to Finite Element Methods
Chapte 1.4 Intoduction to Finite Element Methods Afte eading this chapte, you should e ale to: 1. Undestand the asics of finite element methods using a one-dimensional polem. In the last fifty yeas, the
More informationTitle. Author(s)Y. IMAI; T. TSUJII; S. MOROOKA; K. NOMURA. Issue Date Doc URL. Type. Note. File Information
Title CALCULATION FORULAS OF DESIGN BENDING OENTS ON TH APPLICATION OF THE SAFETY-ARGIN FRO RC STANDARD TO Autho(s)Y. IAI; T. TSUJII; S. OROOKA; K. NOURA Issue Date 013-09-1 Doc URL http://hdl.handle.net/115/538
More informationApplication of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates
Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad
More informationShape optimization of a flywheel
Stuct Multidisc Optim 19, 74 81 Spinge-Velag 000 Shape optimization of a flywheel G.R. Kess Abstact This pape deals with finding the best thickness distibution along the adius of a centally boed flywheel
More informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More information2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum
2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known
More informationDiffusion and Transport. 10. Friction and the Langevin Equation. Langevin Equation. f d. f ext. f () t f () t. Then Newton s second law is ma f f f t.
Diffusion and Tanspot 10. Fiction and the Langevin Equation Now let s elate the phenomena of ownian motion and diffusion to the concept of fiction, i.e., the esistance to movement that the paticle in the
More informationContact impedance of grounded and capacitive electrodes
Abstact Contact impedance of gounded and capacitive electodes Andeas Hödt Institut fü Geophysik und extateestische Physik, TU Baunschweig The contact impedance of electodes detemines how much cuent can
More informationd 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c
Chapte 6 Geneal Relativity 6.1 Towads the Einstein equations Thee ae seveal ways of motivating the Einstein equations. The most natual is pehaps though consideations involving the Equivalence Pinciple.
More informationNuclear size corrections to the energy levels of single-electron atoms
Nuclea size coections to the enegy levels of single-electon atoms Babak Nadii Nii a eseach Institute fo Astonomy and Astophysics of Maagha (IAAM IAN P. O. Box: 554-44. Abstact A study is made of nuclea
More informationIn the previous section we considered problems where the
5.4 Hydodynamically Fully Developed and Themally Developing Lamina Flow In the pevious section we consideed poblems whee the velocity and tempeatue pofile wee fully developed, so that the heat tansfe coefficient
More informationSupplementary material for the paper Platonic Scattering Cancellation for Bending Waves on a Thin Plate. Abstract
Supplementay mateial fo the pape Platonic Scatteing Cancellation fo Bending Waves on a Thin Plate M. Fahat, 1 P.-Y. Chen, 2 H. Bağcı, 1 S. Enoch, 3 S. Guenneau, 3 and A. Alù 2 1 Division of Compute, Electical,
More informationDynamics of Rotating Discs
Dynamics of Rotating Discs Mini Poject Repot Submitted by Subhajit Bhattachaya (0ME1041) Unde the guidance of Pof. Anivan Dasgupta Dept. of Mechanical Engineeing, IIT Khaagpu. Depatment of Mechanical Engineeing,
More informationAnalytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra
Poceedings of the 006 IASME/SEAS Int. Conf. on ate Resouces, Hydaulics & Hydology, Chalkida, Geece, May -3, 006 (pp7-) Analytical Solutions fo Confined Aquifes with non constant Pumping using Compute Algeba
More informationMU+CU+KU=F MU+CU+KU=0
MEEN 67 Handout # MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING ^ Symmetic he motion of a n-dof linea system is descibed by the set of 2 nd ode diffeential equations MU+CU+KU=F t whee U (t) and F
More informationRotor Blade Performance Analysis with Blade Element Momentum Theory
Available online at www.sciencediect.com ScienceDiect Enegy Pocedia 5 (7 ) 3 9 The 8 th Intenational Confeence on Applied Enegy ICAE6 Roto Blade Pefomance Analysis with Blade Element Momentum Theoy Faisal
More information2 Governing Equations
2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,
More informationInternational Journal of Scientific & Engineering Research, Volume 7, Issue 4, April ISSN
Intenational Jounal of Scientific & Engineeing Reseach, Volume 7, Issue 4, Apil-06 67 ISSN 9-558 Reduction at Payload Using Isolation System in a ypical Launch Vehicle Lekshmi Piya S,D. Sundaaajan, D.K
More informationAnalysis of high speed machining center spindle dynamic unit structure performance Yuan guowei
Intenational Confeence on Intelligent Systems Reseach and Mechatonics Engineeing (ISRME 0) Analysis of high speed machining cente spindle dynamic unit stuctue pefomance Yuan guowei Liaoning jidian polytechnic,dan
More informationINFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS. Abstract
INFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS Mohammad Mohammadi, National Cente fo Physical Acoustics, Univesity of Mississippi, MS Caig J. Hicey, National Cente fo Physical Acoustics,
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In Chaptes 2 and 4 we have studied kinematics, i.e., we descibed the motion of objects using paametes such as the position vecto, velocity, and acceleation without any insights
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In chaptes 2 and 4 we have studied kinematics i.e. descibed the motion of objects using paametes such as the position vecto, velocity and acceleation without any insights as to
More informationDymore User s Manual Two- and three dimensional dynamic inflow models
Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1
More informationA Three-Dimensional Magnetic Force Solution Between Axially-Polarized Permanent-Magnet Cylinders for Different Magnetic Arrangements
Poceedings of the 213 Intenational Confeence on echanics, Fluids, Heat, Elasticity Electomagnetic Fields A Thee-Dimensional agnetic Foce Solution Between Axially-Polaied Pemanent-agnet Cylindes fo Diffeent
More informationThe acoustic waves propagation in a cylindrical waveguide with the laminar flow
The acoustic waves popagation in a cylindical waveguide with the lamina flow Alexande Petov 1, and Valentina Rumyantseva, 1 MISIS National Univesity of Science and Technology Bauman Moscow State Technical
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationThermo-Mechanical Model for Wheel Rail Contact using Coupled. Point Contact Elements
IM214 28-3 th July, ambidge, England hemo-mechanical Model fo heel Rail ontact using oupled Point ontact Elements *J. Neuhaus¹ and. Sexto 1 1 hai of Mechatonics and Dynamics, Univesity of Padebon, Pohlweg
More informationANALYSIS OF PRESSURE VARIATION OF FLUID IN AN INFINITE ACTING RESERVOIR
Nigeian Jounal of Technology (NIJOTECH) Vol. 36, No. 1, Januay 2017, pp. 80 86 Copyight Faculty of Engineeing, Univesity of Nigeia, Nsukka, Pint ISSN: 0331-8443, Electonic ISSN: 2467-8821 www.nijotech.com
More informationDonnishJournals
DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationTransformation of the Navier-Stokes Equations in Curvilinear Coordinate Systems with Maple
Global Jounal of Pue and Applied Mathematics. ISSN 0973-1768 Volume 12, Numbe 4 2016, pp. 3315 3325 Reseach India Publications http://www.ipublication.com/gjpam.htm Tansfomation of the Navie-Stokes Equations
More information2.25 Advanced Fluid Mechanics
MIT Depatment of Mechanical Engineeing 2.25 Advanced Fluid Mechanics Poblem 4.27 This poblem is fom Advanced Fluid Mechanics Poblems by A.H. Shapio and A.A. Sonin u(,t) pg Gas Liquid, density Conside a
More informationA new class of exact solutions of the Navier Stokes equations for swirling flows in porous and rotating pipes
Advances in Fluid Mechanics VIII 67 A new class of exact solutions of the Navie Stokes equations fo swiling flows in poous and otating pipes A. Fatsis1, J. Stathaas2, A. Panoutsopoulou3 & N. Vlachakis1
More informationDuality between Statical and Kinematical Engineering Systems
Pape 00, Civil-Comp Ltd., Stiling, Scotland Poceedings of the Sixth Intenational Confeence on Computational Stuctues Technology, B.H.V. Topping and Z. Bittna (Editos), Civil-Comp Pess, Stiling, Scotland.
More informationTwo Dimensional Inertial Flow of a Viscous Fluid in a Corner
Applied Mathematical Sciences, Vol., 207, no. 9, 407-424 HIKARI Ltd, www.m-hikai.com https://doi.og/0.2988/ams.207.62282 Two Dimensional Inetial Flow of a Viscous Fluid in a Cone A. Mahmood and A.M. Siddiqui
More informationLecture 10. Vertical coordinates General vertical coordinate
Lectue 10 Vetical coodinates We have exclusively used height as the vetical coodinate but thee ae altenative vetical coodinates in use in ocean models, most notably the teainfollowing coodinate models
More informationMulti-Objective Optimization Algorithms for Finite Element Model Updating
Multi-Objective Optimization Algoithms fo Finite Element Model Updating E. Ntotsios, C. Papadimitiou Univesity of Thessaly Geece Outline STRUCTURAL IDENTIFICATION USING MEASURED MODAL DATA Weighted Modal
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationComputational Methods of Solid Mechanics. Project report
Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate
More informationSolution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere
Applied Physics eseach; Vol. 9, No. 6; 7 ISSN 96-969 E-ISSN 96-9647 Published by Canadian Cente of Science and Education Solution of a Spheically Symmetic Static Poblem of Geneal elativity fo an Elastic
More informationReview: Electrostatics and Magnetostatics
Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion
More informationSchool of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007
School of Electical and Compute Engineeing, Conell Univesity ECE 303: Electomagnetic Fields and Waves Fall 007 Homewok 8 Due on Oct. 19, 007 by 5:00 PM Reading Assignments: i) Review the lectue notes.
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationCOMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS
Pogess In Electomagnetics Reseach, PIER 73, 93 105, 2007 COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS T.-X. Song, Y.-H. Liu, and J.-M. Xiong School of Mechanical Engineeing
More informationElectromagnetism Physics 15b
lectomagnetism Physics 15b Lectue #20 Dielectics lectic Dipoles Pucell 10.1 10.6 What We Did Last Time Plane wave solutions of Maxwell s equations = 0 sin(k ωt) B = B 0 sin(k ωt) ω = kc, 0 = B, 0 ˆk =
More informationNumerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More informationMathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden
Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty
More informationEstimation of the Correlation Coefficient for a Bivariate Normal Distribution with Missing Data
Kasetsat J. (Nat. Sci. 45 : 736-74 ( Estimation of the Coelation Coefficient fo a Bivaiate Nomal Distibution with Missing Data Juthaphon Sinsomboonthong* ABSTRACT This study poposes an estimato of the
More informationApplied Aerodynamics
Applied Aeodynamics Def: Mach Numbe (M), M a atio of flow velocity to the speed of sound Compessibility Effects Def: eynolds Numbe (e), e ρ c µ atio of inetial foces to viscous foces iscous Effects If
More informationInstitute of Fundamental Sciences-Mathematics,Massey University, Palmerston North,NZ
Poceedings20th NZ Geothemal Wokshop 1998 FLUID FLOW A FLASHING CYCLONE SEPARATOR Institute of Fundamental SciencesMathematics,Massey Univesity, Palmeston Noth,NZ SUMMARY In geothemal steamgatheingsystems
More informationDESIGN OF INTERMEDIATE RING STIFFENERS FOR COLUMN-SUPPORTED CYLINDRICAL STEEL SHELLS
DESIGN OF INTERMEDIATE RING STIFFENERS FOR COLUMN-SUPPORTED CYLINDRICAL STEEL SHELLS Öze Zeybek Reseach Assistant Depatment of Civil Engineeing, Middle East Technical Univesity, Ankaa, Tukey E-mail: ozeybek@metu.edu.t
More informationPhysics 221 Lecture 41 Nonlinear Absorption and Refraction
Physics 221 Lectue 41 Nonlinea Absoption and Refaction Refeences Meye-Aendt, pp. 97-98. Boyd, Nonlinea Optics, 1.4 Yaiv, Optical Waves in Cystals, p. 22 (Table of cystal symmeties) 1. Intoductoy Remaks.
More informationEVALUATION OF S.I.F FOR CRACK EMANATING AT 45 0 ORIENTATION FROM A HOLE IN PRESSURISED CYLINDER USING FEA
EVALUATION OF S.I.F FOR CRAC EMANATING AT 45 0 ORIENTATION FROM A HOLE IN PRESSURISED CYLINDER USING FEA AASH.D.A 1, ANAND.A 2, G.V.GNANENDRA REDDY 3, SUDEV.L.J 4 1,3 Depatment of Mechanical Engineeing,
More informationPearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms
Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two
More informationAalborg Universitet. Load Estimation from Natural input Modal Analysis Aenlle, Manuel López; Brincker, Rune; Canteli, Alfonso Fernández
Aalbog Univesitet Load Estimation fom atual input Modal Analysis Aenlle, Manuel López; Bincke, Rune; Canteli, Alfonso Fenández Published in: Confeence Poceedings Publication date: 005 Document Vesion Publishe's
More information1 Similarity Analysis
ME43A/538A/538B Axisymmetic Tubulent Jet 9 Novembe 28 Similaity Analysis. Intoduction Conside the sketch of an axisymmetic, tubulent jet in Figue. Assume that measuements of the downsteam aveage axial
More informationOn the force and torque on systems of rigid bodies: A remark on an integral formula due to Howe
PHYSICS OF FLUIDS 19, 057108 007 On the foce and toque on systems of igid bodies: A emak on an integal fomula due to Howe Clodoaldo Gotta Ragazzo a Instituto de Matemática e Estatística, Univesidade de
More informationPerturbation to Symmetries and Adiabatic Invariants of Nonholonomic Dynamical System of Relative Motion
Commun. Theo. Phys. Beijing, China) 43 25) pp. 577 581 c Intenational Academic Publishes Vol. 43, No. 4, Apil 15, 25 Petubation to Symmeties and Adiabatic Invaiants of Nonholonomic Dynamical System of
More informationAnalysis and Optimization of a Special Type of Dielectric Loaded Resonant Cavity for Mobile Communication Filters
328 Analysis and Optimization of a Special Type of Dielectic Loaded Resonant Cavity fo Mobile Communication Filtes Haold S. Showes, Banmali S. Rawat *, Syam S. Challa Depatment of Electical and Biomedical
More informationChapter 9 Dynamic stability analysis III Lateral motion (Lectures 33 and 34)
Pof. E.G. Tulapukaa Stability and contol Chapte 9 Dynamic stability analysis Lateal motion (Lectues 33 and 34) Keywods : Lateal dynamic stability - state vaiable fom of equations, chaacteistic equation
More informationFlux. Area Vector. Flux of Electric Field. Gauss s Law
Gauss s Law Flux Flux in Physics is used to two distinct ways. The fist meaning is the ate of flow, such as the amount of wate flowing in a ive, i.e. volume pe unit aea pe unit time. O, fo light, it is
More informationTHE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2
THE LAPLACE EQUATION The Laplace (o potential) equation is the equation whee is the Laplace opeato = 2 x 2 u = 0. in R = 2 x 2 + 2 y 2 in R 2 = 2 x 2 + 2 y 2 + 2 z 2 in R 3 The solutions u of the Laplace
More informationI. CONSTRUCTION OF THE GREEN S FUNCTION
I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,
More informationModeling Fermi Level Effects in Atomistic Simulations
Mat. Res. Soc. Symp. Poc. Vol. 717 Mateials Reseach Society Modeling Femi Level Effects in Atomistic Simulations Zudian Qin and Scott T. Dunham Depatment of Electical Engineeing, Univesity of Washington,
More informationExercise sheet 8 (Modeling large- and small-scale flows) 8.1 Volcanic ash from the Eyjafjallajökull
Execise sheet 8 (Modeling lage- and small-scale flows) last edited June 18, 2018 These lectue notes ae based on textbooks by White [13], Çengel & al.[16], and Munson & al.[18]. Except othewise indicated,
More information11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.
Fall 2007 Qualifie Pat II 12 minute questions 11) A thin, unifom od of mass M is suppoted by two vetical stings, as shown below. Find the tension in the emaining sting immediately afte one of the stings
More informationBIFURCATION ANALYSIS FOR A QUASI-LINEAR CENTRIFUGAL PENDULUM ABSORBER
11 th Intenational Confeence on Vibation Poblems Z. Dimitovová et al. (eds.) Lisbon, Potugal, 9-12 Septembe 213 BIFURCATION ANALYSIS FOR A QUASI-LINEAR CENTRIFUGAL PENDULUM ABSORBER Eugen B. Keme* 1, Mikhail
More informationUncertainty in Operational Modal Analysis of Hydraulic Turbine Components
Intenational Jounal of Fluid Machiney and Systems Vol. 2, No. 4, Octobe-Decembe 2009 Oiginal Pape (Invited) Uncetainty in Opeational Modal Analysis of Hydaulic Tubine Components Matin Gagnon 1, S.-Antoine
More informationA Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation
Mathematical Modelling and Analysis Publishe: Taylo&Fancis and VGTU Volume 22 Numbe 3, May 27, 3 32 http://www.tandfonline.com/tmma https://doi.og/.3846/3926292.27.39329 ISSN: 392-6292 c Vilnius Gediminas
More informationME 425: Aerodynamics
ME 5: Aeodynamics D ABM Toufique Hasan Pofesso Depatment of Mechanical Engineeing, BUET Lectue- 8 Apil 7 teachebuetacbd/toufiquehasan/ toufiquehasan@mebuetacbd ME5: Aeodynamics (Jan 7) Flow ove a stationay
More informationField emission of Electrons from Negatively Charged Cylindrical Particles with Nonlinear Screening in a Dusty Plasma
Reseach & Reviews: Jounal of Pue and Applied Physics Field emission of Electons fom Negatively Chaged Cylindical Paticles with Nonlinea Sceening in a Dusty Plasma Gyan Pakash* Amity School of Engineeing
More informationNew procedure for dynamic structural reanalysis
New pocedue fo dynamic stuctual eanalysis NATASA TRISOVIC, TATJANA LAZOVIC, CASLAV MITROVIC, ALEKSANDAR MARINKOVIC, MIHAILO LAZAREVIC, DRAGOSLAV SUMARAC, ZORAN GOLUBOVIC Faculty of Mechanical Engineeing,
More informationEELE 3331 Electromagnetic I Chapter 4. Electrostatic fields. Islamic University of Gaza Electrical Engineering Department Dr.
EELE 3331 Electomagnetic I Chapte 4 Electostatic fields Islamic Univesity of Gaza Electical Engineeing Depatment D. Talal Skaik 212 1 Electic Potential The Gavitational Analogy Moving an object upwad against
More informationF-IF Logistic Growth Model, Abstract Version
F-IF Logistic Gowth Model, Abstact Vesion Alignments to Content Standads: F-IFB4 Task An impotant example of a model often used in biology o ecology to model population gowth is called the logistic gowth
More informationEM-2. 1 Coulomb s law, electric field, potential field, superposition q. Electric field of a point charge (1)
EM- Coulomb s law, electic field, potential field, supeposition q ' Electic field of a point chage ( ') E( ) kq, whee k / 4 () ' Foce of q on a test chage e at position is ee( ) Electic potential O kq
More informationChapter 13 Gravitation
Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects
More informationCHAPTER IV RADIATION BY SIMPLE ACOUSTIC SOURCE. or by vibratory forces acting directly on the fluid, or by the violent motion of the fluid itself.
CHAPTER IV RADIATION BY SIMPLE ACOUSTIC SOURCE 4.1 POINT SOURCE Sound waves ae geneated by the vibation of any solid body in contact with the fluid medium o by vibatoy foces acting diectly on the fluid,
More informationA Comparison and Contrast of Some Methods for Sample Quartiles
A Compaison and Contast of Some Methods fo Sample Quatiles Anwa H. Joade and aja M. Latif King Fahd Univesity of Petoleum & Mineals ABSTACT A emainde epesentation of the sample size n = 4m ( =, 1, 2, 3)
More informationA Multivariate Normal Law for Turing s Formulae
A Multivaiate Nomal Law fo Tuing s Fomulae Zhiyi Zhang Depatment of Mathematics and Statistics Univesity of Noth Caolina at Chalotte Chalotte, NC 28223 Abstact This pape establishes a sufficient condition
More informationPhysics 2A Chapter 10 - Moment of Inertia Fall 2018
Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.
More informationOn a quantity that is analogous to potential and a theorem that relates to it
Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich
More information