Material point method

Last updated

The material point method (MPM) is a numerical technique used to simulate the behavior of solids, liquids, gases, and any other continuum material. Especially, it is a robust spatial discretization method for simulating multi-phase (solid-fluid-gas) interactions. In the MPM, a continuum body is described by a number of small Lagrangian elements referred to as 'material points'. These material points are surrounded by a background mesh/grid that is used to calculate terms such as the deformation gradient. Unlike other mesh-based methods like the finite element method, finite volume method or finite difference method, the MPM is not a mesh based method and is instead categorized as a meshless/meshfree or continuum-based particle method, examples of which are smoothed particle hydrodynamics and peridynamics. Despite the presence of a background mesh, the MPM does not encounter the drawbacks of mesh-based methods (high deformation tangling, advection errors etc.) which makes it a promising and powerful tool in computational mechanics.

Contents

The MPM was originally proposed, as an extension of a similar method known as FLIP (a further extension of a method called PIC) to computational solid dynamics, in the early 1990 by Professors Deborah L. Sulsky, Zhen Chen and Howard L. Schreyer at University of New Mexico. After this initial development, the MPM has been further developed both in the national labs as well as the University of New Mexico, Oregon State University, University of Utah and more across the US and the world. Recently the number of institutions researching the MPM has been growing with added popularity and awareness coming from various sources such as the MPM's use in the Disney film Frozen .

The algorithm

An MPM simulation consists of the following stages:

(Prior to the time integration phase)

  1. Initialization of grid and material points.
    1. A geometry is discretized into a collection of material points, each with its own material properties and initial conditions (velocity, stress, temperature, etc.)
    2. The grid, being only used to provide a place for gradient calculations is normally made to cover an area large enough to fill the expected extent of computational domain needed for the simulation.

(During the time integration phase - explicit formulation)

  1. Material point quantities are extrapolated to grid nodes.
    1. Material point mass (), momenta (), stresses (), and external forces () are extrapolated to the nodes at the corners of the cells within which the material points reside. This is most commonly done using standard linear shape functions (), the same used in FEM.
    2. The grid use the material point values to create the masses (), velocities (), internal and external force vectors (,) for the nodes:
  2. Equations of motion are solved on the grid.
    1. Newton's 2nd Law is solved to obtain the nodal acceleration ()
    2. New nodal velocities are found ().
  3. Derivative terms are extrapolated back to material points
    1. Material point acceleration (), deformation gradient () (or strain rate () depending on the strain theory used) is extrapolated from the surrounding nodes using similar shape functions to before ().
    2. Variables on the material points: positions, velocities, strains, stresses etc. are then updated with these rates depending on integration scheme of choice and a suitable constitutive model.
  4. Resetting of grid.
    Now that the material points are fully updated at the next time step, the grid is reset to allow for the next time step to begin.

History of PIC/MPM

The PIC was originally conceived to solve problems in fluid dynamics, and developed by Harlow at Los Alamos National Laboratory in 1957. [1] One of the first PIC codes was the Fluid-Implicit Particle (FLIP) program, which was created by Brackbill in 1986 [2] and has been constantly in development ever since. Until the 1990s, the PIC method was used principally in fluid dynamics.

Motivated by the need for better simulating penetration problems in solid dynamics, Sulsky, Chen and Schreyer started in 1993 to reformulate the PIC and develop the MPM, with funding from Sandia National Laboratories. [3] The original MPM was then further extended by Bardenhagen et al.. to include frictional contact, [4] which enabled the simulation of granular flow, [5] and by Nairn to include explicit cracks [6] and crack propagation (known as CRAMP).

Recently, an MPM implementation based on a micro-polar Cosserat continuum [7] has been used to simulate high-shear granular flow, such as silo discharge. MPM's uses were further extended into Geotechnical engineering with the recent development of a quasi-static, implicit MPM solver which provides numerically stable analyses of large-deformation problems in Soil mechanics. [8]

Annual workshops on the use of MPM are held at various locations in the United States. The Fifth MPM Workshop was held at Oregon State University, in Corvallis, OR, on April 2 and 3, 2009.

Applications of PIC/MPM

The uses of the PIC or MPM method can be divided into two broad categories: firstly, there are many applications involving fluid dynamics, plasma physics, magnetohydrodynamics, and multiphase applications. The second category of applications comprises problems in solid mechanics.

Fluid dynamics and multiphase simulations

The PIC method has been used to simulate a wide range of fluid-solid interactions, including sea ice dynamics, [9] penetration of biological soft tissues, [10] fragmentation of gas-filled canisters, [11] dispersion of atmospheric pollutants, [12] multiscale simulations coupling molecular dynamics with MPM, [13] [14] and fluid-membrane interactions. [15] In addition, the PIC-based FLIP code has been applied in magnetohydrodynamics and plasma processing tools, and simulations in astrophysics and free-surface flow. [16]

As a result of a joint effort between UCLA's mathematics department and Walt Disney Animation Studios, MPM was successfully used to simulate snow in the 2013 animated film Frozen . [17] [18] [19]

Solid mechanics

MPM has also been used extensively in solid mechanics, to simulate impact, penetration, collision and rebound, as well as crack propagation. [20] [21] MPM has also become a widely used method within the field of soil mechanics: it has been used to simulate granular flow, quickness test of sensitive clays, [22] landslides, [23] [24] [25] silo discharge, pile driving, fall-cone test, [26] [27] [28] [29] bucket filling, and material failure; and to model soil stress distribution, [30] compaction, and hardening. It is now being used in wood mechanics problems such as simulations of transverse compression on the cellular level including cell wall contact. [31] The work also received the George Marra Award for paper of the year from the Society of Wood Science and Technology. [32]

Classification of PIC/MPM codes

MPM in the context of numerical methods

One subset of numerical methods are Meshfree methods, which are defined as methods for which "a predefined mesh is not necessary, at least in field variable interpolation". Ideally, a meshfree method does not make use of a mesh "throughout the process of solving the problem governed by partial differential equations, on a given arbitrary domain, subject to all kinds of boundary conditions," although existing methods are not ideal and fail in at least one of these respects. Meshless methods, which are also sometimes called particle methods, share a "common feature that the history of state variables is traced at points (particles) which are not connected with any element mesh, the distortion of which is a source of numerical difficulties." As can be seen by these varying interpretations, some scientists consider MPM to be a meshless method, while others do not. All agree, however, that MPM is a particle method.

The Arbitrary Lagrangian Eulerian (ALE) methods form another subset of numerical methods which includes MPM. Purely Lagrangian methods employ a framework in which a space is discretised into initial subvolumes, whose flowpaths are then charted over time. Purely Eulerian methods, on the other hand, employ a framework in which the motion of material is described relative to a mesh that remains fixed in space throughout the calculation. As the name indicates, ALE methods combine Lagrangian and Eulerian frames of reference.

Subclassification of MPM/PIC

PIC methods may be based on either the strong form collocation or a weak form discretisation of the underlying partial differential equation (PDE). Those based on the strong form are properly referred to as finite-volume PIC methods. Those based on the weak form discretisation of PDEs may be called either PIC or MPM.

MPM solvers can model problems in one, two, or three spatial dimensions, and can also model axisymmetric problems. MPM can be implemented to solve either quasi-static or dynamic equations of motion, depending on the type of problem that is to be modeled. Several versions of MPM include Generalized Interpolation Material Point Method [33] ;Convected Particle Domain Interpolation Method; [34] Convected Particle Least Squares Interpolation Method. [35]

The time-integration used for MPM may be either explicit or implicit. The advantage to implicit integration is guaranteed stability, even for large timesteps. On the other hand, explicit integration runs much faster and is easier to implement.

Advantages

Compared to FEM

Unlike FEM, MPM does not require periodical remeshing steps and remapping of state variables, and is therefore better suited to the modeling of large material deformations. In MPM, particles and not the mesh points store all the information on the state of the calculation. Therefore, no numerical error results from the mesh returning to its original position after each calculation cycle, and no remeshing algorithm is required.

The particle basis of MPM allows it to treat crack propagation and other discontinuities better than FEM, which is known to impose the mesh orientation on crack propagation in a material. Also, particle methods are better at handling history-dependent constitutive models.

Compared to pure particle methods

Because in MPM nodes remain fixed on a regular grid, the calculation of gradients is trivial.

In simulations with two or more phases it is rather easy to detect contact between entities, as particles can interact via the grid with other particles in the same body, with other solid bodies, and with fluids.

Disadvantages of MPM

MPM is more expensive in terms of storage than other methods, as MPM makes use of mesh as well as particle data. MPM is more computationally expensive than FEM, as the grid must be reset at the end of each MPM calculation step and reinitialised at the beginning of the following step. Spurious oscillation may occur as particles cross the boundaries of the mesh in MPM, although this effect can be minimized by using generalized interpolation methods (GIMP). In MPM as in FEM, the size and orientation of the mesh can impact the results of a calculation: for example, in MPM, strain localisation is known to be particularly sensitive to mesh refinement. One stability problem in MPM that does not occur in FEM is the cell-crossing errors and null-space errors [36] because the number of integration points (material points) does not remain constant in a cell.

Notes

  1. Johnson, N. L. (1996). "The legacy and future of CFD at Los Alamos". Proceedings of the 1996 Canadian CFD Conference. OSTI   244662.
  2. Brackbill, J. U.; Ruppel, H. M. (1986). "FLIP: A method for adaptively zoned, particle-in-cell calculations of fluid flows in two dimensions". Journal of Computational Physics. 65 (2): 314–343. Bibcode:1986JCoPh..65..314B. doi:10.1016/0021-9991(86)90211-1. ISSN   0021-9991.
  3. Sulsky, D.; Chen, Z.; Schreyer, H. L. (1994). "A particle method for history-dependent materials". Computer Methods in Applied Mechanics and Engineering. 118 (1): 179–196. doi:10.1016/0045-7825(94)90112-0. ISSN   0045-7825.
  4. Bardenhagen, S. G.; Brackbill, J. U.; Sulsky, D. L. (1998). "Shear deformation in granular materials". doi:10.2172/329539. OSTI   329539.{{cite journal}}: Cite journal requires |journal= (help)
  5. Więckowski, Zdzisław; Youn, Sung-Kie; Yeon, Jeoung-Heum (1999). "A particle-in-cell solution to the silo discharging problem". International Journal for Numerical Methods in Engineering. 45 (9): 1203–1225. Bibcode:1999IJNME..45.1203W. doi:10.1002/(SICI)1097-0207(19990730)45:9<1203::AID-NME626>3.0.CO;2-C. ISSN   1097-0207.
  6. Nairn, J. A. (2003). "Material Point Method Calculations with Explicit Cracks". Computer Modeling in Engineering & Sciences. 4 (6): 649–664. doi:10.3970/cmes.2003.004.649.
  7. Coetzee, Corne J. (2004). The modelling of granular flow using the particle-in-cell method (PhD thesis). Stellenbosch : University of Stellenbosch.
  8. Beuth, L., Coetzee, C.J., Bonnier, P. and van den Berg, P. "Formulation and validation of a quasi-static material point method." In 10th International Symposium on Numerical Methods in Geomechanics, 2007.
  9. Wang, R.-X; Ji, S.-Y.; Shen, Hung Tao; Yue, Q.-J. (2005). "Modified PIC method for sea ice dynamics". China Ocean Engineering. 19: 457–468 via ResearchGate.
  10. Ionescu, I., Guilkey, J., Berzins, M., Kirby, R., and Weiss, J. "Computational simulation of penetrating trauma in biological soft tissues using MPM."
  11. Banerjee, Biswajit (2012). "Material point method simulations of fragmenting cylinders". ResearchGate. arXiv: 1201.2439 . Bibcode:2012arXiv1201.2439B . Retrieved 2019-06-18.
  12. Patankar, N. A.; Joseph, D. D. (2001). "Lagrangian numerical simulation of particulate flows". International Journal of Multiphase Flow. 27 (10): 1685–1706. doi:10.1016/S0301-9322(01)00025-8. ISSN   0301-9322.
  13. Lu, H.; Daphalapurkar, N. P.; Wang, B.; Roy, S.; Komanduri, R. (2006). "Multiscale simulation from atomistic to continuum – coupling molecular dynamics (MD) with the material point method (MPM)". Philosophical Magazine. 86 (20): 2971–2994. Bibcode:2006PMag...86.2971L. doi:10.1080/14786430600625578. ISSN   1478-6435. S2CID   137383632.
  14. Ma, Jin (2006). Multiscale Simulation Using the Generalized Interpolation Material Point Method, Discrete Dislocations and Molecular Dynamics (PhD thesis). Oklahoma State University.
  15. York, Allen R.; Sulsky, Deborah; Schreyer, Howard L. (2000). "Fluid–membrane interaction based on the material point method". International Journal for Numerical Methods in Engineering. 48 (6): 901–924. Bibcode:2000IJNME..48..901Y. doi:10.1002/(SICI)1097-0207(20000630)48:6<901::AID-NME910>3.0.CO;2-T. ISSN   1097-0207.
  16. Liu, Wing Kam; Li, Shaofan (2002). "Meshfree and particle methods and their applications". Applied Mechanics Reviews. 55 (1): 1–34. Bibcode:2002ApMRv..55....1L. doi:10.1115/1.1431547. ISSN   0003-6900. S2CID   17197495.
  17. Marquez, Letisia (February 27, 2014). "UCLA's mathematicians bring snow to life for Disney's "Frozen"". UCLA Today. Archived from the original on 10 March 2014. Retrieved 6 March 2014.
  18. Alexey Stomakhin; Craig Schroeder; Lawrence Chai; Joseph Teran; Andrew Selle (August 2013). "A material point method for snow simulation" (PDF). Walt Disney Animation Studios. Archived from the original (PDF) on 24 March 2014. Retrieved 6 March 2014.
  19. "Making of Disney's Frozen: A Material Point Method For Snow Simulation". CG Meetup. November 21, 2013. Retrieved 18 January 2014.
  20. Karuppiah, Venkatesh (2004). Implementation of irregular mesh in MPM for simulation of mixed mode crack opening in tension (Master's thesis). Oklahoma State University.
  21. Daphalapurkar, Nitin P.; Lu, Hongbing; Coker, Demir; Komanduri, Ranga (2007-01-01). "Simulation of dynamic crack growth using the generalized interpolation material point (GIMP) method". International Journal of Fracture. 143 (1): 79–102. doi:10.1007/s10704-007-9051-z. ISSN   1573-2673. S2CID   20013793.
  22. Tran, Quoc-Anh; Solowski, Wojciech; Thakur, Vikas; Karstunen, Minna (2017). "Modelling of the Quickness Test of Sensitive Clays Using the Generalized Interpolation Material Point Method". Landslides in Sensitive Clays. Advances in Natural and Technological Hazards Research. Vol. 46. pp. 323–326. doi:10.1007/978-3-319-56487-6_29. ISBN   978-3-319-56486-9.
  23. Tran, Quoc-Anh; Solowski, Wojciech (2019). "Generalized Interpolation Material Point Method modelling of large deformation problems including strain-rate effects – Application to penetration and progressive failure problems". Computers and Geotechnics. 106 (1): 249–265. Bibcode:2019CGeot.106..249T. doi: 10.1016/j.compgeo.2018.10.020 .
  24. Llano-Serna, Marcelo A.; Farias, Márcio M.; Pedroso, Dorival M. (2016). "An assessment of the material point method for modelling large scale run-out processes in landslides". Landslides. 13 (5): 1057–1066. doi:10.1007/s10346-015-0664-4. ISSN   1612-510X. S2CID   130645666.
  25. Llano Serna, Marcelo Alejandro; Muniz-de Farias, Márcio; Martínez-Carvajal, Hernán Eduardo (2015-12-21). "Numerical modelling of Alto Verde landslide using the material point method". DYNA. 82 (194): 150–159. doi: 10.15446/dyna.v82n194.48179 . ISSN   2346-2183.
  26. Tran, Quoc-Anh; Solowski, Wojciech (2019). "Generalized Interpolation Material Point Method modelling of large deformation problems including strain-rate effects – Application to penetration and progressive failure problems". Computers and Geotechnics. 106 (1): 249–265. Bibcode:2019CGeot.106..249T. doi: 10.1016/j.compgeo.2018.10.020 .
  27. Tran, Quoc-Anh; Solowski, Wojciech; Karstunen, Minna; Korkiala-Tanttua, Leena (2017). "Modelling of Fall-cone Tests with Strain-rate Effects". Procedia Engineering. 175: 293–301. doi: 10.1016/j.proeng.2017.01.029 .
  28. Llano-Serna, M.A.; Farias, M.M.; Pedroso, D.M.; Williams, David J.; Sheng, D. (2016). "Simulations of Fall Cone Test in Soil Mechanics Using the Material Point Method". Applied Mechanics and Materials. 846: 336–341. doi:10.4028/www.scientific.net/AMM.846.336. ISSN   1662-7482. S2CID   113653285.
  29. Llano-Serna, M; Farias, M (2014-06-03), Hicks, Michael; Brinkgreve, Ronald; Rohe, Alexander (eds.), "Use of generalized material point method (GIMP) to simulate shallow wedge penetration", Numerical Methods in Geotechnical Engineering, CRC Press, pp. 259–264, doi:10.1201/b17017-48, ISBN   9781138001466
  30. Llano-Serna, M.A.; Farias, M.M. (2016). "Validación numérica, teórica y experimental del método del punto material para resolver problemas geotécnicos". Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería (in Spanish). 32 (2): 110–115. doi: 10.1016/j.rimni.2015.02.008 . hdl: 2117/167257 .
  31. Nairn, John A. (2007). "Numerical Simulations of Transverse Compression and Densification in Wood". Wood and Fiber Science. 38 (4): 576–591. ISSN   0735-6161.
  32. "Society of Wood Science and Technology: George Marra Award Recipients". 2007. Archived from the original on 2007-09-23. Retrieved 2019-06-18.
  33. Bardenhagen, S. G.; Kober, E. M. (2004). "The Generalized Interpolation Material Point Method". Computer Modeling in Engineering & Sciences. 5: 477–496. doi:10.3970/cmes.2004.005.477.
  34. Sadeghirad, A.; Brannon, R. M.; Burghardt, J. (2011). "A convected particle domain interpolation technique to extend applicability of the material point method for problems involving massive deformations". International Journals for Numerical Methods in Engineering. 86 (12): 1435–1456. Bibcode:2011IJNME..86.1435S. doi:10.1002/nme.3110. S2CID   16715144.
  35. Tran, Quoc-Anh; Solowski, Wojciech; Berzins, Martin; Gulkey, James (2020). "A convected particle least square interpolation material point method". International Journals for Numerical Methods in Engineering. 121 (6): 1068–1100. Bibcode:2020IJNME.121.1068T. doi:10.1002/nme.6257. S2CID   209961739.
  36. Tran, Quoc-Anh; Solowski, Wojciech (2017). "Temporal and null-space filter for the material point method". International Journal for Numerical Methods in Engineering. 120 (3): 328–360. doi: 10.1002/nme.6138 .

Related Research Articles

<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.

<span class="mw-page-title-main">Crash simulation</span> Virtual recreation of a destructive car crash

A crash simulation is a virtual recreation of a destructive crash test of a car or a highway guard rail system using a computer simulation in order to examine the level of safety of the car and its occupants. Crash simulations are used by automakers during computer-aided engineering (CAE) analysis for crashworthiness in the computer-aided design (CAD) process of modelling new cars. During a crash simulation, the kinetic energy, or energy of motion, that a vehicle has before the impact is transformed into deformation energy, mostly by plastic deformation (plasticity) of the car body material, at the end of the impact.

<i>N</i>-body simulation Simulation of a dynamical system of particles

In physics and astronomy, an N-body simulation is a simulation of a dynamical system of particles, usually under the influence of physical forces, such as gravity. N-body simulations are widely used tools in astrophysics, from investigating the dynamics of few-body systems like the Earth-Moon-Sun system to understanding the evolution of the large-scale structure of the universe. In physical cosmology, N-body simulations are used to study processes of non-linear structure formation such as galaxy filaments and galaxy halos from the influence of dark matter. Direct N-body simulations are used to study the dynamical evolution of star clusters.

<span class="mw-page-title-main">Lattice Boltzmann methods</span> Class of computational fluid dynamics methods

The lattice Boltzmann methods (LBM), originated from the lattice gas automata (LGA) method (Hardy-Pomeau-Pazzis and Frisch-Hasslacher-Pomeau models), is a class of computational fluid dynamics (CFD) methods for fluid simulation. Instead of solving the Navier–Stokes equations directly, a fluid density on a lattice is simulated with streaming and collision (relaxation) processes. The method is versatile as the model fluid can straightforwardly be made to mimic common fluid behaviour like vapour/liquid coexistence, and so fluid systems such as liquid droplets can be simulated. Also, fluids in complex environments such as porous media can be straightforwardly simulated, whereas with complex boundaries other CFD methods can be hard to work with.

<span class="mw-page-title-main">Soft-body dynamics</span> Computer graphics simulation of deformable objects

Soft-body dynamics is a field of computer graphics that focuses on visually realistic physical simulations of the motion and properties of deformable objects. The applications are mostly in video games and films. Unlike in simulation of rigid bodies, the shape of soft bodies can change, meaning that the relative distance of two points on the object is not fixed. While the relative distances of points are not fixed, the body is expected to retain its shape to some degree. The scope of soft body dynamics is quite broad, including simulation of soft organic materials such as muscle, fat, hair and vegetation, as well as other deformable materials such as clothing and fabric. Generally, these methods only provide visually plausible emulations rather than accurate scientific/engineering simulations, though there is some crossover with scientific methods, particularly in the case of finite element simulations. Several physics engines currently provide software for soft-body simulation.

<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.

The fast multipole method (FMM) is a numerical technique that was developed to speed up the calculation of long-ranged forces in the n-body problem. It does this by expanding the system Green's function using a multipole expansion, which allows one to group sources that lie close together and treat them as if they are a single source.

<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 fracture mechanics, the energy release rate, , is the rate at which energy is transformed as a material undergoes fracture. Mathematically, the energy release rate is expressed as the decrease in total potential energy per increase in fracture surface area, and is thus expressed in terms of energy per unit area. Various energy balances can be constructed relating the energy released during fracture to the energy of the resulting new surface, as well as other dissipative processes such as plasticity and heat generation. The energy release rate is central to the field of fracture mechanics when solving problems and estimating material properties related to fracture and fatigue.

<span class="mw-page-title-main">Finite element method</span> Numerical method for solving physical or engineering problems

The finite element method (FEM) is a popular method for numerically solving differential equations arising in engineering and mathematical modeling. Typical problem areas of interest include the traditional fields of structural analysis, heat transfer, fluid flow, mass transport, and electromagnetic potential.

In numerical analysis, transfinite interpolation is a means to construct functions over a planar domain in such a way that they match a given function on the boundary. This method is applied in geometric modelling and in the field of finite element method.

In fluid mechanics, meteorology (weather) and oceanography, a trajectory traces the motion of a single point, often called a parcel, in the flow.

Biology Monte Carlo methods (BioMOCA) have been developed at the University of Illinois at Urbana-Champaign to simulate ion transport in an electrolyte environment through ion channels or nano-pores embedded in membranes. It is a 3-D particle-based Monte Carlo simulator for analyzing and studying the ion transport problem in ion channel systems or similar nanopores in wet/biological environments. The system simulated consists of a protein forming an ion channel (or an artificial nanopores like a Carbon Nano Tube, CNT), with a membrane (i.e. lipid bilayer) that separates two ion baths on either side. BioMOCA is based on two methodologies, namely the Boltzmann transport Monte Carlo (BTMC) and particle-particle-particle-mesh (P3M). The first one uses Monte Carlo method to solve the Boltzmann equation, while the later splits the electrostatic forces into short-range and long-range components.

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.

Smoothed finite element methods (S-FEM) are a particular class of numerical simulation algorithms for the simulation of physical phenomena. It was developed by combining meshfree methods with the finite element method. S-FEM are applicable to solid mechanics as well as fluid dynamics problems, although so far they have mainly been applied to the former.

The Kansa method is a computer method used to solve partial differential equations. Its main advantage is it is very easy to understand and program on a computer. It is much less complicated than the finite element method. Another advantage is it works well on multi variable problems. The finite element method is complicated when working with more than 3 space variables and time.

<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.

OceanParcels, “Probably A Really Computationally Efficient Lagrangian Simulator”, is a set of Python classes and methods that is used to track particles like water, plankton and plastics. It uses the output of ocean general circulation model (OGCMs). OceanParcels main goal is to process the increasingly large amounts of data that is governed by OGCM's. The flow dynamics are simulated using Lagrangian modelling and the geophysical fluid dynamics are simulated with Eulerian modelling or provided through experimental data. OceanParcels is dependent on two principles, namely the ability to read external data sets from different formats and customizable kernels to define particle dynamics.

<span class="mw-page-title-main">Momentum mapping format</span> Technique in the Material Point Method (MPM)

Momentum mapping format is a key technique in the Material Point Method (MPM) for transferring physical quantities such as momentum, mass, and stress between a material point and a background grid.