FLIP: A method for adaptively zoned, particle-in-cell calculations of fluid flows in two dimensions

Brackbill, J. U.; Ruppel, H.M.

doi:10.1016/0021-9991(86)90211-1

Cited by 456 publications

(268 citation statements)

References 10 publications

Supporting

Mentioning

268

Contrasting

Order By: Relevance

“…FLIPMHD3D: 2D, Lagrangian, PIC MHD FLIPMHD3D is a computational tool for PIC visco-resistive MHD based on the FLIP-MHD method developed by Brackbill (see Ruppel 1986 andBrackbill 1991). The code is designed to solve the resistive MHD equations (here in Gaussian units),…”

Section: Modeling Approachmentioning

confidence: 99%

Plasma Physical Parameters Along Cme-Driven Shocks. Ii. Observation–simulation Comparison

Bacchini

Susino

Бемпорад

et al. 2015

ApJ

View full text Add to dashboard Cite

In this work, we compare the spatial distribution of the plasma parameters along the 1999 June 11 coronal mass ejection (CME)-driven shock front with the results obtained from a CME-like event simulated with the FLIPMHD3D code, based on the FLIP-MHD particle-in-cell method. The observational data are retrieved from the combination of white-light coronagraphic data (for the upstream values) and the application of the RankineHugoniot equations (for the downstream values). The comparison shows a higher compression ratio X and Alfvénic Mach number M A at the shock nose, and a stronger magnetic field deflection d toward the flanks, in agreement with observations. Then, we compare the spatial distribution of M A with the profiles obtained from the solutions of the shock adiabatic equation relating M A , X, and Bn q (the angle between the upstream magnetic field and the shock front normal) for the special cases of parallel and perpendicular shock, and with a semi-empirical expression for a generically oblique shock. The semi-empirical curve approximates the actual values of M A very well, if the effects of a non-negligible shock thickness sh d and plasma-to magnetic pressure ratio u b are taken into account throughout the computation. Moreover, the simulated shock turns out to be supercritical at the nose and sub-critical at the flanks. Finally, we develop a new one-dimensional Lagrangian ideal MHD method based on the GrAALE code, to simulate the ion-electron temperature decoupling due to the shock transit. Two models are used, a simple solar wind model and a variable-γ model. Both produce results in agreement with observations, the second one being capable of introducing the physics responsible for the additional electron heating due to secondary effects (collisions, Alfvén waves, etc.).

show abstract

Section: Modeling Approachmentioning

confidence: 99%

Plasma Physical Parameters Along Cme-Driven Shocks. Ii. Observation–simulation Comparison

Bacchini

Susino

Бемпорад

et al. 2015

ApJ

View full text Add to dashboard Cite

show abstract

“…The idea of combining fixed meshes with moving particles is not new. It was introduced for convection diffusion problems in [15] and was used in the so-called Particle in Cell method (PIC) [16] and later in its extension called the Material Point Method (MPM) [17]. All these methods use a Finite Element (FE) background mesh.…”

Section: Introductionmentioning

confidence: 99%

Lagrangian versus Eulerian integration errors

Idelsohn

Oñate

Nigro

et al. 2015

Computer Methods in Applied Mechanics and Engineering

View full text Add to dashboard Cite

show abstract

“…The material point method (MPM) was described by Sulsky et al [5,6] as an extension to the FLIP (FluidImplicit Particle) method of Brackbill [7], which itself is an extension of the particle-in-cell (PIC) method of Harlow [8].…”

Section: Materials Point Methodsmentioning

confidence: 99%

Mesoscale numerical modeling of plastic bonded explosives under shock loading

Shang

Zhao

et al. 2015

EPJ Web of Conferences

View full text Add to dashboard Cite

Abstract. Mesoscale responses of plastic bonded explosives under shock loading are investigated using material point method as implemented in the Uintah Computational Framework. The two-dimensional geometrical model which can approximately reflect the mesoscopic structure of plastic bonded explosives was created based on the Voronoi tessellation. Shock loading for the explosive was performed by a piston moving at a constant velocity. For the purpose of investigating the influence of shock strength on the responses of explosives, two different velocities for the piston were used, 200 m/s and 400 m/s, respectively. The simulation results indicate that under shock loading there forms some stress localizations on the grain boundary of explosive. These stress localizations lead to large plastic deformations, and the plastic strain energy transforms to thermal energy immediately, causing temperature to rise rapidly and form some hot spots on grain boundary areas. The comparison between two different piston velocities shows that with increasing shock strength, the distribution of plastic strain and temperature does not have significant change, but their values increase obviously. Namely, the higher the shock strength is, the higher the hot spot temperature will be.

show abstract

FLIP: A method for adaptively zoned, particle-in-cell calculations of fluid flows in two dimensions

Cited by 456 publications

References 10 publications

Plasma Physical Parameters Along Cme-Driven Shocks. Ii. Observation–simulation Comparison

Plasma Physical Parameters Along Cme-Driven Shocks. Ii. Observation–simulation Comparison

Lagrangian versus Eulerian integration errors

Mesoscale numerical modeling of plastic bonded explosives under shock loading

Contact Info

Product

Resources

About