This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
|
|
- Merilyn Joseph
- 6 years ago
- Views:
Transcription
1 This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution and sharing with colleagues. Other uses, including reproduction and distribution, or selling or licensing copies, or posting to personal, institutional or third party websites are prohibited. In most cases authors are permitted to post their version of the article (e.g. in Word or Tex form) to their personal website or institutional repository. Authors requiring further information regarding Elsevier s archiving and manuscript policies are encouraged to visit:
2 Journal of Non-Crystalline Solids 355 (009) Contents lists available at ScienceDirect Journal of Non-Crystalline Solids journal homepage: Finite element analysis of auxetic obstacle deformation and fluid flow in a channel T. Strek a, P. Kedziora b, B. Maruszewski a, *, A.A. Pozniak b, K.V. Tretiakov b, K.W. Wojciechowski b, * a Poznan University of Technology, Institute of Applied Mechanics, ul. Piotrowo 3, Poznan, Poland b Institute of Molecular Physics, Polish Academy of Sciences, ul. M. Smoluchowskiego 1, Poznan, Poland article info abstract Article history: Available online 1 June 009 PACS: 6.0.-x 6.0.dj 4.10.ad 4.11.Fg The influence of the Poisson s ratio of a viscoelastic obstacle on the flow in a two-dimensional channel was studied by a finite element computer simulation. The arbitrary Lagrangian Eulerian (ALE) formulation for deformed mesh was applied to handle the moving boundaries of the computational domain. An anomalous deformation of the obstacle was observed when the Poisson s ratio was negative and close to its minimum value. It was also found that the maximum velocity of the fluid in the system decreased with the decreasing Poisson s ratio. Ó 009 Elsevier B.V. All rights reserved. Keywords: Auxetics Fluid-solid interaction Mechanical properties 1. Introduction Materials for which the Poisson s ratio [8] is negative were first manufactured by Lakes [] and later named auxetics by Evans [6]. Typically, the Poisson s ratio of isotropic materials is positive [8]. Although a negative Poisson s ratio has not been observed in most everyday materials, the classical theory of elasticity states that it is possible for isotropic three-dimensional materials to exhibit the Poisson s ratio within the range of 1 m þ1= [8] whereas two-dimensional isotropic systems can exhibit Poisson s ratios within the range of 1 m þ1 [15]. In anisotropic materials Poisson s ratios can have any positive or negative value [13]. Recently, increasing interest has been shown in the field of materials with a negative Poisson s ratio [11] which are interesting from both an intellectual point of view as well as in terms of potential applications. Studies of simple model systems exhibiting a negative Poisson s ratio or including auxetic parts should help recognize possible unconventional properties of such systems that may result in finding new potential applications of auxetics. The problem studied in this paper concerns a horizontal flow of an incompressible viscous fluid along a channel with an auxetic obstacle formed by a thin plate normal to its walls (Fig. 1). The main aim of our considerations is to find how the obstacle behaves under the flow. For the sake of simplicity, a version of the problem is considered within a plane stress state [3,4,9]. * Corresponding authors. addresses: tomasz.strek@put.poznan.pl (T. Strek), bogdan.maruszewski@put.poznan.pl (B. Maruszewski), kww@ifmpan.poznan.pl (K.W. Wojciechowski). Exact analytical treatment of the problem s mechanical model, even in its simplest version considered in this paper, is very difficult. Fortunately, solid fluid interactions can be modeled by approximate finite element method (FEM) techniques [16]. In particular, the arbitrary Lagrangian Eulerian (ALE) method can be used to follow the FEM mesh geometry deformation resulting from the motion of the computational domain boundaries [4,5]. The technique is very useful for a simulation of solid fluid mechanics. We applied this approach to improve the accuracy of the mathematical description of our problem. This requires a few words of comment. In the special case of the Lagrangian method, the mesh motion follows the motion of the solid material. However, when the material motion is more complicated, as in a fluid flow case, that method is not appropriate. In such a situation the Eulerian method, where the mesh is fixed, is often used even though this method cannot account for the moving boundaries. Therefore, we used the combined ALE technique to describe the fluid flow with detailed tracking of the obstacle boundaries together with the obstacle deformation.. Mathematical model Since the fluid flow along a channel with an obstacle concerns two interacting media, a mathematical description of the mechanics of the entire structure has to consist of two parts. The first part concerning the flow domain is described by an initial-boundary value problem based on the Navier Stokes equation. The second part concerns the solid obstacle which is governed by the viscoelastic equation of motion with proper boundary and initial conditions /$ - see front matter Ó 009 Elsevier B.V. All rights reserved. doi: /j.jnoncrysol
3 1388 T. Strek et al. / Journal of Non-Crystalline Solids 355 (009) The initial condition for (1) is based on the assumption that the fluid at t ¼ 0 is motionless (v ¼ 0) and p ¼ 0. The equation of motion for the solid obstacle q u þ r r ¼ G; where q s is the obstacle density, u is the displacement Þ is the damping term, r is the stress tensor and G is the body force. Fig. 1. The geometry of the studied system. As has been mentioned above, the fluid flow in the channel is described by the Navier Stokes ¼ rpþrsþf; ð1þ q þ v rv r v ¼ 0; where F is the volume force affecting the fluid and v is the fluid velocity, q 0 is the fluid density, p is the pressure, S ¼ g rv þðrvþ T is the extra stress tensor coming from fluid viscosity and g is the dynamic viscosity. In the studied case the force affecting the fluid has been omitted, i.e. F ¼ 0. The boundary conditions resulting from the stated problem are the following: v ¼ 0 at the lower and upper channel surface because there is no slip there, y m in ¼ 4m 0 h ðy 1Þ for y h0; hi, where h is the obstacle height h (inlet parabolic profile condition), ð pi þ SÞn ¼ p 0 n, where p 0 ¼ 0 is the outlet pressure, v ¼ 0 at the clamped obstacle surfaces (no slip condition), v at the deformed obstacle surfaces (u is the vector of the obstacle point). Table 1 The parameters used in simulations. Quantity Symbol Unit Value q 0 kg/m g Pa s 10 3 q s kg/m E kpa 0 m 0.99, 0., 0.3, +0.3 a dm s 1 1 b dk S 10 3 m 0 m/s 1.0 Table Mesh data. Quantity Mesh A Mesh B Mesh C Global number of mesh points Global number of mesh elements Fig.. Total displacement in a deformed auxetic obstacle for Poisson s ratio m ¼ 0:99 with meshes: (A) A, (B) B, and (C) C from Table.
4 T. Strek et al. / Journal of Non-Crystalline Solids 355 (009) Since the material of the obstacle has been assumed to be isotropic, the following constitutive equation writes: r ¼ De ¼ kðr uþi þ le; where I is the unit matrix, e ¼ 1 ðru þðruþt Þ is the strain tensor and k and l denote the Lamé constants as follows: Em k ¼ ð1 þ mþð1 mþ ; l ¼ E ð1 þ mþ ; ð4þ where E is Young s modulus, m is the Poisson s ratio and D is the elasticity matrix, ð3þ 3 1 m m m m 1 m m E m m 1 m D ¼ ð1 þ mþð1 mþ m 0 0 : m m and in the same notation as in (5) the stress and strain read: ð5þ Fig. 3. Total displacement (A) and first principal strain (B) in a deformed auxetic obstacle (m ¼ 0:99) after time t = 0.0 s. Fig. 4. Velocity field of fluid in a deformed domain at the steady state: (A) m ¼ 0:99 (B) m ¼ 0: (C) m ¼ 0:3 (D) m ¼ 0:3 for v 0 = 1 m/s.
5 1390 T. Strek et al. / Journal of Non-Crystalline Solids 355 (009) r ¼ 6 4 r x r y r z s xy s yz s xz 3 5 e x e y e z e ¼ c xy 6 4 c yz c xz 3 : ð6þ 5 Looking at the considered model geometry (Fig. 1) one can note that when H tends to infinity, in each plane x y along z-direction the same mechanical situation may be assumed. This means that all z components of the stress are equal to zero what encourages one to describe the problem within a two-dimensional (D) plane stress state. For finite H, Eqs. () (5) are transformed to that state by using the following averaging operation [9]: f ðx; y; tþ ¼ 1 H Z H 0 f ðx; y; z; tþdz; where f stands for each variable which should be averaged according to the plane stress mode. As the results one has, 3 r x 6 r ¼ 4 r y 5; ð8þ s xy 3 D ¼ E 1 m 0 6 m : ð9þ 1 m m Moreover [9]: ðþ l ¼ l; lk k ¼ k þ l ¼ Em 1 m : ð10þ Neglecting now the body force in (), assuming the Rayleigh damping in the obstacle [10,14] and its small deformations in the form: e ¼ 1 r u þðruþ T ; ð11þ we arrive at the final equation of motion in terms of () [3,4]: q s u r c ru þ þ a dm q ¼ 0; ð1þ where is c the flux matrix: " # E 1 m H 0 c ¼ E 0 H : ð13þ ð1þmþ Two damping coefficients have been introduced in (1): a dm is the mass damping parameter and b dk is the stiffness damping parameter [10,14]. H is the width of the obstacle and channel in z-direction (Fig. 1(A)), then r x y. The boundary conditions for (1) result from the fact that all free obstacle boundaries are loaded by the fluid as follows: n; ð14þ r n ¼ pi þ g rv þðrvþ T where n is the normal vector to the obstacle boundary. This load represents a sum of pressure and viscous forces. The edge load Fig. 5. Total displacement in a deformed obstacle at the steady state: (A) m ¼ 0:99 (B) m ¼ 0: (C) m ¼ 0:3 (D) m ¼ 0:3 for v 0 = 1 m/s.
6 T. Strek et al. / Journal of Non-Crystalline Solids 355 (009) Fig. 6. First principal strain in a deformed obstacle near the fixed boundary: (A) m ¼ 0:99, (B) m ¼ 0:, (C) m ¼ 0:3, (D) m ¼ 0:3 for v 0 = 1 m/s. has been defined as a force per area taking the thickness into account. The initial conditions for (1) read: the displacement u ¼ 0 and there is no initial stress, strain and body load on the obstacle. 3. Numerical simulations Following Fig. 1, the detailed geometry of the problem is the following: the flow channel is h ¼ 100 lm high (y-direction) and l ¼ 300 lm long (x-direction), the flank surfaces (upper and lower) of the obstacle are clamped, the obstacle dimensions: 5 lm width (x-direction), 50 lm height (y-direction), the obstacle is placed 0 lm away from the channel left boundary, the system width is H ¼ 10 m(z-direction). The remaining data have been collected in Table 1. The flow and deformation of the structure were investigated within the time interval between 0 and 0.0 s. The problem was solved using the COMSOL code with the solver DASPK [4,1,] and the direct UMFPACK linear system solver for FEM matrices [4] with 458 Lagrange quadratic triangular mesh elements. It should be noted that the total displacement and first principal strain values in all the considered examples are small comparing to the obstacle width. Since the obstacle deformation is very small we (artificially) scaled up the deformation magnitude on the figures to see the deformation details in the vicinity of the plate s end. To confirm the reliability of the results we performed the finite element analysis with different meshes. The number of mesh points and the number of mesh elements are presented in Table. It can be seen in Fig. (A) (C) that the results are practically indistinguishable for all the three fine meshes. Figs. 5 and 6 show in detail only the upper part of the deformed obstacle. All Figs. 3 6 were drawn based on the FEM calculation results using 1586 mesh elements (mesh B). Fig. 3(A) and (B) present the deformed obstacle after flow for t ¼ 0:0 s. Fig. 3A presents the total displacement in a deformed auxetic (m ¼ 0:99) obstacle. Fig. 3(B) presents the first principal strain of the obstacle. It is easy to show that the time-dependent flow with a constant inlet parabolic velocity profile reaches a steady-state after time t = 0.0 s. Figs. 4 6 present the numerical results for the flow with maximum inlet velocity v 0 ¼ 1 m/s. Fig. 4 presents the steady state velocity field and the deformed obstacles with different Poisson s ratio value (m ¼ 0:99; 0:; 0:3 and 0:3). Fig. 5 presents a total displacement in the deformed obstacle. Fig. 6 presents the first principal strain in the considered obstacle. Let us still consider how Young s modulus values are connected with the elastic properties of the obstacle both for classical and auxetic states. In Fig., we have shown a situation where that modulus increases for the material with the Poisson s ratio m ¼ 0:99.
7 139 T. Strek et al. / Journal of Non-Crystalline Solids 355 (009) Fig.. First principal strain of an auxetic obstacle vs. Young s modulus (A) E = N/m, (B) E =10 5 N/m, (C) E =10 6 N/m, (D) E =10 N/m. 4. Conclusions A finite element analysis of fluid solid interactions was applied to study the deformations of an auxetic obstacle in a D channel flow. It was found that the displacement of the obstacle points coming up to the clamped edges along the plate vanished faster than for the case of a non-auxetic material (cf. Figs., 3(A), 5(A) and (D)). Figs. 3(B) and 6(A) indicate that there is practically no deformation of the auxetic obstacle in the regions directly close to its clamped edges contrary to the situation for the non-auxetic material shown in Fig. 6(D). A surprising feature of the auxetic plate can be seen in Figs., 3(A) and 5(A) where the displacement vector firstly turns oppositely during the auxetic obstacle deformation (growing away from the clamped edges) and then it is in the direction of the fluid flow contrary to the non-auxetic case (Fig. 5(D)), where such a turn was not observed and the displacement vector had always the same direction as the fluid flow. Similar behavior was observed when analyzing the deformation of a three-dimensional plate [1]. This effect is most significant in materials of low values of Young s modulus. Finally, the obstacle with a negative Poisson s ratio influences the fluid velocity field in such a way that the values of velocity in the vicinity of the channel walls are lower compared to the classical (non-auxetic) case (Fig. 4). This means that the probability of a turbulence or even cavitation decreases when an auxetic obstacle is used. Acknowledgement A part of this work was supported by Grant N0 00 3/151 (00-010) of the Ministry of Sciences and Higher Education. References [1] K.E. Brenan, S.L. Campbell, L.R. Petzold, Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations, nd Ed., Elsevier, New York, SIAM, [] P.N. Brown, A.C. Hindmarsh, L.R. Petzold, SIAM J. Sci. Comput. 15 (1994) 146. [3] Cosmol Structure Mechanics Module, FEMLAB 3.1 User s Guide, Comsol, AB, 004. [4] Comsol Multiphysics 3.3 User s Guide, Comsol AB, 006. [5] J. Donea, A. Huerta, J.-Ph. Ponthot, A. Rodrígues-Ferran, in: Erwin Stein, René de Borst, Thomas J.R. Hughes (Eds.), Encyclopedia of Computational Mechanics, John Wiley, 004. [6] K.E. Evans, M.A. Nkansah, I.J. Hutchinson, S.C. Rogers, Nature 353 (1991) 14. [] R. Lakes, Science 35 (198) [8] L.D. Landau, E.M. Lifshits, Theory of Elasticity, Pergamon, London, [9] R.B. Hetnarski, J. Ignaczak, Mathematical Theory of Elasticity, Taylor & Francis, New York, London, 004. [10] J.F. Semblat, J. Sound Vib. 06 (5) (199) 41. [11] Ch. Smith, K.W. Wojciechowski, Phys. Stat. Solidi (b) 45 (008) 486. [1] T. Strek, B. Maruszewski, J. Narojczyk, K.W. Wojciechowski, J. Non-Cryst. Solids 354 (008) 445. [13] T.C.T. Ting, T.Y. Chen, Q. J. Mech. Appl. Math. 58 (005) 3. [14] T. Trombetti, S. Silvestri, J. Sound Vib. 9 (006) 1. [15] K.W. Wojciechowski, Phys. Lett. A13 (1989) 60. [16] O.C. Zienkiewicz, R.L. Taylor, The Finite Element Method, 5th Ed., vols. 1 3, Butterworth-Heinemann, Oxford, 000.
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationANOMALOUS DEFORMATION OF CONSTRAINED AUXETIC SQUARE
Anomalous Rev.Adv.Mater.Sci. deformation 23(2010) of constrained 169-174 auxetic square 169 ANOMALOUS DEFORMATION OF CONSTRAINED AUXETIC SQUARE Artur A. PoŸniak 1, Henryk Kamiñski 2, Przemys³aw Kêdziora
More informationFinite Element Analysis of Dynamic Properties of Thermally Optimal Two-phase Composite Structure
Vibrations in Physical Systems Vol.26 (2014) Finite Element Analysis of Dynamic Properties of Thermally Optimal Two-phase Composite Structure Abstract Maria NIENARTOWICZ Institute of Applied Mechanics,
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationComputer simulation of the Poisson's ratio of soft polydisperse discs at zero temperature
Computer simulation of the Poisson's ratio of soft polydisperse discs at zero temperature JAKUB NAROJCZYK KRZYSZTOF W. WOJCIECHOWSKI Institute of Molecular Physics Polish Academy of Sciences ul. M. Smoluchowskiego
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationForced Response of Plate with Viscoelastic Auxetic Dampers
Vibrations in Physical Systems 2018, 29, 2018003 (1 of 9) Abstract Forced Response of Plate with Viscoelastic Auxetic Dampers Tomasz STREK Poznan University of Technology, Institute of Applied Mechanics
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research education use, including for instruction at the authors institution
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationSimulation of fluid-structure interaction for an elastic cylinder in an axial flow
Advances in Fluid Mechanics IX 151 Simulation of fluid-structure interaction for an elastic cylinder in an axial flow Z. G. Liu 1, Y. Liu 1 & J. Lu 2 1 Department of Mechanical Engineering, Hong Kong Polytechnic
More informationA multiscale framework for lubrication analysis of bearings with textured surface
A multiscale framework for lubrication analysis of bearings with textured surface *Leiming Gao 1), Gregory de Boer 2) and Rob Hewson 3) 1), 3) Aeronautics Department, Imperial College London, London, SW7
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationMODELING AND FEM ANALYSIS OF DYNAMIC PROPERTIES OF THERMALLY OPTIMAL COMPOSITE MATERIALS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationAnalysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements
Acta Mech. Sin. (2016) 32(2):295 300 DOI 10.1007/s10409-015-0534-2 RESEARCH PAPER Analysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements V. H. Carneiro 1
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationNatural frequency analysis of fluid-conveying pipes in the ADINA system
Journal of Physics: Conference Series OPEN ACCESS Natural frequency analysis of fluid-conveying pipes in the ADINA system To cite this article: L Wang et al 2013 J. Phys.: Conf. Ser. 448 012014 View the
More informationPoisson s ratio of rectangular anti-chiral lattices with disorder
Poisson s ratio of rectangular anti-chiral lattices with disorder Artur A. Poźniak 1 Krzysztof W. Wojciechowski 2 1 Institute of Physics, Poznań Universty of Technology, ul. Nieszawska 13A, 60-695 Poznań,
More informationControlling elastic wave with isotropic transformation materials
Controlling elastic wave with isotropic transformation materials Zheng Chang, Jin Hu, a, Gengkai Hu, b, Ran Tao and Yue Wang School of Aerospace Engineering, Beijing Institute of Technology, 0008,Beijing,
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
This is a sample cover image for this issue. The actual cover is not yet available at this time.) This article appeared in a journal published by Elsevier. The attached copy is furnished to the author
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationThe Full-System Approach for Elastohydrodynamic Lubrication
Excerpt from the Proceedings of the COMSOL Conference 009 Milan The Full-System Approach for Elastohydrodynamic Lubrication Nicolas Fillot 1*, Thomas Doki-Thonon 1, Wassim Habchi 1 Université de Lyon,
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationAdd-on unidirectional elastic metamaterial plate cloak
Add-on unidirectional elastic metamaterial plate cloak Min Kyung Lee *a and Yoon Young Kim **a,b a Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro, Gwanak-gu, Seoul,
More informationExact and Numerical Solution of Pure Torsional Shaft
Australian Journal of Basic and Applied Sciences, 4(8): 3043-3052, 2010 ISSN 1991-8178 Exact and Numerical Solution of Pure Torsional Shaft 1 Irsyadi Yani, 2 M.A Hannan, 1 Hassan Basri, and 2 E. Scavino
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sample cover image for this issue. The actual cover is not yet available at this time.) This article appeared in a journal published by Elsevier. The attached copy is furnished to the author
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationOn the displacement of two immiscible Stokes fluids in a 3D Hele-Shaw cell
On the displacement of two immiscible Stokes fluids in a 3D Hele-Shaw cell Gelu Paşa Abstract. In this paper we study the linear stability of the displacement of two incompressible Stokes fluids in a 3D
More informationPlanar auxeticity from various inclusions
Planar auxeticity from various inclusions Artur A. Poźniak 1 Krzysztof W. Wojciechowski 2 Joseph N. Grima 3 Luke Mizzi 4 1 Institute of Physics, Poznań Universty of Technology, Piotrowo 3, 6-695 Poznań
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationMathematical Model of Blood Flow in Carotid Bifurcation
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Mathematical Model of Blood Flow in Carotid Bifurcation E. Muraca *,1, V. Gramigna 1, and G. Fragomeni 1 1 Department of Experimental Medicine
More informationAdvances in Engineering Software
Advances in Engineering Software 41 (21) 712 728 Contents lists available at ScienceDirect Advances in Engineering Software journal homepage: www.elsevier.com/locate/advengsoft The quadratic MITC plate
More informationApplication of pseudo-symmetric technique in dynamic analysis of concrete gravity dams
Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams V. Lotfi Department of Civil and Environmental Engineering, Amirkabir University, Iran Abstract A new approach is
More informationThis file contains the author's copy of:
This file contains the author's copy of: Tasnim, S.H., and Collins, M.R., "Suppressing Natural Convection in a Differentially Heated Square Cavity with an Arc Shaped Baffle", International Communications
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationEarly development drug formulation on a chip: Fabrication of nanoparticles using a microfluidic spray dryer. Supplemental Information
Early development drug formulation on a chip: Fabrication of nanoparticles using a microfluidic spray dryer Julian Thiele, a,b Maike Windbergs, a Adam R. Abate, a Martin Trebbin, b Ho Cheung Shum, a,c
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More informationSound Attenuation by Hearing Aid Earmold Tubing
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Sound Attenuation by Hearing Aid Earmold Tubing Mads J. Herring Jensen Widex A/S Ny Vestergaardsvej 19, DK-3500 Vaerloese, Denmark, mjj@widex.dk
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 FREQUENCY DEPENDENCY AND ANISOTROPY OF THE ELASTIC CONSTANTS OF (NON-)POROUS MATERIALS AND THEIR INFLUENCE ON THE USAGE IN BUILDING
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sample cover image for this issue. The actual cover is not yet available at this time. This article appeared in a journal published by Elsevier. The attached copy is furnished to the author
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationConstitutive Relations
Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sample cover image for this issue. The actual cover is not yet available at this time.) This article appeared in a journal published by Elsevier. The attached copy is furnished to the author
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationFinite Element Model of a complex Glass Forming Process as a Tool for Control Optimization
Excerpt from the Proceedings of the COMSOL Conference 29 Milan Finite Element Model of a complex Glass Forming Process as a Tool for Control Optimization Felix Sawo and Thomas Bernard Fraunhofer Institute
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationTOPOLOGY OPTIMIZATION OF NAVIER-STOKES FLOW IN MICROFLUIDICS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationstatus solidi Computer simulation of heat transfer through a ferrofluid
physica pss status solidi basic solid state physics b Computer simulation of heat transfer through a ferrofluid omasz Strek and Hubert Jopek Institute of Applied Mechanics, Poznan University of echnology,
More informationPeristaltic Pump. Introduction. Model Definition
Peristaltic Pump Introduction In a peristaltic pump, rotating rollers are squeezing a flexible tube. As the pushed down rollers move along the tube, the fluid in the tube follows the motion. The main advantage
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationThe stiffness of plates
The stiffness of plates 1. Introduction The word plate is a collective term for elements in which forces can be transferred in two directions. Floors, walls, bridge slabs and laminates are all plates.
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationComputational Modelling of Vibrations Transmission Loss of Auxetic Lattice Structure
Vibrations in Physical Systems Vol. 27 (2016) Abstract Computational Modelling of Vibrations Transmission Loss of Auxetic Lattice Structure Eligiusz IDCZAK Institute of Applied Mechanics, Poznan University
More informationFluid Mechanics. Jim Emery 9/3/ Bernoulli s Law 2. 2 Torricelli s Law 3. 3 Time to Empty a Tank 3. 4 Viscosity 4. 5 The Acceleration 5
Fluid Mechanics Jim Emery 9/3/2010 Contents 1 Bernoulli s Law 2 2 Torricelli s Law 3 3 Time to Empty a Tank 3 4 Viscosity 4 5 The Acceleration 5 6 The Equation of Motion Due to Cauchy 7 7 The Deformation
More informationNumerical Solutions for the Lévêque Problem of Boundary Layer Mass or Heat Flux
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Numerical Solutions for the Lévêque Problem of Boundary Layer Mass or Heat Flux Ekkehard Holzbecher Weierstrass Institute for Applied
More informationSTRUCTURAL ANALYSIS OF A WESTFALL 2800 MIXER, BETA = 0.8 GFS R1. By Kimbal A. Hall, PE. Submitted to: WESTFALL MANUFACTURING COMPANY
STRUCTURAL ANALYSIS OF A WESTFALL 2800 MIXER, BETA = 0.8 GFS-411519-1R1 By Kimbal A. Hall, PE Submitted to: WESTFALL MANUFACTURING COMPANY OCTOBER 2011 ALDEN RESEARCH LABORATORY, INC. 30 Shrewsbury Street
More informationDEVELOPING A MICRO-SCALE MODEL OF SOIL FREEZING
Proceedings of ALGORITMY 2016 pp. 234 243 DEVELOPING A MICRO-SCALE MODEL OF SOIL FREEZING A. ŽÁK, M. BENEŠ, AND T.H. ILLANGASEKARE Abstract. In this contribution, we analyze thermal and mechanical effects
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationEffect of interfacial dislocations on ferroelectric phase stability and domain morphology in a thin film a phase-field model
JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 4 15 AUGUST 2003 Effect of interfacial dislocations on ferroelectric phase stability and domain morphology in a thin film a phase-field model S. Y. Hu, Y. L.
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationViscosity of magmas containing highly deformable bubbles
Journal of Volcanology and Geothermal Research 105 (2001) 19±24 www.elsevier.nl/locate/jvolgeores Viscosity of magmas containing highly deformable bubbles M. Manga a, *, M. Loewenberg b a Department of
More informationBack Calculation of Rock Mass Modulus using Finite Element Code (COMSOL)
Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL) Amirreza Ghasemi 1. Introduction Deformability is recognized as one of the most important parameters governing the behavior of rock
More informationThe Finite Element Method for Computational Structural Mechanics
The Finite Element Method for Computational Structural Mechanics Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) January 29, 2010 Martin Kronbichler (TDB) FEM for CSM January
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationOptimization of the Gas Flow in a GEM Tracker with COMSOL and TENDIGEM Development. Presented at the 2011 COMSOL Conference
Optimization of the Gas Flow in a GEM Tracker with COMSOL and TENDIGEM Development Francesco Noto Presented at the 2011 COMSOL Conference V. Bellini, E. Cisbani, V. De Smet, F. Librizzi, F. Mammoliti,
More informationINTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA)
INTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA) Title 3 Column (full page) 2 Column What is Finite Element Analysis? 1 Column Half page The Finite Element Method The Finite Element Method (FEM)
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationNiiranen, Jarkko; Khakalo, Sergei; Balobanov, Viacheslav Strain gradient elasticity theories in lattice structure modelling
Powered by TCPDF (www.tcpdf.org) This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Niiranen, Jarkko; Khakalo, Sergei;
More informationFlow-Induced Vibration Analysis of Supported Pipes with a Crack
Flow-Induced Vibration Analysis of Supported Pipes with a Crack Jin-Hyuk Lee 1 *, Samer Masoud Al-Said,3 1 Department of Mechanical Engineering, American University of Sharjah, Sharjah, UAE, Department
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationBack Analysis of Measured Displacements of Tunnels
Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe
More informationAA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations
AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationMultiphysics Analysis of Electromagnetic Flow Valve
Multiphysics Analysis of Electromagnetic Flow Valve Jeffrey S. Crompton *, Kyle C. Koppenhoefer, and Sergei Yushanov AltaSim Technologies, LLC *Corresponding author: 13 E. Wilson Bridge Road, Suite 14,
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationModeling of Humidification in Comsol Multiphysics 4.4
Modeling of Humidification in Comsol Multiphysics 4.4 Indrajit Wadgaonkar *1 and Suresh Arikapudi 1 1 Tata Motors Ltd. Pimpri, Pune, India, 411018. *Corresponding author: Indrajit Wadgaonkar, Tata Motors
More information