FINITE ELEMENT MODELING OF TRANSIENT WAVE PHENOMENA AT
|
|
- Russell Hart
- 5 years ago
- Views:
Transcription
1 FINITE ELEMENT MODELING OF TRANSIENT WAVE PHENOMENA AT SOLIDIFLUID INTERFACES T. Xue, W. Lord, S. Udpa, L. Udpa and M. Mina Department of Electrical and Computer Engineering Iowa State University Ames, IA INTRODUCTION The solid/fluid interface appears in many ultrasonic measurement systems. Models for the system must take account of the interface. Analytical models for wave phenomena at the interface (especially curved interfaces) are either difficult or subject to severe approximation. The finite element method is ideal for this especially when the problem domain is bounded. A survey of this subject has been given by Kalinowski [1]. In this paper, an axisymmetric finite element model is developed for a solid medium and a fluid medium in contact. Displacement is used as the primary variable in the solid media and pressure in the fluid. The scalar pressure in the fluid medium makes the total degrees of freedom less than if displacement is used. The global mass matrix and stiffness matrix are rendered symmetric by introducing a potential variable for the fluid medium [2]. The final finite element equations are solved by the explicit integration approach. FINITE ELEMENT MODELING Solid Medium The governing equation for an isotropic medium with neglected body force and viscous damping is as follows Upon discretization of the problem domain to quadrilateral isoparametric elements, the following matrix equation is derived using either the weighted residual or energy functional approach [3-6] (2) where K is the global stiffness matrix, M is the global mass matrix, F incorporates the traction boundary condition, and U is the displacement vector. (1) Review o/progress in Quantitative Nondestructive Evaluation, Vol. 15 Edited by D.O, Thompson and D.E. Chimenti, Plenum Press, New York,
2 Fluid Medium The governing wave equation for pressure is as follows: n2 1.. v P=-P c 2 (3) where p=pressure and c=l wave velocity. Eq. (3) can be derived from Eq. (1) provided that A. =-11 (4) with a further assumption that U r = O. Then Uz stands for the pressure field. It can be shown that in the pressure analogy, the z component of surface traction becomes the product of the fluid bulk modulus and the normal derivative of the pressure. The finite element equations assume the following form M/p+KfP=f where f = Apc 2 ap / an. Note that from now on p stands for the fluid density. (5) SolidlFluid Coupling At the interface, the force applied to the solid by the fluid is - pa, while the force applied to the fluid by the solid is A(pC)2 Un' then the couple equation is as follows (6) Where A is an area matrix associated with the interfacing nodes. Let p = -q, we have t where G = -J fdt o In a more compact form, (8) Direct Time Integration Using a central difference approximation for the second derivatives and a backward difference approximation for the first derivatives, we have 1 q ( J...D-~M)Uq (J... D 1 M)Uq (9)!l.t2 MUt+flt - F, K + /).t!l.t2 t +!l.t /).t2 t-flt 300
3 O.3cm Front surface interface I --r z c=1500 mj 0=1000 k2/m V1=6300 mj V.=3100 mj O.7cm p=2700 kglm 1 Back surface 2cm Fig. 1. Axiymmetric geometry containing a solid/fluid interface. The mass matrix is diagonalized through a mass lumping technique, then Stability Conditions The explicit central difference scheme is conditionally stable. Then the integration time step is restricted for numerical stability by the following critical value Mer =2/wmax, r" "'max < - \l:g I 1/2 (11) where g is a geometric parameter. For a square element with sides h, the condition is which is valid for the lumped mass matrix. h M::;;- V; (12) RESULTS AND DISCUSSION The basic axisymmetric geometry shown in Fig. 1 is modeled using the finite element procedure presented above. The fluid and solid media are taken as water and aluminum respectively and their mechanical properties are indicated in the figure. A raised cosine function [3-6] with a center frequency of 5 MHz is used as the driving signal, while the source location and structure are varied for a few cases. In each case the source is also symmetric about the z-axis. As the first case, we assume a point source inside the solid medium 0.2 cm from the interface. The driving force of the source is along the axial direction. The wave profiles at different time instants are shown in Fig. 2. Note from now on that in each profile plot the solid side is represented by the z-displacement while the fluid side by the potential which is normalized such that the magnitude is on the same order as the z-displacement in the solid. Fig. 2(a) indicates that the reflection and refraction of the longitudinal wave front have occurred at the interface and the shear wave front from the source has reached the interface. In Fig. 2(b) the reflection and refraction of the shear wave as well as the reflection of the 301
4 (a) (b) (e) 1:=2 f..ls. (e) Fig. 2. Wove pmfil" (0' a point ",,,,,, imido tho 'olid medi"m.(a) '"I ~', (b) ~1.5 1";, and 302
5 (a) (b) /.Ls, and (c)!::4 /.Ls. F;g. 3. Wave profile, fo, a POint -Jil<" 'ou"", on the flu;d,"nace, (aj M. 75 ~', (b J /;3.25 (C) 303
6 (a) (b) Fig. 4. Wave profiles Eo< a 'PhericaJJy foe",,,,, Wave from tho num m'''lium to tho "'lid "'odium, 9_<6" (a) M"" (b) "'3.5 #', and (e) ""'",. 304 (c)
7 (a) (b) Fig. 5. Wave Profiles fo, a 'ph",;oa/ly foou,," wave from 'he fluid me<fiwn '0 the 'Olid medium, 9_>9" (a) "'3.5 ~', (b) "'4 ~', and (0) "'5 ~,. (c) 305
8 longitudinal wave from the back surface are shown. Each reflection associated with the solid medium is accompanied by a mode conversion. Fig. 2(c) records the various wave fronts at a later time. As a second case, we consider a point-like source, i.e., a small finite aperture, on the front surface of the fluid medium. The aperture diameter is 2a=0.01 cm, which corresponds to ka=1t/3. The wave profiles at three different time instants are shown in Fig. 3. It can been seen that after the longitudinal wave front reaches the interface, both reflection and refraction occur. Among the refracted waves, we can identify the longitudinal, shear, head and interface Rayleigh waves. A further point to note is that as the interface Rayleigh wave propagates, it radiates (or leaks) waves (at the Rayleigh angle) back to the fluid medium. By examining the time dependent signals at points along the interface, the interface Rayleigh velocity is estimated to be VR=2940 mis, which is slightly greater than the corresponding Rayleigh wave velocity on a traction free surface. As a third case, assume a spherically focused wave from an aperture with diameter of 1.5 cm and focal length of 2 cm. The wave profiles are shown in Fig. 4. From Fig. 4(a) we see the reflected longitudinal wave and the refracted shear wave. At that time instant a longitudinal wave has not been refracted into the solid medium. This is expected as the maximum angle of incidence is between the critical angles for longitudinal and shear waves (9max =20.56,91=13.77, 9s=28.94 ). An interface Rayleigh wave does not exist as the incident angle is always smaller than the critical Rayleigh angle which is estimated to be 9 R= As the incident wave proceeds as shown in Figs. 4(b) and (c), the incident angle becomes smaller and then a focused longitudinal wave is transmitted into the solid medium. It is defocused after passing the focal point. As a final case, consider a spherically focused wave similar to the previous wave except that the focal length is 1 cm. The corresponding wave profiles are shown in Fig. 5. An essential difference is that the maximum angle of incidence (9max=36.87 ) is greater than all the three critical angles. As the incident wave proceeds, the interface Rayleigh wave, shear wave, longitudinal wave and head wave gradually appear. The head wave on the fluid side radiated from the interface Rayleigh wave is clearly indicated especially when the interface wave passes through the axis of symmetry. CONCLUSIONS Transient wave propagation phenomena associated with solid/fluid interfaces have been presented using the finite element method. Various wave profiles are predicted based on different source conditions. In this study, the longitudinal wave velocity in the fluid medium is smaller than the Rayleigh wave velocity for the solid medium. The results indicate that an interface Rayleigh wave is generated when an incident wave from the fluid medium matches the critical Rayleigh angle, and a propagating interface Rayleigh wave radiates waves back into the fluid medium. Although only planar interfaces are considered here, the extension to arbitrarily curved interfaces is trivial for the finite element method with the expense of some extra computer resources. REFERENCES 1. A. 1. Kalinowski, in Shock and Vibration Computer Programs, Reviews and Summaries, Edited by W. Pikley and B. Pikley, (1975), pp G. C. Everstine, J. Sound and Vibration, 79, (1981), pp R. Ludwig. and W. Lord, IEEE Trans. UFFC, 35, (1988), pp R. Ludwig. and W. Lord, IEEE Trans. UFFC, 36, (1989), pp W. Lord, R. Ludwig and Z. You, J. NDE, 9, (1990), pp Z. You, M. Lust, R. Ludwig and W. Lord, IEEE Trans. UFFC 38, (1991), pp
ULTRASONIC INSPECT ABILITY MODELS FOR JET ENGINE FORGINGS
ULTRASONIC INSPECT ABILITY MODELS FOR JET ENGINE FORGINGS INTRODUCTION T. A. Gray Center for Nondestructive Evaluation Iowa State University Ames, IA 50011 Ultrasonic inspections of axially symmetric forgings,
More informationModeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics
Modeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics Anthony L. Bonomo *1 and Marcia J. Isakson 1 1 Applied
More informationValidation of aspects of BeamTool
Vol.19 No.05 (May 2014) - The e-journal of Nondestructive Testing - ISSN 1435-4934 www.ndt.net/?id=15673 Validation of aspects of BeamTool E. GINZEL 1, M. MATHESON 2, P. CYR 2, B. BROWN 2 1 Materials Research
More informationInspection of Spar-Core Bond in Helicopter Rotor Blades Using Finite Element Analysis
Inspection of Spar-Core Bond in Helicopter Rotor Blades Using Finite Element Analysis Sunil Kishore Chakrapani* a,b, Daniel J. Barnard a, and Vinay Dayal a,b a Center for NDE, Iowa State University, Ames,
More informationRA Y -MODELING FOR COMPUTER SIMULATION OF ULTRASONIC TESTING
RA Y -MODELING FOR COMPUTER SIMULATION OF ULTRASONIC TESTING INTRODUCTION T. Furukawa and K. Date Department of Materials Processing, Faculty of Engineering Tohoku University Sendai, Miyagi, 9'iYJ-77,
More informationRECONSTRUCTION OF A PISTON TRANSDUCER BEAM USING MULTI-GAUSSIAN
RECONSTRUCTION OF A PISTON TRANSDUCER BEAM USING MULTI-GAUSSIAN BEAMS (MGB) AND ITS APPLICATIONS INTRODUCTION A. Minachi*, F. 1. Margetan** and R. B. Thompson** * Southwest Research Institute San Antonio,
More informationULTRASONIC TRANSDUCER RADIATION THROUGH A CURVED FLUID-SOLID
ULTRASONIC TRANSDUCER RADIATION THROUGH A CURVED FLUID-SOLID INTERFACE INTRODUCTION Terence P. Lerch and Lester W. Schmerr Center for NDE and the Department of Aerospace Engineering and Engineering Mechanics
More informationWall thickness= Inlet: Prescribed mass flux. All lengths in meters kg/m, E Pa, 0.3,
Problem description Problem 30: Analysis of fluid-structure interaction within a pipe constriction It is desired to analyze the flow and structural response within the following pipe constriction: 1 1
More informationAPPLICATION OF ULTRASONIC BEAM MODELING TO PHASED ARRAY
APPLICATION OF ULTRASONIC BEAM MODELING TO PHASED ARRAY TESTING OF COMPLEX GEOMETRY COMPONENTS INTRODUCflON O. Roy, S. Mahaut, M. Serre, Commissariat 11 l'energie Atomique, CEAlCEREM, CE Saclay France
More informationVALIDITY OF THE GAUSS-HERMITE ßEAM MODEL IN AN ANISOTROPIC
VALIDITY OF THE GAUSS-HERMITE ßEAM MODEL IN AN ANISOTROPIC LAYERED MEDIUM: COMPARISON TO THE FINITE ELEMENT METHOD A. Minachi, Z. You, R.ß. Thompson, and W. Lord Center for Nondestructive Evaluation and
More informationULTRASONIC WAVE PROPAGATION THROUGH NOZZLES AND PIPES WITH
ULTRASONIC WAVE PROPAGATION THROUGH NOZZLES AND PIPES WITH CLADDINGS AROUND THEIR INNER WALLS INTRODUCTION A. Minachi and R. B. Thompson Center for NDE Iowa State University Ames, Iowa 5001 J The inner
More informationVALIDATION 3D RAY TRACING FOR UT EXAMINATION OF THE NOZZLE
VALIDATION 3D RAY TRACING FOR UT EXAMINATION OF THE NOZZLE INNER BLEND REGION INTRODUCTION Douglas MacDonald Greg Selby EPRI NDE Center Charlotte, NC 28221 Mathew Koshy Weidlinger Associates, Inc. Los
More informationS. J. Wormley, B. P. Newberry, M. S. Hughes D. K. Hsu, and D. O. Thompson Center for NDE Iowa State University Ames, Iowa 50011
APPLICATION OF GAUSS-HERMITE BEAM MODEL TO THE DESIGN OF ULTRASONIC PROBES S. J. Wormley, B. P. Newberry, M. S. Hughes D. K. Hsu, and D. O. Thompson Center for NDE Iowa State University Ames, Iowa 50011
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationOptimization to Reduce Automobile Cabin Noise
EngOpt 2008 - International Conference on Engineering Optimization Rio de Janeiro, Brazil, 01-05 June 2008. Optimization to Reduce Automobile Cabin Noise Harold Thomas, Dilip Mandal, and Narayanan Pagaldipti
More informationCOMPUTATIONALLY EFFICIENT RAY TRACING ALGORITHM FOR SIMULATION OF TRANSDUCER FIELDS IN ANISOTROPIC MATERIALS
Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation NDE 2011, December 8-10, 2011 COMPUTATIONALLY EFFICIENT RAY TRACING ALGORITHM FOR SIMULATION OF TRANSDUCER FIELDS IN ANISOTROPIC
More informationSHEAR WAVE WEDGE FOR LASER ULTRASONICS INTRODUCTION
SHEAR WAVE WEDGE FOR LASER ULTRASONICS INTRODUCTION R. Daniel Costleyl, Vimal Shah2, Krishnan Balasubramaniam2, and Jagdish Singhl 1 Diagnostic Instrumentation and Analysis Laboratories 2 Department of
More informationVOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL
VOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL The following is a description of the model setups and input/output parameters for benchmarking analytical volcanic deformation models
More informationCHAPTER 4 CFD AND FEA ANALYSIS OF DEEP DRAWING PROCESS
54 CHAPTER 4 CFD AND FEA ANALYSIS OF DEEP DRAWING PROCESS 4.1 INTRODUCTION In Fluid assisted deep drawing process the punch moves in the fluid chamber, the pressure is generated in the fluid. This fluid
More informationAcoustic Wave Crack Detection: A First Principles Approach
Acoustic Wave Crack Detection: A First Principles Approach Roger W. Pryor, Ph.D. Pryor Knowledge Systems, Inc. * 4918 Malibu Drive, Bloomfield Hills, MI, 48302, rwpryor@pksez1.com Abstract: Mechanical
More informationSIMULATION OF ULTRASONIC EXAMINATIONS USING SURFACE ACOUSTIC WAVES
SIMULATION OF ULTASONIC EXAMINATIONS USING SUFACE ACOUSTIC WAVES N. Leymarie 1, A. Lhémery 1, P. Calmon 1 S. Chatillon 1 and. Coulette 2 1 Commissariat à l Énergie Atomique, CEA-Saclay, Gif-sur-Yvette
More informationNUMERICAL SIMULATION OF SELF-FOCUSING OF AN ARRAY ON. Center for Quality Engineering and Failure Prevention Northwestern University Evanston, IL 60208
NUMERICAL SIMULATION OF SELF-FOCUSING OF AN ARRAY ON AN INTERIOR CRACK INTRODUCTION Ming Zhang and J. D. Achenbach Center for Quality Engineering and Failure Prevention Northwestern University Evanston,
More informationModeling Submerged Structures Loaded by Underwater Explosions with ABAQUS/Explicit
Modeling Submerged Structures Loaded by Underwater Explosions with ABAQUS/Explicit David B. Woyak ABAQUS Solutions Northeast, LLC Abstract: Finite element analysis can be used to predict the transient
More informationSaurabh GUPTA and Prabhu RAJAGOPAL *
8 th International Symposium on NDT in Aerospace, November 3-5, 2016 More info about this article: http://www.ndt.net/?id=20609 Interaction of Fundamental Symmetric Lamb Mode with Delaminations in Composite
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationFinite Element Modeling and Multiphysics Simulation of Air Coupled Ultrasonic with Time Domain Analysis
More Info at Open Access Database www.ndt.net/?id=15194 Finite Element Modeling and Multiphysics Simulation of Air Coupled Ultrasonic with Time Domain Analysis Bikash Ghose 1, a, Krishnan Balasubramaniam
More informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationAdvanced Image Reconstruction Methods for Photoacoustic Tomography
Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic
More informationMonochromatic Transfer Matrix method for acoustic field simulation thorough media boundaries
Available online at www.sciencedirect.com Physics Physics Procedia 3 (2010) 00 (2009) 883 890 000 000 www.elsevier.com/locate/procedia International Congress on Ultrasonics, Universidad de Santiago de
More informationAnalysis of Fluid-Structure Interaction Effects of Liquid-Filled Container under Drop Testing
Kasetsart J. (Nat. Sci.) 42 : 165-176 (2008) Analysis of Fluid-Structure Interaction Effects of Liquid-Filled Container under Drop Testing Chakrit Suvanjumrat*, Tumrong Puttapitukporn and Satjarthip Thusneyapan
More informationECNDT Poster 7
ECNDT 2006 - Poster 7 A Ray Technique to Calculate the Multiple Reflections and Leaky Wave Propagation from a Single or Multilayer Plate for the Detection of Critical Disbonds in Layered Media J. SADLER,
More informationGuidelines for proper use of Plate elements
Guidelines for proper use of Plate elements In structural analysis using finite element method, the analysis model is created by dividing the entire structure into finite elements. This procedure is known
More informationNUMERICAL MODELING OF DISPERSION PHENOMENA IN THICK PLATE
NUMERICAL MODELING OF DISPERSION PHENOMENA IN THICK PLATE Petr Hora, Jiří Michálek Abstract This paper report on a technique for the analysis of propagating multimode signals. The method involves a -D
More informationThe geometry of reflection and refraction Wave conversion and reflection coefficient
7.2.2 Reflection and Refraction The geometry of reflection and refraction Wave conversion and reflection coefficient The geometry of reflection and refraction A wave incident on a boundary separating two
More informationCALCULATING SNELL S LAW
CALCULATING SNELL S LAW What is Snell s Law? Snell s Law defines angular relationships of sound waves crossing an interface. The calculations for Snell s Law will come into play when we are conducting
More informationCHAPTER-10 DYNAMIC SIMULATION USING LS-DYNA
DYNAMIC SIMULATION USING LS-DYNA CHAPTER-10 10.1 Introduction In the past few decades, the Finite Element Method (FEM) has been developed into a key indispensable technology in the modeling and simulation
More informationTABLE OF CONTENTS SECTION 2 BACKGROUND AND LITERATURE REVIEW... 3 SECTION 3 WAVE REFLECTION AND TRANSMISSION IN RODS Introduction...
TABLE OF CONTENTS SECTION 1 INTRODUCTION... 1 1.1 Introduction... 1 1.2 Objectives... 1 1.3 Report organization... 2 SECTION 2 BACKGROUND AND LITERATURE REVIEW... 3 2.1 Introduction... 3 2.2 Wave propagation
More informationExample 24 Spring-back
Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static
More informationDYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION
DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION 7 DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION In this tutorial the influence of a vibrating source on its surrounding soil is studied.
More informationHomework Set 3 Due Thursday, 07/14
Homework Set 3 Due Thursday, 07/14 Problem 1 A room contains two parallel wall mirrors, on opposite walls 5 meters apart. The mirrors are 8 meters long. Suppose that one person stands in a doorway, in
More informationContinued Fraction Absorbing Boundary Conditions for Transient Elastic Wave Propagation Modeling
Continued Fraction Absorbing Boundary Conditions for Transient Elastic Wave Propagation Modeling Md Anwar Zahid and Murthy N. Guddati 1 Summary This paper presents a novel absorbing boundary condition
More informationModule 1: Introduction to Finite Element Analysis. Lecture 4: Steps in Finite Element Analysis
25 Module 1: Introduction to Finite Element Analysis Lecture 4: Steps in Finite Element Analysis 1.4.1 Loading Conditions There are multiple loading conditions which may be applied to a system. The load
More informationA Verification Study of ABAQUS AC3D8R Elements for Acoustic Wave Propagation
A Verification Study of ABAQUS AC3D8R Elements for Acoustic Wave Propagation by Michael Robert Hubenthal A Project Submitted to the Graduate Faculty of Rensselaer Polytechnic Institute in Partial Fulfillment
More informationCOMPUTER AIDED ENGINEERING. Part-1
COMPUTER AIDED ENGINEERING Course no. 7962 Finite Element Modelling and Simulation Finite Element Modelling and Simulation Part-1 Modeling & Simulation System A system exists and operates in time and space.
More informationTime-Domain Dynamic Analysis of Helical Gears with Reduced Housing Model
2013-01-1898 Published 05/13/2013 Copyright 2013 SAE International doi:10.4271/2013-01-1898 saeaero.saejournals.org Time-Domain Dynamic Analysis of Helical Gears with Reduced Housing Model Vijaya Kumar
More informationNumerical Calculations of Stability of Spherical Shells
Mechanics and Mechanical Engineering Vol. 14, No. 2 (2010) 325 337 c Technical University of Lodz Numerical Calculations of Stability of Spherical Shells Tadeusz Niezgodziński Department of Dynamics Technical
More informationMEEN 3360 Engineering Design and Simulation Spring 2016 Homework #1 This homework is due Tuesday, February 2, From your book and other
MEEN 3360 Engineering Design and Simulation Spring 2016 Homework #1 This homework is due Tuesday, February 2, 2016 1.0 From your book and other sources, write a paragraph explaining CFD and finite element
More information2: Static analysis of a plate
2: Static analysis of a plate Topics covered Project description Using SolidWorks Simulation interface Linear static analysis with solid elements Finding reaction forces Controlling discretization errors
More informationNOISE PROPAGATION FROM VIBRATING STRUCTURES
NOISE PROPAGATION FROM VIBRATING STRUCTURES Abstract R. Helfrich, M. Spriegel (INTES GmbH, Germany) Noise and noise exposure are becoming more important in product development due to environmental legislation.
More informationFederal Institute for Materials Research and Testing (BAM), Unter den Eichen 87, Berlin, Germany
Jannis Bulling 1, Jens Prager 1, Fabian Krome 1 1 Federal Institute for Materials Research and Testing (BAM), Unter den Eichen 87, 12205 Berlin, Germany Abstract: This paper addresses the computation of
More informationTutorial 2. Modeling Periodic Flow and Heat Transfer
Tutorial 2. Modeling Periodic Flow and Heat Transfer Introduction: Many industrial applications, such as steam generation in a boiler or air cooling in the coil of an air conditioner, can be modeled as
More informationInfluence of geometric imperfections on tapered roller bearings life and performance
Influence of geometric imperfections on tapered roller bearings life and performance Rodríguez R a, Calvo S a, Nadal I b and Santo Domingo S c a Computational Simulation Centre, Instituto Tecnológico de
More informationChapter 11 Application of High Speed Imaging in Particle Dynamics Study with Explosives
Chapter 11 Application of High Speed Imaging in Particle Dynamics Study with Explosives Elena Jacobs and Vilem Petr Abstract High-speed imaging can provide invaluable information during any explosive testing.
More informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationLECTURE 37: Ray model of light and Snell's law
Lectures Page 1 Select LEARNING OBJECTIVES: LECTURE 37: Ray model of light and Snell's law Understand when the ray model of light is applicable. Be able to apply Snell's Law of Refraction to any system.
More informationSimulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg
Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second
More informationLIGO Scissors Table Static Test and Analysis Results
LIGO-T980125-00-D HYTEC-TN-LIGO-31 LIGO Scissors Table Static Test and Analysis Results Eric Swensen and Franz Biehl August 30, 1998 Abstract Static structural tests were conducted on the LIGO scissors
More informationUsing MSC.Nastran for Explicit FEM Simulations
3. LS-DYNA Anwenderforum, Bamberg 2004 CAE / IT III Using MSC.Nastran for Explicit FEM Simulations Patrick Doelfs, Dr. Ingo Neubauer MSC.Software GmbH, D-81829 München, Patrick.Doelfs@mscsoftware.com Abstract:
More informationLASER GENERATION OF RAYLEIGH AND LAMB WAVES FOR ULTRASONIC
LASER GENERATION OF RAYLEIGH AND LAMB WAVES FOR ULTRASONIC NONDESTRUCTIVE TESTING R. Daniel Costley, Jr. and Yves H. Berthelot Woodruff School of Mechanical Engineering Georgia Institute of Technology
More information1. Particle Scattering. Cogito ergo sum, i.e. Je pense, donc je suis. - René Descartes
1. Particle Scattering Cogito ergo sum, i.e. Je pense, donc je suis. - René Descartes Generally gas and particles do not scatter isotropically. The phase function, scattering efficiency, and single scattering
More informationSimulation of Flow Development in a Pipe
Tutorial 4. Simulation of Flow Development in a Pipe Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D turbulent fluid flow in a pipe. The pipe networks are common
More informationReal-time full matrix capture with auto-focussing of known geometry through dual layered media
Real-time full matrix capture with auto-focussing of known geometry through dual layered media Mark Sutcliffe, Miles Weston, Ben Dutton and Ian Cooper TWI NDT Validation Centre Heol Cefn Gwrgan, Margam,
More informationThe Dynamic Response of an Euler-Bernoulli Beam on an Elastic Foundation by Finite Element Analysis using the Exact Stiffness Matrix
Journal of Physics: Conference Series The Dynamic Response of an Euler-Bernoulli Beam on an Elastic Foundation by Finite Element Analysis using the Exact Stiffness Matrix To cite this article: Jeong Soo
More informationSAMCEF for ROTORS. Chapter 3.2: Rotor modeling. This document is the property of SAMTECH S.A. MEF A, Page 1
SAMCEF for ROTORS Chapter 3.2: Rotor modeling This document is the property of SAMTECH S.A. MEF 101-03-2-A, Page 1 Table of contents Introduction Introduction 1D Model 2D Model 3D Model 1D Models: Beam-Spring-
More informationGENERAL AUTOMATED FLAW DETECTION SCHEME FOR NDE X-RAY IMAGES
GENERAL AUTOMATED FLAW DETECTION SCHEME FOR NDE X-RAY IMAGES Karl W. Ulmer and John P. Basart Center for Nondestructive Evaluation Department of Electrical and Computer Engineering Iowa State University
More informationAbout the Author. Acknowledgements
About the Author Dr. Paul Kurowski obtained his MSc and PhD in Applied Mechanics from Warsaw Technical University. He completed postdoctoral work at Kyoto University. Dr. Kurowski is an Assistant Professor
More informationNDT OF SPECIMEN OF COMPLEX GEOMETRY USING ULTRASONIC ADAPTIVE
NDT OF SPECIMEN OF COMPLEX GEOMETRY USING ULTRASONIC ADAPTIVE TECHNIQUES - THE F.A.U.S.T. SYSTEM INTRODUCTION O. Roy, S. Mahaut, M. Serre Commissariat a I'Energie Atomique CEAlCEREM, CE Saclay France Phased
More informationSimulating Underbelly Blast Events using Abaqus/Explicit - CEL
U.S. Army Research, Development and Engineering Command Simulating Underbelly Blast Events using Abaqus/Explicit - CEL J. Jablonski, P. Carlucci (U.S. Army ARDEC) R. Thyagarajan (U.S. Army TARDEC) B. Nandi,
More informationIntroduction to FEM Modeling
Total Analysis Solution for Multi-disciplinary Optimum Design Apoorv Sharma midas NFX CAE Consultant 1 1. Introduction 2. Element Types 3. Sample Exercise: 1D Modeling 4. Meshing Tools 5. Loads and Boundary
More informationMixed Mode Fracture of Through Cracks In Nuclear Reactor Steam Generator Helical Coil Tube
Journal of Materials Science & Surface Engineering Vol. 3 (4), 2015, pp 298-302 Contents lists available at http://www.jmsse.org/ Journal of Materials Science & Surface Engineering Mixed Mode Fracture
More informationTable of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Contents.
Table of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Preface page 2 Part I Fluid Dynamics and Waves 7 1 Elements of fluid dynamics 9 1.1
More informationSAFT-Reconstruction in ultrasonic immersion technique using phased array transducers
SAFT-Reconstruction in ultrasonic immersion technique using phased array transducers J. Kitze, J. Prager, R. Boehm, U. Völz, H.-J. Montag Federal Institute for Materials Research and Testing Berlin, Division
More informationSimulation of Overhead Crane Wire Ropes Utilizing LS-DYNA
Simulation of Overhead Crane Wire Ropes Utilizing LS-DYNA Andrew Smyth, P.E. LPI, Inc., New York, NY, USA Abstract Overhead crane wire ropes utilized within manufacturing plants are subject to extensive
More informationAuto Injector Syringe. A Fluent Dynamic Mesh 1DOF Tutorial
Auto Injector Syringe A Fluent Dynamic Mesh 1DOF Tutorial 1 2015 ANSYS, Inc. June 26, 2015 Prerequisites This tutorial is written with the assumption that You have attended the Introduction to ANSYS Fluent
More informationRay Optics I. Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex
Phys 531 Lecture 8 20 September 2005 Ray Optics I Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex Today shift gears, start applying
More informationUsing the Eulerian Multiphase Model for Granular Flow
Tutorial 21. Using the Eulerian Multiphase Model for Granular Flow Introduction Mixing tanks are used to maintain solid particles or droplets of heavy fluids in suspension. Mixing may be required to enhance
More informationDEVELOPMENT OF A NUMERICAL MODEL FOR SIMULATIONS OF SPLIT HOPKINSON PRESSURE BAR
DEVELOPMENT OF A NUMERICAL MODEL FOR SIMULATIONS OF SPLIT HOPKINSON PRESSURE BAR Afdhal 1, Annisa Jusuf 1, Muhammad Agus Kariem 2 and Leonardo Gunawan 1 1 Lightweight Structures Research Group, Faculty
More informationDevelopment of Compressible-Incompressible Link to Efficiently Model Bubble Dynamics near Floating Body Chao-Tsung Hsiao 1, Georges L.
Development of Compressible-Incompressible Link to Efficiently Model Bubble Dynamics near Floating Body Chao-Tsung Hsiao 1, Georges L. Chahine 1 DYNAFLOW, INC. 1061-J Iron Bridge Road, Jessup, MD 0794,
More informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
More informationModeling Unsteady Compressible Flow
Tutorial 4. Modeling Unsteady Compressible Flow Introduction In this tutorial, FLUENT s density-based implicit solver is used to predict the timedependent flow through a two-dimensional nozzle. As an initial
More informationTerminal Falling Velocity of a Sand Grain
Terminal Falling Velocity of a Sand Grain Introduction The first stop for polluted water entering a water work is normally a large tank, where large particles are left to settle. More generally, gravity
More informationCloth Simulation. COMP 768 Presentation Zhen Wei
Cloth Simulation COMP 768 Presentation Zhen Wei Outline Motivation and Application Cloth Simulation Methods Physically-based Cloth Simulation Overview Development References 2 Motivation Movies Games VR
More informationNumerical solution of Rayleigh-Lamb frequency equation for real, imaginary and complex wavenumbers
Numerical solution of Rayleigh-Lamb frequency equation for real, imaginary and complex wavenumbers Michal Šofer 1,*, Petr Ferfecki 1,2, Pavel Šofer 3 1 Department of Applied Mechanics, Faculty of Mechanical
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
More informationLS-DYNA s Linear Solver Development Phase 2: Linear Solution Sequence
LS-DYNA s Linear Solver Development Phase 2: Linear Solution Sequence Allen T. Li 1, Zhe Cui 2, Yun Huang 2 1 Ford Motor Company 2 Livermore Software Technology Corporation Abstract This paper continues
More informationChallenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure
Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure In the final year of his engineering degree course a student was introduced to finite element analysis and conducted an assessment
More informationAnalysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions
2010 NDIA GROUND VEHICLE SYSTEMS ENGINEERING AND TECHNOLOGY SYMPOSIUM MODELING & SIMULATION, TESTING AND VALIDATION (MSTV) MINI-SYMPOSIUM AUGUST 17-19 DEARBORN, MICHIGAN Analysis Comparison between CFD
More informationSimulation of NDT Inspection in 3D Elastic Waveguide Involving Arbitrary Defect
19 th World Conference on Non-Destructive Testing 2016 Simulation of NDT Inspection in 3D Elastic Waveguide Involving Arbitrary Defect Vahan BARONIAN 1, Karim JEZZINE 2 1 CEA LIST, Gif-sur-Yvette, France
More informationAcoustic Source Location in Vehicle Cabins and Free-field with Nearfield Acoustical Holography via Acoustic Arrays.
Acoustic Source Location in Vehicle Cabins and Free-field with Nearfield Acoustical Holography via Acoustic Arrays. D.L.Hallman S.M. Dumbacher B.W. Libbey J.S. Bolton D.L. Brown M.J. Lally Herrick Laboratories
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationSolution Recording and Playback: Vortex Shedding
STAR-CCM+ User Guide 6663 Solution Recording and Playback: Vortex Shedding This tutorial demonstrates how to use the solution recording and playback module for capturing the results of transient phenomena.
More informationCoupling of surface roughness to the performance of computer-generated holograms
Coupling of surface roughness to the performance of computer-generated holograms Ping Zhou* and Jim Burge College of Optical Sciences, University of Arizona, Tucson, Arizona 85721, USA *Corresponding author:
More informationULTRASONIC TESTING AND FLAW CHARACTERIZATION. Alex KARPELSON Kinectrics Inc., Toronto, Canada
ULTRASONIC TESTING AND FLAW CHARACTERIZATION Alex KARPELSON Kinectrics Inc., Toronto, Canada 1. Introduction Ultrasonic Testing (UT) is a commonly used inspection method. Various techniques are employed
More informationOptimised corrections for finite-difference modelling in two dimensions
Optimized corrections for 2D FD modelling Optimised corrections for finite-difference modelling in two dimensions Peter M. Manning and Gary F. Margrave ABSTRACT Finite-difference two-dimensional correction
More informationAcknowledgements. Prof. Dan Negrut Prof. Darryl Thelen Prof. Michael Zinn. SBEL Colleagues: Hammad Mazar, Toby Heyn, Manoj Kumar
Philipp Hahn Acknowledgements Prof. Dan Negrut Prof. Darryl Thelen Prof. Michael Zinn SBEL Colleagues: Hammad Mazar, Toby Heyn, Manoj Kumar 2 Outline Motivation Lumped Mass Model Model properties Simulation
More informationGUIDED WAVE ANIMATION USING SEMI-ANALYTICAL FINITE ELEMENT METHOD T. Hayashi Nagoya Institute of Technology, Nagoya, Japan
GUIDED WAVE ANIMATION USING SEMI-ANALYTICAL FINITE ELEMENT METHOD T. Hayashi Nagoya Institute of Technology, Nagoya, Japan Abstract: Many experimental and theoretical studies have shown the great potential
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More information2. Fundamental Equations of Elastodynamic Finite Integration Technique
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16455 Elastic Wave Modeling in Complex Geometries
More informationPhysics 11. Unit 8 Geometric Optics Part 1
Physics 11 Unit 8 Geometric Optics Part 1 1.Review of waves In the previous section, we have investigated the nature and behaviors of waves in general. We know that all waves possess the following characteristics:
More information