Finite pointset method

Last updated

In applied mathematics, the name finite pointset method is a general approach for the numerical solution of problems in continuum mechanics, such as the simulation of fluid flows. In this approach (often abbreviated as FPM) the medium is represented by a finite set of points, each endowed with the relevant local properties of the medium such as density, velocity, pressure, and temperature. [1]

Contents

The sampling points can move with the medium, as in the Lagrangian approach to fluid dynamics or they may be fixed in space while the medium flows through them, as in the Eulerian approach. A mixed Lagrangian-Eulerian approach may also be used. The Lagrangian approach is also known (especially in the computer graphics field) as particle method.

Finite pointset methods are meshfree methods and therefore are easily adapted to domains with complex and/or time-evolving geometries and moving phase boundaries (such as a liquid splashing into a container, or the blowing of a glass bottle) without the software complexity that would be required to handle those features with topological data structures. They can be useful in non-linear problems involving viscous fluids, heat and mass transfer, linear and non-linear elastic or plastic deformations, etc.

Description

In the simplest implementations, the finite point set is stored as an unstructured list of points in the medium. In the Lagrangian approach the points move with the medium, and points may be added or deleted in order to maintain a prescribed sampling density. The point density is usually prescribed by a smoothing length defined locally. In the Eulerian approach the points are fixed in space, but new points may be added where there is need for increased accuracy. So, in both approaches the nearest neighbors of a point are not fixed, and are determined again at each time step.

Advantages

This method has various advantages over grid-based techniques; for example, it can handle fluid domains, which change naturally, whereas grid based techniques require additional computational effort. The finite points have to completely cover the whole flow domain, i.e. the point cloud has to fulfill certain quality criteria (finite points are not allowed to form “holes” which means finite points have to find sufficiently numerous neighbours; also, finite points are not allowed to cluster; etc.).

The finite point cloud is a geometrical basis, which allows for a numerical formulation making FPM a general finite difference idea applied to continuum mechanics. That especially means, if the point reduced to a regular cubic point grid, then FPM would reduce to a classical finite difference method. The idea of general finite differences also means that FPM is not based on a weak formulation like Galerkin's approach. Rather, FPM is a strong formulation which models differential equations by direct approximation of the occurring differential operators. The method used is a moving least squares idea which was especially developed for FPM.

History

In order to overcome the disadvantages of the classical methods many approaches have been developed to simulate such flows. [2] [3] [4] [5] [6] [7] A classical grid free Lagrangian method is Smoothed Particle Hydrodynamics (SPH), which was originally introduced to solve problems in astrophysics. [8] [9]

It has since been extended to simulate the compressible Euler equations in fluid dynamics and applied to a wide range of problems. [10] [11] [12] The method has also been extended to simulate inviscid incompressible free surface flows. [13] . The implementation of the boundary conditions is the main problem of the SPH method.

Another approach for solving fluid dynamic equations in a grid free framework is the moving least squares or least squares method. [1] [14] [15] [16] [17] [7] With this approach boundary conditions can be implemented in a natural way just by placing the finite points on boundaries and prescribing boundary conditions on them. [15] The robustness of this method is shown by the simulation results in the field of airbag deployment in car industry. Here, the membrane (or boundary) of the airbag changes very rapidly in time and takes a quite complicated shape (Kuhnert et al. 2000).

Tiwari et al. (2003) performed simulations of incompressible flows as the limit of the compressible Navier–Stokes equations with some stiff equation of state. [18] This approach was first used in Monaghan (1992) to simulate incompressible free surface flows by SPH. The incompressible limit is obtained by choosing a very large speed of sound in the equation of state such that the Mach number becomes small. However, the large value of the speed of sound restricts the time step to be very small due to the CFL-condition. [10]

The projection method of Chorin is a widely used approach to solve problems governed by the incompressible Navier–Stokes equation in a grid based structure. [19] In Tiwari et al. (2001), this method has been applied to a grid free framework with the help of the weighted least squares method. The scheme gives accurate results for the incompressible Navier–Stokes equations. The occurring Poisson equation for the pressure field is solved by a grid free method. It has been shown that the Poisson equation can be solved accurately by this approach for any boundary conditions. The Poisson solver can be adapted to the weighted least squares approximation procedure with the condition that the Poisson equation and the boundary condition must be satisfied on each finite point. This is a local iteration procedure. [17]

Software

Related Research Articles

<span class="mw-page-title-main">Computational fluid dynamics</span> Analysis and solving of problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

<span class="mw-page-title-main">Particle-in-cell</span> Mathematical technique used to solve a certain class of partial differential equations

In plasma physics, the particle-in-cell (PIC) method refers to a technique used to solve a certain class of partial differential equations. In this method, individual particles in a Lagrangian frame are tracked in continuous phase space, whereas moments of the distribution such as densities and currents are computed simultaneously on Eulerian (stationary) mesh points.

<span class="mw-page-title-main">Smoothed-particle hydrodynamics</span> Method of hydrodynamics simulation

Smoothed-particle hydrodynamics (SPH) is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and fluid flows. It was developed by Gingold and Monaghan and Lucy in 1977, initially for astrophysical problems. It has been used in many fields of research, including astrophysics, ballistics, volcanology, and oceanography. It is a meshfree Lagrangian method, and the resolution of the method can easily be adjusted with respect to variables such as density.

In numerical analysis, adaptive mesh refinement (AMR) is a method of adapting the accuracy of a solution within certain sensitive or turbulent regions of simulation, dynamically and during the time the solution is being calculated. When solutions are calculated numerically, they are often limited to predetermined quantified grids as in the Cartesian plane which constitute the computational grid, or 'mesh'. Many problems in numerical analysis, however, do not require a uniform precision in the numerical grids used for graph plotting or computational simulation, and would be better suited if specific areas of graphs which needed precision could be refined in quantification only in the regions requiring the added precision. Adaptive mesh refinement provides such a dynamic programming environment for adapting the precision of the numerical computation based on the requirements of a computation problem in specific areas of multi-dimensional graphs which need precision while leaving the other regions of the multi-dimensional graphs at lower levels of precision and resolution.

Fluid mechanics is the branch of physics concerned with the mechanics of fluids and the forces on them. It has applications in a wide range of disciplines, including mechanical, aerospace, civil, chemical, and biomedical engineering, as well as geophysics, oceanography, meteorology, astrophysics, and biology.

<span class="mw-page-title-main">Fluid–structure interaction</span>

Fluid–structure interaction (FSI) is the interaction of some movable or deformable structure with an internal or surrounding fluid flow. Fluid–structure interactions can be stable or oscillatory. In oscillatory interactions, the strain induced in the solid structure causes it to move such that the source of strain is reduced, and the structure returns to its former state only for the process to repeat.

<span class="mw-page-title-main">Meshfree methods</span> Methods in numerical analysis not requiring knowledge of neighboring points

In the field of numerical analysis, meshfree methods are those that do not require connection between nodes of the simulation domain, i.e. a mesh, but are rather based on interaction of each node with all its neighbors. As a consequence, original extensive properties such as mass or kinetic energy are no longer assigned to mesh elements but rather to the single nodes. Meshfree methods enable the simulation of some otherwise difficult types of problems, at the cost of extra computing time and programming effort. The absence of a mesh allows Lagrangian simulations, in which the nodes can move according to the velocity field.

<span class="mw-page-title-main">Bram van Leer</span> Dutch mathematician

Bram van Leer is Arthur B. Modine Emeritus Professor of aerospace engineering at the University of Michigan, in Ann Arbor. He specializes in Computational fluid dynamics (CFD), fluid dynamics, and numerical analysis. His most influential work lies in CFD, a field he helped modernize from 1970 onwards. An appraisal of his early work has been given by C. Hirsch (1979)

Stokesian dynamics is a solution technique for the Langevin equation, which is the relevant form of Newton's 2nd law for a Brownian particle. The method treats the suspended particles in a discrete sense while the continuum approximation remains valid for the surrounding fluid, i.e., the suspended particles are generally assumed to be significantly larger than the molecules of the solvent. The particles then interact through hydrodynamic forces transmitted via the continuum fluid, and when the particle Reynolds number is small, these forces are determined through the linear Stokes equations. In addition, the method can also resolve non-hydrodynamic forces, such as Brownian forces, arising from the fluctuating motion of the fluid, and interparticle or external forces. Stokesian Dynamics can thus be applied to a variety of problems, including sedimentation, diffusion and rheology, and it aims to provide the same level of understanding for multiphase particulate systems as molecular dynamics does for statistical properties of matter. For rigid particles of radius suspended in an incompressible Newtonian fluid of viscosity and density , the motion of the fluid is governed by the Navier–Stokes equations, while the motion of the particles is described by the coupled equation of motion:

<span class="mw-page-title-main">Fluid animation</span> Computer graphics techniques for generating realistic animations of fluids

Fluid animation refers to computer graphics techniques for generating realistic animations of fluids such as water and smoke. Fluid animations are typically focused on emulating the qualitative visual behavior of a fluid, with less emphasis placed on rigorously correct physical results, although they often still rely on approximate solutions to the Euler equations or Navier–Stokes equations that govern real fluid physics. Fluid animation can be performed with different levels of complexity, ranging from time-consuming, high-quality animations for films, or visual effects, to simple and fast animations for real-time animations like computer games.

<span class="mw-page-title-main">Volume of fluid method</span> Free-surface modelling technique

In computational fluid dynamics, the volume of fluid (VOF) method is a family of free-surface modelling techniques, i.e. numerical techniques for tracking and locating the free surface. They belong to the class of Eulerian methods which are characterized by a mesh that is either stationary or is moving in a certain prescribed manner to accommodate the evolving shape of the interface. As such, VOF methods are advection schemes capturing the shape and position of the interface, but are not standalone flow solving algorithms. The Navier–Stokes equations describing the motion of the flow have to be solved separately.

In computational fluid dynamics, the immersed boundary method originally referred to an approach developed by Charles Peskin in 1972 to simulate fluid-structure (fiber) interactions. Treating the coupling of the structure deformations and the fluid flow poses a number of challenging problems for numerical simulations. In the immersed boundary method the fluid is represented in an Eulerian coordinate system and the structure is represented in Lagrangian coordinates. For Newtonian fluids governed by the Navier–Stokes equations, the fluid equations are

The moving particle semi-implicit (MPS) method is a computational method for the simulation of incompressible free surface flows. It is a macroscopic, deterministic particle method developed by Koshizuka and Oka (1996).

The multiphase particle-in-cell method (MP-PIC) is a numerical method for modeling particle-fluid and particle-particle interactions in a computational fluid dynamics (CFD) calculation. The MP-PIC method achieves greater stability than its particle-in-cell predecessor by simultaneously treating the solid particles as computational particles and as a continuum. In the MP-PIC approach, the particle properties are mapped from the Lagrangian coordinates to an Eulerian grid through the use of interpolation functions. After evaluation of the continuum derivative terms, the particle properties are mapped back to the individual particles. This method has proven to be stable in dense particle flows, computationally efficient, and physically accurate. This has allowed the MP-PIC method to be used as particle-flow solver for the simulation of industrial-scale chemical processes involving particle-fluid flows.

The viscous vortex domains (VVD) method is a mesh-free method of computational fluid dynamics for directly numerically solving 2D Navier-Stokes equations in Lagrange coordinates. It doesn't implement any turbulence model and free of arbitrary parameters. The main idea of this method is to present vorticity field with discrete regions (domains), which travel with diffusive velocity relatively to fluid and conserve their circulation. The same approach was used in Diffusion Velocity method of Ogami and Akamatsu, but VVD uses other discrete formulas

<span class="mw-page-title-main">Gerris (software)</span> Computer Software

Gerris is computer software in the field of computational fluid dynamics (CFD). Gerris was released as free and open-source software, subject to the requirements of the GNU General Public License (GPL), version 2 or any later.

Thomas Yizhao Hou is the Charles Lee Powell Professor of Applied and Computational Mathematics in the Department of Computing and Mathematical Sciences at the California Institute of Technology. He is known for his work in numerical analysis and mathematical analysis.

The finite point method (FPM) is a meshfree method for solving partial differential equations (PDEs) on scattered distributions of points. The FPM was proposed in the mid-nineties in, and with the purpose to facilitate the solution of problems involving complex geometries, free surfaces, moving boundaries and adaptive refinement. Since then, the FPM has evolved considerably, showing satisfactory accuracy and capabilities to deal with different fluid and solid mechanics problems.

<span class="mw-page-title-main">FEATool Multiphysics</span>

FEATool Multiphysics is a physics, finite element analysis (FEA), and partial differential equation (PDE) simulation toolbox. FEATool Multiphysics features the ability to model fully coupled heat transfer, fluid dynamics, chemical engineering, structural mechanics, fluid-structure interaction (FSI), electromagnetics, as well as user-defined and custom PDE problems in 1D, 2D (axisymmetry), or 3D, all within a graphical user interface (GUI) or optionally as script files. FEATool has been employed and used in academic research, teaching, and industrial engineering simulation contexts.

<span class="mw-page-title-main">Numerical modeling (geology)</span> Technique to solve geological problems by computational simulation

In geology, numerical modeling is a widely applied technique to tackle complex geological problems by computational simulation of geological scenarios.

References

  1. 1 2 T., Belytschko; Y., Krongauz; M., Flemming; D., Organ; S., Liu W. K. (1996), "Smoothing and accelerated computations in the element-free Galerkin method", Journal of Computational and Applied Mathematics , 74 (1–2): 111–126, doi:10.1016/0377-0427(96)00020-9
  2. ↑ P., Hansbo (1992), "The characteristic streamline diffusion method for the time-dependent incompressible Navier-Stokes equations", Computer Methods in Applied Mechanics and Engineering, 99 (2–3): 171–186, Bibcode:1992CMAME..99..171H, doi:10.1016/0045-7825(92)90039-M
  3. ↑ Harlow, Francis H.; E., Welch J. (1965), "Numerical study of large amplitude free surface motions", Physics of Fluids , 8 (5): 842–851, doi:10.1063/1.1761784
  4. ↑ J., Kelecy F.; H., Pletcher R. (1997), "The development of free surface capturing approach for multi dimensional free surface flows in closed containers", Journal of Computational Physics , 138 (2): 939–980, Bibcode:1997JCoPh.138..939K, doi:10.1006/jcph.1997.5847
  5. ↑ B., Kothe D.; C., Mjolsness R. (1992), "RIPPLE: A new model for incompressible flows with free surfaces", AIAA Journal , 30 (11): 2694–2700, Bibcode:1992AIAAJ..30.2694K, doi:10.2514/3.11286
  6. ↑ V., Maronnier; M., Picasso; J., Rappaz (1999), "Numerical simulation of free surface flows", Journal of Computational Physics , 155 (2): 439–455, Bibcode:1999JCoPh.155..439M, doi:10.1006/jcph.1999.6346
  7. 1 2 S., Tiwari; S., Manservisi (2000), Modeling incompressible Navier-Stokes flows by LSQ-SPH, Berichte des Fraunhofer ITWM, Kaiserslautern, Germany
  8. ↑ B., Lucy L. (1977), "A numerical approach to the testing of the fission hypothesis", Astronomical Journal, 82: 1013–1024, Bibcode:1977AJ.....82.1013L, doi:10.1086/112164
  9. ↑ A., Gingold R.; J., Monaghan J. (1977), "Smoothed particle hydrodynamics: theory and application to non-spherical stars", Monthly Notices of the Royal Astronomical Society , 181 (3): 375–389, doi: 10.1093/mnras/181.3.375
  10. 1 2 J., Monaghan J. (1992), "Smoothed particle hydrodynamics", Annual Review of Astronomy and Astrophysics , 30: 543–574, Bibcode:1992ARA&A..30..543M, doi:10.1146/annurev.aa.30.090192.002551
  11. ↑ J., Monaghan J.; A., Gingold R. (1983), "Shock Simulation by particle method SPH", Journal of Computational Physics , 52 (2): 374–389, Bibcode:1983JCoPh..52..374M, doi:10.1016/0021-9991(83)90036-0
  12. ↑ P., Morris J.; J., Fox P.; Y., Zhu (1997), "Modeling Low Reynolds Number Incompressible Flows Using SPH", Journal of Computational Physics , 136 (1): 214–226, Bibcode:1997JCoPh.136..214M, doi:10.1006/jcph.1997.5776
  13. ↑ J., Monaghan J. (1994), "Simulating free surface flows with SPH", Journal of Computational Physics , 110 (2): 399–406, Bibcode:1994JCoPh.110..399M, doi:10.1006/jcph.1994.1034
  14. ↑ A., Dilts G. (1996), Moving least squares particle hydrodynamics. I: consistency and stability, Hydrodynamics methods group report, Los Alamos National Laboratory, doi:10.1002/(SICI)1097-0207(19990320)44:8<1115::AID-NME547>3.0.CO;2-L
  15. 1 2 J., Kuhnert (1999), General smoothed particle hydrodynamics, Ph.D. thesis, Kaiserslautern University, Germany
  16. ↑ J., Kuhnert (2000), An upwind finite pointset method for compressible Euler and Navier-Stokes equations, preprint, ITWM, Kaiserslautern, Germany
  17. 1 2 S., Tiwari S.; J., Kuhnert J. (2001), Grid free method for solving Poisson equation, Berichte des Fraunhofer ITWM, Kaiserslautern, Germany, ISSN   1434-9973
  18. ↑ S., Tiwari; J., Kuhnert (2003), "Particle method for simulations of free surface flows", in Hou, Thomas Y.; Tadmor, Eitan (eds.), Hyperbolic Problems: Theory, Numerics, Applications: Proceedings of the Ninth International Conference on Hyperbolic Problems held in CalTech, Pasadena, March 25–29, 2002, doi:10.1007/978-3-642-55711-8, ISBN   978-3-642-55711-8
  19. ↑ A., Chorin A. (1968), "Numerical solution of the Navier-Stokes equations", Journal of Mathematics of Computation, 22 (104): 745–762, doi:10.2307/2004575, JSTOR   2004575
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.