Improvement of Calculation Stability for Slow Fluid Flow Analysis Using Particle Method *
|
|
- Loreen Reynolds
- 6 years ago
- Views:
Transcription
1 Materials Transactions, Vol. 58, No. 3 (2017) pp. 479 to Japan Foundry Engineering Society Improvement of Calculation Stability for Slow Fluid Flow Analysis Using Particle Method * Naoya Hirata and Koichi Anzai Department of Metallurgy, Graduate School of Engineering, Tohoku University, Sendai , Japan MPS (Moving particle semi-implicit) method is one of the popular particle methods. It is well known that calculation instability occurs especially for a slow fluid flow analysis, and the instability appears as unnatural pressure and velocity oscillation. This is because the MPS method adopts predictor-corrector method, and tentative position of particles are not corrected enough at the corrector phase. Therefore an improved multiple relaxation method is proposed in this study to improve the calculation stability. Pressure and velocity correction is conducted multiply in the correction phase, and also adjusting coefficient is adopted in the source term in the Poisson equation for the pressure. As a result, the proposed method improved the pressure and velocity distribution instability. The method was applied to complex shaped casting with a shoulder. The conventional method showed unnatural high pressure and the fluid hopped, whereas smooth and natural pressure field could be obtained using the proposed method. [doi: /matertrans.f-m ] (Received October 21, 2016; Accepted December 2, 2016; Published February 25, 2017) Keywords: flow, simulation, particle method, moving particle semi-implicit (MPS) method 1. Introduction Recently, particle methods, which are based on a fully Lagrangian method, have been developing to calculate complex phenomena occurring in the casting processes 1 4). Nowadays, the typical particle methods for computational continuum dynamics are usually categorized into two methods, SPH (Smoothed Particle Hydrodynamics) 2) method and MPS (Moving Particle Semi-implicit) 3,4) method. Discrete objects are used as calculation elements in the particle methods, and they can be placed and moved freely in the space. This feature allows the particle methods to simulate the heat and mass transfer phenomena that are observed in the casting process more easily and directly than other methods that use the calculation lattice. However, the particle methods have still some problems. Calculation stability is one of the significant problems in the particle methods. Casting process involves a variety of phenomena which occur simultaneously from the pouring to the solidification. Because it is difficult to expect what kind of phenomena will occur at where and at which timing, higher calculation stability is required for the integrated simulation using particle methods. It is known that the oscillation of velocity and pressure occurs in the flow simulation, and the oscillation decreases the calculation stability especially for a slow flow 5). Hirata et al. have proposed a stable and rapid calculation method for the slow flow by ignoring inertia force and adjusting the gravitational force after the fluid flow can be assumed to be almost stopped 5). However, the procedure requires precise estimation of the flow stop time. The original flow calculation by the particle method 3) is unstable for the slow fluid flow. Therefore the flow calculation without any improvement will decrease the calculation stability when the flow becomes slow during the casting processes. Recently, Hirata et al. reported the multiple relaxation method which improves calculation stability and speed 6) in the case of slow * This Paper was Originally Published in Japanese in J. JFS. 87 (2015) fluid flow. However, the oscillation of pressure and velocity are still not negligible. Therefore further improvement in the accuracy is expected. In this study, we tried further improvement of the calculation stability and accuracy of the particle method in the case of slow fluid flow by improving the multiple relaxation method. 2. Numerical Method We adopted MPS method to combine flow and solidification simulations. In the MPS method, the calculations such as gradient operation are evaluated by calculating an interaction among the surrounding particles. The interaction is calculated with the particles by weighted averaging within a certain distance, which is referred to as the kernel size. A summary of the proposed calculation algorithms is shown in this section. 2.1 Weight function and particle interaction models A weight function w is used in the MPS method to calculate the interaction of particles by weighted averaging. Equation (1) was used in this study because of its higher stability for the flow calculation 7). w(r, c k r 0 ) r2 2 7πc k c k2 r + 6r3 (0 r < 0.5c 2 0 c k3 r 3 k r 0 ) 0 = πc 2r 3 (0.5c k r 0 r < c k r 0 ) k c k r 0 0 (c k r 0 < r) Here, r is a distance between particles i and j. c k r 0 is the kernel size which represents interaction calculation range. The MPS method uses particle size r 0 which has the same significance as the lattice size in the Eulerian methods such as FDM (finite difference method). c k is the kernel size coefficient and usually varies between 2 and 4. c k = 2.1 was used in this study as a commonly used value 3,7). Particle number density n i of the particle i is used in the particle method to calculate the interaction between particle i (1)
2 480 N. Hirata and K. Anzai and the surrounding particles. n i is the sum of the weight function of the particles surrounding the particle i, and is expressed as follows. n i = i j w r j r i, c k r 0 (2) Here, r i and r j are the position vectors of particle i and j, respectively. Assuming that the fluid is incompressible and considering that the particle number density is directly related to the fluid density, n i is equal to n 0 in the case that a particle has no surface particle around it. And we can use n 0 for the mass conservation condition in the incompressible fluid flow analysis using the MPS method. Particle interaction models are used to describe the differential operators in the MPS method. If ϕ i and ϕ j are arbitrary scalars at positions r i and r j, then the particle interaction models for the differential operators at the particle i can be expressed as follows. φ i = d (φ i φ j ) (r j r i ) w r j r i, c k r 0 (3) n i i j r j r i r j r i 2 φ i = 2d (φ i φ j ) n w r 2 j r i, c k r 0 (4) i i j r j r i Equation (3) is a gradient model, and eq. (4) is a Laplacian model. Suffixes i and j represent the assigned numbers of particles. The d is the number of space dimensions. 2.2 Flow simulation The governing equations for incompressible flow calculation by the MPS method consist of the continuity equation and the Navier-Stokes equation as follows. Dρ Dt = 0, Du Dt = 1 ρ p + ν 2 u + f (5) Here, ρ indicates density (kg/m 3 ); u, velocity vector (m/s); p, pressure (Pa); ν, kinematic viscosity (m 2 /s); and f (m/s 2 ), the body force acceleration vector including gravity. The lefthand side of the second equation, D/Dt, is the Lagrangian derivative, having the same significance as temporal differentiation. Flow calculation by the MPS method is based on the predictor-corrector method similar to that in the case of the FDM or other Eulerian methods. The major differences between the MPS method and the FDM are the formulation of the Poisson equation of pressure and the calculation of particle position. In the MPS method, the tentative particle position is calculated using tentative velocity in the prediction phase. Therefore, the correction in position is calculated along with the correction in velocity in the correction phase. The Poisson equation for pressure in the MPS method is described using the tentative particle number density n* and n 0 as follows 1). 2 p k+1 = ρ t 2 n n 0 n 0 (6) Here, n* is the tentative particle number density calculated using the tentative particle positions r* and eq. (2). In the MPS method, the pressure is calculated so as to make the tentative particle number density n* to the constant value n 0, whereas making the divergence of velocity to zero in the FDM. In the flow simulation with a free surface, the pressure at the free surface is assumed to be zero. The criteria to determine the free surface particles is described as following equation because the particle number density decrease at the free surface. n i < β (7) n 0 β = 0.97 is usually used in the MPS method with a different weight function to the eq. (1) 3). In this study, β = was used for eq. (1), which shows the most natural free surface particle criteria. The time increment Δt is calculated by using Courant number C n using the following equation. t = C n r 0 u max Here, u max is the maximum velocity of the particles. It is known that the use of small Δt sometimes worsens the calculation stability. This is because the source term of the Poisson equation includes Δt 2 in its denominator, and too small Δt produces unnatural high pressure. In this study, constant value C n = 0.1 was used. 3. Improvement of Calculation Stability for Slow Fluid Flow The authors have reported that the multiple relaxation method improves calculation stability and speed in the previous report 6). However unnatural pressure oscillation and unnatural movement of the surface particle are still at issues. Pressure oscillation should be avoided as far as possible to predict a feeding behavior to the solidification shrinkage precisely. Therefore an improved multiple relaxation method was proposed in this study to avoid unnatural pressure oscillation. 3.1 Improved multiple relaxation method The Poisson equation of the pressure in the MPS method uses the source term based on the particle number density n i, and the tentative value n* is expected to be corrected to the constant value n 0. Therefore, if an excess concentration of particles occurs in the prediction phase, excess positional correction occurs because of excess high pressure. Such excess correction is one of the reasons of velocity and pressure oscillation occurred in the particle method calculation. The multiple relaxation method improves the stability and accuracy of flow simulation using the particle method by calculating the Poisson equation more than once, so as to the particle number density n i become closer to the constant value n 0 6). However, a possibility of excess correction of the position is still remain even by using multiple relaxations; sometimes the excess correction will let particles move a distance more than expected, and the oscillation of velocity and position will occur during the multiple relaxation procedures. Therefore, an artificial adjustment coefficient of the source term in the Poisson equation c poi was introduced to reduce the amount of positional correction in one relaxation process to settle the particles at (8)
3 Improvement of Calculation Stability for Slow Fluid Flow Analysis Using Particle Method 481 Table 1 Physical properties for simulation. Casting Materials Pure Al Density, kg/m Kinematic Viscosity, m 2 /s 10 5 Table 2 Calculation condition. Δt max (s) MR (Multiple relaxation) c poi Fig. 1 Flow chart for the flow simulation. Fig. 3 Pressure distribution at 0.1 s. Fig. 2 Calculation model. the expected position without excess movement. The flowchart for the flow simulation with the improved multiple relaxation method is shown in the Fig. 1. The Poisson equation is solved more than once in the multiple relaxation method. Accordingly the velocity and position are also corrected more than once during one time step. The improved method uses the adjustment coefficient c poi in the source term in the Poisson equation for the pressure. 3.2 Calculation model Two-dimensional flow calculation is conducted for still fluid as shown in the Fig. 2. The width of casting (fluid) is 40 mm and the height is 80 mm. Particle size is 2 mm; thus, we use 800 particles for the casting, 384 particles for the mold. The calculation is conducted for 2 seconds. The maximum time increment Δt max, the number of multiple relaxations (MR) and the adjustment coefficient in the source term of the Poisson equation c poi, and the pressure value at three points were measured. The measuring points were at the center in the horizontal directions and the depth from the casting surface were h = 30, mm. The physical properties and the calculation conditions are shown in Table 1 and 2. The relationship between Δt max and the number of multiple relaxation is determined to optimize the stability and the efficiency of calculation referring the past research 6). Four c poi values (c poi = 1.0, 0.5, 0.2, 0.1) were used for each MR and Δt max conditions. 3.3 Results Pressure Figure 3 shows a pressure distribution at 0.1 s. From now on, the calculation conditions are expressed using the combination of the number of multiple relaxation and c poi ; for example, in the case of two times of multiple relaxation and c poi = 0.5, the condition is expressed as MR An ideal hydrostatic pressure at the bottom of the casting is 1960 Pa. Figure 3(a) is the result in the case of the basic MPS method without multiple relaxations and c poi = 1.0. The pressure distribution was non-uniform, and the maximum value exceeded 5000 Pa. Such non-uniform distribution was not settled during 2 s calculation, and the particles continued to flow in
4 482 N. Hirata and K. Anzai whole unnaturally. Figure 3(b) is the result in the case of two times multiple relaxations with c poi = 1.0. The pressure at the bottom was around 2000 Pa which is similar to the ideal value, however remains non-uniform distribution. Figure 3(c) is the results in the case of five times multiple relaxations and c poi = 0.5. In this case, significant improvement was found compared with the former two results. We could obtain almost similar results under the conditions shown in the Table 2, except for MR1-1.0 and MR2-1.0 as shown in the Fig. 3(a) and (b). Next, an average pressure value during 2 s at h = 70 mm and a standard deviation of the pressure are shown in the Fig. 4. Black dots are the average value and the error bars correspond to the standard deviation. The ideal pressure value at h = 70 mm is 1715 Pa. The average values indicated by the black dots agreed well with the ideal value in each condition. On the other hand, the standard deviation of MR1-1.0 and MR2-1.0 were larger than the other conditions, showing a similar tendency in the Fig. 3. As a result, the ideal pressure distribution could be achieved in the case that the standard deviation was less than 200 Pa in this study. For further discussion, the pressure change in time under typical six conditions are shown in Fig. 5. The figures show the pressure oscillation at h = 30, 50, 70 mm. Figure 5(a) to Fig. 4 Pressure at h = 0.07 m. Fig. 5 Pressure oscillation.
5 Improvement of Calculation Stability for Slow Fluid Flow Analysis Using Particle Method 483 (c) show the results in the case of c poi = 1.0. The results show that the multiple relaxation reduced the pressure oscillation, and the effect was more obvious if the number of multiple relaxations is higher. Figure 5(d) to (f) are the case of c poi = 0.5. Figure 5(d) shows that the pressure oscillation could not be reduced enough without multiple relaxations. From the Fig. 5(e) and (f), the combination of the multiple relaxations and adjustment coefficient c poi significantly reduces the pressure oscillation. A relationship between calculation conditions and the fluid behavior was summarized in Fig. 6. The circle indicates the conditions showing the excellent stability of still fluid with less pressure oscillation. Cross mark means poor stability conditions which show an unnatural rapid increase of pressure, and the calculation resulted in failure with a divergence of the velocity field. Diamond shows the conditions which do not result in failure, however, unnatural behavior such as hopping of fluid occurs. Overall, a use of multiple relaxations and less c poi conditions improved calculation stability. Whereas in the case of 2 times multiple relaxations, too small c poi could reduce the stability. This is because 2 times multiple relaxations is not enough for pressure, velocity and positional correction in the correction phase if c poi is too small Velocity Figure 7 shows a variation of the maximum velocity of particles. Figure 7(a) and (b) are the results in the case of MR1 and MR5, respectively. The results in the case of MR2 and MR10 were similar to the case of MR5. In the case without multiple relaxations shown in the Fig. 7(a), the use of smaller c poi decrease the maximum velocity, whereas the values were larger than any case with multiple relaxations shown in the Fig. 7(b). However, the maximum velocity was over 0.05 m/s and the average velocity was over 0.02 m/s even using the multiple relaxations. Therefore, if the maximum or average velocity is expected to be under such values, some artificial operation will be useful; such as ignoring the inertia force of fluid particles after pouring as reported by Hirata et al. 5) Pressure stability in complex shaped casting Finally, we discuss the pressure stability in the complex shaped casting. In this section, three-dimensional flow calculation was conducted for a cylindrical casting with a shoulder. Fig. 7 Maximum velocity of particles. Fig. 6 Effect of multiple relaxation and source term adjustment on the pressure stability of the flow simulation for gently placed fluid. Circle indicates excellent stability, diamond shows good stability but large oscillation of pressure occurs, cross means poor stability. Fig. 8 Pressure distribution in the casting with shoulder.
6 484 N. Hirata and K. Anzai The height of the casting was 180 mm, and the upper half of the casting was 25 mm and the lower half was 40 mm in diameter, respectively. The physical properties used in the calculation are shown in the Table 1, and the fluid was poured from the upper side with a rate of 100 g/s and pouring finish within 4 s. The results are shown in the Fig. 8. Figure 8(a) is the result of MR1-1.0 at 4.5 s, showing unnatural pressure distribution at the bottom of the casting, and hopping of the casting particles. At the last moment, the fluid flow was very slow. However the pressure in the casting was unnaturally high under the shoulder. In such shapes with shoulders, positional correction in the correction phase tends to be insufficient, and the pressure under the shoulder will be accumulated unnaturally. As a result, unnatural hopping occurs, and the calculation will result in failure. In the Fig. 8(b), we can find that the proposed method improves the calculation stability. The calculated pressure value at the bottom agreed well with the theoretical value, 4410 Pa. It is possible to improve calculation stability by using smaller Δt. However the pressure or velocity oscillation problem will remain, and also, longer calculation time will be required because of higher maximum velocity. The series of results show that the multiple relaxation method and the use of adjustment coefficient c poi are effective to reduce unnatural oscillation of the pressure and the velocity, and to improve calculation stability for slow fluid flow, even for the complex shaped casting with the shoulder. Improved multiple relaxation method can improve the unnatural velocity decay and calculation stability for faster fluid flow. The results of the faster flow calculation are discussed in the next report 8). 4. Conclusions The improved multiple relaxation method was proposed to improve calculation stability for slow fluid flow based on the MPS method; one of the particle methods. The multiple relaxation method is combined with the adjustment coefficient of the source term in the Poisson equation of pressure. As a result, the pressure oscillation was significantly reduced by the proposed method, and the calculation stability was also improved. The proposed method can calculate stably even for the complex shaped casting with the shoulder. Acknowledgement The authors wish to thank the young researchers support program 2012 of JFS for its support. REFERENCES 1) S.Koshizuka: Suuchi Ryutai Rikigaku, (Baifukan, 1997). 2) S. Koshizuka: Ryushihou Simulation, (Baifukan, 2008). 3) S. Koshizuka: Ryushihou, (Maruzen, 2005). 4) S. Koshizuka:, Ryushihou Nyumon, (Maruzen, 2014). 5) N. Hirata and K. Anzai: Mater. Trans. 52 (2011) ) N. Hirata and K. Anzai: J.JFS 86 (2014) ) N. Hirata and K. Anzai: J.JFS 83 (2011) ) N. Hirata and K. Anzai: J.JFS 88 (2016)
Math background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationFluid-soil multiphase flow simulation by an SPH-DEM coupled method
Fluid-soil multiphase flow simulation by an SPH-DEM coupled method *Kensuke Harasaki 1) and Mitsuteru Asai 2) 1), 2) Department of Civil and Structural Engineering, Kyushu University, 744 Motooka, Nishi-ku,
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 informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
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 informationMulti-physics CFD simulation of three-phase flow with MPS method
APCOM & ISCM 11-14 th December, 2013, Singapore Abstract Multi-physics CFD simulation of three-phase flow with MPS method *Ryouhei Takahashi¹, Makoto Yamamoto 2 and Hiroshi Kitada 1 1 CMS Corporation,
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationEffect of Initial Geometry on Surface Flatness after Arc Welding Analyzed with MPS Method
Journal of Mechanics Engineering and Automation 5 (2015) 63-67 doi: 10.17265/2159-5275/2015.02.001 D DAVID PUBLISHING Effect of Initial Geometry on Surface Flatness after Arc Welding Analyzed with MPS
More informationLecture 3: 1. Lecture 3.
Lecture 3: 1 Lecture 3. Lecture 3: 2 Plan for today Summary of the key points of the last lecture. Review of vector and tensor products : the dot product (or inner product ) and the cross product (or vector
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationSmoothed Particle Hydrodynamics (SPH) Huamin Wang
Smoothed Particle Hydrodynamics (SPH) Huamin Wang Fluid Representation Fluid is represented using a set of particles. Particle i has position x i, velocity v i, and mass m i. Animation Example Using 10
More information1/3/2011. This course discusses the physical laws that govern atmosphere/ocean motions.
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationDan s Morris s Notes on Stable Fluids (Jos Stam, SIGGRAPH 1999)
Dan s Morris s Notes on Stable Fluids (Jos Stam, SIGGRAPH 1999) This is intended to be a detailed by fairly-low-math explanation of Stam s Stable Fluids, one of the key papers in a recent series of advances
More informationSmooth Particle Hydrodynamic (SPH) Presented by: Omid Ghasemi Fare Nina Zabihi XU Zhao Miao Zhang Sheng Zhi EGEE 520
Smooth Particle Hydrodynamic (SPH) Presented by: Omid Ghasemi Fare Nina Zabihi XU Zhao Miao Zhang Sheng Zhi EGEE 520 OUTLINE Ø Introduction and Historical Perspective: Ø General Principles: Ø Governing
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 informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationPhysics-Based Animation
CSCI 5980/8980: Special Topics in Computer Science Physics-Based Animation 13 Fluid simulation with grids October 20, 2015 Today Presentation schedule Fluid simulation with grids Course feedback survey
More informationLaminar Boundary Layers. Answers to problem sheet 1: Navier-Stokes equations
Laminar Boundary Layers Answers to problem sheet 1: Navier-Stokes equations The Navier Stokes equations for d, incompressible flow are + v ρ t + u + v v ρ t + u v + v v = 1 = p + µ u + u = p ρg + µ v +
More informationSIMULATION OF THE RESITIVE FORCES ACTING ON THE BUCKET OF WHEEL LOADER BY USE OF DEM
SIMULATION OF THE RESITIVE FORCES ACTING ON THE BUCKET OF WHEEL LOADER BY USE OF DEM Hiroshi TAKAHASHI Department of Geoscience and Technology Graduate School of Engineering Tohoku University Sendai 980-8579,
More informationRational derivation of the Boussinesq approximation
Rational derivation of the Boussinesq approximation Kiyoshi Maruyama Department of Earth and Ocean Sciences, National Defense Academy, Yokosuka, Kanagawa 239-8686, Japan February 22, 2019 Abstract This
More informationGS388 Handout: Radial density distribution via the Adams-Williamson equation 1
GS388 Handout: Radial density distribution via the Adams-Williamson equation 1 TABLE OF CONTENTS ADIABATIC COMPRESSION: THE ADAMS WILLIAMSON EQUATION...1 EFFECT OF NON-ADIABATIC TEMPERATURE GRADIENT...3
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationDYNAMICS OF LIQUEFIED SEDIMENT FLOW. Advances in Natural and Technological Hazards Research Vol. 19
DYNAMICS OF LIQUEFIED SEDIMENT FLOW Advances in Natural and Technological Hazards Research Vol. 9 THE DYNAMICS OF LIQUEFIED SEDIMENT FLOW UNDERGOING PROGRESSIVE SOLIDIFICATION S. SASSA Disaster Prevention
More informationd v 2 v = d v d t i n where "in" and "rot" denote the inertial (absolute) and rotating frames. Equation of motion F =
Governing equations of fluid dynamics under the influence of Earth rotation (Navier-Stokes Equations in rotating frame) Recap: From kinematic consideration, d v i n d t i n = d v rot d t r o t 2 v rot
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationDEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationNumerical Simulation of the Hagemann Entrainment Experiments
CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationAn improved MPS method for numerical simulations of convective heat transfer problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2006; 51:31 47 Published online 15 November 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/d.1106
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 information1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed in an inertial system = rate of change of Ua following the motion in an inertial
More informationSmoothed Particle Hydrodynamics (SPH) 4. May 2012
Smoothed Particle Hydrodynamics (SPH) 4. May 2012 Calculating density SPH density estimator Weighted summation over nearby particles: ρ(r) = N neigh b=1 m bw (r r b, h) W weight function with dimension
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationDaniel J. Jacob, Models of Atmospheric Transport and Chemistry, 2007.
1 0. CHEMICAL TRACER MODELS: AN INTRODUCTION Concentrations of chemicals in the atmosphere are affected by four general types of processes: transport, chemistry, emissions, and deposition. 3-D numerical
More informationSPH Molecules - a model of granular materials
SPH Molecules - a model of granular materials Tatiana Capone DITS, Univeristy of Roma (la Sapienza) Roma, Italy Jules Kajtar School of Mathematical Sciences Monash University Vic. 3800, Australia Joe Monaghan
More informationFluid Mechanics. du dy
FLUID MECHANICS Technical English - I 1 th week Fluid Mechanics FLUID STATICS FLUID DYNAMICS Fluid Statics or Hydrostatics is the study of fluids at rest. The main equation required for this is Newton's
More informationJournal of Fluid Science and Technology
Science and Technology Effect of Molecular Diffusivities on Countergradient Scalar Transfer in a Strong Stable Stratified Flow (Study on the Linear and Nonlinear Processes by using RDT) Kouji NAGATA, Takashi
More informationThe effect of natural convection on solidification in tall tapered feeders
ANZIAM J. 44 (E) ppc496 C511, 2003 C496 The effect of natural convection on solidification in tall tapered feeders C. H. Li D. R. Jenkins (Received 30 September 2002) Abstract Tall tapered feeders (ttfs)
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationEffect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process
, pp. 844 850 Effect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process Baokuan LI and Fumitaka TSUKIHASHI 1) Department of Thermal Engineering, The School
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More informationParticle-based Fluids
Particle-based Fluids Particle Fluids Spatial Discretization Fluid is discretized using particles 3 Particles = Molecules? Particle approaches: Molecular Dynamics: relates each particle to one molecule
More informationEquivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme
Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,
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 information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationConservation of Mass Conservation of Energy Scaling Analysis. ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationCSCI1950V Project 4 : Smoothed Particle Hydrodynamics
CSCI1950V Project 4 : Smoothed Particle Hydrodynamics Due Date : Midnight, Friday March 23 1 Background For this project you will implement a uid simulation using Smoothed Particle Hydrodynamics (SPH).
More informationSink particle accretion test
Sink particle accretion test David A. Hubber & Stefanie Walch 1 Objectives Simulate spherically-symmetric Bondi accretion onto a sink particle for an isothermal gas. Calculate the accretion rate onto a
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationNumerical simulation of a small-scale snow avalanche tests using non-newtonian SPH model
A2, Vol. 70, No. 2 Vol. 17, I_681-I_690, 2014. Numerical simulation of a small-scale snow avalanche tests using non-newtonian SPH model Ahmed M. ABDELRAZEK 1, Ichiro KIMURA 2, and Yasuyuki SHIMIZU 3 1
More informationEXTENDED FREE SURFACE FLOW MODEL BASED ON THE LATTICE BOLTZMANN APPROACH
METALLURGY AND FOUNDRY ENGINEERING Vol. 36, 2010, No. 2 Micha³ Szucki*, Józef S. Suchy***, Pawe³ ak*, Janusz Lelito**, Beata Gracz* EXTENDED FREE SURFACE FLOW MODEL BASED ON THE LATTICE BOLTZMANN APPROACH
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationSimulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls
Fluid Structure Interaction V 85 Simulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls K. Fujita Osaka City University,
More informationNumerical simulation of landslide impulsive waves by incompressible smoothed particle hydrodynamics
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2008; 56:209 232 Published online 30 May 2007 in Wiley InterScience (www.interscience.wiley.com)..1526 Numerical simulation
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationLiquid Sloshing in a Rotating, Laterally Oscillating Cylindrical Container
Universal Journal of Mechanical Engineering 5(3): 97-101, 2017 DOI: 10.13189/ujme.2017.050304 http://www.hrpub.org Liquid Sloshing in a Rotating, Laterally Oscillating Cylindrical Container Yusuke Saito,
More informationHeat Transfer Benchmark Problems Verification of Finite Volume Particle (FVP) Method-based Code
PROCEEDING OF 3 RD INTERNATIONAL CONFERENCE ON RESEARCH, IMPLEMENTATION AND EDUCATION OF MATHEMATICS AND SCIENCE YOGYAKARTA, 16 17 MAY 2016 Heat Transfer Benchmark Problems Verification of Finite Volume
More informationTreecodes for Cosmology Thomas Quinn University of Washington N-Body Shop
Treecodes for Cosmology Thomas Quinn University of Washington N-Body Shop Outline Motivation Multipole Expansions Tree Algorithms Periodic Boundaries Time integration Gravitational Softening SPH Parallel
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationFluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows
Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 17 Laminar and Turbulent flows Welcome back to the video course on fluid mechanics. In
More information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationSummary We develop an unconditionally stable explicit particle CFD scheme: Boltzmann Particle Hydrodynamics (BPH)
Summary Steps in BPH Space is divided into Cartesian cells Finite number of particles ~10 6-10 7 Particles fly freely between t n and t n+1 Mass, momentum and total energy are conserved Relax into a (LTE)
More informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationEffect of Carrier Gas Flow Behavior on Performance of Separation by Using Ultrasonic Atomization
Effect of Carrier Gas Flow Behavior on Performance of Separation by Using Ultrasonic Atomization Yoshiyuki Bando 1, Keiji Yasuda 1, Akira Matsuoka 1 and Yasuhito Kawase 2 1. Department of Chemical Engineering,
More informationChapter 4: Fundamental Forces
Chapter 4: Fundamental Forces Newton s Second Law: F=ma In atmospheric science it is typical to consider the force per unit mass acting on the atmosphere: Force mass = a In order to understand atmospheric
More informationThe Johns Hopkins Turbulence Databases (JHTDB)
The Johns Hopkins Turbulence Databases (JHTDB) HOMOGENEOUS BUOYANCY DRIVEN TURBULENCE DATA SET Data provenance: D. Livescu 1 Database Ingest and Web Services: C. Canada 1, K. Kalin 2, R. Burns 2 & IDIES
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationIntroduction to Fluid Dynamics
Introduction to Fluid Dynamics Roger K. Smith Skript - auf englisch! Umsonst im Internet http://www.meteo.physik.uni-muenchen.de Wählen: Lehre Manuskripte Download User Name: meteo Password: download Aim
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
More informationPredictor-Corrector Finite-Difference Lattice Boltzmann Schemes
Applied Mathematical Sciences, Vol. 9, 2015, no. 84, 4191-4199 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.5138 Predictor-Corrector Finite-Difference Lattice Boltzmann Schemes G. V.
More informationESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria
ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr DeCaria References: An Introduction to Dynamic Meteorology, Holton MOMENTUM EQUATIONS The momentum equations governing the ocean or atmosphere
More informationINDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226
INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence
More informationKELVIN-HELMHOLTZ INSTABILITY BY SPH
II International Conference on Particle-based Methods - Fundamentals and Applications PARTICLES 2011 E. Oñate and D.R.J. Owen (Eds) KELVIN-HELMHOLTZ INSTABILITY BY SPH M. S. Shadloo, M. Yildiz Faculty
More informationA fundamental study of the flow past a circular cylinder using Abaqus/CFD
A fundamental study of the flow past a circular cylinder using Abaqus/CFD Masami Sato, and Takaya Kobayashi Mechanical Design & Analysis Corporation Abstract: The latest release of Abaqus version 6.10
More informationFluid Mechanics Abdusselam Altunkaynak
Fluid Mechanics Abdusselam Altunkaynak 1. Unit systems 1.1 Introduction Natural events are independent on units. The unit to be used in a certain variable is related to the advantage that we get from it.
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationMOVING PARTICLE SEMI-IMPLICIT METHOD: FULLY LAGRANGIAN ANALYSIS OF INCOMPRESSIBLE FLOWS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 ECCOMAS MOVING PARTICLE SEMI-IMPLICIT METHOD: FULLY LAGRANGIAN ANALYSIS OF INCOMPRESSIBLE
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationFundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
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 information