A local adaptive discretization algorithm for Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Spreng, Fabian; Schnabel, Dirk; Mueller, Alexandra; Eberhard, Peter
2014-06-01
In this paper, an extension to the Smoothed Particle Hydrodynamics (SPH) method is proposed that allows for an adaptation of the discretization level of a simulated continuum at runtime. By combining a local adaptive refinement technique with a newly developed coarsening algorithm, one is able to improve the accuracy of the simulation results while reducing the required computational cost at the same time. For this purpose, the number of particles is, on the one hand, adaptively increased in critical areas of a simulation model. Typically, these are areas that show a relatively low particle density and high gradients in stress or temperature. On the other hand, the number of SPH particles is decreased for domains with a high particle density and low gradients. Besides a brief introduction to the basic principle of the SPH discretization method, the extensions to the original formulation providing such a local adaptive refinement and coarsening of the modeled structure are presented in this paper. After having introduced its theoretical background, the applicability of the enhanced formulation, as well as the benefit gained from the adaptive model discretization, is demonstrated in the context of four different simulation scenarios focusing on solid continua. While presenting the results found for these examples, several properties of the proposed adaptive technique are discussed, e.g. the conservation of momentum as well as the existing correlation between the chosen refinement and coarsening patterns and the observed quality of the results.
Adaptive particle refinement and derefinement applied to the smoothed particle hydrodynamics method
NASA Astrophysics Data System (ADS)
Barcarolo, D. A.; Le Touzé, D.; Oger, G.; de Vuyst, F.
2014-09-01
SPH simulations are usually performed with a uniform particle distribution. New techniques have been recently proposed to enable the use of spatially varying particle distributions, which encouraged the development of automatic adaptivity and particle refinement/derefinement algorithms. All these efforts resulted in very interesting and promising procedures leading to more efficient and faster SPH simulations. In this article, a family of particle refinement techniques is reviewed and a new derefinement technique is proposed and validated through several test cases involving both free-surface and viscous flows. Besides, this new procedure allows higher resolutions in the regions requiring increased accuracy. Moreover, several levels of refinement can be used with this new technique, as often encountered in adaptive mesh refinement techniques in mesh-based methods.
NASA Astrophysics Data System (ADS)
García-Senz, Domingo; Cabezón, Rubén M.; Escartín, José A.; Ebinger, Kevin
2014-10-01
Context. The smoothed-particle hydrodynamics (SPH) technique is a numerical method for solving gas-dynamical problems. It has been applied to simulate the evolution of a wide variety of astrophysical systems. The method has a second-order accuracy, with a resolution that is usually much higher in the compressed regions than in the diluted zones of the fluid. Aims: We propose and check a method to balance and equalize the resolution of SPH between high- and low-density regions. This method relies on the versatility of a family of interpolators called sinc kernels, which allows increasing the interpolation quality by varying only a single parameter (the exponent of the sinc function). Methods: The proposed method was checked and validated through a number of numerical tests, from standard one-dimensional Riemann problems in shock tubes, to multidimensional simulations of explosions, hydrodynamic instabilities, and the collapse of a Sun-like polytrope. Results: The analysis of the hydrodynamical simulations suggests that the scheme devised to equalize the accuracy improves the treatment of the post-shock regions and, in general, of the rarefacted zones of fluids while causing no harm to the growth of hydrodynamic instabilities. The method is robust and easy to implement with a low computational overload. It conserves mass, energy, and momentum and reduces to the standard SPH scheme in regions of the fluid that have smooth density gradients.
NASA Astrophysics Data System (ADS)
Power, C.; Read, J. I.; Hobbs, A.
2014-06-01
We simulate cosmological galaxy cluster formation using three different approaches to solving the equations of non-radiative hydrodynamics - classic smoothed particle hydrodynamics (SPH), novel SPH with a higher order dissipation switch (SPHS), and an adaptive mesh refinement (AMR) method. Comparing spherically averaged entropy profiles, we find that SPHS and AMR approaches result in a well-defined entropy core that converges rapidly with increasing mass and force resolution. In contrast, the central entropy profile in the SPH approach is sensitive to the cluster's assembly history and shows poor numerical convergence. We trace this disagreement to the known artificial surface tension in SPH that appears at phase boundaries. Varying systematically numerical dissipation in SPHS, we study the contributions of numerical and physical dissipation to the entropy core and argue that numerical dissipation is required to ensure single-valued fluid quantities in converging flows. However, provided it occurs only at the resolution limit and does not propagate errors to larger scales, its effect is benign - there is no requirement to build `sub-grid' models of unresolved turbulence for galaxy cluster simulations. We conclude that entropy cores in non-radiative galaxy cluster simulations are physical, resulting from entropy generation in shocked gas during cluster assembly.
NASA Astrophysics Data System (ADS)
Rasia, Elena; Lau, Erwin T.; Borgani, Stefano; Nagai, Daisuke; Dolag, Klaus; Avestruz, Camille; Granato, Gian Luigi; Mazzotta, Pasquale; Murante, Giuseppe; Nelson, Kaylea; Ragone-Figueroa, Cinthia
2014-08-01
Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10%-30%. The largest bias originates from both violation of hydrostatic equilibrium (HE) and an additional temperature bias caused by inhomogeneities in the X-ray-emitting intracluster medium (ICM). To elucidate this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle hydrodynamics (SPH) or adaptive-mesh refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in nonradiative simulations, whereas it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with that observed in the central regions of clusters. Around R 500 the temperature inhomogeneities of the SPH simulations can generate twice the typical HE mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high-sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii.
Rasia, Elena; Lau, Erwin T.; Nagai, Daisuke; Avestruz, Camille; Borgani, Stefano; Dolag, Klaus; Granato, Gian Luigi; Murante, Giuseppe; Ragone-Figueroa, Cinthia; Mazzotta, Pasquale; Nelson, Kaylea
2014-08-20
Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10%-30%. The largest bias originates from both violation of hydrostatic equilibrium (HE) and an additional temperature bias caused by inhomogeneities in the X-ray-emitting intracluster medium (ICM). To elucidate this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle hydrodynamics (SPH) or adaptive-mesh refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in nonradiative simulations, whereas it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with that observed in the central regions of clusters. Around R {sub 500} the temperature inhomogeneities of the SPH simulations can generate twice the typical HE mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high-sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii.
Nonequilibrium flows with smooth particle applied mechanics
Kum, O.
1995-07-01
Smooth particle methods are relatively new methods for simulating solid and fluid flows through they have a 20-year history of solving complex hydrodynamic problems in astrophysics, such as colliding planets and stars, for which correct answers are unknown. The results presented in this thesis evaluate the adaptability or fitness of the method for typical hydrocode production problems. For finite hydrodynamic systems, boundary conditions are important. A reflective boundary condition with image particles is a good way to prevent a density anomaly at the boundary and to keep the fluxes continuous there. Boundary values of temperature and velocity can be separately controlled. The gradient algorithm, based on differentiating the smooth particle expression for (u{rho}) and (T{rho}), does not show numerical instabilities for the stress tensor and heat flux vector quantities which require second derivatives in space when Fourier`s heat-flow law and Newton`s viscous force law are used. Smooth particle methods show an interesting parallel linking to them to molecular dynamics. For the inviscid Euler equation, with an isentropic ideal gas equation of state, the smooth particle algorithm generates trajectories isomorphic to those generated by molecular dynamics. The shear moduli were evaluated based on molecular dynamics calculations for the three weighting functions, B spline, Lucy, and Cusp functions. The accuracy and applicability of the methods were estimated by comparing a set of smooth particle Rayleigh-Benard problems, all in the laminar regime, to corresponding highly-accurate grid-based numerical solutions of continuum equations. Both transient and stationary smooth particle solutions reproduce the grid-based data with velocity errors on the order of 5%. The smooth particle method still provides robust solutions at high Rayleigh number where grid-based methods fails.
Some cautionary remarks about smoothed particle hydrodynamics
NASA Technical Reports Server (NTRS)
Hernquist, Lars
1993-01-01
Potential difficulties with smoothed particle hydrodynamics are discussed. In particular, empirical tests are used to demonstrate that the errors resulting from the use of variable smoothing can be much larger than commonly believed. Fortunately, however, these errors, which are normally small, do not appear to promote instability on small scales, such as fragmentation in self-gravitating fluids. Still, while SPH remains a useful tool for many problems of astrophysical interest, a rigorous formulation of it, which is adaptive but still satisfies conservation properties, is clearly wanting.
Progress in smooth particle hydrodynamics
Wingate, C.A.; Dilts, G.A.; Mandell, D.A.; Crotzer, L.A.; Knapp, C.E.
1998-07-01
Smooth Particle Hydrodynamics (SPH) is a meshless, Lagrangian numerical method for hydrodynamics calculations where calculational elements are fuzzy particles which move according to the hydrodynamic equations of motion. Each particle carries local values of density, temperature, pressure and other hydrodynamic parameters. A major advantage of SPH is that it is meshless, thus large deformation calculations can be easily done with no connectivity complications. Interface positions are known and there are no problems with advecting quantities through a mesh that typical Eulerian codes have. These underlying SPH features make fracture physics easy and natural and in fact, much of the applications work revolves around simulating fracture. Debris particles from impacts can be easily transported across large voids with SPH. While SPH has considerable promise, there are some problems inherent in the technique that have so far limited its usefulness. The most serious problem is the well known instability in tension leading to particle clumping and numerical fracture. Another problem is that the SPH interpolation is only correct when particles are uniformly spaced a half particle apart leading to incorrect strain rates, accelerations and other quantities for general particle distributions. SPH calculations are also sensitive to particle locations. The standard artificial viscosity treatment in SPH leads to spurious viscosity in shear flows. This paper will demonstrate solutions for these problems that they and others have been developing. The most promising is to replace the SPH interpolant with the moving least squares (MLS) interpolant invented by Lancaster and Salkauskas in 1981. SPH and MLS are closely related with MLS being essentially SPH with corrected particle volumes. When formulated correctly, JLS is conservative, stable in both compression and tension, does not have the SPH boundary problems and is not sensitive to particle placement. The other approach to
Generating Optimal Initial Conditions for Smoothed Particle Hydrodynamics Simulations
NASA Astrophysics Data System (ADS)
Diehl, S.; Rockefeller, G.; Fryer, C. L.; Riethmiller, D.; Statler, T. S.
2015-12-01
We review existing smoothed particle hydrodynamics setup methods and outline their advantages, limitations, and drawbacks. We present a new method for constructing initial conditions for smoothed particle hydrodynamics simulations, which may also be of interest for N-body simulations, and demonstrate this method on a number of applications. This new method is inspired by adaptive binning techniques using weighted Voronoi tessellations. Particles are placed and iteratively moved based on their proximity to neighbouring particles and the desired spatial resolution. This new method can satisfy arbitrarily complex spatial resolution requirements.
Numerical Convergence In Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Zhu, Qirong; Hernquist, Lars; Li, Yuexing
2015-02-01
We study the convergence properties of smoothed particle hydrodynamics (SPH) using numerical tests and simple analytic considerations. Our analysis shows that formal numerical convergence is possible in SPH only in the joint limit N → ∞, h → 0, and Nnb → ∞, where N is the total number of particles, h is the smoothing length, and Nnb is the number of neighbor particles within the smoothing volume used to compute smoothed estimates. Previous work has generally assumed that the conditions N → ∞ and h → 0 are sufficient to achieve convergence, while holding Nnb fixed. We demonstrate that if Nnb is held fixed as the resolution is increased, there will be a residual source of error that does not vanish as N → ∞ and h → 0. Formal numerical convergence in SPH is possible only if Nnb is increased systematically as the resolution is improved. Using analytic arguments, we derive an optimal compromise scaling for Nnb by requiring that this source of error balance that present in the smoothing procedure. For typical choices of the smoothing kernel, we find Nnb vpropN 0.5. This means that if SPH is to be used as a numerically convergent method, the required computational cost does not scale with particle number as O(N), but rather as O(N 1 + δ), where δ ≈ 0.5, with a weak dependence on the form of the smoothing kernel.
Smoothed particle hydrodynamics with smoothed pseudo-density
NASA Astrophysics Data System (ADS)
Yamamoto, Satoko; Saitoh, Takayuki R.; Makino, Junichiro
2015-06-01
In this paper, we present a new formulation of smoothed particle hydrodynamics (SPH), which, unlike the standard SPH (SSPH), is well behaved at the contact discontinuity. The SSPH scheme cannot handle discontinuities in density (e.g., the contact discontinuity and the free surface), because it requires that the density of fluid is positive and continuous everywhere. Thus there is inconsistency in the formulation of the SSPH scheme at discontinuities of the fluid density. To solve this problem, we introduce a new quantity associated with particles and the "density" of that quantity. This "density" evolves through the usual continuity equation with an additional artificial diffusion term, in order to guarantee the continuity of the "density." We use this "density," or pseudo-density, instead of the mass density, to formulate our SPH scheme. We call our new method SPH with smoothed pseudo-density, and we show that it is physically consistent and can handle discontinuities quite well.
Conservative Smoothing on an Adaptive Quadrilateral Grid
NASA Astrophysics Data System (ADS)
Sun, M.; Takayama, K.
1999-03-01
The Lax-Wendroff scheme can be freed of spurious oscillations by introducing conservative smoothing. In this paper the approach is first tested in 1-D modeling equations and then extended to multidimensional flows by the finite volume method. The scheme is discretized by a space-splitting method on an adaptive quadrilateral grid. The artificial viscosity coefficients in the conservative smoothing step are specially designed to capture slipstreams and vortices. Algorithms are programmed using a vectorizable data structure, under which not only the flow solver but also the adaptation procedure is well vectorized. The good resolution and high efficiency of the approach are demonstrated in calculating both unsteady and steady compressible flows with either weak or strong shock waves.
NUMERICAL CONVERGENCE IN SMOOTHED PARTICLE HYDRODYNAMICS
Zhu, Qirong; Li, Yuexing; Hernquist, Lars
2015-02-10
We study the convergence properties of smoothed particle hydrodynamics (SPH) using numerical tests and simple analytic considerations. Our analysis shows that formal numerical convergence is possible in SPH only in the joint limit N → ∞, h → 0, and N{sub nb} → ∞, where N is the total number of particles, h is the smoothing length, and N{sub nb} is the number of neighbor particles within the smoothing volume used to compute smoothed estimates. Previous work has generally assumed that the conditions N → ∞ and h → 0 are sufficient to achieve convergence, while holding N{sub nb} fixed. We demonstrate that if N{sub nb} is held fixed as the resolution is increased, there will be a residual source of error that does not vanish as N → ∞ and h → 0. Formal numerical convergence in SPH is possible only if N{sub nb} is increased systematically as the resolution is improved. Using analytic arguments, we derive an optimal compromise scaling for N{sub nb} by requiring that this source of error balance that present in the smoothing procedure. For typical choices of the smoothing kernel, we find N{sub nb} ∝N {sup 0.5}. This means that if SPH is to be used as a numerically convergent method, the required computational cost does not scale with particle number as O(N), but rather as O(N {sup 1} {sup +} {sup δ}), where δ ≈ 0.5, with a weak dependence on the form of the smoothing kernel.
SPHGR: Smoothed-Particle Hydrodynamics Galaxy Reduction
NASA Astrophysics Data System (ADS)
Thompson, Robert
2015-02-01
SPHGR (Smoothed-Particle Hydrodynamics Galaxy Reduction) is a python based open-source framework for analyzing smoothed-particle hydrodynamic simulations. Its basic form can run a baryonic group finder to identify galaxies and a halo finder to identify dark matter halos; it can also assign said galaxies to their respective halos, calculate halo & galaxy global properties, and iterate through previous time steps to identify the most-massive progenitors of each halo and galaxy. Data about each individual halo and galaxy is collated and easy to access. SPHGR supports a wide range of simulations types including N-body, full cosmological volumes, and zoom-in runs. Support for multiple SPH code outputs is provided by pyGadgetReader (ascl:1411.001), mainly Gadget (ascl:0003.001) and TIPSY (ascl:1111.015).
Particle splitting in smoothed particle hydrodynamics based on Voronoi diagram
NASA Astrophysics Data System (ADS)
Chiaki, Gen; Yoshida, Naoki
2015-08-01
We present a novel method for particle splitting in smoothed particle hydrodynamics simulations. Our method utilizes the Voronoi diagram for a given particle set to determine the position of fine daughter particles. We perform several test simulations to compare our method with a conventional splitting method in which the daughter particles are placed isotropically over the local smoothing length. We show that, with our method, the density deviation after splitting is reduced by a factor of about 2 compared with the conventional method. Splitting would smooth out the anisotropic density structure if the daughters are distributed isotropically, but our scheme allows the daughter particles to trace the original density distribution with length-scales of the mean separation of their parent. We apply the particle splitting to simulations of the primordial gas cloud collapse. The thermal evolution is accurately followed to the hydrogen number density of 1012 cm-3. With the effective mass resolution of ˜10-4 M⊙ after the multistep particle splitting, the protostellar disc structure is well resolved. We conclude that the method offers an efficient way to simulate the evolution of an interstellar gas and the formation of stars.
Smoothed particle hydrodynamics with GRAPE-1A
NASA Technical Reports Server (NTRS)
Umemura, Masayuki; Fukushige, Toshiyuki; Makino, Junichiro; Ebisuzaki, Toshikazu; Sugimoto, Daiichiro; Turner, Edwin L.; Loeb, Abraham
1993-01-01
We describe the implementation of a smoothed particle hydrodynamics (SPH) scheme using GRAPE-1A, a special-purpose processor used for gravitational N-body simulations. The GRAPE-1A calculates the gravitational force exerted on a particle from all other particles in a system, while simultaneously making a list of the nearest neighbors of the particle. It is found that GRAPE-1A accelerates SPH calculations by direct summation by about two orders of magnitudes for a ten thousand-particle simulation. The effective speed is 80 Mflops, which is about 30 percent of the peak speed of GRAPE-1A. Also, in order to investigate the accuracy of GRAPE-SPH, some test simulations were executed. We found that the force and position errors are smaller than those due to representing a fluid by a finite number of particles. The total energy and momentum were conserved within 0.2-0.4 percent and 2-5 x 10 exp -5, respectively, in simulations with several thousand particles. We conclude that GRAPE-SPH is quite effective and sufficiently accurate for self-gravitating hydrodynamics.
Impact modeling with Smooth Particle Hydrodynamics
Stellingwerf, R.F.; Wingate, C.A.
1993-07-01
Smooth Particle Hydrodynamics (SPH) can be used to model hypervelocity impact phenomena via the addition of a strength of materials treatment. SPH is the only technique that can model such problems efficiently due to the combination of 3-dimensional geometry, large translations of material, large deformations, and large void fractions for most problems of interest. This makes SPH an ideal candidate for modeling of asteroid impact, spacecraft shield modeling, and planetary accretion. In this paper we describe the derivation of the strength equations in SPH, show several basic code tests, and present several impact test cases with experimental comparisons.
Workshop on advances in smooth particle hydrodynamics
Wingate, C.A.; Miller, W.A.
1993-12-31
This proceedings contains viewgraphs presented at the 1993 workshop held at Los Alamos National Laboratory. Discussed topics include: negative stress, reactive flow calculations, interface problems, boundaries and interfaces, energy conservation in viscous flows, linked penetration calculations, stability and consistency of the SPH method, instabilities, wall heating and conservative smoothing, tensors, tidal disruption of stars, breaking the 10,000,000 particle limit, modelling relativistic collapse, SPH without H, relativistic KSPH avoidance of velocity based kernels, tidal compression and disruption of stars near a supermassive rotation black hole, and finally relativistic SPH viscosity and energy.
Conduction Modelling Using Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Cleary, Paul W.; Monaghan, Joseph J.
1999-01-01
Heat transfer is very important in many industrial and geophysical problems. Because these problems often have complicated fluid dynamics, there are advantages in solving them using Lagrangian methods like smoothed particle hydrodynamics (SPH). Since SPH particles become disordered, the second derivative terms may be estimated poorly, especially when materials with different properties are adjacent. In this paper we show how a simple alteration to the standard SPH formulation ensures continuity of heat flux across discontinuities in material properties. A set of rules is formulated for the construction of isothermal boundaries leading to accurate conduction solutions. A method for accurate prediction of heat fluxes through isothermal boundaries is also given. The accuracy of the SPH conduction solutions is demonstrated through a sequence of test problems of increasing complexity.
Multiscale modeling with smoothed dissipative particle dynamics.
Kulkarni, Pandurang M; Fu, Chia-Chun; Shell, M Scott; Leal, L Gary
2013-06-21
In this work, we consider two issues related to the use of Smoothed Dissipative Particle Dynamics (SDPD) as an intermediate mesoscale model in a multiscale scheme for solution of flow problems when there are local parts of a macroscopic domain that require molecular resolution. The first is to demonstrate that SDPD with different levels of resolution can accurately represent the fluid properties from the continuum scale all the way to the molecular scale. Specifically, while the thermodynamic quantities such as temperature, pressure, and average density remain scale-invariant, we demonstrate that the dynamic properties are quantitatively consistent with an all-atom Lennard-Jones reference system when the SDPD resolution approaches the atomistic scale. This supports the idea that SDPD can serve as a natural bridge between molecular and continuum descriptions. In the second part, a simple multiscale methodology is proposed within the SDPD framework that allows several levels of resolution within a single domain. Each particle is characterized by a unique physical length scale called the smoothing length, which is inversely related to the local number density and can change on-the-fly. This multiscale methodology is shown to accurately reproduce fluid properties for the simple problem of steady and transient shear flow. PMID:23802949
A Particle-Particle Collision Model for Smoothed Profile Method
NASA Astrophysics Data System (ADS)
Mohaghegh, Fazlolah; Mousel, John; Udaykumar, H. S.
2014-11-01
Smoothed Profile Method (SPM) is a type of continuous forcing approach that adds the particles to the fluid using a forcing. The fluid-structure interaction is through a diffuse interface which avoids sudden transition from solid to fluid. The SPM simulation as a monolithic approach uses an indicator function field in the whole domain based on the distance from each particle's boundary where the possible particle-particle interaction can occur. A soft sphere potential based on the indicator function field has been defined to add an artificial pressure to the flow pressure in the potential overlapping regions. Thus, a repulsion force is obtained to avoid overlapping. Study of two particles which impulsively start moving in an initially uniform flow shows that the particle in the wake of the other one will have less acceleration leading to frequent collisions. Various Reynolds numbers and initial distances have been chosen to test the robustness of the method. Study of Drafting-Kissing Tumbling of two cylindrical particles shows a deviation from the benchmarks due to lack of rotation modeling. The method is shown to be accurate enough for simulating particle-particle collision and can easily be extended for particle-wall modeling and for non-spherical particles.
A Smoothed Particle Hydrodynamics approach for poroelasticity
NASA Astrophysics Data System (ADS)
Osorno, Maria; Steeb, Holger
2016-04-01
Within the framework of the SHynergie project we look to investigate hydraulic fracturing and crack evolving in poroelastic media. We model biphasic media assuming incompressible solid grain and incompressible pore liquid. Modeling evolving fractures and fracture networks in elastic and poroelastic media by mesh-based numerical approaches, like X-FEM, is especially in 3-dim a challenging task. Therefore, we propose a meshless particle method for fractured media based on the Smoothed Particle Hydrodynamics (SPH) approach. SPH is a meshless Lagrangian method highly suitable for the simulation of large deformations including free surfaces and/or interfaces. Within the SPH method, the computational domain is discretized with particles, avoiding the computational expenses of meshing. Our SPH solution is implemented in a parallel computational framework, which allows to simulate large domains more representative of the scale of our study cases. Our implementation is carefully validated against classical mesh-based approaches and compared with classical solutions for consolidation problems. Furthermore, we discuss fracture initiation and propagation in poroelastic rocks at the reservoir scale.
An analysis of smoothed particle hydrodynamics
Swegle, J.W.; Attaway, S.W.; Heinstein, M.W.; Mello, F.J.; Hicks, D.L.
1994-03-01
SPH (Smoothed Particle Hydrodynamics) is a gridless Lagrangian technique which is appealing as a possible alternative to numerical techniques currently used to analyze high deformation impulsive loading events. In the present study, the SPH algorithm has been subjected to detailed testing and analysis to determine its applicability in the field of solid dynamics. An important result of the work is a rigorous von Neumann stability analysis which provides a simple criterion for the stability or instability of the method in terms of the stress state and the second derivative of the kernel function. Instability, which typically occurs only for solids in tension, results not from the numerical time integration algorithm, but because the SPH algorithm creates an effective stress with a negative modulus. The analysis provides insight into possible methods for removing the instability. Also, SPH has been coupled into the transient dynamics finite element code PRONTO, and a weighted residual derivation of the SPH equations has been obtained.
Computational brittle fracture using smooth particle hydrodynamics
Mandell, D.A.; Wingate, C.A.; Schwalbe, L.A.
1996-10-01
We are developing statistically based, brittle-fracture models and are implementing them into hydrocodes that can be used for designing systems with components of ceramics, glass, and/or other brittle materials. Because of the advantages it has simulating fracture, we are working primarily with the smooth particle hydrodynamics code SPBM. We describe a new brittle fracture model that we have implemented into SPBM. To illustrate the code`s current capability, we have simulated a number of experiments. We discuss three of these simulations in this paper. The first experiment consists of a brittle steel sphere impacting a plate. The experimental sphere fragment patterns are compared to the calculations. The second experiment is a steel flyer plate in which the recovered steel target crack patterns are compared to the calculated crack patterns. We also briefly describe a simulation of a tungsten rod impacting a heavily confined alumina target, which has been recently reported on in detail.
Chemical mixing in smoothed particle hydrodynamics simulations
NASA Astrophysics Data System (ADS)
Greif, Thomas H.; Glover, Simon C. O.; Bromm, Volker; Klessen, Ralf S.
2009-02-01
We introduce a simple and efficient algorithm for diffusion in smoothed particle hydrodynamics (SPH) simulations and apply it to the problem of chemical mixing. Based on the concept of turbulent diffusion, we link the diffusivity of a pollutant to the local physical conditions and can thus resolve mixing in space and time. We apply our prescription to the evolution of an idealized supernova remnant and find that we can model the distribution of heavy elements without having to explicitly resolve hydrodynamic instabilities in the post-shock gas. Instead, the dispersal of the pollutant is implicitly modelled through its dependence on the local velocity dispersion. Our method can thus be used in any SPH simulation that investigates chemical mixing but lacks the necessary resolution on small scales. Potential applications include the enrichment of the interstellar medium in present-day galaxies, as well as the intergalactic medium at high redshifts.
An implicit Smooth Particle Hydrodynamic code
Charles E. Knapp
2000-04-01
An implicit version of the Smooth Particle Hydrodynamic (SPH) code SPHINX has been written and is working. In conjunction with the SPHINX code the new implicit code models fluids and solids under a wide range of conditions. SPH codes are Lagrangian, meshless and use particles to model the fluids and solids. The implicit code makes use of the Krylov iterative techniques for solving large linear-systems and a Newton-Raphson method for non-linear corrections. It uses numerical derivatives to construct the Jacobian matrix. It uses sparse techniques to save on memory storage and to reduce the amount of computation. It is believed that this is the first implicit SPH code to use Newton-Krylov techniques, and is also the first implicit SPH code to model solids. A description of SPH and the techniques used in the implicit code are presented. Then, the results of a number of tests cases are discussed, which include a shock tube problem, a Rayleigh-Taylor problem, a breaking dam problem, and a single jet of gas problem. The results are shown to be in very good agreement with analytic solutions, experimental results, and the explicit SPHINX code. In the case of the single jet of gas case it has been demonstrated that the implicit code can do a problem in much shorter time than the explicit code. The problem was, however, very unphysical, but it does demonstrate the potential of the implicit code. It is a first step toward a useful implicit SPH code.
Simulating Ice Particle Melting using Smooth Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Kuo, Kwo-Sen; Pelissier, Craig
2015-04-01
To measure precipitation from space requires an accurate estimation of the collective scattering properties of particles suspended in a precipitating column. It is well known that the complicated and typically unknowable shapes of the solid precipitation particles cause much uncertainty in the retrievals involving such particles. This remote-sensing problem becomes even more difficult with the "melting layer" containing partially melted ice particles, where both the geometric shape and liquid-solid fraction of the hydrometeors are variables.. For the scattering properties of these particles depend not only on their shapes, but also their melt-water fraction,and the spatial distribution of liquid and ice within. To obtain an accurate estimation thus requires a set of "realistic" particle geometries and a method to determine the melt-water distribution at various stages in the melting process. Once this is achieved, a suitable method can be used to compute the scattering properties. In previous work, the growth of a set of astoundingly realistic ice particles has been simulated using the "Snowfake" algorithm of Gravner and Griffeath. To simulate the melting process of these particles, the method of Smooth Particle Hydrodynamics (SPH) is used. SPH is a mesh-less particle-based approach where kinematic and thermal dynamics is controlled entirely through two-body interactions between neighboring SPH particles. An important property of SPH is that the interaction at boundaries between air/ice/water is implicitly taken care of. This is crucial for this work since those boundaries are complex and vary throughout the melting process. We present the SPH implementation and a simulation, using highly parallel Graphic Processing Units (GPUs), with ~1 million SPH particles to represent one of the generated ice particle geometries. We plan to use this method, especially its parallelized version, to simulate the melting of all the "Snowfake" particles (~10,000 of them) in our
NASA Astrophysics Data System (ADS)
Stevens, Mark R.; Gutchess, Dan; Checka, Neal; Snorrason, Magnús
2006-05-01
Image exploitation algorithms for Intelligence, Surveillance and Reconnaissance (ISR) and weapon systems are extremely sensitive to differences between the operating conditions (OCs) under which they are trained and the extended operating conditions (EOCs) in which the fielded algorithms are tested. As an example, terrain type is an important OC for the problem of tracking hostile vehicles from an airborne camera. A system designed to track cars driving on highways and on major city streets would probably not do well in the EOC of parking lots because of the very different dynamics. In this paper, we present a system we call ALPS for Adaptive Learning in Particle Systems. ALPS takes as input a sequence of video images and produces labeled tracks. The system detects moving targets and tracks those targets across multiple frames using a multiple hypothesis tracker (MHT) tightly coupled with a particle filter. This tracker exploits the strengths of traditional MHT based tracking algorithms by directly incorporating tree-based hypothesis considerations into the particle filter update and resampling steps. We demonstrate results in a parking lot domain tracking objects through occlusions and object interactions.
An investigation of particles suspension using smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Pazouki, Arman; Negrut, Dan
2013-11-01
This contribution outlines a method for the direct numerical simulation of rigid body suspensions in a Lagrangian-Lagrangian framework using extended Smoothed Particle Hydrodynamics (XSPH) method. The dynamics of the arbitrarily shaped rigid bodies is fully resolved via Boundary Condition Enforcing (BCE) markers and updated according to the general Newton-Euler equations of motion. The simulation tool, refered to herien as Chrono::Fluid, relies on a parallel implementation that runs on Graphics Processing Unit (GPU) cards. The simulation results obtained for transient Poiseuille flow, migration of cylinder and sphere in Poiseuille flow, and distribution of particles at different cross sections of the laminar flow of dilute suspension were respectively within 0.1%, 1%, and 5% confidence interval of analytical and experimental results reported in the literature. It was shown that at low Reynolds number, Re = O(1), the radial migration (a) behaves non-monotonically as the particles relative distance (distance over diameter) increases from zero to two; and (b) decreases as the particle skewness and size increases. The scaling of Chrono::Fluid was demonstrated in conjunction with a suspension dynamics analysis in which the number of ellipsoids went up to 3e4. Financial support was provided in part by National Science Foundation grant NSF CMMI-084044.
Smoothed-particle hydrodynamics and nonequilibrium molecular dynamics
Hoover, W. G.; Hoover, C. G.
1993-08-01
Gingold, Lucy, and Monaghan invented a grid-free version of continuum mechanics ``smoothed-particle hydrodynamics,`` in 1977. It is a likely contributor to ``hybrid`` simulations combining atomistic and continuum simulations. We describe applications of this particle-based continuum technique from the closely-related standpoint of nonequilibrium molecular dynamics. We compare chaotic Lyapunov spectra for atomistic solids and fluids with those which characterize a two-dimensional smoothed-particle fluid system.
Guided Adaptive Image Smoothing via Directional Anisotropic Structure Measurement.
Zang, Yu; Huang, Hua; Zhang, Lei
2015-09-01
Image smoothing prefers a good metric to identify dominant structures from textures adaptive of intensity contrast. In this paper, we drop on a novel directional anisotropic structure measurement (DASM) toward adaptive image smoothing. With observations on psychological perception regarding anisotropy, non-periodicity and local directionality, DASM can well characterize structures and textures independent on their contrast scales. By using such measurement as constraint, we design a guided adaptive image smoothing scheme by improving extrema localization and envelopes construction in a structure-aware manner. Our approach can well suppresses the staircase-like artifacts and blur of structures that appear in previous methods, which better suits structure-preserving image smoothing task. The algorithm is performed on a space-filling curve as the reduced domain, so it is very fast and much easy to implement in practice. We make comprehensive comparisons with previous state-of-the-art methods for a variety of applications. Experimental results demonstrate the merit using our DASM as metric to identify structures, and the effectiveness and efficiency of our adaptive image smoothing approach to produce commendable results. PMID:26357284
Application of smoothed particle hydrodynamics method in aerodynamics
NASA Astrophysics Data System (ADS)
Cortina, Miguel
2014-11-01
Smoothed Particle Hydrodynamics (SPH) is a meshless Lagrangian method in which the domain is represented by particles. Each particle is assigned properties such as mass, pressure, density, temperature, and velocity. These properties are then evaluated at the particle positions using a smoothing kernel that integrates over the values of the surrounding particles. In the present study the SPH method is first used to obtain numerical solutions for fluid flows over a cylinder and then we are going to apply the same principle over an airfoil obstacle.
Likelihood Methods for Adaptive Filtering and Smoothing. Technical Report #455.
ERIC Educational Resources Information Center
Butler, Ronald W.
The dynamic linear model or Kalman filtering model provides a useful methodology for predicting the past, present, and future states of a dynamic system, such as an object in motion or an economic or social indicator that is changing systematically with time. Recursive likelihood methods for adaptive Kalman filtering and smoothing are developed.…
Smoothed Particle Hydrodynamics with Time Varying, Piecewise Constant Smoothing Length Profiles
NASA Astrophysics Data System (ADS)
Børve, S.; Omang, M.; Trulsen, J.
2000-12-01
Smoothed Particle Hydrodynamics (SPH) has proven to be a very useful numerical tool in studying a number of widely different astrophysical problems. Still, used on many other types of problems the method faces problems concerning efficiency and accuracy compared to that of modern grid-based methods. Essential to efficiency is maintaining a near-optimal particle distribution and smoothing length profile that reflects the physics of the problem. This means, directing computer resources towards those regions and time intervals where the action is taking place and not being wasted where nothing is happening. In the literature researchers have tried to achieve these goals by combining the Lagrangian nature of the SPH method with a smoothing length profile varying smoothly in space and time. To make the SPH method better suited for accurately describing a wider range of problems, a scheme containing two novel features is proposed. First, the scheme assumes a piecewise constant smoothing length profile. To avoid substantial errors near steps in the smoothing length profile, alternative forms of the SPH equations of motion is used. Secondly, a predictive attitude towards optimizing the particle distribution is introduced by activating a mass, momentum and internal energy conservation regularization process at intervals. The main challenge faced by the scheme has been to put the newly optimized smoothing length profile into use without severely altering the underlying physics. To achieve this, the entire set of particles is redefined in the process. The basic ideas behind this scheme is briefly described. Finally, the results from several hydrodynamical and magnetohydrodynamical tests in one and two dimensions are presented. This work is funded by the Research Council of Norway.
Rapid Structured Volume Grid Smoothing and Adaption Technique
NASA Technical Reports Server (NTRS)
Alter, Stephen J.
2006-01-01
A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reductions in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.
Rapid Structured Volume Grid Smoothing and Adaption Technique
NASA Technical Reports Server (NTRS)
Alter, Stephen J.
2004-01-01
A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reduction in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.
Particle systems for adaptive, isotropic meshing of CAD models
Levine, Joshua A.; Whitaker, Ross T.
2012-01-01
We present a particle-based approach for generating adaptive triangular surface and tetrahedral volume meshes from computer-aided design models. Input shapes are treated as a collection of smooth, parametric surface patches that can meet non-smoothly on boundaries. Our approach uses a hierarchical sampling scheme that places particles on features in order of increasing dimensionality. These particles reach a good distribution by minimizing an energy computed in 3D world space, with movements occurring in the parametric space of each surface patch. Rather than using a pre-computed measure of feature size, our system automatically adapts to both curvature as well as a notion of topological separation. It also enforces a measure of smoothness on these constraints to construct a sizing field that acts as a proxy to piecewise-smooth feature size. We evaluate our technique with comparisons against other popular triangular meshing techniques for this domain. PMID:23162181
Adaptive response of pulmonary arterial smooth muscle to length change.
Syyong, Harley; Cheung, Christine; Solomon, Dennis; Seow, Chun Y; Kuo, Kuo H
2008-04-01
Hypervasoconstriction is associated with pulmonary hypertension and dysfunction of the pulmonary arterial smooth muscle (PASM) is implicated. However, relatively little is known about the mechanical properties of PASM. Recent advances in our understanding of plastic adaptation in smooth muscle may shed light on the disease mechanism. In this study, we determined whether PASM is capable of adapting to length changes (especially shortening) and regain its contractile force. We examined the time course of length adaptation in PASM in response to step changes in length and to length oscillations mimicking the periodic stretches due to pulsatile arterial pressure. Rings from sheep pulmonary artery were mounted on myograph and stimulated using electrical field stimulation (12-16 s, 20 V, 60 Hz). The length-force relationship was determined at L(ref) to 0.6 L(ref), where L(ref) was a reference length close to the in situ length of PASM. The response to length oscillations was determined at L(ref), after the muscle was subjected to length oscillation of various amplitudes for 200 s at 1.5 Hz. Release (or stretch) of resting PASM from L(ref) to 0.6 (and vice versa) was followed by a significant force recovery (73 and 63%, respectively), characteristic of length adaptation. All recoveries of force followed a monoexponential time course. Length oscillations with amplitudes ranging from 5 to 20% L(ref) caused no significant change in force generation in subsequent contractions. It is concluded that, like many smooth muscles, PASM possesses substantial capability to adapt to changes in length. Under pathological conditions, this could contribute to hypervasoconstriction in pulmonary hypertension. PMID:18218913
Smoothed Particle Hydrodynamics Model for Reactive Transport and Mineral Precipitation
Tartakovsky, Alexandre M.; Scheibe, Timothy D.; Redden, George; Meakin, Paul; Fang, Yilin
2006-06-30
A new Lagrangian particle model based on smoothed particle hydrodynamics was used to simulate pore scale precipitation reactions. The side-by-side injection of reacting solutions into two halves of a two-dimensional granular porous medium was simulated. Precipitation on grain surfaces occurred along a narrow zone in the middle of the domain, where the reacting solutes mixed to generate a supersaturated reaction product. The numerical simulations qualitatively reproduced the behavior observed in related laboratory experiments.
SPH (smoothed particle hydrodynamics) simulations of hypervelocity impacts
Cloutman, L.D.
1991-01-24
The smoothed particle hydrodynamics (SPH) method has been used to simulate several cases of hypervelocity impact in an exploratory study to determine the suitability of the method for such problems. The calculations compare favorably with experimental results and with other numerical simulations. We discuss the requirements that must be satisfied for SPH to produce accurate simulations of such problems. 18 refs., 9 figs.
Water Flow Simulation using Smoothed Particle Hydrodynamics (SPH)
NASA Technical Reports Server (NTRS)
Vu, Bruce; Berg, Jared; Harris, Michael F.
2014-01-01
Simulation of water flow from the rainbird nozzles has been accomplished using the Smoothed Particle Hydrodynamics (SPH). The advantage of using SPH is that no meshing is required, thus the grid quality is no longer an issue and accuracy can be improved.
Launch Environment Water Flow Simulations Using Smoothed Particle Hydrodynamics
NASA Technical Reports Server (NTRS)
Vu, Bruce T.; Berg, Jared J.; Harris, Michael F.; Crespo, Alejandro C.
2015-01-01
This paper describes the use of Smoothed Particle Hydrodynamics (SPH) to simulate the water flow from the rainbird nozzle system used in the sound suppression system during pad abort and nominal launch. The simulations help determine if water from rainbird nozzles will impinge on the rocket nozzles and other sensitive ground support elements.
Visual enhancement of unmixed multispectral imagery using adaptive smoothing
Lemeshewsky, G.P.
2004-01-01
Adaptive smoothing (AS) has been previously proposed as a method to smooth uniform regions of an image, retain contrast edges, and enhance edge boundaries. The method is an implementation of the anisotropic diffusion process which results in a gray scale image. This paper discusses modifications to the AS method for application to multi-band data which results in a color segmented image. The process was used to visually enhance the three most distinct abundance fraction images produced by the Lagrange constraint neural network learning-based unmixing of Landsat 7 Enhanced Thematic Mapper Plus multispectral sensor data. A mutual information-based method was applied to select the three most distinct fraction images for subsequent visualization as a red, green, and blue composite. A reported image restoration technique (partial restoration) was applied to the multispectral data to reduce unmixing error, although evaluation of the performance of this technique was beyond the scope of this paper. The modified smoothing process resulted in a color segmented image with homogeneous regions separated by sharpened, coregistered multiband edges. There was improved class separation with the segmented image, which has importance to subsequent operations involving data classification.
A method of smoothed particle hydrodynamics using spheroidal kernels
NASA Technical Reports Server (NTRS)
Fulbright, Michael S.; Benz, Willy; Davies, Melvyn B.
1995-01-01
We present a new method of three-dimensional smoothed particle hydrodynamics (SPH) designed to model systems dominated by deformation along a preferential axis. These systems cause severe problems for SPH codes using spherical kernels, which are best suited for modeling systems which retain rough spherical symmetry. Our method allows the smoothing length in the direction of the deformation to evolve independently of the smoothing length in the perpendicular plane, resulting in a kernel with a spheroidal shape. As a result the spatial resolution in the direction of deformation is significantly improved. As a test case we present the one-dimensional homologous collapse of a zero-temperature, uniform-density cloud, which serves to demonstrate the advantages of spheroidal kernels. We also present new results on the problem of the tidal disruption of a star by a massive black hole.
Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
NASA Astrophysics Data System (ADS)
Dehnen, Walter; Aly, Hossam
2012-09-01
The numerical convergence of smoothed particle hydrodynamics (SPH) can be severely restricted by random force errors induced by particle disorder, especially in shear flows, which are ubiquitous in astrophysics. The increase in the number NH of neighbours when switching to more extended smoothing kernels at fixed resolution (using an appropriate definition for the SPH resolution scale) is insufficient to combat these errors. Consequently, trading resolution for better convergence is necessary, but for traditional smoothing kernels this option is limited by the pairing (or clumping) instability. Therefore, we investigate the suitability of the Wendland functions as smoothing kernels and compare them with the traditional B-splines. Linear stability analysis in three dimensions and test simulations demonstrate that the Wendland kernels avoid the pairing instability for all NH, despite having vanishing derivative at the origin (disproving traditional ideas about the origin of this instability; instead, we uncover a relation with the kernel Fourier transform and give an explanation in terms of the SPH density estimator). The Wendland kernels are computationally more convenient than the higher order B-splines, allowing large NH and hence better numerical convergence (note that computational costs rise sublinear with NH). Our analysis also shows that at low NH the quartic spline kernel with NH ≈ 60 obtains much better convergence than the standard cubic spline.
Using the symmetries in 3D Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Cruz-Pérez, Juan P.; González, José A.
2011-10-01
We apply the cartoon SPH technique to solve numerically the Newtonian Euler's equations in scenarios with spherical symmetry with the standard Smoothed Particle Hydrodynamics method, using the equations written in cartesian coordinates. This implementation allow us to increase the resolution of the simulations in order to obtain accurate results. We test it using the Shock tube and the isothermal collapse of a perfect fluid in spherical symmetry.
Hyperbolic Divergence Cleaning Method for Godunov Smoothed Particle Magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Iwasaki, K.; Inutsuka, S.-I.
2013-04-01
In this paper, we implement a divergence cleaning method into Godunov smoothed particle magnetohydrodynamics (GSPM). In the GSPM, to describe MHD shocks accurately, a Riemann solver is applied to the SPH method instead of artificial viscosity and resistivity that have been used in previous works. We confirmed that the divergence cleaning method reduces divergence errors significantly. The performance of the method is demonstrated in the numerical simulations of a strongly magnetized gas and bipolar outflow from the first core.
Numerical simulations of glass impacts using smooth particle hydrodynamics
Mandell, D.A.; Wingate, C.A.
1996-05-01
As part of a program to develop advanced hydrocode design tools, we have implemented a brittle fracture model for glass into the SPHINX smooth particle hydrodynamics code. We have evaluated this model and the code by predicting data from one-dimensional flyer plate impacts into glass. Since fractured glass properties, which are needed in the model, are not available, we did sensitivity studies of these properties, as well as sensitivity studies to determine the number of particles needed in the calculations. The numerical results are in good agreement with the data. {copyright} {ital 1996 American Institute of Physics.}
Numerical simulations of glass impacts using smooth particle hydrodynamics
Mandell, D.A.; Wingate, C.A.
1995-07-01
As part of a program to develop advanced hydrocode design tools, we have implemented a brittle fracture model for glass into the SPHINX smooth particle hydrodynamics code. We have evaluated this model and the code by predicting data from one-dimensional flyer plate impacts into glass. Since fractured glass properties, which are needed in the model, are not available, we did sensitivity studies of these properties, as well as sensitivity studies to determine the number of particles needed in the calculations. The numerical results are in good agreement with the data.
SPLASH: An Interactive Visualization Tool for Smoothed Particle Hydrodynamics Simulations
NASA Astrophysics Data System (ADS)
Price, Daniel J.
2011-03-01
SPLASH (formerly SUPERSPHPLOT) is a visualization tool for output from (astrophysical) simulations using the Smoothed Particle Hydrodynamics (SPH) method in one, two and three dimensions. It is written in Fortran 90 and utilises the PGPLOT graphics subroutine library to do the actual plotting. It is based around a command-line menu structure but utilises the interactive capabilities of PGPLOT to manipulate data interactively in the plotting window. SPLASH is a fully interactive program; visualizations can be changed rapidly at the touch of a button (e.g. zooming, rotating, shifting cross section positions etc). Data is read directly from the code dump format giving rapid access to results and the visualization is advanced forwards and backwards through timesteps by single keystrokes. SPLASH uses the SPH kernel to render plots of not only density but other physical quantities, giving a smooth representation of the data.
Modeling nanoscale hydrodynamics by smoothed dissipative particle dynamics
NASA Astrophysics Data System (ADS)
Lei, Huan; Mundy, Christopher J.; Schenter, Gregory K.; Voulgarakis, Nikolaos K.
2015-05-01
Thermal fluctuation and hydrophobicity are two hallmarks of fluid hydrodynamics on the nano-scale. It is a challenge to consistently couple the small length and time scale phenomena associated with molecular interaction with larger scale phenomena. The development of this consistency is the essence of mesoscale science. In this study, we use a nanoscale fluid model based on smoothed dissipative particle dynamics that accounts for the phenomena associated with density fluctuations and hydrophobicity. We show consistency in the fluctuation spectrum across scales. In doing so, it is necessary to account for finite fluid particle size. Furthermore, we demonstrate that the present model can capture the void probability and solvation free energy of nonpolar hard particles of different sizes. The present fluid model is well suited for an understanding of emergent phenomena in nano-scale fluid systems.
Multiscale simulation of ideal mixtures using smoothed dissipative particle dynamics.
Petsev, Nikolai D; Leal, L Gary; Shell, M Scott
2016-02-28
Smoothed dissipative particle dynamics (SDPD) [P. Español and M. Revenga, Phys. Rev. E 67, 026705 (2003)] is a thermodynamically consistent particle-based continuum hydrodynamics solver that features scale-dependent thermal fluctuations. We obtain a new formulation of this stochastic method for ideal two-component mixtures through a discretization of the advection-diffusion equation with thermal noise in the concentration field. The resulting multicomponent approach is consistent with the interpretation of the SDPD particles as moving volumes of fluid and reproduces the correct fluctuations and diffusion dynamics. Subsequently, we provide a general multiscale multicomponent SDPD framework for simulations of molecularly miscible systems spanning length scales from nanometers to the non-fluctuating continuum limit. This approach reproduces appropriate equilibrium properties and is validated with simulation of simple one-dimensional diffusion across multiple length scales. PMID:26931689
Multiscale simulation of ideal mixtures using smoothed dissipative particle dynamics
NASA Astrophysics Data System (ADS)
Petsev, Nikolai D.; Leal, L. Gary; Shell, M. Scott
2016-02-01
Smoothed dissipative particle dynamics (SDPD) [P. Español and M. Revenga, Phys. Rev. E 67, 026705 (2003)] is a thermodynamically consistent particle-based continuum hydrodynamics solver that features scale-dependent thermal fluctuations. We obtain a new formulation of this stochastic method for ideal two-component mixtures through a discretization of the advection-diffusion equation with thermal noise in the concentration field. The resulting multicomponent approach is consistent with the interpretation of the SDPD particles as moving volumes of fluid and reproduces the correct fluctuations and diffusion dynamics. Subsequently, we provide a general multiscale multicomponent SDPD framework for simulations of molecularly miscible systems spanning length scales from nanometers to the non-fluctuating continuum limit. This approach reproduces appropriate equilibrium properties and is validated with simulation of simple one-dimensional diffusion across multiple length scales.
Modeling nanoscale hydrodynamics by smoothed dissipative particle dynamics
Lei, Huan; Mundy, Christopher J.; Schenter, Gregory K.; Voulgarakis, Nikolaos
2015-05-21
Thermal fluctuation and hydrophobicity are two hallmarks of fluid hydrodynamics on the nano-scale. It is a challenge to consistently couple the small length and time scale phenomena associated with molecular interaction with larger scale phenomena. The development of this consistency is the essence of mesoscale science. In this study, we develop a nanoscale fluid model based on smoothed dissipative particle dynamics that accounts for the phenomena of associated with density fluctuations and hydrophobicity. We show consistency in the fluctuation spectrum across scales. In doing so, it is necessary to account for finite fluid particle size. Furthermore, we demonstrate that the present model can capture of the void probability and solvation free energy of apolar particles of different sizes. The present fluid model is well suited for a understanding emergent phenomena in nano-scale fluid systems.
Numerical solution of the nonlinear Schrödinger equation using smoothed-particle hydrodynamics.
Mocz, Philip; Succi, Sauro
2015-05-01
We formulate a smoothed-particle hydrodynamics numerical method, traditionally used for the Euler equations for fluid dynamics in the context of astrophysical simulations, to solve the nonlinear Schrödinger equation in the Madelung formulation. The probability density of the wave function is discretized into moving particles, whose properties are smoothed by a kernel function. The traditional fluid pressure is replaced by a quantum pressure tensor, for which a robust discretization is found. We demonstrate our numerical method on a variety of numerical test problems involving the simple harmonic oscillator, soliton-soliton collision, Bose-Einstein condensates, collapsing singularities, and dark matter halos governed by the Gross-Pitaevskii-Poisson equation. Our method is conservative, applicable to unbounded domains, and is automatically adaptive in its resolution, making it well suited to study problems with collapsing solutions. PMID:26066276
Numerical solution of the nonlinear Schrödinger equation using smoothed-particle hydrodynamics
NASA Astrophysics Data System (ADS)
Mocz, Philip; Succi, Sauro
2015-05-01
We formulate a smoothed-particle hydrodynamics numerical method, traditionally used for the Euler equations for fluid dynamics in the context of astrophysical simulations, to solve the nonlinear Schrödinger equation in the Madelung formulation. The probability density of the wave function is discretized into moving particles, whose properties are smoothed by a kernel function. The traditional fluid pressure is replaced by a quantum pressure tensor, for which a robust discretization is found. We demonstrate our numerical method on a variety of numerical test problems involving the simple harmonic oscillator, soliton-soliton collision, Bose-Einstein condensates, collapsing singularities, and dark matter halos governed by the Gross-Pitaevskii-Poisson equation. Our method is conservative, applicable to unbounded domains, and is automatically adaptive in its resolution, making it well suited to study problems with collapsing solutions.
Smoothed dissipative particle dynamics with angular momentum conservation
Müller, Kathrin Fedosov, Dmitry A. Gompper, Gerhard
2015-01-15
Smoothed dissipative particle dynamics (SDPD) combines two popular mesoscopic techniques, the smoothed particle hydrodynamics and dissipative particle dynamics (DPD) methods, and can be considered as an improved dissipative particle dynamics approach. Despite several advantages of the SDPD method over the conventional DPD model, the original formulation of SDPD by Español and Revenga (2003) [9], lacks angular momentum conservation, leading to unphysical results for problems where the conservation of angular momentum is essential. To overcome this limitation, we extend the SDPD method by introducing a particle spin variable such that local and global angular momentum conservation is restored. The new SDPD formulation (SDPD+a) is directly derived from the Navier–Stokes equation for fluids with spin, while thermal fluctuations are incorporated similarly to the DPD method. We test the new SDPD method and demonstrate that it properly reproduces fluid transport coefficients. Also, SDPD with angular momentum conservation is validated using two problems: (i) the Taylor–Couette flow with two immiscible fluids and (ii) a tank-treading vesicle in shear flow with a viscosity contrast between inner and outer fluids. For both problems, the new SDPD method leads to simulation predictions in agreement with the corresponding analytical theories, while the original SDPD method fails to capture properly physical characteristics of the systems due to violation of angular momentum conservation. In conclusion, the extended SDPD method with angular momentum conservation provides a new approach to tackle fluid problems such as multiphase flows and vesicle/cell suspensions, where the conservation of angular momentum is essential.
Smoothed dissipative particle dynamics with angular momentum conservation
NASA Astrophysics Data System (ADS)
Müller, Kathrin; Fedosov, Dmitry A.; Gompper, Gerhard
2015-01-01
Smoothed dissipative particle dynamics (SDPD) combines two popular mesoscopic techniques, the smoothed particle hydrodynamics and dissipative particle dynamics (DPD) methods, and can be considered as an improved dissipative particle dynamics approach. Despite several advantages of the SDPD method over the conventional DPD model, the original formulation of SDPD by Español and Revenga (2003) [9], lacks angular momentum conservation, leading to unphysical results for problems where the conservation of angular momentum is essential. To overcome this limitation, we extend the SDPD method by introducing a particle spin variable such that local and global angular momentum conservation is restored. The new SDPD formulation (SDPD+a) is directly derived from the Navier-Stokes equation for fluids with spin, while thermal fluctuations are incorporated similarly to the DPD method. We test the new SDPD method and demonstrate that it properly reproduces fluid transport coefficients. Also, SDPD with angular momentum conservation is validated using two problems: (i) the Taylor-Couette flow with two immiscible fluids and (ii) a tank-treading vesicle in shear flow with a viscosity contrast between inner and outer fluids. For both problems, the new SDPD method leads to simulation predictions in agreement with the corresponding analytical theories, while the original SDPD method fails to capture properly physical characteristics of the systems due to violation of angular momentum conservation. In conclusion, the extended SDPD method with angular momentum conservation provides a new approach to tackle fluid problems such as multiphase flows and vesicle/cell suspensions, where the conservation of angular momentum is essential.
Tensor classification of structure in smoothed particle hydrodynamics density fields
NASA Astrophysics Data System (ADS)
Forgan, Duncan; Bonnell, Ian; Lucas, William; Rice, Ken
2016-04-01
As hydrodynamic simulations increase in scale and resolution, identifying structures with non-trivial geometries or regions of general interest becomes increasingly challenging. There is a growing need for algorithms that identify a variety of different features in a simulation without requiring a `by eye' search. We present tensor classification as such a technique for smoothed particle hydrodynamics (SPH). These methods have already been used to great effect in N-Body cosmological simulations, which require smoothing defined as an input free parameter. We show that tensor classification successfully identifies a wide range of structures in SPH density fields using its native smoothing, removing a free parameter from the analysis and preventing the need for tessellation of the density field, as required by some classification algorithms. As examples, we show that tensor classification using the tidal tensor and the velocity shear tensor successfully identifies filaments, shells and sheet structures in giant molecular cloud simulations, as well as spiral arms in discs. The relationship between structures identified using different tensors illustrates how different forces compete and co-operate to produce the observed density field. We therefore advocate the use of multiple tensors to classify structure in SPH simulations, to shed light on the interplay of multiple physical processes.
SPHS: smoothed particle hydrodynamics with a higher order dissipation switch
NASA Astrophysics Data System (ADS)
Read, J. I.; Hayfield, T.
2012-06-01
We present a novel implementation of smoothed particle hydrodynamics that uses the spatial derivative of the velocity divergence as a higher order dissipation switch. Our switch - which is second order accurate - detects flow convergence before it occurs. If particle trajectories are going to cross, we switch on the usual SPH artificial viscosity, as well as conservative dissipation in all advected fluid quantities (e.g. the entropy). The viscosity and dissipation terms (that are numerical errors) are designed to ensure that all fluid quantities remain single valued as particles approach one another, to respect conservation laws, and to vanish on a given physical scale as the resolution is increased. SPHS alleviates a number of known problems with 'classic' SPH, successfully resolving mixing, and recovering numerical convergence with increasing resolution. An additional key advantage is that - treating the particle mass similarly to the entropy - we are able to use multimass particles, giving significantly improved control over the refinement strategy. We present a wide range of code tests including the Sod shock tube, Sedov-Taylor blast wave, Kelvin-Helmholtz Instability, the 'blob test' and some convergence tests. Our method performs well on all tests, giving good agreement with analytic expectations.
Simulating Brittle Fracture of Rocks using Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Das, Rajarshi; Cleary, Paul W.
2009-05-01
Numerical modelling can assist in understanding and predicting complex fracture processes. Smoothed Particle Hydrodynamics (SPH) is a particle-based Lagrangian method that is particularly suited to the analysis of fracture due to its capacity to model large deformation and to track free surfaces generated. A damage model is used to predict the fracture of elastic solids. The damage parameter represents the volume-averaged micro-fracture of the volume of material represented by an SPH particle. Evolution of damage is predicted using the strain history of each particle. Damage inhibits the transmission of tensile stress between particles, and once it reaches unity, the interface becomes unable to transmit tensile stress, resulting in a macro-crack. Connected macro-cracks lead to complete fragmentation. In this paper, we explore the ability of an SPH-based damage model to predict brittle fracture of rocks during impact. Rock shape is found to have considerable influence on the fracture process, the fragment sizes, the energy dissipation during impact, and the post-fracture motion of the fragments.
Water pipe flow simulation using improved virtual particles on smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Ting, E. S.; Yeak, S. H.
2014-12-01
Smoothed Particle Hydrodynamics (SPH) is a meshless method used widely to solve problems such as fluid flows. Due to its meshless property, it is ideal to solve problems on complex geometry. In this paper, boundary treatment were implied for the rectangular pipe flow simulations using SPH. The repulsive force is applied to the boundary particles along with the improved virtual particles on different geometry alignment. The water flow is solved using incompressible SPH and will be examined throughout the simulation. Results from this simulation will be compared with single layered virtual particles. Based on the result of the study, it is found that the improved virtual particles is more accurate and stable.
Smoothed particle hydrodynamics simulations of gas and dust mixtures
NASA Astrophysics Data System (ADS)
Booth, R. A.; Sijacki, D.; Clarke, C. J.
2015-10-01
We present a `two-fluid' implementation of dust in smoothed particle hydrodynamics (SPH) in the test particle limit. The scheme is able to handle both short and long stopping times and reproduces the short friction time limit, which is not properly handled in other implementations. We apply novel tests to verify its accuracy and limitations, including multidimensional tests that have not been previously applied to the drag-coupled dust problem and which are particularly relevant to self-gravitating protoplanetary discs. Our tests demonstrate several key requirements for accurate simulations of gas-dust mixtures. First, in standard SPH particle jitter can degrade the dust solution, even when the gas density is well reproduced. The use of integral gradients, a Wendland kernel and a large number of neighbours can control this, albeit at a greater computational cost. Secondly, when it is necessary to limit the artificial viscosity we recommend using the Cullen & Dehnen switch, since the alternative, using α ˜ 0.1, can generate a large velocity noise up to σv ≲ 0.3cs in the dust particles. Thirdly, we find that an accurate dust density estimate requires >400 neighbours, since, unlike the gas, the dust particles do not feel regularization forces. This density noise applies to all particle-based two-fluid implementations of dust, irrespective of the hydro solver and could lead to numerically induced fragmentation. Although our tests show accurate dusty gas simulations are possible, care must be taken to minimize the contribution from numerical noise.
SPHRAY: A Smoothed Particle Hydrodynamics Ray Tracer for Radiative Transfer
NASA Astrophysics Data System (ADS)
Altay, Gabriel; Croft, Rupert A. C.; Pelupessy, Inti
2011-03-01
SPHRAY, a Smoothed Particle Hydrodynamics (SPH) ray tracer, is designed to solve the 3D, time dependent, radiative transfer (RT) equations for arbitrary density fields. The SPH nature of SPHRAY makes the incorporation of separate hydrodynamics and gravity solvers very natural. SPHRAY relies on a Monte Carlo (MC) ray tracing scheme that does not interpolate the SPH particles onto a grid but instead integrates directly through the SPH kernels. Given initial conditions and a description of the sources of ionizing radiation, the code will calculate the non-equilibrium ionization state (HI, HII, HeI, HeII, HeIII, e) and temperature (internal energy/entropy) of each SPH particle. The sources of radiation can include point like objects, diffuse recombination radiation, and a background field from outside the computational volume. The MC ray tracing implementation allows for the quick introduction of new physics and is parallelization friendly. A quick Axis Aligned Bounding Box (AABB) test taken from computer graphics applications allows for the acceleration of the raytracing component. We present the algorithms used in SPHRAY and verify the code by performing all the test problems detailed in the recent Radiative Transfer Comparison Project of Iliev et. al. The Fortran 90 source code for SPHRAY and example SPH density fields are made available online.
Integral approximations to classical diffusion and smoothed particle hydrodynamics
Du, Q.; Lehoucq, Richard B.; Tartakovsky, Alexandre M.
2015-04-01
The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary. The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. An immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.
Application of Smooth-Particle Hydrodynamics in Metal Machining
NASA Astrophysics Data System (ADS)
Zahedi, Abolfazl; Li, Simin; Roy, Anish; Babitsky, Vladimir; Silberschmidt, Vadim V.
2012-08-01
The finite element (FE) method has been extensively used to model complex cutting processes. However, due to large strains in a process zone, leading to increased element distortions, such simulations are confronted with numerical difficulties. Smooth-particle hydrodynamics (SPH) is a mesh-free computational method, which has been used to simulate multi-body problems. In this paper we present a 3D hybrid modelling approach for orthogonal micro-machining of a copper single crystal with the use of SPH and continuum FE. The model is implemented in a commercial FE software ABAQUS/Explicit. The study is used to gain insight into the effects of crystallographic anisotropy on the machining response of f.c.c. cubic metals.
An hourglass control algorithm for Lagrangian Smooth Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Ganzenmüller, Georg C.
2015-04-01
This paper presents a stabilization scheme which addresses the rank-deficiency problem in meshless collocation methods for solid mechanics. Specifically, Smooth-Particle Hydrodynamics (SPH) in the Total Lagrangian formalism is considered. This method is rank-deficient in the sense that the SPH approximation of the deformation gradient is not unique with respect to the positions of the integration points. The non-uniqueness can result in the formation of zero-energy modes. If undetected, these modes can grow and completely dominate the solution. Here, an algorithm is introduced, which effectively suppresses these modes in a fashion similar to hour-glass control mechanisms in Finite-Element methods. Simulations utilizing this control algorithm result exhibit much improved stability, accuracy, and error convergence properties. In contrast to an alternative method which eliminates zero-energy modes, namely the use of additional integration points, the here presented algorithm is easy to implement and computationally very efficient.
Integral approximations to classical diffusion and smoothed particle hydrodynamics
Du, Qiang; Lehoucq, R. B.; Tartakovsky, A. M.
2014-12-31
The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary.more » The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. As a result, an immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.« less
Simulations of dolphin kick swimming using smoothed particle hydrodynamics.
Cohen, Raymond C Z; Cleary, Paul W; Mason, Bruce R
2012-06-01
In competitive human swimming the submerged dolphin kick stroke (underwater undulatory swimming) is utilized after dives and turns. The optimal dolphin kick has a balance between minimizing drag and maximizing thrust while also minimizing the physical exertion required of the swimmer. In this study laser scans of athletes are used to provide realistic swimmer geometries in a single anatomical pose. These are rigged and animated to closely match side-on video footage. Smoothed Particle Hydrodynamics (SPH) fluid simulations are performed to evaluate variants of this swimming stroke technique. This computational approach provides full temporal and spatial information about the flow moving around the deforming swimmer model. The effects of changes in ankle flexibility and stroke frequency are investigated through a parametric study. The results suggest that the net streamwise force on the swimmer is relatively insensitive to ankle flexibility but is strongly dependent on kick frequency. PMID:21840077
Integral approximations to classical diffusion and smoothed particle hydrodynamics
Du, Qiang; Lehoucq, R. B.; Tartakovsky, A. M.
2014-12-31
The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary. The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. As a result, an immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.
Modeling of cast systems using smoothed-particle hydrodynamics
NASA Astrophysics Data System (ADS)
Cleary, Paul; Prakash, Mahesh; Ha, Joseph; Sinnott, Matthew; Nguyen, Thang; Grandfield, John
2004-03-01
To understand and control the filling process for metals in high-pressure die casting and ingot casting, researchers have used new flow-simulation software for the modeling of mold filling. Smoothed-particle hydrodynamics (SPH) is a non-conventional computational fluid dynamics method that has been successfully applied to these problems. Due to its mesh-free nature, it can handle complex splashing free surface flows and the differential motion of multiple solid-casting equipment components relatively easily. The ability of SPH to predict the detailed filling patterns of real large-scale automotive die castings is demonstrated in this study, and the use of SPH simulation for wheel shape optimization in ingot casting based on minimizing oxide generation while increasing the throughput is also presented.
Protostellar jets and magnetised turbulence with smoothed particle magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Tricco, Terrence
2016-01-01
Magnetic fields are an integral component of the formation of stars. During my thesis work, I built new methods to model magnetic fields in smoothed particle magnetohydrodynamics which enforce the divergence-free constraint on the magnetic field and reduce numerical dissipation of the magnetic field. Using these methods, we have performed simulations of isolated protostar formation, studying the production of jets and outflows of material and their effect on transporting angular momentum away from the protostar and reducing the efficiency of star formation. A major code comparison project on the small-scale turbulent dynamo amplification of magnetic fields was performed, using conditions representative of molecular clouds, the formation site of stars. The results were compared against results from grid-based methods, finding excellent agreement on their statistics and qualitative behaviour. I will outline the numerical methods developed, and present the results from our protostar and molecular cloud simulations.
NASA Astrophysics Data System (ADS)
Kordilla, J.; Tartakovsky, A. M.; Geyer, T.
2014-12-01
Unsaturated flow in fractured porous media exhibits highly complex flow dynamics and a wide range of intermittent flow processes. Especially in wide aperture fractures, flow processes may be dominated by gravitational instead of capillary forces leading to a deviation from the classical volume effective approaches (Richard's equation, Van Genuchten type relationships). The existence of various flow modes such as droplets, rivulets, turbulent and adsorbed films is well known, however, their spatial and temporal distribution within fracture networks is still an open question partially due to the lack of appropriate modeling tools. With our work we want to gain a deeper understanding of the underlying flow and transport dynamics in unsaturated fractured media in order to support the development of more refined upscaled methods, applicable on catchment scales. We present pore- and fracture-scale flow simulations obtained with a Smoothed Particle Hydrodynamics (SPH) model. The model allows to simulate free-surface flow dynamics including the effect of surface tension for a wide range of wetting conditions. Several empirical and semi-analytical solutions are used to verify the model. We show that our results satisfy the empirical scaling laws for droplet velocity and critical contact angle. Due to the efficient generation of surface tension via particle-particle interaction forces the dynamic wetting of surfaces as well as the velocity enhancement of droplets on saturated surfaces can readily be obtained. Furthermore, we study the effect of surface roughness on droplet velocities. Lastly, we present flow and transport simulations in the presence of an adjacent porous matrix in order to investigate its influence on the fracture surface flow dynamics and transport across the matrix-fracture interface.
Smoothed particle hydrodynamics modelling for failure in metals
NASA Astrophysics Data System (ADS)
Strand, Russell K.
It is generally regarded to be a difficult task to model multiple fractures leading to fragmentation in metals subjected to high strain rates using numerical methods. Meshless methods such as Smoothed Particle Hydrodynamics (SPH) are well suited to the application of fracture mechanics, since they are not prone to the problems associated with mesh tangling. This research demonstrates and validates a numerical inter-particle fracture model for the initiation, growth and subsequent failure in metals at high strain rate, applicable within a Total Lagrangian SPH scheme. Total Lagrangian SPH performs calculations in the reference state of a material and therefore the neighbourhoods remain fixed throughout the computation; this allows the inter-particle bonds to be stored and tracked as material history parameters. Swegle (2000) showed that the SPH momentum equation can be rearranged in terms of a particle-particle interaction area. By reducing this area to zero via an inter-particle damage parameter, the principles of continuum damage mechanics can be observed without the need for an effective stress term, held at the individual particles.. This research makes use of the Cochran-Banner damage growth model which has been updated for 3D damage and makes the appropriate modifications for inter-particle damage growth. The fracture model was tested on simulations of a 1D flyer plate impact test and the results were compared to experimental data. Some limited modelling was also conducted in 2 and 3 dimensions and promising results were observed. Research was also performed into the mesh sensitivity of the explosively driven Mock- Holt experiment. 3D simulations using the Eulerian SPH formulation were conducted and the best results were observed with a radial packing arrangement. An in-depth assessment of the Monaghan repulsive force correction was also conducted in attempt to eliminate the presence of the SPH tensile instability and stabilise the available Eulerian SPH code
Hydrodynamic simulations with the Godunov smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Murante, G.; Borgani, S.; Brunino, R.; Cha, S.-H.
2011-10-01
We present results based on an implementation of the Godunov smoothed particle hydrodynamics (GSPH), originally developed by Inutsuka, in the GADGET-3 hydrodynamic code. We first review the derivation of the GSPH discretization of the equations of moment and energy conservation, starting from the convolution of these equations with the interpolating kernel. The two most important aspects of the numerical implementation of these equations are (a) the appearance of fluid velocity and pressure obtained from the solution of the Riemann problem between each pair of particles, and (b) the absence of an artificial viscosity term. We carry out three different controlled hydrodynamical three-dimensional tests, namely the Sod shock tube, the development of Kelvin-Helmholtz instabilities in a shear-flow test and the 'blob' test describing the evolution of a cold cloud moving against a hot wind. The results of our tests confirm and extend in a number of aspects those recently obtained by Cha, Inutsuka & Nayakshin: (i) GSPH provides a much improved description of contact discontinuities, with respect to smoothed particle hydrodynamics (SPH), thus avoiding the appearance of spurious pressure forces; (ii) GSPH is able to follow the development of gas-dynamical instabilities, such as the Kevin-Helmholtz and the Rayleigh-Taylor ones; (iii) as a result, GSPH describes the development of curl structures in the shear-flow test and the dissolution of the cold cloud in the 'blob' test. Besides comparing the results of GSPH with those from standard SPH implementations, we also discuss in detail the effect on the performances of GSPH of changing different aspects of its implementation: choice of the number of neighbours, accuracy of the interpolation procedure to locate the interface between two fluid elements (particles) for the solution of the Riemann problem, order of the reconstruction for the assignment of variables at the interface, choice of the limiter to prevent oscillations of
Length adaptation of smooth muscle contractile filaments in response to sustained activation.
Stålhand, Jonas; Holzapfel, Gerhard A
2016-05-21
Airway and bladder smooth muscles are known to undergo length adaptation under sustained contraction. This adaptation process entails a remodelling of the intracellular actin and myosin filaments which shifts the peak of the active force-length curve towards the current length. Smooth muscles are therefore able to generate the maximum force over a wide range of lengths. In contrast, length adaptation of vascular smooth muscle has attracted very little attention and only a handful of studies have been reported. Although their results are conflicting on the existence of a length adaptation process in vascular smooth muscle, it seems that, at least, peripheral arteries and arterioles undergo such adaptation. This is of interest since peripheral vessels are responsible for pressure regulation, and a length adaptation will affect the function of the cardiovascular system. It has, e.g., been suggested that the inward remodelling of resistance vessels associated with hypertension disorders may be related to smooth muscle adaptation. In this study we develop a continuum mechanical model for vascular smooth muscle length adaptation by assuming that the muscle cells remodel the actomyosin network such that the peak of the active stress-stretch curve is shifted towards the operating point. The model is specialised to hamster cheek pouch arterioles and the simulated response to stepwise length changes under contraction. The results show that the model is able to recover the salient features of length adaptation reported in the literature. PMID:26925813
Hybrid molecular-continuum simulations using smoothed dissipative particle dynamics
Petsev, Nikolai D.; Leal, L. Gary; Shell, M. Scott
2015-01-28
We present a new multiscale simulation methodology for coupling a region with atomistic detail simulated via molecular dynamics (MD) to a numerical solution of the fluctuating Navier-Stokes equations obtained from smoothed dissipative particle dynamics (SDPD). In this approach, chemical potential gradients emerge due to differences in resolution within the total system and are reduced by introducing a pairwise thermodynamic force inside the buffer region between the two domains where particles change from MD to SDPD types. When combined with a multi-resolution SDPD approach, such as the one proposed by Kulkarni et al. [J. Chem. Phys. 138, 234105 (2013)], this method makes it possible to systematically couple atomistic models to arbitrarily coarse continuum domains modeled as SDPD fluids with varying resolution. We test this technique by showing that it correctly reproduces thermodynamic properties across the entire simulation domain for a simple Lennard-Jones fluid. Furthermore, we demonstrate that this approach is also suitable for non-equilibrium problems by applying it to simulations of the start up of shear flow. The robustness of the method is illustrated with two different flow scenarios in which shear forces act in directions parallel and perpendicular to the interface separating the continuum and atomistic domains. In both cases, we obtain the correct transient velocity profile. We also perform a triple-scale shear flow simulation where we include two SDPD regions with different resolutions in addition to a MD domain, illustrating the feasibility of a three-scale coupling.
Multi-phase shock simulations with smoothed particle hydrodynamics (SPH)
NASA Astrophysics Data System (ADS)
Omang, M. G.; Trulsen, J. K.
2014-09-01
In this paper we present an approach to the implementation of a multi-phase description in the numerical Smoothed Particle Hydrodynamics method. The work is based on previous work, but has been modified to suit the applications of interest, in this case shock propagation through dusty gases. Theoretical models for multi-phase systems rely on the introduction of a number of terms describing the interaction between the different phases; drag and heat exchange are two examples. These terms contain parameters, the value of many of which must be determined empirically. We present results on the effect of changing values of some of the important parameters and compare our results to experimental and numerical results published in the literature. Our numerical results generally agree well with published results, taking uncertainties concerning accuracy in existing experimental data and details in the choice of parameters for numerical results into consideration. In particular, we find that a reduction in dust particle size is an efficient way of increasing shock retardation for a given dust loading.
Simulations of reactive transport and precipitation with smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Tartakovsky, Alexandre M.; Meakin, Paul; Scheibe, Timothy D.; Eichler West, Rogene M.
2007-03-01
A numerical model based on smoothed particle hydrodynamics (SPH) was developed for reactive transport and mineral precipitation in fractured and porous materials. Because of its Lagrangian particle nature, SPH has several advantages for modeling Navier-Stokes flow and reactive transport including: (1) in a Lagrangian framework there is no non-linear term in the momentum conservation equation, so that accurate solutions can be obtained for momentum dominated flows and; (2) complicated physical and chemical processes such as surface growth due to precipitation/dissolution and chemical reactions are easy to implement. In addition, SPH simulations explicitly conserve mass and linear momentum. The SPH solution of the diffusion equation with fixed and moving reactive solid-fluid boundaries was compared with analytical solutions, Lattice Boltzmann [Q. Kang, D. Zhang, P. Lichtner, I. Tsimpanogiannis, Lattice Boltzmann model for crystal growth from supersaturated solution, Geophysical Research Letters, 31 (2004) L21604] simulations and diffusion limited aggregation (DLA) [P. Meakin, Fractals, scaling and far from equilibrium. Cambridge University Press, Cambridge, UK, 1998] model simulations. To illustrate the capabilities of the model, coupled three-dimensional flow, reactive transport and precipitation in a fracture aperture with a complex geometry were simulated.
Hybrid molecular-continuum simulations using smoothed dissipative particle dynamics
NASA Astrophysics Data System (ADS)
Petsev, Nikolai D.; Leal, L. Gary; Shell, M. Scott
2015-01-01
We present a new multiscale simulation methodology for coupling a region with atomistic detail simulated via molecular dynamics (MD) to a numerical solution of the fluctuating Navier-Stokes equations obtained from smoothed dissipative particle dynamics (SDPD). In this approach, chemical potential gradients emerge due to differences in resolution within the total system and are reduced by introducing a pairwise thermodynamic force inside the buffer region between the two domains where particles change from MD to SDPD types. When combined with a multi-resolution SDPD approach, such as the one proposed by Kulkarni et al. [J. Chem. Phys. 138, 234105 (2013)], this method makes it possible to systematically couple atomistic models to arbitrarily coarse continuum domains modeled as SDPD fluids with varying resolution. We test this technique by showing that it correctly reproduces thermodynamic properties across the entire simulation domain for a simple Lennard-Jones fluid. Furthermore, we demonstrate that this approach is also suitable for non-equilibrium problems by applying it to simulations of the start up of shear flow. The robustness of the method is illustrated with two different flow scenarios in which shear forces act in directions parallel and perpendicular to the interface separating the continuum and atomistic domains. In both cases, we obtain the correct transient velocity profile. We also perform a triple-scale shear flow simulation where we include two SDPD regions with different resolutions in addition to a MD domain, illustrating the feasibility of a three-scale coupling.
Kordilla, Jannes; Tartakovsky, Alexandre M.; Geyer, Tobias
2013-09-01
Flow on fracture surfaces has been identified by many authors as an important flow process in unsaturated fractured rock formations. Given the complexity of flow dynamics on such small scales, robust numerical methods have to be employed in order to capture the highly dynamic interfaces and flow intermittency. In this work we present microscale free-surface flow simulations using a three-dimensional multiphase Smoothed Particle Hydrodynamics (SPH) code. Pairwise solid-fluid and fluid-fluid interaction forces are used to control the wetting behavior and cover a wide range of static and transient contact angles as well as Reynolds numbers encountered in droplet flow on rock surfaces. We validate our model via comparison with existing empirical and semi-analyical solutions for droplet flow. We use the model to investigate the occurence of adsorbed trailing films of droplets under various flow conditions and its importance for the flow dynamics when films and droplets coexist. We show that flow velocities are higher on prewetted surfaces covered by a thin film which is qualitatively attributed to the enhanced dynamic wetting and dewetting at the trailing and advancing contact line.
Smoothed Particle Inference: A Kilo-Parametric Method for X-ray Galaxy Cluster Modeling
Peterson, John R.; Marshall, P.J.; Andersson, K.; /Stockholm U. /SLAC
2005-08-05
We propose an ambitious new method that models the intracluster medium in clusters of galaxies as a set of X-ray emitting smoothed particles of plasma. Each smoothed particle is described by a handful of parameters including temperature, location, size, and elemental abundances. Hundreds to thousands of these particles are used to construct a model cluster of galaxies, with the appropriate complexity estimated from the data quality. This model is then compared iteratively with X-ray data in the form of adaptively binned photon lists via a two-sample likelihood statistic and iterated via Markov Chain Monte Carlo. The complex cluster model is propagated through the X-ray instrument response using direct sampling Monte Carlo methods. Using this approach the method can reproduce many of the features observed in the X-ray emission in a less assumption-dependent way that traditional analyses, and it allows for a more detailed characterization of the density, temperature, and metal abundance structure of clusters. Multi-instrument X-ray analyses and simultaneous X-ray, Sunyaev-Zeldovich (SZ), and lensing analyses are a straight-forward extension of this methodology. Significant challenges still exist in understanding the degeneracy in these models and the statistical noise induced by the complexity of the models.
Dusty gas with one fluid in smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Laibe, Guillaume; Price, Daniel J.
2014-05-01
In a companion paper we have shown how the equations describing gas and dust as two fluids coupled by a drag term can be re-formulated to describe the system as a single-fluid mixture. Here, we present a numerical implementation of the one-fluid dusty gas algorithm using smoothed particle hydrodynamics (SPH). The algorithm preserves the conservation properties of the SPH formalism. In particular, the total gas and dust mass, momentum, angular momentum and energy are all exactly conserved. Shock viscosity and conductivity terms are generalized to handle the two-phase mixture accordingly. The algorithm is benchmarked against a comprehensive suit of problems: DUSTYBOX, DUSTYWAVE, DUSTYSHOCK and DUSTYOSCILL, each of them addressing different properties of the method. We compare the performance of the one-fluid algorithm to the standard two-fluid approach. The one-fluid algorithm is found to solve both of the fundamental limitations of the two-fluid algorithm: it is no longer possible to concentrate dust below the resolution of the gas (they have the same resolution by definition), and the spatial resolution criterion h < csts, required in two-fluid codes to avoid over-damping of kinetic energy, is unnecessary. Implicit time-stepping is straightforward. As a result, the algorithm is up to ten billion times more efficient for 3D simulations of small grains. Additional benefits include the use of half as many particles, a single kernel and fewer SPH interpolations. The only limitation is that it does not capture multi-streaming of dust in the limit of zero coupling, suggesting that in this case a hybrid approach may be required.
Metal Diffusion in Smoothed Particle Hydrodynamics Simulations of Dwarf Galaxies
NASA Astrophysics Data System (ADS)
Williamson, David; Martel, Hugo; Kawata, Daisuke
2016-05-01
We perform a series of smoothed particle hydrodynamics simulations of isolated dwarf galaxies to compare different metal mixing models. In particular, we examine the role of diffusion in the production of enriched outflows and in determining the metallicity distributions of gas and stars. We investigate different diffusion strengths by changing the pre-factor of the diffusion coefficient, by varying how the diffusion coefficient is calculated from the local velocity distribution, and by varying whether the speed of sound is included as a velocity term. Stronger diffusion produces a tighter [O/Fe]–[Fe/H] distribution in the gas and cuts off the gas metallicity distribution function at lower metallicities. Diffusion suppresses the formation of low-metallicity stars, even with weak diffusion, and also strips metals from enriched outflows. This produces a remarkably tight correlation between “metal mass-loading” (mean metal outflow rate divided by mean metal production rate) and the strength of diffusion, even when the diffusion coefficient is calculated in different ways. The effectiveness of outflows at removing metals from dwarf galaxies and the metal distribution of the gas is thus dependent on the strength of diffusion. By contrast, we show that the metallicities of stars are not strongly dependent on the strength of diffusion, provided that some diffusion is present.
Simulating Magnetized Laboratory Plasmas with Smoothed Particle Hydrodynamics
Johnson, Jeffrey N.
2009-01-01
The creation of plasmas in the laboratory continues to generate excitement in the physics community. Despite the best efforts of the intrepid plasma diagnostics community, the dynamics of these plasmas remains a difficult challenge to both the theorist and the experimentalist. This dissertation describes the simulation of strongly magnetized laboratory plasmas with Smoothed Particle Hydrodynamics (SPH), a method born of astrophysics but gaining broad support in the engineering community. We describe the mathematical formulation that best characterizes a strongly magnetized plasma under our circumstances of interest, and we review the SPH method and its application to astrophysical plasmas based on research by Phillips [1], Buerve [2], and Price and Monaghan [3]. Some modifications and extensions to this method are necessary to simulate terrestrial plasmas, such as a treatment of magnetic diffusion based on work by Brookshaw [4] and by Atluri [5]; we describe these changes as we turn our attention toward laboratory experiments. Test problems that verify the method are provided throughout the discussion. Finally, we apply our method to the compression of a magnetized plasma performed by the Compact Toroid Injection eXperiment (CTIX) [6] and show that the experimental results support our computed predictions.
Metal Diffusion in Smoothed Particle Hydrodynamics Simulations of Dwarf Galaxies
NASA Astrophysics Data System (ADS)
Williamson, David; Martel, Hugo; Kawata, Daisuke
2016-05-01
We perform a series of smoothed particle hydrodynamics simulations of isolated dwarf galaxies to compare different metal mixing models. In particular, we examine the role of diffusion in the production of enriched outflows and in determining the metallicity distributions of gas and stars. We investigate different diffusion strengths by changing the pre-factor of the diffusion coefficient, by varying how the diffusion coefficient is calculated from the local velocity distribution, and by varying whether the speed of sound is included as a velocity term. Stronger diffusion produces a tighter [O/Fe]-[Fe/H] distribution in the gas and cuts off the gas metallicity distribution function at lower metallicities. Diffusion suppresses the formation of low-metallicity stars, even with weak diffusion, and also strips metals from enriched outflows. This produces a remarkably tight correlation between “metal mass-loading” (mean metal outflow rate divided by mean metal production rate) and the strength of diffusion, even when the diffusion coefficient is calculated in different ways. The effectiveness of outflows at removing metals from dwarf galaxies and the metal distribution of the gas is thus dependent on the strength of diffusion. By contrast, we show that the metallicities of stars are not strongly dependent on the strength of diffusion, provided that some diffusion is present.
The small-scale turbulent dynamo in smoothed particle magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Tricco, T. S.; Price, D. J.; Federrath, C.
2016-05-01
Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small- scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of starforming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning method may be used to achieve ∇ • B = 0 to machine precision, albeit at significant computational expense.
SPHGal: smoothed particle hydrodynamics with improved accuracy for galaxy simulations
NASA Astrophysics Data System (ADS)
Hu, Chia-Yu; Naab, Thorsten; Walch, Stefanie; Moster, Benjamin P.; Oser, Ludwig
2014-09-01
We present the smoothed particle hydrodynamics (SPH) implementation SPHGal, which combines some recently proposed improvements in GADGET. This includes a pressure-entropy formulation with a Wendland kernel, a higher order estimate of velocity gradients, a modified artificial viscosity switch with a modified strong limiter, and artificial conduction of thermal energy. With a series of idealized hydrodynamic tests, we show that the pressure-entropy formulation is ideal for resolving fluid mixing at contact discontinuities but performs conspicuously worse at strong shocks due to the large entropy discontinuities. Including artificial conduction at shocks greatly improves the results. In simulations of Milky Way like disc galaxies a feedback-induced instability develops if too much artificial viscosity is introduced. Our modified artificial viscosity scheme prevents this instability and shows efficient shock capturing capability. We also investigate the star formation rate and the galactic outflow. The star formation rates vary slightly for different SPH schemes while the mass loading is sensitive to the SPH scheme and significantly reduced in our favoured implementation. We compare the accretion behaviour of the hot halo gas. The formation of cold blobs, an artefact of simple SPH implementations, can be eliminated efficiently with proper fluid mixing, either by conduction and/or by using a pressure-entropy formulation.
Smoothed Particle Hydrodynamics for water wave propagation in a channel
NASA Astrophysics Data System (ADS)
Omidvar, Pourya; Norouzi, Hossein; Zarghami, Ahad
2015-01-01
In this paper, Smoothed Particle Hydrodynamics (SPH) is used to simulate the propagation of waves in an intermediate depth water channel. The major advantage of using SPH is that no special treatment of the free surface is required, which is advantageous for simulating highly nonlinear flows with possible wave breaking. The SPH method has an option of different formulations with their own advantages and drawbacks to be implemented. Here, we apply the classical and Arbitrary Lagrange-Euler (ALE) formulation for wave propagation in a water channel. The classical SPH should come with an artificial viscosity which stabilizes the numerical algorithm and increases the accuracy. Here, we will show that the use of classical SPH with an artificial viscosity may cause the waves in the channel to decay. On the other hand, we will show that using the ALE-SPH algorithm with a Riemann solver is more stable, and in addition to producing the pressure fields with much less numerical noise, the waves propagate in the channel without dissipation.
The giant impact simulations with density independent smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Hosono, Natsuki; Saitoh, Takayuki R.; Makino, Junichiro; Genda, Hidenori; Ida, Shigeru
2016-06-01
At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity.
Smoothed aggregation adaptive spectral element-based algebraic multigrid
Energy Science and Technology Software Center (ESTSC)
2015-01-20
SAAMGE provides parallel methods for building multilevel hierarchies and solvers that can be used for elliptic equations with highly heterogeneous coefficients. Additionally, hierarchy adaptation is implemented allowing solving multiple problems with close coefficients without rebuilding the hierarchy.
AN ADAPTIVE PARTICLE-MESH GRAVITY SOLVER FOR ENZO
Passy, Jean-Claude; Bryan, Greg L.
2014-11-01
We describe and implement an adaptive particle-mesh algorithm to solve the Poisson equation for grid-based hydrodynamics codes with nested grids. The algorithm is implemented and extensively tested within the astrophysical code Enzo against the multigrid solver available by default. We find that while both algorithms show similar accuracy for smooth mass distributions, the adaptive particle-mesh algorithm is more accurate for the case of point masses, and is generally less noisy. We also demonstrate that the two-body problem can be solved accurately in a configuration with nested grids. In addition, we discuss the effect of subcycling, and demonstrate that evolving all the levels with the same timestep yields even greater precision.
Investigating the global collapse of filaments using smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Clarke, S. D.; Whitworth, A. P.
2015-05-01
We use smoothed particle hydrodynamic simulations of cold, uniform density, self-gravitating filaments, to investigate their longitudinal collapse time-scales; these time-scales are important because they determine the time available for a filament to fragment into cores. A filament is initially characterized by its line-mass, μO, its radius, RO (or equivalently its density ρ O= μ O/π RO^2), and its aspect ratio, AO (≡ZO/RO, where ZO is its half-length). The gas is only allowed to contract longitudinally, i.e. parallel to the symmetry axis of the filament (the z-axis). Pon et al. (2012) have considered the global dynamics of such filaments analytically. They conclude that short filaments (AO ≲ 5) collapse along the z-axis more-or-less homologously, on a time-scale tHOM ˜ 0.44 AO (GρO)-1/2; in contrast, longer filaments (AO ≳ 5) undergo end-dominated collapse, i.e. two dense clumps form at the ends of the filament and converge on the centre sweeping up mass as they go, on a time-scale t_{END} ˜ 0.98 AO^{1/2} (Gρ O)^{-1/2}. Our simulations do not corroborate these predictions. First, for all AO ≳ 2, the collapse time satisfies a single equation t_{COL}˜ (0.49+0.26AO)(Gρ O)^{-1/2}, which for large AO is much longer than the Pon et al. prediction. Secondly, for all AO ≳ 2, the collapse is end-dominated. Thirdly, before being swept up, the gas immediately ahead of an end-clump is actually accelerated outwards by the gravitational attraction of the approaching clump, resulting in a significant ram pressure. For high aspect ratio filaments, the end-clumps approach an asymptotic inward speed, due to the fact that they are doing work both accelerating and compressing the gas they sweep up. Pon et al. appear to have neglected the outward acceleration and its consequences.
Ham, Bumsub; Min, Dongbo; Sohn, Kwanghoon
2013-03-01
Anisotropic diffusion has been known to be closely related to adaptive smoothing and discretized in a similar manner. This paper revisits a fundamental relationship between two approaches. It is shown that adaptive smoothing and anisotropic diffusion have different theoretical backgrounds by exploring their characteristics with the perspective of normalization, evolution step size, and energy flow. Based on this principle, adaptive smoothing is derived from a second order partial differential equation (PDE), not a conventional anisotropic diffusion, via the coupling of Fick's law with a generalized continuity equation where a "source" or "sink" exists, which has not been extensively exploited. We show that the source or sink is closely related to the asymmetry of energy flow as well as the normalization term of adaptive smoothing. It enables us to analyze behaviors of adaptive smoothing, such as the maximum principle and stability with a perspective of a PDE. Ultimately, this relationship provides new insights into application-specific filtering algorithm design. By modeling the source or sink in the PDE, we introduce two specific diffusion filters, the robust anisotropic diffusion and the robust coherence enhancing diffusion, as novel instantiations which are more robust against the outliers than the conventional filters. PMID:23193236
NASA Astrophysics Data System (ADS)
Klapp, Jaime; di G Sigalotti, Leonardo; Troconis, Jorge; Sira, Eloy; Pena, Franklin; ININ-IVIC Team; Cinvestav-UAM-A Team
2014-11-01
We study numerically liquid-vapor phase separation in two-dimensional, nonisothermal, van der Waals (vdW) liquid drops using the method of Smoothed Particle Hydrodynamics (SPH). In contrast to previous SPH simulations of drop formation, our approach is fully adaptive and follows the diffuse interface model for a single-component fluid, where a reversible, capillary (Korteweg) force is added to the equations of motion to model the rapid but smooth transition of physical quantities through the interface separating the bulk phases. Surface tension arises naturally from the cohesive part of the vdW equation of state and the capillary forces. The drop models all start from a square-shaped liquid and spinodal decomposition is investigated for a range of initial densities and temperatures. The simulations predict the formation of stable, subcritical liquid drops with a vapor atmosphere, with the densities and temperatures of coexisting liquid and vapor in the vdW phase diagram closely matching the binodal curve. We find that the values of surface tension, as determined from the Young-Laplace equation, are in good agreement with the results of independent numerical simulations and experimental data. The models also predict the increase of the vapor pressure with temperature and the fitting to the numerical data reproduces very well the Clausius-Clapeyron relation, thus allowing for the calculation of the vaporization pressure for this vdW fluid. Cinvestav-Abacus.
NASA Astrophysics Data System (ADS)
Sigalotti, Leonardo Di G.; Troconis, Jorge; Sira, Eloy; Peña-Polo, Franklin; Klapp, Jaime
2014-07-01
We study numerically liquid-vapor phase separation in two-dimensional, nonisothermal, van der Waals (vdW) liquid drops using the method of smoothed particle hydrodynamics (SPH). In contrast to previous SPH simulations of drop formation, our approach is fully adaptive and follows the diffuse-interface model for a single-component fluid, where a reversible, capillary (Korteweg) force is added to the equations of motion to model the rapid but smooth transition of physical quantities through the interface separating the bulk phases. Surface tension arises naturally from the cohesive part of the vdW equation of state and the capillary forces. The drop models all start from a square-shaped liquid and spinodal decomposition is investigated for a range of initial densities and temperatures. The simulations predict the formation of stable, subcritical liquid drops with a vapor atmosphere, with the densities and temperatures of coexisting liquid and vapor in the vdW phase diagram closely matching the binodal curve. We find that the values of surface tension, as determined from the Young-Laplace equation, are in good agreement with the results of independent numerical simulations and experimental data. The models also predict the increase of the vapor pressure with temperature and the fitting to the numerical data reproduces very well the Clausius-Clapeyron relation, thus allowing for the calculation of the vaporization pressure for this vdW fluid.
Sigalotti, Leonardo Di G; Troconis, Jorge; Sira, Eloy; Peña-Polo, Franklin; Klapp, Jaime
2014-07-01
We study numerically liquid-vapor phase separation in two-dimensional, nonisothermal, van der Waals (vdW) liquid drops using the method of smoothed particle hydrodynamics (SPH). In contrast to previous SPH simulations of drop formation, our approach is fully adaptive and follows the diffuse-interface model for a single-component fluid, where a reversible, capillary (Korteweg) force is added to the equations of motion to model the rapid but smooth transition of physical quantities through the interface separating the bulk phases. Surface tension arises naturally from the cohesive part of the vdW equation of state and the capillary forces. The drop models all start from a square-shaped liquid and spinodal decomposition is investigated for a range of initial densities and temperatures. The simulations predict the formation of stable, subcritical liquid drops with a vapor atmosphere, with the densities and temperatures of coexisting liquid and vapor in the vdW phase diagram closely matching the binodal curve. We find that the values of surface tension, as determined from the Young-Laplace equation, are in good agreement with the results of independent numerical simulations and experimental data. The models also predict the increase of the vapor pressure with temperature and the fitting to the numerical data reproduces very well the Clausius-Clapeyron relation, thus allowing for the calculation of the vaporization pressure for this vdW fluid. PMID:25122383
Unit Commitment by Adaptive Particle Swarm Optimization
NASA Astrophysics Data System (ADS)
Saber, Ahmed Yousuf; Senjyu, Tomonobu; Miyagi, Tsukasa; Urasaki, Naomitsu; Funabashi, Toshihisa
This paper presents an Adaptive Particle Swarm Optimization (APSO) for Unit Commitment (UC) problem. APSO reliably and accurately tracks a continuously changing solution. By analyzing the social model of standard PSO for the UC problem of variable size and load demand, adaptive criteria are applied on PSO parameters and the global best particle (knowledge) based on the diversity of fitness. In this proposed method, PSO parameters are automatically adjusted using Gaussian modification. To increase the knowledge, the global best particle is updated instead of a fixed one in each generation. To avoid the method to be frozen, idle particles are reset. The real velocity is digitized (0/1) by a logistic function for binary UC. Finally, the benchmark data and methods are used to show the effectiveness of the proposed method.
A new adaptive exponential smoothing method for non-stationary time series with level shifts
NASA Astrophysics Data System (ADS)
Monfared, Mohammad Ali Saniee; Ghandali, Razieh; Esmaeili, Maryam
2014-07-01
Simple exponential smoothing (SES) methods are the most commonly used methods in forecasting and time series analysis. However, they are generally insensitive to non-stationary structural events such as level shifts, ramp shifts, and spikes or impulses. Similar to that of outliers in stationary time series, these non-stationary events will lead to increased level of errors in the forecasting process. This paper generalizes the SES method into a new adaptive method called revised simple exponential smoothing (RSES), as an alternative method to recognize non-stationary level shifts in the time series. We show that the new method improves the accuracy of the forecasting process. This is done by controlling the number of observations and the smoothing parameter in an adaptive approach, and in accordance with the laws of statistical control limits and the Bayes rule of conditioning. We use a numerical example to show how the new RSES method outperforms its traditional counterpart, SES.
Molelcular dynamics, smooth particle applied mechanics and Clausius` inequality
Hoover, W.G.
1996-03-01
Recent developments in molecular dynamics furnish new interconnections among three classical fields: particle mechanics, continuum mechanics, and thermodynamics. The resulting links clarify the importance of Lyapunov instability to irreversibility.
High energy scattering of Dirac particles on smooth potentials
NASA Astrophysics Data System (ADS)
Han, Nguyen Suan; Dung, Le Anh; Xuan, Nguyen Nhu; Thang, Vu Toan
2016-08-01
The derivation of the Glauber type representation for the high energy scattering amplitude of particles of spin 1/2 is given within the framework of the Dirac equation in the Foldy-Wouthuysen (FW) representation and two-component formalism. The differential cross-sections on the Yukawa and Gaussian potentials are also considered and discussed.
Bifurcating Particle Swarms in Smooth-Walled Fractures
NASA Astrophysics Data System (ADS)
Pyrak-Nolte, L. J.; Sun, H.
2010-12-01
Particle swarms can occur naturally or from industrial processes where small liquid drops containing thousands to millions of micron-size to colloidal-size particles are released over time from seepage or leaks into fractured rock. The behavior of these particle swarms as they fall under gravity are affected by particle interactions as well as interactions with the walls of the fractures. In this paper, we present experimental results on the effect of fractures on the cohesiveness of the swarm and the formation of bifurcation structures as they fall under gravity and interact with the fracture walls. A transparent cubic sample (100 mm x 100 mm x 100 mm) containing a synthetic fracture with uniform aperture distributions was optically imaged to quantify the effect of confinement within fractures on particle swarm formation, swarm velocity, and swarm geometry. A fracture with a uniform aperture distribution was fabricated from two polished rectangular prisms of acrylic. A series of experiments were performed to determine how swarm movement and geometry are affected as the walls of the fracture are brought closer together from 50 mm to 1 mm. During the experiments, the fracture was fully saturated with water. We created the swarms using two different particle sizes in dilute suspension (~ 1.0% by mass). The particles were 3 micron diameter fluorescent polymer beads and 25 micron diameter soda-lime glass beads. Experiments were performed using swarms that ranged in size from 5 µl to 60 µl. The swarm behavior was imaged using an optical fluorescent imaging system composed of a CCD camera illuminated by a 100 mW diode-pumped doubled YAG laser. As a swarm falls in an open-tank of water, it forms a torroidal shape that is stable as long as no ambient or background currents exist in the water tank. When a swarm is released into a fracture with an aperture less than 5 mm, the swarm forms the torroidal shape but it is distorted because of the presence of the walls. The
ADVANCING THE ION BEAM THIN FILM PLANARIZATION PROCESS FOR THE SMOOTHING OF SUBSTRATE PARTICLES
Mirkarimi, P B; Spiller, E; Baker, S L; Robinson, J C; Stearns, D G; Liddle, J A; Salmassi, F; Liang, T; Stivers, A R
2004-10-19
For a number of technologies small substrate contaminants are undesirable, and for one technology in particular, extreme ultraviolet lithography (EUVL), they can be a very serious issue. We have demonstrated that the Ion Beam Thin Film Planarization Process, a coating process designed to planarize substrate asperities, can be extended to smooth {approx}70 nm and {approx}80 nm diameter particles on EUVL reticle substrates to a height of {approx}0.5 nm, which will render them noncritical in an EUVL printing process. We demonstrate this smoothing process using controlled nanoscale substrate particles and lines fabricated with an e-beam lithography process. The above smoothing process was also modified to yield an excellent reflectance/wavelength uniformity and a good EUV reflectivity for the multilayer, which is required for EUVL reticles. Cross-sectional TEM on a smoothed substrate line defect shows excellent agreement with results obtained from our multilayer growth model.
3D Smoothed Particle Hydrodynamics Models of Betelgeuse's Bow Shock
NASA Astrophysics Data System (ADS)
Mohamed, S.; Mackey, J.; Langer, N.
2013-05-01
Betelgeuse, the bright red supergiant (RSG) in Orion, is a runaway star. Its supersonic motion through the interstellar medium has resulted in the formation of a bow shock, a cometary structure pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. We show that the bow shock morphology depends substantially on the growth timescale for Rayleigh-Taylor versus Kelvin-Helmholtz instabilities. We discuss our models in light of the recent Herschel, GALEX and VLA observations. If the mass in the bow shock shell is low (~few × 10-3 M⊙), as seems to be implied by the AKARI and Herschel observations, then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular, smooth bow shock shell is consistent with this conclusion. We further discuss the implications of our results, in particular, the possibility that Betelgeuse may have only recently entered the RSG phase.
Kelvin-Helmholtz instabilities with Godunov smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Cha, Seung-Hoon; Inutsuka, Shu-Ichiro; Nayakshin, Sergei
2010-04-01
Numerical simulations for the non-linear development of Kelvin-Helmholtz instability in two different density layers have been performed with the particle-based method (Godunov SPH) developed by Inutsuka. The Godunov SPH can describe the Kelvin-Helmholtz instability even with a high-density contrast, while the standard SPH shows the absence of the instability across a density gradient. The interaction of a dense blob with a hot ambient medium has been performed also. The Godunov SPH describes the formation and evolution of the fingers due to the combinations of Rayleigh-Taylor, Richtmyer-Meshkov and Kelvin-Helmholtz instabilities. The blob test result coincides well with the results of the grid-based codes. An inaccurate handling of a density gradient in the standard SPH has been pointed out as the direct reason of the absence of the instabilities. An unphysical force happens at the density gradient even in a pressure equilibrium, and repulses particles from the initial density discontinuity. Therefore, the initial perturbation damps, and a gap form at the discontinuity. The unphysical force has been studied in terms of the consistency of a numerical scheme. Contrary to the standard SPH, the momentum equation of the Godunov SPH does not use the particle approximation, and has been derived from the kernel convolution or a new Lagrangian function. The new Lagrangian function used in the Godunov SPH is more analogous to the real Lagrangian function for continuum. The momentum equation of the Godunov SPH has much better linear consistency, so the unphysical force is greatly reduced compared to the standard SPH in a high density contrast.
Volume conservation issues in incompressible smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Nair, Prapanch; Tomar, Gaurav
2015-09-01
A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates.
NASA Astrophysics Data System (ADS)
Esmaili Sikarudi, M. A.; Nikseresht, A. H.
2016-01-01
Smoothed particle hydrodynamics is a robust Lagrangian particle method which is widely used in various applications, from astrophysics to hydrodynamics and heat conduction. It has intrinsic capabilities for simulating large deformation, composites, multiphysics events, and multiphase fluid flows. It is vital to use reliable boundary conditions when boundary value problems like heat conduction or Poisson equation for incompressible flows are solved. Since smoothed particle hydrodynamics is not a boundary fitted grids method, implementation of boundary conditions can be problematic. Many methods have been proposed for enhancing the accuracy of implementation of boundary conditions. In the present study a new approach for facilitating the implementation of Robin and Neumann boundary conditions is proposed and proven to give accurate results. Also there is no need to use complicated preprocessing as in virtual particle method. The new method is compared to an equivalent one dimensional moving least square scheme and it is shown that the present method is less sensitive to particle disorder.
Kapilevich, L V; Zaĭtseva, T N; Nosarev, A V; D'iakova, E Iu; Petlina, Z R; Ogorodova, L M; Ageev, B G; Magaeva, A A; Itin, V I; Terekhova, O G; Medvedev, M A
2012-02-01
Contractile responses of airways segments of porpoises inhaling nanopowder CoFe2O4 were stidued by means of a mechanographic method. Inhalation of the nanosize particles of CoFe2O4 in vivo and in vitro testing the nanomaterial on isolated smooth muscles led to potentiation histaminergic, cholinergic contractile activity in airways of porpoises and to strengthening of adrenergic relaxing answers. Nanosize particles vary amplitude of hyperpotassium reductions in smooth muscle segments of airways similarly to the effect of depolymerizing drug colchicine. PMID:22650066
Deformation of Soft Tissue and Force Feedback Using the Smoothed Particle Hydrodynamics
Liu, Xuemei; Wang, Ruiyi; Li, Yunhua; Song, Dongdong
2015-01-01
We study the deformation and haptic feedback of soft tissue in virtual surgery based on a liver model by using a force feedback device named PHANTOM OMNI developed by SensAble Company in USA. Although a significant amount of research efforts have been dedicated to simulating the behaviors of soft tissue and implementing force feedback, it is still a challenging problem. This paper introduces a kind of meshfree method for deformation simulation of soft tissue and force computation based on viscoelastic mechanical model and smoothed particle hydrodynamics (SPH). Firstly, viscoelastic model can present the mechanical characteristics of soft tissue which greatly promotes the realism. Secondly, SPH has features of meshless technique and self-adaption, which supply higher precision than methods based on meshes for force feedback computation. Finally, a SPH method based on dynamic interaction area is proposed to improve the real time performance of simulation. The results reveal that SPH methodology is suitable for simulating soft tissue deformation and force feedback calculation, and SPH based on dynamic local interaction area has a higher computational efficiency significantly compared with usual SPH. Our algorithm has a bright prospect in the area of virtual surgery. PMID:26417380
Deformation of Soft Tissue and Force Feedback Using the Smoothed Particle Hydrodynamics.
Liu, Xuemei; Wang, Ruiyi; Li, Yunhua; Song, Dongdong
2015-01-01
We study the deformation and haptic feedback of soft tissue in virtual surgery based on a liver model by using a force feedback device named PHANTOM OMNI developed by SensAble Company in USA. Although a significant amount of research efforts have been dedicated to simulating the behaviors of soft tissue and implementing force feedback, it is still a challenging problem. This paper introduces a kind of meshfree method for deformation simulation of soft tissue and force computation based on viscoelastic mechanical model and smoothed particle hydrodynamics (SPH). Firstly, viscoelastic model can present the mechanical characteristics of soft tissue which greatly promotes the realism. Secondly, SPH has features of meshless technique and self-adaption, which supply higher precision than methods based on meshes for force feedback computation. Finally, a SPH method based on dynamic interaction area is proposed to improve the real time performance of simulation. The results reveal that SPH methodology is suitable for simulating soft tissue deformation and force feedback calculation, and SPH based on dynamic local interaction area has a higher computational efficiency significantly compared with usual SPH. Our algorithm has a bright prospect in the area of virtual surgery. PMID:26417380
Smooth particle hydrodynamics: theory and application to the origin of the moon
Benz, W.
1986-01-01
The origin of the moon is modeled by the so-called smooth particle hydrodynamics (SPH) method (Lucy, 1977, Monaghan 1985) which substitutes to the fluid a finite set of extended particles, the hydrodynamics equations reduce to the equation of motion of individual particles. These equations of motion differ only from the standard gravitational N-body problem insofar that pressure gradients and viscosity terms have to be added to the gradient of the potential to derive the forces between the particles. The numerical tools developed for ''classical'' N-body problems can therefore be readily applied to solve 3 dimensional hydroynamical problems. 12 refs., 1 fig.
Stability of bumps in piecewise smooth neural fields with nonlinear adaptation
NASA Astrophysics Data System (ADS)
Kilpatrick, Zachary P.; Bressloff, Paul C.
2010-06-01
We study the linear stability of stationary bumps in piecewise smooth neural fields with local negative feedback in the form of synaptic depression or spike frequency adaptation. The continuum dynamics is described in terms of a nonlocal integrodifferential equation, in which the integral kernel represents the spatial distribution of synaptic weights between populations of neurons whose mean firing rate is taken to be a Heaviside function of local activity. Discontinuities in the adaptation variable associated with a bump solution means that bump stability cannot be analyzed by constructing the Evans function for a network with a sigmoidal gain function and then taking the high-gain limit. In the case of synaptic depression, we show that linear stability can be formulated in terms of solutions to a system of pseudo-linear equations. We thus establish that sufficiently strong synaptic depression can destabilize a bump that is stable in the absence of depression. These instabilities are dominated by shift perturbations that evolve into traveling pulses. In the case of spike frequency adaptation, we show that for a wide class of perturbations the activity and adaptation variables decouple in the linear regime, thus allowing us to explicitly determine stability in terms of the spectrum of a smooth linear operator. We find that bumps are always unstable with respect to this class of perturbations, and destabilization of a bump can result in either a traveling pulse or a spatially localized breather.
NASA Astrophysics Data System (ADS)
Nassauer, Benjamin; Liedke, Thomas; Kuna, Meinhard
2016-03-01
In the present paper, the direct coupling of a discrete element method (DEM) with polyhedral particles and smoothed particle hydrodynamics (SPH) is presented. The two simulation techniques are fully coupled in both ways through interaction forces between the solid DEM particles and the fluid SPH particles. Thus this simulation method provides the possibility to simulate the individual movement of polyhedral, sharp-edged particles as well as the flow field around these particles in fluid-saturated granular matter which occurs in many technical processes e.g. wire sawing, grinding or lapping. The coupled method is exemplified and validated by the simulation of a particle in a shear flow, which shows good agreement with analytical solutions.
A Fast Variational Method for the Construction of Resolution Adaptive C-Smooth Molecular Surfaces.
Bajaj, Chandrajit L; Xu, Guoliang; Zhang, Qin
2009-05-01
We present a variational approach to smooth molecular (proteins, nucleic acids) surface constructions, starting from atomic coordinates, as available from the protein and nucleic-acid data banks. Molecular dynamics (MD) simulations traditionally used in understanding protein and nucleic-acid folding processes, are based on molecular force fields, and require smooth models of these molecular surfaces. To accelerate MD simulations, a popular methodology is to employ coarse grained molecular models, which represent clusters of atoms with similar physical properties by psuedo- atoms, resulting in coarser resolution molecular surfaces. We consider generation of these mixed-resolution or adaptive molecular surfaces. Our approach starts from deriving a general form second order geometric partial differential equation in the level-set formulation, by minimizing a first order energy functional which additionally includes a regularization term to minimize the occurrence of chemically infeasible molecular surface pockets or tunnel-like artifacts. To achieve even higher computational efficiency, a fast cubic B-spline C(2) interpolation algorithm is also utilized. A narrow band, tri-cubic B-spline level-set method is then used to provide C(2) smooth and resolution adaptive molecular surfaces. PMID:19802355
Ijpma, Gijs; Al-Jumaily, Ahmed M; Cairns, Simeon P; Sieck, Gary C
2011-09-01
Length adaptation in airway smooth muscle (ASM) is attributed to reorganization of the cytoskeleton, and in particular the contractile elements. However, a constantly changing lung volume with tidal breathing (hence changing ASM length) is likely to restrict full adaptation of ASM for force generation. There is likely to be continuous length adaptation of ASM between states of incomplete or partial length adaption. We propose a new model that assimilates findings on myosin filament polymerization/depolymerization, partial length adaptation, isometric force, and shortening velocity to describe this continuous length adaptation process. In this model, the ASM adapts to an optimal force-generating capacity in a repeating cycle of events. Initially the myosin filament, shortened by prior length changes, associates with two longer actin filaments. The actin filaments are located adjacent to the myosin filaments, such that all myosin heads overlap with actin to permit maximal cross-bridge cycling. Since in this model the actin filaments are usually longer than myosin filaments, the excess length of the actin filament is located randomly with respect to the myosin filament. Once activated, the myosin filament elongates by polymerization along the actin filaments, with the growth limited by the overlap of the actin filaments. During relaxation, the myosin filaments dissociate from the actin filaments, and then the cycle repeats. This process causes a gradual adaptation of force and instantaneous adaptation of shortening velocity. Good agreement is found between model simulations and the experimental data depicting the relationship between force development, myosin filament density, or shortening velocity and length. PMID:21659490
A relativistic smoothed particle hydrodynamics method tested with the shock tube
NASA Astrophysics Data System (ADS)
Mann, Patrick J.
1991-12-01
The smoothed particle hydrodynamics method is applied to an ADM 3 + 1 formulation of the equations for relativistic fluid flow. In particular the one-dimensional shock tube is addressed. Three codes are described. The first is a straightforward extension of classic SPH, while the other two are modifications which allow for time-dependent smoothing lengths. The first of these modifications approximates the internal energy density, while the second approximates the total energy density. Two smoothing forms are tested: an artificial viscosity and the direct method of A.J. Baker [Finite Element Computation Fluid Mechanics (Hemisphere, New York, 1983)]. The results indicate that the classic SPH code with particle-particle based artificial viscosity is reasonably accurate and very consistent. It gives quite sharp edges and flat plateaus, but the velocity plateau is significantly overestimated, and an oscillation can appear in the rarefaction wave. The modified versions with Baker smoothing procedure better results for moderate initial conditions, but begin to show spikes when the initial density jump is large. Generally the results are comparable to simple finite element and finite difference methods.
Modeling of liquid-vapor phase change using smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Das, A. K.; Das, P. K.
2015-12-01
A model has been proposed based on smoothed particle hydrodynamics to describe gas liquid phase change. Pseudo particles of zero mass are initially placed to locate the interface. Mass generated due to phase change is assigned to the pseudo particles and their positions are updated at intervals to track the mobility of the interface. The developed algorithm has been used to simulate vapor formation around solid spheres both in the absence of gravity and in the normal gravitational field. Finally, bubble growth over a hot horizontal surface due to boiling has been simulated. Simulated results showed good matching with the reported literature.
Sigalotti, Leonardo Di G; Troconis, Jorge; Sira, Eloy; Peña-Polo, Franklin; Klapp, Jaime
2015-07-01
The rapid evaporation and explosive boiling of a van der Waals (vdW) liquid drop in microgravity is simulated numerically in two-space dimensions using the method of smoothed particle hydrodynamics. The numerical approach is fully adaptive and incorporates the effects of surface tension, latent heat, mass transfer across the interface, and liquid-vapor interface dynamics. Thermocapillary forces are modeled by coupling the hydrodynamics to a diffuse-interface description of the liquid-vapor interface. The models start from a nonequilibrium square-shaped liquid of varying density and temperature. For a fixed density, the drop temperature is increased gradually to predict the point separating normal boiling at subcritical heating from explosive boiling at the superheat limit for this vdW fluid. At subcritical heating, spontaneous evaporation produces stable drops floating in a vapor atmosphere, while at near-critical heating, a bubble is nucleated inside the drop, which then collapses upon itself, leaving a smaller equilibrated drop embedded in its own vapor. At the superheat limit, unstable bubble growth leads to either fragmentation or violent disruption of the liquid layer into small secondary drops, depending on the liquid density. At higher superheats, explosive boiling occurs for all densities. The experimentally observed wrinkling of the bubble surface driven by rapid evaporation followed by a Rayleigh-Taylor instability of the thin liquid layer and the linear growth of the bubble radius with time are reproduced by the simulations. The predicted superheat limit (T(s)≈0.96) is close to the theoretically derived value of T(s)=1 at zero ambient pressure for this vdW fluid. PMID:26274283
NASA Astrophysics Data System (ADS)
Sigalotti, Leonardo Di G.; Troconis, Jorge; Sira, Eloy; Peña-Polo, Franklin; Klapp, Jaime
2015-07-01
The rapid evaporation and explosive boiling of a van der Waals (vdW) liquid drop in microgravity is simulated numerically in two-space dimensions using the method of smoothed particle hydrodynamics. The numerical approach is fully adaptive and incorporates the effects of surface tension, latent heat, mass transfer across the interface, and liquid-vapor interface dynamics. Thermocapillary forces are modeled by coupling the hydrodynamics to a diffuse-interface description of the liquid-vapor interface. The models start from a nonequilibrium square-shaped liquid of varying density and temperature. For a fixed density, the drop temperature is increased gradually to predict the point separating normal boiling at subcritical heating from explosive boiling at the superheat limit for this vdW fluid. At subcritical heating, spontaneous evaporation produces stable drops floating in a vapor atmosphere, while at near-critical heating, a bubble is nucleated inside the drop, which then collapses upon itself, leaving a smaller equilibrated drop embedded in its own vapor. At the superheat limit, unstable bubble growth leads to either fragmentation or violent disruption of the liquid layer into small secondary drops, depending on the liquid density. At higher superheats, explosive boiling occurs for all densities. The experimentally observed wrinkling of the bubble surface driven by rapid evaporation followed by a Rayleigh-Taylor instability of the thin liquid layer and the linear growth of the bubble radius with time are reproduced by the simulations. The predicted superheat limit (Ts≈0.96 ) is close to the theoretically derived value of Ts=1 at zero ambient pressure for this vdW fluid.
On the dynamics of nonlinear, unsteady landslide flow within the smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Khvostova, O.; Averbukh, E.
2012-04-01
In the present study the idea of landslide modeling by particle method is described. Smoothed particle hydrodynamics was invented in 1977 by Leon Lucy and independently by Bob Gingold and Joe Monaghan [1]. It was used for astrophysics phenomena's simulation. Later it was adapted for hydrodynamics, gas dynamics and solid body problems. Landslides can be caused by the influence of different factors. Landslides occur when the angle of inclination of the slope of the slope or if the slope is burdened with loose material. A landslide flow is a thin homogeneous layer of nearly incompressible fluid. It is considered that at the initial moment shifted part of a ground mass is splitting and turning into liquid of several layers which then is streaming down along the slope. The landslide flow motion is described with the Navie-Stocks set of equations: D→u-= - 1\\upsidedownBigTriangle P + μ \\upsidedownBigTriangle →u + g Dt ρ (1) D-ρ = 0, Dt (2) where u is velocity vector, t is time, ρ is a flow density, P is a pressure, μ is a viscosity coefficient, g is gravity. Continuum discretization by finite number of lagrangian particles is the main idea of SPH [2,3]. Particles moves with the flow and arbitrary connectivity is allowed. Therefore, SPH does not need a grid to calculate spatial derivatives. For any field A(r), involved in equation (1), e.g. pressure, density, viscosity etc., we consider an approximation with a finite function: A(r) = ∫ω A (r')W (r- r',h)dr' (3) where A is a desired field, r is a radius-vector, W is an interpolating kernel. The free boundary condition problem is discussed. Finding the particles on a free surface is described. Also the surface tension force defining is shown. Described method is implemented and mathematical modeling of landslide flows motion along slope is simulated. Different types of slopes are considered: with constant and variable steepness, long and wide. Wave-breaking effects near the wall are shown. Findings are analyzed
Smoothed particle hydrodynamics simulation of shear-induced powder migration in injection moulding.
Kauzlarić, David; Pastewka, Lars; Meyer, Hagen; Heldele, Richard; Schulz, Michael; Weber, Oxana; Piotter, Volker; Hausselt, Jürgen; Greiner, Andreas; Korvink, Jan G
2011-06-13
We present the application of the smoothed particle hydrodynamics (SPH) discretization scheme to Phillips' model for shear-induced particle migration in concentrated suspensions. This model provides an evolution equation for the scalar mean volume fraction of idealized spherical solid particles of equal diameter which is discretized by the SPH formalism. In order to obtain a discrete evolution equation with exact conservation properties we treat in fact the occupied volume of the solid particles as the degree of freedom for the fluid particles. We present simulation results in two- and three-dimensional channel flow. The two-dimensional results serve as a verification by a comparison to analytic solutions. The three-dimensional results are used for a comparison with experimental measurements obtained from computer tomography of injection moulded ceramic microparts. We observe the best agreement of measurements with snapshots of the transient simulation for a ratio D(c)/D(η)=0.1 of the two model parameters. PMID:21536579
Cen, Guanjun; Zeng, Xianru; Long, Xiuzhen; Wei, Dewei; Gao, Xuyuan; Zeng, Tao
2015-01-01
In insects, the frequency distribution of the measurements of sclerotized body parts is generally used to classify larval instars and is characterized by a multimodal overlap between instar stages. Nonparametric methods with fixed bandwidths, such as histograms, have significant limitations when used to fit this type of distribution, making it difficult to identify divisions between instars. Fixed bandwidths have also been chosen somewhat subjectively in the past, which is another problem. In this study, we describe an adaptive kernel smoothing method to differentiate instars based on discontinuities in the growth rates of sclerotized insect body parts. From Brooks’ rule, we derived a new standard for assessing the quality of instar classification and a bandwidth selector that more accurately reflects the distributed character of specific variables. We used this method to classify the larvae of Austrosimulium tillyardianum (Diptera: Simuliidae) based on five different measurements. Based on head capsule width and head capsule length, the larvae were separated into nine instars. Based on head capsule postoccipital width and mandible length, the larvae were separated into 8 instars and 10 instars, respectively. No reasonable solution was found for antennal segment 3 length. Separation of the larvae into nine instars using head capsule width or head capsule length was most robust and agreed with Crosby’s growth rule. By strengthening the distributed character of the separation variable through the use of variable bandwidths, the adaptive kernel smoothing method could identify divisions between instars more effectively and accurately than previous methods. PMID:26546689
NASA Astrophysics Data System (ADS)
Danilewicz, Andrzej; Sikora, Zbigniew
2015-02-01
A theoretical base of SPH method, including the governing equations, discussion of importance of the smoothing function length, contact formulation, boundary treatment and finally utilization in hydrocode simulations are presented. An application of SPH to a real case of large penetrations (crater creating) into the soil caused by falling mass in Dynamic Replacement Method is discussed. An influence of particles spacing on method accuracy is presented. An example calculated by LS-DYNA software is discussed. Chronological development of Smooth Particle Hydrodynamics is presented. Theoretical basics of SPH method stability and consistency in SPH formulation, artificial viscosity and boundary treatment are discussed. Time integration techniques with stability conditions, SPH+FEM coupling, constitutive equation and equation of state (EOS) are presented as well.
NASA Astrophysics Data System (ADS)
Tartakovsky, Alexandre M.; Panchenko, Alexander
2016-01-01
We present a novel formulation of the Pairwise Force Smoothed Particle Hydrodynamics (PF-SPH) model and use it to simulate two- and three-phase flows in bounded domains. In the PF-SPH model, the Navier-Stokes equations are discretized with the Smoothed Particle Hydrodynamics (SPH) method, and the Young-Laplace boundary condition at the fluid-fluid interface and the Young boundary condition at the fluid-fluid-solid interface are replaced with pairwise forces added into the Navier-Stokes equations. We derive a relationship between the parameters in the pairwise forces and the surface tension and static contact angle. Next, we demonstrate the model's accuracy under static and dynamic conditions. Finally, we use the Pf-SPH model to simulate three phase flow in a porous medium.
NASA Astrophysics Data System (ADS)
Sandalski, Stou
Smooth particle hydrodynamics is an efficient method for modeling the dynamics of fluids. It is commonly used to simulate astrophysical processes such as binary mergers. We present a newly developed GPU accelerated smooth particle hydrodynamics code for astrophysical simulations. The code is named
Numerical Studies on the Explosive Welding by Smoothed Particle Hydrodynamics (SPH)
NASA Astrophysics Data System (ADS)
Tanaka, Katsumi
2007-06-01
A particular characteristic of an explosively produced weld is that the profile of the weld interface often has a regular wavy appearance. An effect of detached shock wave and jetting on the metal interface of explosive welding has been shown by SPH (Smoothed particle hydrodynamics). Numerical results show wavy interface which is observed in several experiments. High speed jet between interface and Karman vortex after oblique impact of a flyer plate to a parent plate were major mechanism of explosive welding.
MODELING OF FRICTION STIR WELDING (FSW) PROCESS USING SMOOTH PARTICLE HYDRODYNAMICS (SPH)
Tartakovsky, Alexandre M.; Grant, Glenn J.; Sun, Xin; Khaleel, Mohammad A.
2006-06-14
We present a novel modeling approach to simulate FSW process that may have significant advantages over current traditional finite element or finite difference based methods. The proposed model is based on Smoothed Particle Hydrodynamics (SPH) method, a fully Lagrangian particle method that can simulate the dynamics of interfaces, large material deformations, void formations and material's strain and temperature history without employing complex tracking schemes. Two- and three-dimensional simulations for different tool designs are presented. Preliminary numerical results are in good qualitative agreement with experimental observations.
Crack Propagation in Bi-Material System via Pseudo-Spring Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Chakraborty, Sukanta; Shaw, Amit
2014-05-01
A Smoothed Particles Hydrodynamics (SPH) based framework with material constitutive model is developed to simulate crack initiation and propagation in a bi-material system. An efficient immediate neighbor interaction is formulated by connecting neighbors through pseudo-springs. A damage evolution law defines degradation of the inter-neighbor spring forces and corresponding reduced interaction is introduced in mass, momentum, and energy-conserving particle collocation. The proposed technique is validated through a simple test on a pre-notched bi-material system producing a conformal crack path.
Cleary, Paul W; Prakash, Mahesh
2004-09-15
Particle-based simulation methods, such as the discrete-element method and smoothed particle hydrodynamics, have specific advantages in modelling complex three-dimensional (3D) environmental fluid and particulate flows. The theory of both these methods and their relative advantages compared with traditional methods will be discussed. Examples of 3D flows on realistic topography illustrate the environmental application of these methods. These include the flooding of a river valley as a result of a dam collapse, coastal inundation by a tsunami, volcanic lava flow and landslides. Issues related to validation and quality data availability are also discussed. PMID:15306427
Multi-resolution flow simulations by smoothed particle hydrodynamics via domain decomposition
NASA Astrophysics Data System (ADS)
Bian, Xin; Li, Zhen; Tang, Yu-Hang; Karniadakis, George
2015-11-01
We present a methodology to concurrently couple particle-based methods via a domain decomposition (DD) technique for simulating viscous flows. In particular, we select two resolutions of the smoothed particle hydrodynamics (SPH) method as demonstration. Within the DD framework, a simulation domain is decomposed into two (or more) overlapping sub-domains, each of which has an individual particle scale determined by the local flow physics. Consistency of the two sub-domains is achieved in the overlap region by matching the two independent simulations based on Lagrangian interpolation of state variables and fluxes. The domain decomposition based SPH method (DD-SPH) employs different spatial and temporal resolutions, and hence, each sub-domain has its own smoothing length and time step. As a consequence, particle refinement and de-refinement are performed asynchronously according to individual time advancement of each sub-domain. The proposed strategy avoids SPH force interactions between different resolutions on purpose, so that coupling, in principle, can go beyond SPH - SPH, and may allow SPH to be coupled with other mesoscopic or microscopic particle methods. The DD-SPH method is validated first for a transient Couette flow, where simulation results base. US DOE Collaboratory on Mathematics for Mesoscopic Modeling of Materials (CM4).
Hendrickson, B.; Plimpton, S.; Attaway, S.; Swegle, J.
1996-09-01
Transient dynamics simulations are commonly used to model phenomena such as car crashes, underwater explosions, and the response of shipping containers to high-speed impacts. Physical objects in such a simulation are typically represented by Lagrangian meshes because the meshes can move and deform with the objects as they undergo stress. Fluids (gasoline, water) or fluid-like materials (earth) in the simulation can be modeled using the techniques of smoothed particle hydrodynamics. Implementing a hybrid mesh/particle model on a massively parallel computer poses several difficult challenges. One challenge is to simultaneously parallelize and load-balance both the mesh and particle portions of the computation. A second challenge is to efficiently detect the contacts that occur within the deforming mesh and between mesh elements and particles as the simulation proceeds. These contacts impart forces to the mesh elements and particles which must be computed at each timestep to accurately capture the physics of interest. In this paper we describe new parallel algorithms for smoothed particle hydrodynamics and contact detection which turn out to have several key features in common. Additionally, we describe how to join the new algorithms with traditional parallel finite element techniques to create an integrated particle/mesh transient dynamics simulation. Our approach to this problem differs from previous work in that we use three different parallel decompositions, a static one for the finite element analysis and dynamic ones for particles and for contact detection. We have implemented our ideas in a parallel version of the transient dynamics code PRONTO-3D and present results for the code running on a large Intel Paragon.
Xiao, Zhu; Havyarimana, Vincent; Li, Tong; Wang, Dong
2016-01-01
In this paper, a novel nonlinear framework of smoothing method, non-Gaussian delayed particle smoother (nGDPS), is proposed, which enables vehicle state estimation (VSE) with high accuracy taking into account the non-Gaussianity of the measurement and process noises. Within the proposed method, the multivariate Student’s t-distribution is adopted in order to compute the probability distribution function (PDF) related to the process and measurement noises, which are assumed to be non-Gaussian distributed. A computation approach based on Ensemble Kalman Filter (EnKF) is designed to cope with the mean and the covariance matrix of the proposal non-Gaussian distribution. A delayed Gibbs sampling algorithm, which incorporates smoothing of the sampled trajectories over a fixed-delay, is proposed to deal with the sample degeneracy of particles. The performance is investigated based on the real-world data, which is collected by low-cost on-board vehicle sensors. The comparison study based on the real-world experiments and the statistical analysis demonstrates that the proposed nGDPS has significant improvement on the vehicle state accuracy and outperforms the existing filtering and smoothing methods. PMID:27187405
Xiao, Zhu; Havyarimana, Vincent; Li, Tong; Wang, Dong
2016-01-01
In this paper, a novel nonlinear framework of smoothing method, non-Gaussian delayed particle smoother (nGDPS), is proposed, which enables vehicle state estimation (VSE) with high accuracy taking into account the non-Gaussianity of the measurement and process noises. Within the proposed method, the multivariate Student's t-distribution is adopted in order to compute the probability distribution function (PDF) related to the process and measurement noises, which are assumed to be non-Gaussian distributed. A computation approach based on Ensemble Kalman Filter (EnKF) is designed to cope with the mean and the covariance matrix of the proposal non-Gaussian distribution. A delayed Gibbs sampling algorithm, which incorporates smoothing of the sampled trajectories over a fixed-delay, is proposed to deal with the sample degeneracy of particles. The performance is investigated based on the real-world data, which is collected by low-cost on-board vehicle sensors. The comparison study based on the real-world experiments and the statistical analysis demonstrates that the proposed nGDPS has significant improvement on the vehicle state accuracy and outperforms the existing filtering and smoothing methods. PMID:27187405
NASA Astrophysics Data System (ADS)
Kang, S.; Suh, Y. K.
2011-02-01
The so-called smoothed profile method, originally suggested by Nakayama and Yamamoto and further improved by Luo et al. in 2005 and 2009, respectively, is an efficient numerical solver for fluid-structure interaction problems, which represents the particles by a certain smoothed profile on a fixed grid and constructs some form of body force added into the momentum (Navier-Stokes) equation by ensuring the rigidity of particles. For numerical simulations, the method first advances the flow and pressure fields by integrating the momentum equation except the body-force (momentum impulse) term in time and next updates them by separately taking temporal integration of the body-force term, thus requiring one more Poisson-equation solver for the extra pressure field due to the rigidity of particles to ensure the divergence-free constraint of the total velocity field. In the present study, we propose a simplified version of the smoothed profile method or the one-stage method, which combines the two stages of velocity update (temporal integration) into one to eliminate the necessity for the additional solver and, thus, significantly save the computational cost. To validate the proposed one-stage method, we perform the so-called direct numerical simulations on the two-dimensional motion of multiple inertialess paramagnetic particles in a nonmagnetic fluid subjected to an external uniform magnetic field and compare their results with the existing benchmark solutions. For the validation, we develop the finite-volume version of the direct simulation method by employing the proposed one-stage method. Comparison shows that the proposed one-stage method is very accurate and efficient in direct simulations of such magnetic particulate flows.
Biomechanical effects of environmental and engineered particles on human airway smooth muscle cells
Berntsen, P.; Park, C. Y.; Rothen-Rutishauser, B.; Tsuda, A.; Sager, T. M.; Molina, R. M.; Donaghey, T. C.; Alencar, A. M.; Kasahara, D. I.; Ericsson, T.; Millet, E. J.; Swenson, J.; Tschumperlin, D. J.; Butler, J. P.; Brain, J. D.; Fredberg, J. J.; Gehr, P.; Zhou, E. H.
2010-01-01
The past decade has seen significant increases in combustion-generated ambient particles, which contain a nanosized fraction (less than 100 nm), and even greater increases have occurred in engineered nanoparticles (NPs) propelled by the booming nanotechnology industry. Although inhalation of these particulates has become a public health concern, human health effects and mechanisms of action for NPs are not well understood. Focusing on the human airway smooth muscle cell, here we show that the cellular mechanical function is altered by particulate exposure in a manner that is dependent upon particle material, size and dose. We used Alamar Blue assay to measure cell viability and optical magnetic twisting cytometry to measure cell stiffness and agonist-induced contractility. The eight particle species fell into four categories, based on their respective effect on cell viability and on mechanical function. Cell viability was impaired and cell contractility was decreased by (i) zinc oxide (40–100 nm and less than 44 μm) and copper(II) oxide (less than 50 nm); cell contractility was decreased by (ii) fluorescent polystyrene spheres (40 nm), increased by (iii) welding fumes and unchanged by (iv) diesel exhaust particles, titanium dioxide (25 nm) and copper(II) oxide (less than 5 μm), although in none of these cases was cell viability impaired. Treatment with hydrogen peroxide up to 500 μM did not alter viability or cell mechanics, suggesting that the particle effects are unlikely to be mediated by particle-generated reactive oxygen species. Our results highlight the susceptibility of cellular mechanical function to particulate exposures and suggest that direct exposure of the airway smooth muscle cells to particulates may initiate or aggravate respiratory diseases. PMID:20356875
Biomechanical effects of environmental and engineered particles on human airway smooth muscle cells.
Berntsen, P; Park, C Y; Rothen-Rutishauser, B; Tsuda, A; Sager, T M; Molina, R M; Donaghey, T C; Alencar, A M; Kasahara, D I; Ericsson, T; Millet, E J; Swenson, J; Tschumperlin, D J; Butler, J P; Brain, J D; Fredberg, J J; Gehr, P; Zhou, E H
2010-06-01
The past decade has seen significant increases in combustion-generated ambient particles, which contain a nanosized fraction (less than 100 nm), and even greater increases have occurred in engineered nanoparticles (NPs) propelled by the booming nanotechnology industry. Although inhalation of these particulates has become a public health concern, human health effects and mechanisms of action for NPs are not well understood. Focusing on the human airway smooth muscle cell, here we show that the cellular mechanical function is altered by particulate exposure in a manner that is dependent upon particle material, size and dose. We used Alamar Blue assay to measure cell viability and optical magnetic twisting cytometry to measure cell stiffness and agonist-induced contractility. The eight particle species fell into four categories, based on their respective effect on cell viability and on mechanical function. Cell viability was impaired and cell contractility was decreased by (i) zinc oxide (40-100 nm and less than 44 microm) and copper(II) oxide (less than 50 nm); cell contractility was decreased by (ii) fluorescent polystyrene spheres (40 nm), increased by (iii) welding fumes and unchanged by (iv) diesel exhaust particles, titanium dioxide (25 nm) and copper(II) oxide (less than 5 microm), although in none of these cases was cell viability impaired. Treatment with hydrogen peroxide up to 500 microM did not alter viability or cell mechanics, suggesting that the particle effects are unlikely to be mediated by particle-generated reactive oxygen species. Our results highlight the susceptibility of cellular mechanical function to particulate exposures and suggest that direct exposure of the airway smooth muscle cells to particulates may initiate or aggravate respiratory diseases. PMID:20356875
NASA Astrophysics Data System (ADS)
Kordilla, J.; Tartakovsky, A. M.; Pan, W.; Shigorina, E.; Noffz, T.; Geyer, T.
2015-12-01
Unsaturated flow in fractured porous media exhibits highly complex flow dynamics and a wide range of intermittent flow processes. Especially in wide aperture fractures, flow processes may be dominated by gravitational instead of capillary forces leading to a deviation from the classical volume effective approaches (Richard's equation, Van Genuchten type relationships). The existence of various flow modes such as droplets, rivulets, turbulent and adsorbed films is well known, however, their spatial and temporal distribution within fracture networks is still an open question partially due to the lack of appropriate modeling tools. With our work we want to gain a deeper understanding of the underlying flow and transport dynamics in unsaturated fractured media in order to support the development of more refined upscaled methods, applicable on catchment scales. We present fracture-scale flow simulations obtained with a parallelized Smoothed Particle Hydrodynamics (SPH) model. The model allows us to simulate free-surface flow dynamics including the effect of surface tension for a wide range of wetting conditions in smooth and rough fractures. Due to the highly efficient generation of surface tension via particle-particle interaction forces the dynamic wetting of surfaces can readily be obtained. We validated the model via empirical and semi-analytical solutions and conducted laboratory-scale percolation experiments of unsaturated flow through synthetic fracture systems. The setup allows us to obtain travel time distributions and identify characteristic flow mode distributions on wide aperture fractures intercepted by horizontal fracture elements.
NASA Astrophysics Data System (ADS)
Juan-Mian, Lei; Xue-Ying, Peng
2016-02-01
Kernel gradient free-smoothed particle hydrodynamics (KGF-SPH) is a modified smoothed particle hydrodynamics (SPH) method which has higher precision than the conventional SPH. However, the Laplacian in KGF-SPH is approximated by the two-pass model which increases computational cost. A new kind of discretization scheme for the Laplacian is proposed in this paper, then a method with higher precision and better stability, called Improved KGF-SPH, is developed by modifying KGF-SPH with this new Laplacian model. One-dimensional (1D) and two-dimensional (2D) heat conduction problems are used to test the precision and stability of the Improved KGF-SPH. The numerical results demonstrate that the Improved KGF-SPH is more accurate than SPH, and stabler than KGF-SPH. Natural convection in a closed square cavity at different Rayleigh numbers are modeled by the Improved KGF-SPH with shifting particle position, and the Improved KGF-SPH results are presented in comparison with those of SPH and finite volume method (FVM). The numerical results demonstrate that the Improved KGF-SPH is a more accurate method to study and model the heat transfer problems.
Somasundaram, Deepak S; Trabia, Mohamed; O'Toole, Brendan; Hixson, Robert S
2014-01-23
This paper describes our work to characterize the variables affecting the smoothed particle hydrodynamics (SPH) method in the LS-DYNA package for simulating high-velocity flyer plate impact experiments. LS-DYNA simulations are compared with one-dimensional experimental data of an oxygen-free high-conductivity (OFHC) copper flyer plate impacting another plate of the same material. The comparison is made by measuring the velocity of a point on the back surface of the impact plate using the velocity interferometer system for any reflector (VISAR) technique.
The Ultraviolet View of Multi-Spin Galaxies: Insight from Smooth Particle Hydrodynamic Simulations
NASA Astrophysics Data System (ADS)
Bettoni, D.; Mazzei, P.; Marino, A.; Rampazzo, R.; Galletta, G.; Buson, L. M.
2014-05-01
The UV images of the Galaxy Evolution Explorer satellite revealed that about (30±3)% of early-type galaxies show UV emission indicating a rejuvenation episode. In early-type galaxies with multi-spin components this percentage increases at 50%. We present here the characteristics of this sample and our smooth particle hydrodynamic simulations with chemo-photometric implementation that provide dynamical and morphological information together with the spectral energy distribution at each evolutionary stage. We show our match of the global properties of two early-type galaxies, NGC 3626 and NGC 5173. For these galaxies we can trace their evolutionary path.
Coupled discrete element and smoothed particle hydrodynamics simulations of the die filling process
NASA Astrophysics Data System (ADS)
Breinlinger, Thomas; Kraft, Torsten
2015-08-01
Die filling is an important part of the powder compaction process chain, where defects in the final part can be introduced—or prevented. Simulation of this process is therefore a goal for many part producers and has been studied by some researchers already. In this work, we focus on the influence of the surrounding air on the powder flow. We demonstrate the implementing and coupling of the discrete element method for the granular powder and the smoothed particle hydrodynamics method for the gas flow. Application of the method to the die filling process is demonstrated.
Simulation of wave mitigation by coastal vegetation using smoothed particle hydrodynamics method
NASA Astrophysics Data System (ADS)
Iryanto; Gunawan, P. H.
2016-02-01
Vegetation in coastal area lead to wave mitigation has been studied by some researchers recently. The effect of vegetation forest in coastal area is minimizing the negative impact of wave propagation. In order to describe the effect of vegetation resistance into the water flow, the modified model of framework smoothed hydrodynamics particle has been constructed. In the Lagrangian framework, the Darcy, Manning, and laminar viscosity resistances are added. The effect of each resistances is given in some results of numerical simulations. Simulation of wave mitigation on sloping beach is also given.
DualSPHysics: Open-source parallel CFD solver based on Smoothed Particle Hydrodynamics (SPH)
NASA Astrophysics Data System (ADS)
Crespo, A. J. C.; Domínguez, J. M.; Rogers, B. D.; Gómez-Gesteira, M.; Longshaw, S.; Canelas, R.; Vacondio, R.; Barreiro, A.; García-Feal, O.
2015-02-01
DualSPHysics is a hardware accelerated Smoothed Particle Hydrodynamics code developed to solve free-surface flow problems. DualSPHysics is an open-source code developed and released under the terms of GNU General Public License (GPLv3). Along with the source code, a complete documentation that makes easy the compilation and execution of the source files is also distributed. The code has been shown to be efficient and reliable. The parallel power computing of Graphics Computing Units (GPUs) is used to accelerate DualSPHysics by up to two orders of magnitude compared to the performance of the serial version.
NASA Astrophysics Data System (ADS)
Pedretti, Daniele; Fernàndez-Garcia, Daniel
2013-09-01
Particle tracking methods to simulate solute transport deal with the issue of having to reconstruct smooth concentrations from a limited number of particles. This is an error-prone process that typically leads to large fluctuations in the determined late-time behavior of breakthrough curves (BTCs). Kernel density estimators (KDE) can be used to automatically reconstruct smooth BTCs from a small number of particles. The kernel approach incorporates the uncertainty associated with subsampling a large population by equipping each particle with a probability density function. Two broad classes of KDE methods can be distinguished depending on the parametrization of this function: global and adaptive methods. This paper shows that each method is likely to estimate a specific portion of the BTCs. Although global methods offer a valid approach to estimate early-time behavior and peak of BTCs, they exhibit important fluctuations at the tails where fewer particles exist. In contrast, locally adaptive methods improve tail estimation while oversmoothing both early-time and peak concentrations. Therefore a new method is proposed combining the strength of both KDE approaches. The proposed approach is universal and only needs one parameter (α) which slightly depends on the shape of the BTCs. Results show that, for the tested cases, heavily-tailed BTCs are properly reconstructed with α ≈ 0.5 .
NASA Astrophysics Data System (ADS)
Yang, Xiufeng; Liu, Moubin; Peng, Shiliu
2014-12-01
This paper presents a smoothed particle hydrodynamics (SPH) and element bending group (EBG) coupling method for modeling the interaction of flexible fibers with moving viscous fluids. SPH is a well-developed mesh-free particle method for simulating viscous fluid flows. EBG is also a particle method for modeling flexible bodies. The interaction of flexible fibers with moving viscous fluids is rendered through the interaction of EBG particles for flexible fiber and SPH particles for fluid. In numerical simulation, flexible fibers of different lengths are immersed in a moving viscous fluid driven by a body force. The drag force on the fiber obtained from SPH-EBG simulation agrees well with experimental observations. It is shown that the flexible fiber demonstrates three typical bending modes, including the U-shaped mode, the flapping mode, and the closed mode, and that the flexible fiber experiences a drag reduction due to its reconfiguration by bending. It is also found that the U 4/3 drag scaling law for a flexible fiber is only valid for the U-shaped mode, but not valid for the flapping and closed modes. The results indicate that the reconfiguration of a flexible fiber is caused by the fluid force acting on it, while vortex shedding is of importance in the translations of bending modes.
Yang, Xiufeng; Liu, Moubin; Peng, Shiliu
2014-12-01
This paper presents a smoothed particle hydrodynamics (SPH) and element bending group (EBG) coupling method for modeling the interaction of flexible fibers with moving viscous fluids. SPH is a well-developed mesh-free particle method for simulating viscous fluid flows. EBG is also a particle method for modeling flexible bodies. The interaction of flexible fibers with moving viscous fluids is rendered through the interaction of EBG particles for flexible fiber and SPH particles for fluid. In numerical simulation, flexible fibers of different lengths are immersed in a moving viscous fluid driven by a body force. The drag force on the fiber obtained from SPH-EBG simulation agrees well with experimental observations. It is shown that the flexible fiber demonstrates three typical bending modes, including the U-shaped mode, the flapping mode, and the closed mode, and that the flexible fiber experiences a drag reduction due to its reconfiguration by bending. It is also found that the U4/3 drag scaling law for a flexible fiber is only valid for the U-shaped mode, but not valid for the flapping and closed modes. The results indicate that the reconfiguration of a flexible fiber is caused by the fluid force acting on it, while vortex shedding is of importance in the translations of bending modes. PMID:25615191
NASA Astrophysics Data System (ADS)
Vázquez-Quesada, Adolfo; Bian, Xin; Ellero, Marco
2016-04-01
A three-dimensional model for a suspension of rigid spherical particles in a Newtonian fluid is presented. The solvent is modeled with smoothed particle hydrodynamics method, which takes into account exactly the long-range multi-body hydrodynamic interactions between suspended spheres. Short-range lubrication forces which are necessary to simulate concentrated suspensions, are introduced pair-wisely based on the analytical solution of Stokes equations for approaching/departing objects. Given that lubrication is singular at vanishing solid particle separations, an implicit splitting integration scheme is used to obtain accurate results and at the same time to avoid prohibitively small simulation time steps. Hydrodynamic interactions between solid particles, at both long-range and short-range limits, are verified against theory in the case of two approaching spheres in a quiescent medium and under bulk shear flow, where good agreements are obtained. Finally, numerical results for the suspension viscosity of a many-particle system are shown and compared with analytical solutions available in the dilute and semi-dilute case as well as with previous numerical results obtained in the concentrated limit.
Mandell, D.A.; Wingate, C.A.
1994-08-01
The design of many military devices involves numerical predictions of the material strength and fracture of brittle materials. The materials of interest include ceramics, that are used in armor packages; glass that is used in truck and jeep windshields and in helicopters; and rock and concrete that are used in underground bunkers. As part of a program to develop advanced hydrocode design tools, the authors have implemented a brittle fracture model for glass into the SPHINX smooth particle hydrodynamics code. The authors have evaluated this model and the code by predicting data from one-dimensional flyer plate impacts into glass, and data from tungsten rods impacting glass. Since fractured glass properties, which are needed in the model, are not available, the authors did sensitivity studies of these properties, as well as sensitivity studies to determine the number of particles needed in the calculations. The numerical results are in good agreement with the data.
Mandell, D.A.; Wingate, C.A.; Stellingwwerf, R.F.
1995-12-31
The design of many devices involves numerical predictions of the material strength and fracture of brittle materials. The materials of interest include ceramics that are used in armor packages; glass that is used in windshields; and rock and concrete that are used in oil wells. As part of a program to develop advanced hydrocode design tools, the authors have implemented a brittle fracture model for glass into the SPHINX smooth particle hydrodynamics code. The authors have evaluated this model and the code by predicting data from tungsten rods impacting glass. Since fractured glass properties, which are needed in the model, are not available, they did sensitivity studies of these properties, as well as sensitivity studies to determine the number of particles needed in the calculations. The numerical results are in good agreement with the data.
NASA Astrophysics Data System (ADS)
Lei, Huan; Baker, Nathan A.; Wu, Lei; Schenter, Gregory K.; Mundy, Christopher J.; Tartakovsky, Alexandre M.
2016-08-01
Thermal fluctuations cause perturbations of fluid-fluid interfaces and highly nonlinear hydrodynamics in multiphase flows. In this work, we develop a multiphase smoothed dissipative particle dynamics (SDPD) model. This model accounts for both bulk hydrodynamics and interfacial fluctuations. Interfacial surface tension is modeled by imposing a pairwise force between SDPD particles. We show that the relationship between the model parameters and surface tension, previously derived under the assumption of zero thermal fluctuation, is accurate for fluid systems at low temperature but overestimates the surface tension for intermediate and large thermal fluctuations. To analyze the effect of thermal fluctuations on surface tension, we construct a coarse-grained Euler lattice model based on the mean field theory and derive a semianalytical formula to directly relate the surface tension to model parameters for a wide range of temperatures and model resolutions. We demonstrate that the present method correctly models dynamic processes, such as bubble coalescence and capillary spectra across the interface.
Three-dimensional model of a plasma railgun using smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Jackson, Lloyd M.
Pulsed plasma accelerators are utilized for in-space propulsion and drivers for inertial fusion concepts. Theoretical models are necessary to assist in diagnostic analysis and for developing scaling laws. SPFMax is a new 3D code which uses smoothed particle hydrodynamics (SPH) to simulate fluid flow, and has been designed specifically for modeling plasmas produced by these pulsed devices. A set of gasdynamic test cases were established and utilized to verify the accuracy of SPFMax for modeling the gas dynamics in a railgun. The free expansion confirmed that the gas expands supersonically without exceeding the predicted maximum value. With the square wave test, SPFMax advected the waves with floating point accuracy. Shocks, expansion waves, and contact surfaces were resolved in both 1D and 3D tests with a relatively low number of particles.
Adaptive fuzzy control with smooth inverse for nonlinear systems preceded by non-symmetric dead-zone
NASA Astrophysics Data System (ADS)
Wang, Xingjian; Wang, Shaoping
2016-07-01
In this study, the adaptive output feedback control problem of a class of nonlinear systems preceded by non-symmetric dead-zone is considered. To cope with the possible control signal chattering phenomenon which is caused by non-smooth dead-zone inverse, a new smooth inverse is proposed for non-symmetric dead-zone compensation. For the systematic design procedure of the adaptive fuzzy control algorithm, we combine the backstepping technique and small-gain approach. The Takagi-Sugeno fuzzy logic systems are used to approximate unknown system nonlinearities. The closed-loop stability is studied by using small gain theorem and the closed-loop system is proved to be semi-globally uniformly ultimately bounded. Simulation results indicate that, compared to the algorithm with the non-smooth inverse, the proposed control strategy can achieve better tracking performance and the chattering phenomenon can be avoided effectively.
Luo, Xiongbiao; Wan, Ying; He, Xiangjian; Mori, Kensaku
2015-08-01
This paper proposes an observation-driven adaptive differential evolution algorithm that fuses bronchoscopic video sequences, electromagnetic sensor measurements, and computed tomography images for accurate and smooth bronchoscope three-dimensional motion tracking. Currently an electromagnetic tracker with a position sensor fixed at the bronchoscope tip is commonly used to estimate bronchoscope movements. The large tracking error from directly using sensor measurements, which may be deteriorated heavily by patient respiratory motion and the magnetic field distortion of the tracker, limits clinical applications. How to effectively use sensor measurements for precise and stable bronchoscope electromagnetic tracking remains challenging. We here exploit an observation-driven adaptive differential evolution framework to address such a challenge and boost the tracking accuracy and smoothness. In our framework, two advantageous points are distinguished from other adaptive differential evolution methods: (1) the current observation including sensor measurements and bronchoscopic video images is used in the mutation equation and the fitness computation, respectively and (2) the mutation factor and the crossover rate are determined adaptively on the basis of the current image observation. The experimental results demonstrate that our framework provides much more accurate and smooth bronchoscope tracking than the state-of-the-art methods. Our approach reduces the tracking error from 3.96 to 2.89 mm, improves the tracking smoothness from 4.08 to 1.62 mm, and increases the visual quality from 0.707 to 0.741. PMID:25660001
Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method
NASA Astrophysics Data System (ADS)
Zhang, Mingyu
2010-09-01
The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503-1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710-1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.
Keller, Amy C.; Knaub, Leslie A.; McClatchey, P. Mason; Connon, Chelsea A.; Bouchard, Ron; Miller, Matthew W.; Geary, Kate E.; Walker, Lori A.; Klemm, Dwight J.; Reusch, Jane E. B.
2016-01-01
Diabetes affects more than 330 million people worldwide and causes elevated cardiovascular disease risk. Mitochondria are critical for vascular function, generate cellular reactive oxygen species (ROS), and are perturbed by diabetes, representing a novel target for therapeutics. We hypothesized that adaptive mitochondrial plasticity in response to nutrient stress would be impaired in diabetes cellular physiology via a nitric oxide synthase- (NOS-) mediated decrease in mitochondrial function. Primary smooth muscle cells (SMCs) from aorta of the nonobese, insulin resistant rat diabetes model Goto-Kakizaki (GK) and the Wistar control rat were exposed to high glucose (25 mM). At baseline, significantly greater nitric oxide evolution, ROS production, and respiratory control ratio (RCR) were observed in GK SMCs. Upon exposure to high glucose, expression of phosphorylated eNOS, uncoupled respiration, and expression of mitochondrial complexes I, II, III, and V were significantly decreased in GK SMCs (p < 0.05). Mitochondrial superoxide increased with high glucose in Wistar SMCs (p < 0.05) with no change in the GK beyond elevated baseline concentrations. Baseline comparisons show persistent metabolic perturbations in a diabetes phenotype. Overall, nutrient stress in GK SMCs caused a persistent decline in eNOS and mitochondrial function and disrupted mitochondrial plasticity, illustrating eNOS and mitochondria as potential therapeutic targets. PMID:27034743
Flight test results from a supercritical mission adaptive wing with smooth variable camber
NASA Technical Reports Server (NTRS)
Powers, Sheryll Goecke; Webb, Lannie D.; Friend, Edward L.; Lokos, William A.
1992-01-01
The mission adaptive wing (MAW) consisted of leading- and trailing-edge variable-camber surfaces that could be deflected in flight to provide a near-ideal wing camber shape for any flight condition. These surfaces featured smooth, flexible upper surfaces and fully enclosed lower surfaces, distinguishing them from conventional flaps that have discontinuous surfaces and exposed or semiexposed mechanisms. Camber shape was controlled by either a manual or automatic flight control system. The wing and aircraft were extensively instrumented to evaluate the local flow characteristics and the total aircraft performance. This paper discusses the interrelationships between the wing pressure, buffet, boundary-layer and flight deflection measurement system analyses and describes the flight maneuvers used to obtain the data. The results are for a wing sweep of 26 deg, a Mach number of 0.85, leading and trailing-edge cambers (delta(sub LE/TE)) of 0/2 and 5/10, and angles of attack from 3.0 deg to 14.0 deg. For the well-behaved flow of the delta(sub LE/TE) = 0/2 camber, a typical cruise camber shape, the local and global data are in good agreement with respect to the flow properties of the wing. For the delta(sub LE/TE) = 5/10 camber, a maneuvering camber shape, the local and global data have similar trends and conclusions, but not the clear-cut agreement observed for cruise camber.
Adaptive Optimal Kernel Smooth-Windowed Wigner-Ville Distribution for Digital Communication Signal
NASA Astrophysics Data System (ADS)
Tan, Jo Lynn; Sha'ameri, Ahmad Zuribin
2009-12-01
Time-frequency distributions (TFDs) are powerful tools to represent the energy content of time-varying signal in both time and frequency domains simultaneously but they suffer from interference due to cross-terms. Various methods have been described to remove these cross-terms and they are typically signal-dependent. Thus, there is no single TFD with a fixed window or kernel that can produce accurate time-frequency representation (TFR) for all types of signals. In this paper, a globally adaptive optimal kernel smooth-windowed Wigner-Ville distribution (AOK-SWWVD) is designed for digital modulation signals such as ASK, FSK, and M-ary FSK, where its separable kernel is determined automatically from the input signal, without prior knowledge of the signal. This optimum kernel is capable of removing the cross-terms and maintaining accurate time-frequency representation at SNR as low as 0 dB. It is shown that this system is comparable to the system with prior knowledge of the signal.
Keller, Amy C; Knaub, Leslie A; McClatchey, P Mason; Connon, Chelsea A; Bouchard, Ron; Miller, Matthew W; Geary, Kate E; Walker, Lori A; Klemm, Dwight J; Reusch, Jane E B
2016-01-01
Diabetes affects more than 330 million people worldwide and causes elevated cardiovascular disease risk. Mitochondria are critical for vascular function, generate cellular reactive oxygen species (ROS), and are perturbed by diabetes, representing a novel target for therapeutics. We hypothesized that adaptive mitochondrial plasticity in response to nutrient stress would be impaired in diabetes cellular physiology via a nitric oxide synthase- (NOS-) mediated decrease in mitochondrial function. Primary smooth muscle cells (SMCs) from aorta of the nonobese, insulin resistant rat diabetes model Goto-Kakizaki (GK) and the Wistar control rat were exposed to high glucose (25 mM). At baseline, significantly greater nitric oxide evolution, ROS production, and respiratory control ratio (RCR) were observed in GK SMCs. Upon exposure to high glucose, expression of phosphorylated eNOS, uncoupled respiration, and expression of mitochondrial complexes I, II, III, and V were significantly decreased in GK SMCs (p < 0.05). Mitochondrial superoxide increased with high glucose in Wistar SMCs (p < 0.05) with no change in the GK beyond elevated baseline concentrations. Baseline comparisons show persistent metabolic perturbations in a diabetes phenotype. Overall, nutrient stress in GK SMCs caused a persistent decline in eNOS and mitochondrial function and disrupted mitochondrial plasticity, illustrating eNOS and mitochondria as potential therapeutic targets. PMID:27034743
Xiao, Jianyuan; Liu, Jian; Qin, Hong; Plasma Physics Laboratory, Princeton University, Princeton, New Jersey 08543 ; Yu, Zhi
2013-10-15
Smoothing functions are commonly used to reduce numerical noise arising from coarse sampling of particles in particle-in-cell (PIC) plasma simulations. When applying smoothing functions to symplectic algorithms, the conservation of symplectic structure should be guaranteed to preserve good conservation properties. In this paper, we show how to construct a variational multi-symplectic PIC algorithm with smoothing functions for the Vlasov-Maxwell system. The conservation of the multi-symplectic structure and the reduction of numerical noise make this algorithm specifically suitable for simulating long-term dynamics of plasmas, such as those in the steady-state operation or long-pulse discharge of a super-conducting tokamak. The algorithm has been implemented in a 6D large scale PIC code. Numerical examples are given to demonstrate the good conservation properties of the multi-symplectic algorithm and the reduction of the noise due to the application of smoothing function.
Multi-resolution flow simulations by smoothed particle hydrodynamics via domain decomposition
NASA Astrophysics Data System (ADS)
Bian, Xin; Li, Zhen; Karniadakis, George Em
2015-09-01
We present a methodology to concurrently couple particle-based methods via a domain decomposition (DD) technique for simulating viscous flows. In particular, we select two resolutions of the smoothed particle hydrodynamics (SPH) method as demonstration. Within the DD framework, a simulation domain is decomposed into two (or more) overlapping sub-domains, each of which has an individual particle scale determined by the local flow physics. Consistency of the two sub-domains is achieved in the overlap region by matching the two independent simulations based on Lagrangian interpolation of state variables and fluxes. The domain decomposition based SPH method (DD-SPH) employs different spatial and temporal resolutions, and hence, each sub-domain has its own smoothing length and time step. As a consequence, particle refinement and de-refinement are performed asynchronously according to individual time advancement of each sub-domain. The proposed strategy avoids SPH force interactions between different resolutions on purpose, so that coupling, in principle, can go beyond SPH-SPH, and may allow SPH to be coupled with other mesoscopic or microscopic particle methods. The DD-SPH method is validated first for a transient Couette flow, where simulation results based on proper coupling of spatial-temporal scales agree well with analytical solutions. In particular, we find that the size of the overlap region should be at least rc,1 + 2rc,2, where rc,1 and rc,2 are cut off radii in the two sub-domains with rc,1 ≤rc,2. Subsequently, a perturbation wave is considered traveling either parallel or perpendicular to the hybrid interface. Compressibility is significant if transient behavior at short sonic-time-scale is relevant, while the fluid can be treated as quasi-incompressible at sufficiently long time scale. To this end, we propose a coupling of density fields from the two sub-domains. Finally, a steady Wannier flow is simulated, where a rotating cylinder is placed next to a
Ganzenmüller, Georg C.; Hiermaier, Stefan; Steinhauser, Martin O.
2012-01-01
We propose a thermodynamically consistent and energy-conserving temperature coupling scheme between the atomistic and the continuum domain. The coupling scheme links the two domains using the DPDE (Dissipative Particle Dynamics at constant Energy) thermostat and is designed to handle strong temperature gradients across the atomistic/continuum domain interface. The fundamentally different definitions of temperature in the continuum and atomistic domain – internal energy and heat capacity versus particle velocity – are accounted for in a straightforward and conceptually intuitive way by the DPDE thermostat. We verify the here-proposed scheme using a fluid, which is simultaneously represented as a continuum using Smooth Particle Hydrodynamics, and as an atomistically resolved liquid using Molecular Dynamics. In the case of equilibrium contact between both domains, we show that the correct microscopic equilibrium properties of the atomistic fluid are obtained. As an example of a strong non-equilibrium situation, we consider the propagation of a steady shock-wave from the continuum domain into the atomistic domain, and show that the coupling scheme conserves both energy and shock-wave dynamics. To demonstrate the applicability of our scheme to real systems, we consider shock loading of a phospholipid bilayer immersed in water in a multi-scale simulation, an interesting topic of biological relevance. PMID:23300586
NASA Astrophysics Data System (ADS)
Lind, S. J.; Stansby, P. K.; Rogers, B. D.
2016-03-01
A new two-phase incompressible-compressible Smoothed Particle Hydrodynamics (SPH) method has been developed where the interface is discontinuous in density. This is applied to water-air problems with a large density difference. The incompressible phase requires surface pressure from the compressible phase and the compressible phase requires surface velocity from the incompressible phase. Compressible SPH is used for the air phase (with the isothermal stiffened ideal gas equation of state for low Mach numbers) and divergence-free (projection based) incompressible SPH is used for the water phase, with the addition of Fickian shifting to produce sufficiently homogeneous particle distributions to enable stable, accurate, converged solutions without noise in the pressure field. Shifting is a purely numerical particle regularisation device. The interface remains a true material discontinuity at a high density ratio with continuous pressure and velocity at the interface. This approach with the physics of compressibility and incompressibility represented is novel within SPH and is validated against semi-analytical results for a two-phase elongating and oscillating water drop, analytical results for low amplitude inviscid standing waves, the Kelvin-Helmholtz instability, and a dam break problem with high interface distortion and impact on a vertical wall where experimental and other numerical results are available.
Hoover, W.G.
1999-02-01
The eventual entropy increase of an ideal gas undergoing free expansion, {Delta}S=k&hthinsp;ln(V{sub fin}/V{sub 0}), requires a {open_quotes}coarse-grained{close_quotes} hydrodynamic description because Gibbs{close_quote} fine-grained entropy is unchanged in such a process. Smooth particle applied mechanics (SPAM) is well suited to the simulation and study of such problems because the particles in SPAM simulations can be of any size, from microscopic to macroscopic. SPAM furnishes a natural interpolation, or bridge, linking microscopic molecular dynamics to macroscopic continuum mechanics. We analyze particle-based simulations of ideal-gas free expansions from both the microscopic and the macroscopic points of view, comparing several dynamical estimates for the time development of the system entropy. Most of the entropy increase occurs rapidly, within a single sound traversal time. A local comoving version of turbulent hydrodynamics provides the most useful viewpoint for describing flows of this kind. {copyright} {ital 1999} {ital The American Physical Society}
Crespo, Alejandro C.; Dominguez, Jose M.; Barreiro, Anxo; Gómez-Gesteira, Moncho; Rogers, Benedict D.
2011-01-01
Smoothed Particle Hydrodynamics (SPH) is a numerical method commonly used in Computational Fluid Dynamics (CFD) to simulate complex free-surface flows. Simulations with this mesh-free particle method far exceed the capacity of a single processor. In this paper, as part of a dual-functioning code for either central processing units (CPUs) or Graphics Processor Units (GPUs), a parallelisation using GPUs is presented. The GPU parallelisation technique uses the Compute Unified Device Architecture (CUDA) of nVidia devices. Simulations with more than one million particles on a single GPU card exhibit speedups of up to two orders of magnitude over using a single-core CPU. It is demonstrated that the code achieves different speedups with different CUDA-enabled GPUs. The numerical behaviour of the SPH code is validated with a standard benchmark test case of dam break flow impacting on an obstacle where good agreement with the experimental results is observed. Both the achieved speed-ups and the quantitative agreement with experiments suggest that CUDA-based GPU programming can be used in SPH methods with efficiency and reliability. PMID:21695185
Smoothed Particle Hydrodynamics Stochastic Model for Flow and Transport in Porous Media
Tartakovsky, Alexandre M.; Tartakovsky, Daniel M.; Meakin, Paul
2008-11-03
A meso-scale stochastic Lagrangian particle model was developed and used to simulate conservative and reactive transport in porous media. In the stochastic model, the fluid flow in a porous continuum is governed by a combination of a Langevin equation and continuity equation. Pore-scale velocity fluctuations, the source of hydrodynamic dispersion, are represented by the white noise. A smoothed particle hydrodynamics method was used to solve the governing equations. Changes in the properties of the fluid particles (e.g., the solute concentration) are governed by the advection-diffusion equation. The separate treatment of advective and diffusive mixing in the stochastic transport model is more realistic than the classical advection-dispersion theory, which uses a single effective diffusion coefficient (the dispersion coefficient) to describe both types of mixing leading to over-prediction of mixing induced effective reaction rates. The stochastic model predicts much lower reaction product concentrations in mixing induced reactions. In addition, the dispersion theory predicts more stable fronts (with a higher effective fractal dimension) than the stochastic model during the growth of Rayleigh-Taylor instabilities.
A synthetic 21-cm Galactic Plane Survey of a smoothed particle hydrodynamics galaxy simulation
NASA Astrophysics Data System (ADS)
Douglas, Kevin A.; Acreman, David M.; Dobbs, Clare L.; Brunt, Christopher M.
2010-09-01
We have created synthetic neutral hydrogen (HI) Galactic Plane Survey data cubes covering 90° <= l <= 180°, using a model spiral galaxy from smoothed particle hydrodynamics (SPH) simulations and the radiative transfer code TORUS. The density, temperature and other physical parameters are fed from the SPH simulation into TORUS, where the HI emissivity and opacity are calculated before the 21-cm line emission profile is determined. Our main focus is the observation of outer Galaxy `Perseus arm' HI, with a view to tracing atomic gas as it encounters shock motions as it enters a spiral arm interface, an early step in the formation of molecular clouds. The observation of HI self-absorption features at these shock sites (in both real observations and our synthetic data) allows us to investigate further the connection between cold atomic gas and the onset of molecular cloud formation.
Survey of shock-wave structures of smooth-particle granular flows.
Padgett, D A; Mazzoleni, A P; Faw, S D
2015-12-01
We show the effects of simulated supersonic granular flow made up of smooth particles passing over two prototypical bodies: a wedge and a disk. We describe a way of computationally identifying shock wave locations in granular flows and tabulate the shock wave locations for flow over wedges and disks. We quantify the shock structure in terms of oblique shock angle for wedge impediments and shock standoff distance for disk impediments. We vary granular flow parameters including upstream volume fraction, average upstream velocity, granular temperature, and the collision coefficient of restitution. Both wedges and disks have been used in the aerospace community as prototypical impediments to flowing air in order to investigate the fundamentally different shock structures emanating from sharp and blunt bodies, and we present these results in order to increase the understanding of the fundamental behavior of supersonic granular flow. PMID:26764684
An incompressible smoothed particle hydrodynamics method for the motion of rigid bodies in fluids
NASA Astrophysics Data System (ADS)
Tofighi, N.; Ozbulut, M.; Rahmat, A.; Feng, J. J.; Yildiz, M.
2015-09-01
A two-dimensional incompressible smoothed particle hydrodynamics scheme is presented for simulation of rigid bodies moving through Newtonian fluids. The scheme relies on combined usage of the rigidity constraints and the viscous penalty method to simulate rigid body motion. Different viscosity ratios and interpolation schemes are tested by simulating a rigid disc descending in quiescent medium. A viscosity ratio of 100 coupled with weighted harmonic averaging scheme has been found to provide satisfactory results. The performance of the resulting scheme is systematically tested for cases with linear motion, rotational motion and their combination. The test cases include sedimentation of a single and a pair of circular discs, sedimentation of an elliptic disc and migration and rotation of a circular disc in linear shear flow. Comparison with previous results at various Reynolds numbers indicates that the proposed method captures the motion of rigid bodies driven by flow or external body forces accurately.
NASA Astrophysics Data System (ADS)
Yang, G.; Han, X.; Hu, D. A.
2015-11-01
Modified cylindrical smoothed particle hydrodynamics (MCSPH) approximation equations are derived for hydrodynamics with material strength in axisymmetric cylindrical coordinates. The momentum equation and internal energy equation are represented to be in the axisymmetric form. The MCSPH approximation equations are applied to simulate the process of explosively driven metallic tubes, which includes strong shock waves, large deformations and large inhomogeneities, etc. The meshless and Lagrangian character of the MCSPH method offers the advantages in treating the difficulties embodied in these physical phenomena. Two test cases, the cylinder test and the metallic tube driven by two head-on colliding detonation waves, are presented. Numerical simulation results show that the new form of the MCSPH method can predict the detonation process of high explosives and the expansion process of metallic tubes accurately and robustly.
NASA Technical Reports Server (NTRS)
Libersky, Larry; Allahdadi, Firooz A.; Carney, Theodore C.
1992-01-01
Analysis of interaction occurring between space debris and orbiting structures is of great interest to the planning and survivability of space assets. Computer simulation of the impact events using hydrodynamic codes can provide some understanding of the processes but the problems involved with this fundamental approach are formidable. First, any realistic simulation is necessarily three-dimensional, e.g., the impact and breakup of a satellite. Second, the thickness of important components such as satellite skins or bumper shields are small with respect to the dimension of the structure as a whole, presenting severe zoning problems for codes. Thirdly, the debris cloud produced by the primary impact will yield many secondary impacts which will contribute to the damage and possible breakup of the structure. The problem was approached by choosing a relatively new computational technique that has virtues peculiar to space impacts. The method is called Smoothed Particle Hydrodynamics.
Two-fluid dust and gas mixtures in smoothed particle hydrodynamics: a semi-implicit approach
NASA Astrophysics Data System (ADS)
Lorén-Aguilar, Pablo; Bate, Matthew R.
2014-09-01
A method to avoid the explicit time integration of small dust grains in the two-fluid gas/dust smoothed particle hydrodynamics (SPH) approach is proposed. By assuming a very simple exponential decay model for the relative velocity between the gas and dust components, all the effective characteristics of the drag force can be reproduced. A series of tests has been performed to compare the accuracy of the method with analytical and explicit integration results. We find that the method performs well on a wide range of tests, and can provide large speed-ups over explicit integration when the dust stopping time is small. We have also found that the method is much less dissipative than conventional explicit or implicit two-fluid SPH approaches when modelling dusty shocks.
Simulation of a ceramic impact experiment using the SPHINX smooth particle hydrodynamics code
Mandell, D.A.; Wingate, C.A.; Schwalbe, L.A.
1996-08-01
We are developing statistically based, brittle-fracture models and are implementing them into hydrocodes that can be used for designing systems with components of ceramics, glass, and/or other brittle materials. Because of the advantages it has simulating fracture, we are working primarily with the smooth particle hydrodynamics code SPHINX. We describe a new brittle fracture model that we have implemented into SPHINX, and we discuss how the model differs from others. To illustrate the code`s current capability, we simulate an experiment in which a tungsten rod strikes a target of heavily confined ceramic. Simulations in 3D at relatively coarse resolution yield poor results. However, 2D plane-strain approximations to the test produce crack patterns that are strikingly similar to the data, although the fracture model needs further refinement to match some of the finer details. We conclude with an outline of plans for continuing research and development.
Scattering and bound states of spinless particles in a mixed vector-scalar smooth step potential
Garcia, M.G.; Castro, A.S. de
2009-11-15
Scattering and bound states for a spinless particle in the background of a kink-like smooth step potential, added with a scalar uniform background, are considered with a general mixing of vector and scalar Lorentz structures. The problem is mapped into the Schroedinger-like equation with an effective Rosen-Morse potential. It is shown that the scalar uniform background present subtle and trick effects for the scattering states and reveals itself a high-handed element for formation of bound states. In that process, it is shown that the problem of solving a differential equation for the eigenenergies is transmuted into the simpler and more efficient problem of solving an irrational algebraic equation.
Hoover, Wm G; Hoover, Carol G; Merritt, Elizabeth C
2004-01-01
Smooth-particle applied mechanics (SPAM) provides several approaches to approximate solutions of the continuum equations for both fluids and solids. Though many of the usual formulations conserve mass, (linear) momentum, and energy, the angular momentum is typically not conserved by SPAM. A second difficulty with the usual formulations is that tensile stress states often exhibit an exponentially fast high-frequency short-wavelength instability, "tensile instability." We discuss these twin defects of SPAM and illustrate them for a rotating elastic body. We formulate ways to conserve angular momentum while at the same time delaying the symptoms of tensile instability for many sound-traversal times. These ideas should prove useful in more general situations. PMID:14995750
Coupling of Smoothed Particle Hydrodynamics with Finite Volume method for free-surface flows
NASA Astrophysics Data System (ADS)
Marrone, S.; Di Mascio, A.; Le Touzé, D.
2016-04-01
A new algorithm for the solution of free surface flows with large front deformation and fragmentation is presented. The algorithm is obtained by coupling a classical Finite Volume (FV) approach, that discretizes the Navier-Stokes equations on a block structured Eulerian grid, with an approach based on the Smoothed Particle Hydrodynamics (SPH) method, implemented in a Lagrangian framework. The coupling procedure is formulated in such a way that each solver is applied in the region where its intrinsic characteristics can be exploited in the most efficient and accurate way: the FV solver is used to resolve the bulk flow and the wall regions, whereas the SPH solver is implemented in the free surface region to capture details of the front evolution. The reported results clearly prove that the combined use of the two solvers is convenient from the point of view of both accuracy and computing time.
Smoothed Particle Hydrodynamics pore-scale simulations of unstable immiscible flow in porous media
Bandara, Dunusinghe Mudiyanselage Uditha C.; Tartakovsky, Alexandre M.; Oostrom, Martinus; Palmer, Bruce J.; Grate, Jay W.; Zhang, Changyong
2013-12-01
We have conducted a series of high-resolution numerical experiments using the Pair-Wise Force Smoothed Particle Hydrodynamics (PF-SPH) multiphase flow model. First, we derived analytical expressions relating parameters in the PF-SPH model to the surface tension and static contact angle. Next, we used the model to study viscous fingering, capillary fingering, and stable displacement of immiscible fluids in porous media for a wide range of capillary numbers and viscosity ratios. We demonstrated that the steady state saturation profiles and the boundaries of viscous fingering, capillary fingering, and stable displacement regions compare favorably with micromodel laboratory experimental results. For displacing fluid with low viscosity, we observed that the displacement pattern changes from viscous fingering to stable displacement with increasing injection rate. When a high viscosity fluid is injected, transition behavior from capillary fingering to stable displacement occurred as the flow rate was increased. These observation also agree with the results of the micromodel laboratory experiments.
Particle Swarm Based Collective Searching Model for Adaptive Environment
Cui, Xiaohui; Patton, Robert M; Potok, Thomas E; Treadwell, Jim N
2007-01-01
This report presents a pilot study of an integration of particle swarm algorithm, social knowledge adaptation and multi-agent approaches for modeling the collective search behavior of self-organized groups in an adaptive environment. The objective of this research is to apply the particle swarm metaphor as a model of social group adaptation for the dynamic environment and to provide insight and understanding of social group knowledge discovering and strategic searching. A new adaptive environment model, which dynamically reacts to the group collective searching behaviors, is proposed in this research. The simulations in the research indicate that effective communication between groups is not the necessary requirement for whole self-organized groups to achieve the efficient collective searching behavior in the adaptive environment.
Particle Swarm Based Collective Searching Model for Adaptive Environment
Cui, Xiaohui; Patton, Robert M; Potok, Thomas E; Treadwell, Jim N
2008-01-01
This report presents a pilot study of an integration of particle swarm algorithm, social knowledge adaptation and multi-agent approaches for modeling the collective search behavior of self-organized groups in an adaptive environment. The objective of this research is to apply the particle swarm metaphor as a model of social group adaptation for the dynamic environment and to provide insight and understanding of social group knowledge discovering and strategic searching. A new adaptive environment model, which dynamically reacts to the group collective searching behaviors, is proposed in this research. The simulations in the research indicate that effective communication between groups is not the necessary requirement for whole self-organized groups to achieve the efficient collective searching behavior in the adaptive environment.
Dowie, Jackson; Ansell, Thomas K; Noble, Peter B; Donovan, Graham M
2016-01-01
Length adaptation is a phenomenon observed in airway smooth muscle (ASM) wherein over time there is a shift in the length-tension curve. There is potential for length adaptation to play an important role in airway constriction and airway hyper-responsiveness in asthma. Recent results by Ansell et al., 2015 (JAP 2014 10.1152/japplphysiol.00724.2014) have cast doubt on this role by testing for length adaptation using an intact airway preparation, rather than strips of ASM. Using this technique they found no evidence for length adaptation in intact airways. Here we attempt to resolve this apparent discrepancy by constructing a minimal mathematical model of the intact airway, including ASM which follows the classic length-tension curve and undergoes length adaptation. This allows us to show that (1) no evidence of length adaptation should be expected in large, cartilaginous, intact airways; (2) even in highly compliant peripheral airways, or at more compliant regions of the pressure-volume curve of large airways, the effect of length adaptation would be modest and at best marginally detectable in intact airways; (3) the key parameters which control the appearance of length adaptation in intact airways are airway compliance and the relaxation timescale. The results of this mathematical simulation suggest that length adaptation observed at the level of the isolated ASM may not clearly manifest in the normal intact airway. PMID:26376002
A smooth dissipative particle dynamics method for domains with arbitrary-geometry solid boundaries
NASA Astrophysics Data System (ADS)
Gatsonis, Nikolaos A.; Potami, Raffaele; Yang, Jun
2014-01-01
A smooth dissipative particle dynamics method with dynamic virtual particle allocation (SDPD-DV) for modeling and simulation of mesoscopic fluids in wall-bounded domains is presented. The physical domain in SDPD-DV may contain external and internal solid boundaries of arbitrary geometries, periodic inlets and outlets, and the fluid region. The SDPD-DV method is realized with fluid particles, boundary particles, and dynamically allocated virtual particles. The internal or external solid boundaries of the domain can be of arbitrary geometry and are discretized with a surface grid. These boundaries are represented by boundary particles with assigned properties. The fluid domain is discretized with fluid particles of constant mass and variable volume. Conservative and dissipative force models due to virtual particles exerted on a fluid particle in the proximity of a solid boundary supplement the original SDPD formulation. The dynamic virtual particle allocation approach provides the density and the forces due to virtual particles. The integration of the SDPD equations is accomplished with a velocity-Verlet algorithm for the momentum and a Runge-Kutta for the entropy equation. The velocity integrator is supplemented by a bounce-forward algorithm in cases where the virtual particle force model is not able to prevent particle penetration. For the incompressible isothermal systems considered in this work, the pressure of a fluid particle is obtained by an artificial compressibility formulation for liquids and the ideal gas law for gases. The self-diffusion coefficient is obtained by an implementation of the generalized Einstein and the Green-Kubo relations. Field properties are obtained by sampling SDPD-DV outputs on a post-processing grid that allows harnessing the particle information on desired spatiotemporal scales. The SDPD-DV method is verified and validated with simulations in bounded and periodic domains that cover the hydrodynamic and mesoscopic regimes for
Brown, C; Adcock, A; Azevedo, S; Liebman, J; Bond, E
2010-12-28
Some diagnostics at the National Ignition Facility (NIF), including the Gamma Reaction History (GRH) diagnostic, require multiple channels of data to achieve the required dynamic range. These channels need to be stitched together into a single time series, and they may have non-uniform and redundant time samples. We chose to apply the popular cubic smoothing spline technique to our stitching problem because we needed a general non-parametric method. We adapted one of the algorithms in the literature, by Hutchinson and deHoog, to our needs. The modified algorithm and the resulting code perform a cubic smoothing spline fit to multiple data channels with redundant time samples and missing data points. The data channels can have different, time-varying, zero-mean white noise characteristics. The method we employ automatically determines an optimal smoothing level by minimizing the Generalized Cross Validation (GCV) score. In order to automatically validate the smoothing level selection, the Weighted Sum-Squared Residual (WSSR) and zero-mean tests are performed on the residuals. Further, confidence intervals, both analytical and Monte Carlo, are also calculated. In this paper, we describe the derivation of our cubic smoothing spline algorithm. We outline the algorithm and test it with simulated and experimental data.
Mayoral-Villa, Estela; Alvarado-Rodríguez, Carlos E; Klapp, Jaime; Gómez-Gesteira, Moncho; Sigalotti, Leonardo Di G
2016-04-01
A smoothed particle hydrodynamics (SPH) model is presented for simulating the decay chain transport of radionuclides in confined aqueous solutions. The SPH formulation is based on the open-source parallel code DualSPHysics extended to solve the advective-diffusion equation for the evolution of the concentration field coupled to the fluid-dynamic equations, including the effects of radioactive decay of the tracer contaminants. The performance of the method is demonstrated for environmental engineering problems dealing with the transport of contaminants in still and flowing water. The results from a series of benchmark test calculations are described in two- and three-space dimensions, where the advection, diffusion, and radioactive decay modes are tested separately and in combined form. The accuracy of the present SPH transport model is shown by direct comparison with the analytical solutions and results from other SPH approaches. For a given problem, convergence of the SPH solution is seen to increase with decreasing particle size and spacing. PMID:26921532
NASA Astrophysics Data System (ADS)
Mayoral-Villa, Estela; Alvarado-Rodríguez, Carlos E.; Klapp, Jaime; Gómez-Gesteira, Moncho; Di G. Sigalotti, Leonardo
2016-04-01
A smoothed particle hydrodynamics (SPH) model is presented for simulating the decay chain transport of radionuclides in confined aqueous solutions. The SPH formulation is based on the open-source parallel code DualSPHysics extended to solve the advective-diffusion equation for the evolution of the concentration field coupled to the fluid-dynamic equations, including the effects of radioactive decay of the tracer contaminants. The performance of the method is demonstrated for environmental engineering problems dealing with the transport of contaminants in still and flowing water. The results from a series of benchmark test calculations are described in two- and three-space dimensions, where the advection, diffusion, and radioactive decay modes are tested separately and in combined form. The accuracy of the present SPH transport model is shown by direct comparison with the analytical solutions and results from other SPH approaches. For a given problem, convergence of the SPH solution is seen to increase with decreasing particle size and spacing.
Pu, Jaan Hui; Shao, Songdong
2012-01-01
This research paper presents an incompressible smoothed particle hydrodynamics (ISPH) technique to investigate a regular wave overtopping on the coastal structure of different types. The SPH method is a mesh-free particle modeling approach that can efficiently treat the large deformation of free surface. The incompressible SPH approach employs a true hydrodynamic formulation to solve the fluid pressure that has less pressure fluctuations. The generation of flow turbulence during the wave breaking and overtopping is modeled by a subparticle scale (SPS) turbulence model. Here the ISPH model is used to investigate the wave overtopping over a coastal structure with and without the porous material. The computations disclosed the features of flow velocity, turbulence, and pressure distributions for different structure types and indicated that the existence of a layer of porous material can effectively reduce the wave impact pressure and overtopping rate. The proposed numerical model is expected to provide a promising practical tool to investigate the complicated wave-structure interactions. PMID:22919291
Developing a weakly compressible smoothed particle hydrodynamics model for biological flows
NASA Astrophysics Data System (ADS)
Vasyliv, Yaroslav; Alexeev, Alexander
2014-11-01
Smoothed Particle Hydrodynamics (SPH) is a meshless particle method originally developed for astrophysics applications in 1977. Over the years, limitations of the original formulations have been addressed by different groups to extend the domain of SPH application. In biologically relevant internal flows, two of the several challenges still facing SPH are 1) treatment of inlet, outlet, and no slip boundary conditions and 2) treatment of second derivatives present in the viscous terms. In this work, we develop a 2D weakly compressible SPH (WCSPH) for simulating viscous internal flows which incorporates some of the recent advancements made by groups in the above two areas. The method is validated against several analytical and experimental benchmark solutions for both steady and unsteady laminar flows. In particular, the 2013 U.S. Food and Drug Administration benchmark test case for medical devices - steady forward flow through a nozzle with a sudden contraction and conical diffuser - is simulated for different Reynolds numbers in the laminar region and results are validated against the published experimental and CFD datasets. Support from the National Science Foundation Graduate Research Fellowship Program (NSF GRFP) is gratefully acknowledged.
NASA Astrophysics Data System (ADS)
Smith, Scott C.; Houser, Janet L.; Centrella, Joan M.
1996-02-01
We have carried out three-dimensional numerical simulations of the dynamical bar instability in a rotating star and the resulting gravitational radiation using both an Eulerian code written in cylindrical coordinates and a smooth particle hydrodynamics (SPH) code. The star is modeled initially as a polytrope with index n = 3/2 and Trot|W| ≍ 0.30, where Trot is the rotational kinetic energy and |W| is the gravitational potential energy. In both codes the gravitational field is purely Newtonian, and the gravitational radiation is calculated in the quadrupole approximation. We have run three simulations with the Eulerian code, varying the number of angular zones and the treatment of the boundary between the star and the vacuum. Using the SPH code we did seven runs, varying the number of particles, the artificial viscosity, and the type of initial model. We compare the growth rate and rotation speed of the bar, the mass, and angular momentum distributions, and the gravitational radiation quantities. We highlight the successes and difficulties of both methods and make suggestions for future improvements.
A smooth dissipative particle dynamics method for domains with arbitrary-geometry solid boundaries
NASA Astrophysics Data System (ADS)
Gatsonis, Nikolaos A.; Potami, Raffaele; Yang, Jun
2014-01-01
A smooth dissipative particle dynamics method with dynamic virtual particle allocation (SDPD-DV) for modeling and simulation of mesoscopic fluids in wall-bounded domains is presented. The physical domain in SDPD-DV may contain external and internal solid boundaries of arbitrary geometries, periodic inlets and outlets, and the fluid region. The SDPD-DV method is realized with fluid particles, boundary particles, and dynamically allocated virtual particles. The internal or external solid boundaries of the domain can be of arbitrary geometry and are discretized with a surface grid. These boundaries are represented by boundary particles with assigned properties. The fluid domain is discretized with fluid particles of constant mass and variable volume. Conservative and dissipative force models due to virtual particles exerted on a fluid particle in the proximity of a solid boundary supplement the original SDPD formulation. The dynamic virtual particle allocation approach provides the density and the forces due to virtual particles. The integration of the SDPD equations is accomplished with a velocity-Verlet algorithm for the momentum and a Runge-Kutta for the entropy equation. The velocity integrator is supplemented by a bounce-forward algorithm in cases where the virtual particle force model is not able to prevent particle penetration. For the incompressible isothermal systems considered in this work, the pressure of a fluid particle is obtained by an artificial compressibility formulation for liquids and the ideal gas law for gases. The self-diffusion coefficient is obtained by an implementation of the generalized Einstein and the Green-Kubo relations. Field properties are obtained by sampling SDPD-DV outputs on a post-processing grid that allows harnessing the particle information on desired spatiotemporal scales. The SDPD-DV method is verified and validated with simulations in bounded and periodic domains that cover the hydrodynamic and mesoscopic regimes for
Plasma cell adaptation to enhance particle acceleration
Ragheb, M. S.
2008-06-15
A plasma study is performed in order to construct a cell for plasma acceleration purpose. As well, a multicell design is introduced for the injection of beam driver application. The suggested idea is experimentally demonstrated for two plasma cell configuration. The preformed plasma is obtained by a symmetrically driven capacitive audio frequency discharge. It is featured by its moderate pressure of 0.1-0.2 Torr, low consumption power of 130 W maximum, low discharge voltage and frequency up to 950 V and 20 kHz, respectively, and high plasma density from 10{sup 11} to 10{sup 15} cm{sup -3}. The electron temperature obtained by Langmuir double probe varies from 1 up to 16 eV. It is observed that the increases of the discharge voltage and frequency enlarge the plasma parameters to their maximum values. The plasma cell filled with different gases demonstrates that the Ar and He gases manifest the highest ionization efficiency exceeding 100% at 950 V and 20 kHz. The formed plasma is cold; its density is uniform and stable along the positive column for long competitive lifetime. Showing that it follows the conditions to enhance particle acceleration and in conjunction with its periphery devices form a plasma cell that could be extended to serve this purpose. Demonstrating that an injected electron beam into the extended preformed plasma could follow, to long distance, a continuous trajectory of uniform density. Such plasma generated by H{sub 2} or Ar gases is suggested to be used, respectively, for low-density or higher density beam driver.
Kordilla, Jannes; Pan, Wenxiao; Tartakovsky, Alexandre M.
2014-12-14
We propose a novel Smoothed Particle Hydrodynamics (SPH) discretization of the fully-coupled Landau-Lifshitz-Navier-Stokes (LLNS) and advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations are found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for the coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study the formation of the so-called giant fluctuations of the front between light and heavy fluids with and without gravity, where the light fluid lays on the top of the heavy fluid. We find that the power spectra of the simulated concentration field is in good agreement with the experiments and analytical solutions. In the absence of gravity the the power spectra decays as the power -4 of the wave number except for small wave numbers which diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations resulting in the much weaker dependence of the power spectra on the wave number. Finally the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlying a light fluid. The front dynamics is shown to agree well with the analytical solutions.
Simulation of film and droplet flow on wide aperture fracture using Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Kordilla, J.; Tartakovsky, A.; Geyer, T.
2012-04-01
Fractured media provide rapid flow pathways for water percolating through the unsaturated zone. Film flow has been widely acknowledged as a major rapid flow process with average velocities of 3x10-7m/s (Tokunaga, 1997). Further flow regimes such as droplets, rivulets and falling films may reach much higher velocities while coexisting with films (Ghezzehei, 2004). In order to establish a unified description of multiphase flow at such small scales simulation approaches have to be able to deal with the highly dynamical interfaces and reproduce the physical behavior dominated by capillary, surface tension and gravitational forces. In this work we show simulations of free-surface flow on inclined fracture surfaces using a Smoothed Particle Hydrodynamics (SPH) model (Tartakovsky, 2005). The three-dimensional Lagrangian code employs an interpolation kernel in order to solve the Navier-Stokes equation at an arbitrary set of points (particles). Pairwise fluid-fluid and solid-fluid interaction forces are used to simulate a wide range of wetting conditions and Reynolds numbers encountered in laboratory experiments. Model results are verified with empirical and semianalytical solutions. Contact angles of droplets in a critical state, i.e. at the verge of movement, are compared with laboratory experiments reported in literature. Transient droplet dynamics are shown to reproduce the linear scaling proposed by Podgorski (2001). Depending on Reynolds number and static contact angles droplets leave behind trailing films. In order to investigate the influence of adsorbed films on droplet flow surfaces are prewetted with a thin film and simulations repeated. The results indicate a strong dependence of droplet flow dynamics on the existence of adsorbed films with droplet velocities being tripled under certain conditions. Despite their relatively slow velocities, adsorbed films are shown to be an essential part of unsaturated droplet flow dynamics as they enhance the wetting and
NASA Astrophysics Data System (ADS)
Kordilla, Jannes; Pan, Wenxiao; Tartakovsky, Alexandre
2014-12-01
We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called "giant fluctuations" of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power -4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.
Kordilla, Jannes; Pan, Wenxiao Tartakovsky, Alexandre
2014-12-14
We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called “giant fluctuations” of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power −4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.
Smoothed particle hydrodynamics non-Newtonian model for ice-sheet and ice-shelf dynamics
Pan, W.; Tartakovsky, A. M.; Monaghan, J. J.
2013-06-01
Mathematical modeling of ice sheets is complicated by the non-linearity of the governing equations and boundary conditions. Standard grid-based methods require complex front tracking techniques and have limited capability to handle large material deformations and abrupt changes in bottom topography. As a consequence, numerical methods are usually restricted to shallow ice sheet and ice shelf approximations. We propose a new smoothed particle hydrodynamics (SPH) non-Newtonian model for coupled ice sheet and ice shelf dynamics. SPH, a fully Lagrangian particle method, is highly scalable and its Lagrangian nature and meshless discretization are well suited to the simulation of free surface flows, large material deformation, and material fragmentation. In this paper, SPH is used to study 3D ice sheet/ice shelf behavior, and the dynamics of the grounding line. The steady state position of the grounding line obtained from SPH simulations is in good agreement with laboratory observations for a wide range of simulated bedrock slopes, and density ratios, similar to those of ice and sea water. The numerical accuracy of the SPH algorithm is verif;ed by simulating Poiseuille flow, plane shear flow with free surface and the propagation of a blob of ice along a horizontal surface. In the laboratory experiment, the ice was represented with a viscous Newtonian fluid. In the present work, however, the ice is modeled as both viscous Newtonian fluid and non-Newtonian fluid, such that the effect of non-Newtonian rheology on the dynamics of grounding line was examined. The non-Newtonian constitutive relation is prescribed to be Glen’s law for the creep of polycrystalline ice. A V-shaped bedrock ramp is further introduced to model the real geometry of bedrock slope.
A Smoothed Particle Hydrodynamics Model for Ice Sheet and Ice Shelf Dynamics
Pan, Wenxiao; Tartakovsky, Alexandre M.; Monaghan, Joseph J.
2012-02-08
Mathematical modeling of ice sheets is complicated by the non-linearity of the governing equations and boundary conditions. Standard grid-based methods require complex front tracking techniques and have limited capability to handle large material deformations and abrupt changes in bottom topography. As a consequence, numerical methods are usually restricted to shallow ice sheet and ice shelf approximations. We propose a new smoothed particle hydrodynamics (SPH) model for coupled ice sheet and ice shelf dynamics. SPH is a fully Lagrangian particle method. It is highly scalable and its Lagrangian nature and meshless discretization are well suited to the simulation of free surface flows, large material deformation, and material fragmentation. In this paper SPH is used to study ice sheet/ice shelf behavior, and the dynamics of the grounding line. The steady state position of the grounding line obtained from the SPH simulations is in good agreement with laboratory observations for a wide range of simulated bedrock slopes, and density ratios similar to those of ice and sea water. The numerical accuracy of the SPH algorithm is further verified by simulating the plane shear flow of two immiscible fluids and the propagation of a highly viscous blob of fluid along a horizontal surface. In the experiment, the ice was represented with a viscous newtonian fluid. For consistency, in the described SPH model the ice is also modeled as a viscous newtonian fluid. Typically, ice sheets are modeled as a non-Newtonian fluid, accounting for the changes in the mechanical properties of ice. Implementation of a non-Newtonian rheology in the SPH model is the subject of our ongoing research.
NASA Astrophysics Data System (ADS)
Khorasanizade, Sh.; Sousa, J. M. M.
2016-03-01
A Segmented Boundary Algorithm (SBA) is proposed to deal with complex boundaries and moving bodies in Smoothed Particle Hydrodynamics (SPH). Boundaries are formed in this algorithm with chains of lines obtained from the decomposition of two-dimensional objects, based on simple line geometry. Various two-dimensional, viscous fluid flow cases have been studied here using a truly incompressible SPH method with the aim of assessing the capabilities of the SBA. Firstly, the flow over a stationary circular cylinder in a plane channel was analyzed at steady and unsteady regimes, for a single value of blockage ratio. Subsequently, the flow produced by a moving circular cylinder with a prescribed acceleration inside a plane channel was investigated as well. Next, the simulation of the flow generated by the impulsive start of a flat plate, again inside a plane channel, has been carried out. This was followed by the study of confined sedimentation of an elliptic body subjected to gravity, for various density ratios. The set of test cases was completed with the simulation of periodic flow around a sunflower-shaped object. Extensive comparisons of the results obtained here with published data have demonstrated the accuracy and effectiveness of the proposed algorithms, namely in cases involving complex geometries and moving bodies.
The effect of density estimation on the conservativeness in Smoothed Particle Hydrodynamics
NASA Astrophysics Data System (ADS)
Suresh, Pranav; Kumar, S. S. Prasanna; Patnaik, B. S. V.
2015-11-01
Smoothed Particle Hydrodynamics (SPH) is a popular mesh-free method for solving a wide range of problems that involve interfaces. In SPH, the Lagrangian nature of the method enables mass conservation to be naturally satisfied. However, satisfying the conservation of momentum and energy are indeed formulation dependent. One major aspect of ensuring conservativeness comes from the density estimation. There are two distinct types of density estimation approaches, namely continuity density approach and summation density approach. Both approaches are indeed popular with single and multi-phase flow communities. In the present study, we assess the role of density evaluation on the conservativeness, using several representative numerical examples. In particular, we have simulated the Rayleigh-Taylor instability problem, Non-Boussinesq lock exchange problem, bubble rise in water column etc. Although for shorter time scales of simulation, both methods have similar conservative properties, we observe that for longer time scales, summation-density approach is better. For free surface detection and normal vector computations, efficient computational procedures have been devised.
NASA Astrophysics Data System (ADS)
Dyachkov, Sergey; Parshikov, Anatoly; Zhakhovsky, Vasily
2015-06-01
The machining of materials produces regular micrometer-sized surface perturbations. The microscopic cumulative jets can be generated from such surface under shock loading. It is a problem to trace space-time evolution of such jets with good enough resolution in experimental conditions. Comparative simulations by molecular dynamics (MD) and smoothed-particle hydrodynamics (SPH) methods, using an equation of state consistent with the employed interatomic potential, can shed of light on details of jet formation. The realistic experimental samples can be directly simulated by SPH method, while the linear size of a MD sample is restricted by the order of 100 nm. To compare the SPH and MD simulations the MD results must to be scaled to micrometer-sized samples. We demonstrate that the scaling provides the similar jet velocity profiles and mass distributions obtained by both methods. Furthermore, the simulated results agree well with the experimental observations with Copper and Tin. The effect of surface tension, which guides evolution of nanoscale-sized jet shape, may lead to discrepancies between MD and SPH simulations, especially for weak shocks and small surface perturbations.
NASA Astrophysics Data System (ADS)
Nishiura, Daisuke; Furuichi, Mikito; Sakaguchi, Hide
2015-09-01
The computational performance of a smoothed particle hydrodynamics (SPH) simulation is investigated for three types of current shared-memory parallel computer devices: many integrated core (MIC) processors, graphics processing units (GPUs), and multi-core CPUs. We are especially interested in efficient shared-memory allocation methods for each chipset, because the efficient data access patterns differ between compute unified device architecture (CUDA) programming for GPUs and OpenMP programming for MIC processors and multi-core CPUs. We first introduce several parallel implementation techniques for the SPH code, and then examine these on our target computer architectures to determine the most effective algorithms for each processor unit. In addition, we evaluate the effective computing performance and power efficiency of the SPH simulation on each architecture, as these are critical metrics for overall performance in a multi-device environment. In our benchmark test, the GPU is found to produce the best arithmetic performance as a standalone device unit, and gives the most efficient power consumption. The multi-core CPU obtains the most effective computing performance. The computational speed of the MIC processor on Xeon Phi approached that of two Xeon CPUs. This indicates that using MICs is an attractive choice for existing SPH codes on multi-core CPUs parallelized by OpenMP, as it gains computational acceleration without the need for significant changes to the source code.
NASA Astrophysics Data System (ADS)
Dyachkov, S. A.; Parshikov, A. N.; Zhakhovsky, V. V.
2015-11-01
Experimental methods of observation of early stage of shock-induced ejecta from metal surface with micrometer-sized perturbations are still limited in terms of following a complete sequence of processes having microscale dimensions and nanoscale times. Therefore, simulations by the smoothed particle hydrodynamics (SPH) and molecular dynamics (MD) methods can shed of light on details of micro-jet evolution. The size of simulated sample is too restricted in MD, but the simulations with large enough number of atoms can be scaled well to the sizes of realistic samples. To validate such scaling the comparative MD and SPH simulations of tin samples are performed. SPH simulation takes the realistic experimental sizes, while MD uses the proportionally scaled sizes of samples. It is shown that the velocity and mass distributions along the jets simulated by MD and SPH are in a good agreement. The observed difference in velocity of spikes between MD and experiments can be partially explained by a profound effect of surface tension on jets ejected from the small-scale samples.
A smooth particle hydrodynamics code to model collisions between solid, self-gravitating objects
NASA Astrophysics Data System (ADS)
Schäfer, C.; Riecker, S.; Maindl, T. I.; Speith, R.; Scherrer, S.; Kley, W.
2016-05-01
Context. Modern graphics processing units (GPUs) lead to a major increase in the performance of the computation of astrophysical simulations. Owing to the different nature of GPU architecture compared to traditional central processing units (CPUs) such as x86 architecture, existing numerical codes cannot be easily migrated to run on GPU. Here, we present a new implementation of the numerical method smooth particle hydrodynamics (SPH) using CUDA and the first astrophysical application of the new code: the collision between Ceres-sized objects. Aims: The new code allows for a tremendous increase in speed of astrophysical simulations with SPH and self-gravity at low costs for new hardware. Methods: We have implemented the SPH equations to model gas, liquids and elastic, and plastic solid bodies and added a fragmentation model for brittle materials. Self-gravity may be optionally included in the simulations and is treated by the use of a Barnes-Hut tree. Results: We find an impressive performance gain using NVIDIA consumer devices compared to our existing OpenMP code. The new code is freely available to the community upon request. If you are interested in our CUDA SPH code miluphCUDA, please write an email to Christoph Schäfer. miluphCUDA is the CUDA port of miluph. miluph is pronounced [maßl2v]. We do not support the use of the code for military purposes.
How smooth are particle trajectories in a ΛCDM Universe?
NASA Astrophysics Data System (ADS)
Rampf, Cornelius; Villone, Barbara; Frisch, Uriel
2015-09-01
It is shown here that in a flat, cold dark matter (CDM)-dominated Universe with positive cosmological constant (Λ), modelled in terms of a Newtonian and collisionless fluid, particle trajectories are analytical in time (representable by a convergent Taylor series) until at least a finite time after decoupling. The time variable used for this statement is the cosmic scale factor, i.e. the `a-time', and not the cosmic time. For this, a Lagrangian-coordinate formulation of the Euler-Poisson equations is employed, originally used by Cauchy for 3D incompressible flow. Temporal analyticity for ΛCDM is found to be a consequence of novel explicit all-order recursion relations for the a-time Taylor coefficients of the Lagrangian displacement field, from which we derive the convergence of the a-time Taylor series. A lower bound for the a-time where analyticity is guaranteed and shell-crossing is ruled out is obtained, whose value depends only on Λ and on the initial spatial smoothness of the density field. The largest time interval is achieved when Λ vanishes, i.e. for an Einstein-de Sitter universe. Analyticity holds also if, instead of the a-time, one uses the linear structure growth D-time, but no simple recursion relations are then obtained. The analyticity result also holds when a curvature term is included in the Friedmann equation for the background, but inclusion of a radiation term arising from the primordial era spoils analyticity.
NASA Astrophysics Data System (ADS)
Laibe, Guillaume; Price, Daniel J.
2012-03-01
In a companion paper, we have presented an algorithm for simulating two-fluid gas and dust mixtures in smoothed particle hydrodynamics (SPH). In this paper, we develop an implicit timestepping method that preserves the exact conservation of the both linear and angular momenta in the underlying SPH algorithm, but unlike previous schemes, allows the iterations to converge to arbitrary accuracy and is suited to the treatment of non-linear drag regimes. The algorithm presented in Paper I is also extended to deal with realistic astrophysical drag regimes, including both linear and non-linear Epstein and Stokes drag. The scheme is benchmarked against the test suite presented in Paper I, including (i) the analytic solutions of the DUSTYBOX problem and (ii) solutions of the DUSTYWAVE, DUSTYSHOCK, DUSTYSEDOV and DUSTYDISC obtained with explicit timestepping. We find that the implicit method is 1-10 times faster than the explicit temporal integration when the ratio r between the timestep and the drag stopping time is 1 ≲r≲ 1000.
NASA Astrophysics Data System (ADS)
Braun, Anika; Cuomo, Sabatino; Wang, Xueliang; Zhang, Luqing
2016-04-01
Debris flows and landslide dams are a major natural hazard causing high socioeconomic risk in inhabited mountainous areas. This is also true for vast parts of southwestern China, which are highly prone to slope failures due to several factors, such as a humid climate with high precipitation in the summer months, geological predisposing factors with highly weathered sedimentary rocks and a high seismicity. Not only do the landslides and flooding related to landslide dams threaten residents, buildings and transportation structures, but also do flooding conditions endanger power supply, which relies in this region partly on hydropower. In order to assess the potential of landslides to block rivers, it is crucial to understand which factors influence possible run-out distances and how they can be quantified. In the study we are presenting a numerical modeling analysis for a particular case of a complex landslide in Ningnan county, southwestern China, which transformed into a debris flow and blocked the river and the major road leading through the valley after heavy rainfall. For this purpose a quasi-3D Smooth Particle Hydrodynamics (SPH) model was implemented that can account for geotechnical slope parameters, run-out distance, velocities and deposition heights. A digital terrain model and the geometry information of the landslide were used as input data in order to estimate the relevant geotechnical parameters by back-analysis. The results of the analysis can be used for the prediction of run-out distances for future events at this site and other similar sites.
On the feasibility of using smoothed particle hydrodynamics for underwater explosion calculations
Swegle, J.W.; Attaway, S.W.
1995-02-01
SPH (Smoothed Particle Hydrodynamics) is a gridless Lagrangian technique which is appealing as a possible alternative to numerical techniques currently used to analyze high deformation impulsive loading events. In the present study, the SPH algorithm has been subjected to detailed testing and analysis to determine the feasibility of using PRONTO/SPH for the analysis of various types of underwater explosion problems involving fluid-structure and shock-structure interactions. Of particular interest are effects of bubble formation and collapse and the permanent deformation of thin walled structures due to these loadings. These are exceptionally difficult problems to model. Past attempts with various types of codes have not been satisfactory. Coupling SPH into the finite element code PRONTO represents a new approach to the problem. Results show that the method is well-suited for transmission of loads from underwater explosions to nearby structures, but the calculation of late time effects due to acceleration of gravity and bubble buoyancy will require additional development, and possibly coupling with implicit or incompressible methods.
A smooth particle hydrodynamics code to model collisions between solid, self-gravitating objects
NASA Astrophysics Data System (ADS)
Schäfer, C.; Riecker, S.; Maindl, T. I.; Speith, R.; Scherrer, S.; Kley, W.
2016-04-01
Context. Modern graphics processing units (GPUs) lead to a major increase in the performance of the computation of astrophysical simulations. Owing to the different nature of GPU architecture compared to traditional central processing units (CPUs) such as x86 architecture, existing numerical codes cannot be easily migrated to run on GPU. Here, we present a new implementation of the numerical method smooth particle hydrodynamics (SPH) using CUDA and the first astrophysical application of the new code: the collision between Ceres-sized objects. Aims: The new code allows for a tremendous increase in speed of astrophysical simulations with SPH and self-gravity at low costs for new hardware. Methods: We have implemented the SPH equations to model gas, liquids and elastic, and plastic solid bodies and added a fragmentation model for brittle materials. Self-gravity may be optionally included in the simulations and is treated by the use of a Barnes-Hut tree. Results: We find an impressive performance gain using NVIDIA consumer devices compared to our existing OpenMP code. The new code is freely available to the community upon request. If you are interested in our CUDA SPH code miluphCUDA, please write an email to Christoph Schäfer. miluphCUDA is the CUDA port of miluph. miluph is pronounced [maßl2v]. We do not support the use of the code for military purposes.
Geometrical on-the-fly shock detection in smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Beck, A. M.; Dolag, K.; Donnert, J. M. F.
2016-05-01
We present an on-the-fly geometrical approach for shock detection and Mach number calculation in simulations employing smoothed particle hydrodynamics (SPH). We utilize pressure gradients to select shock candidates and define up- and downstream positions. We obtain hydrodynamical states in the up- and downstream regimes with a series of normal and inverted kernel weightings parallel and perpendicular to the shock normals. Our on-the-fly geometrical Mach detector incorporates well within the SPH formalism and has low computational cost. We implement our Mach detector into the simulation code GADGET and alongside many SPH improvements. We test our shock finder in a sequence of shock tube tests with successively increasing Mach numbers exceeding by far the typical values inside galaxy clusters. For all shocks, we resolve the shocks well and the correct Mach numbers are assigned. An application to a strong magnetized shock tube gives stable results in full magnetohydrodynamic setups. We simulate a merger of two idealized galaxy clusters and study the shock front. Shock structures within the merging clusters as well as the cluster shock are well captured by our algorithm and assigned correct Mach numbers.
MASS TRANSFER IN BINARY STARS USING SMOOTHED PARTICLE HYDRODYNAMICS. II. ECCENTRIC BINARIES
Lajoie, Charles-Philippe; Sills, Alison E-mail: asills@mcmaster.ca
2011-01-10
Despite numerous efforts to better understand binary star evolution, some aspects of it remain poorly constrained. In particular, the evolution of eccentric binaries has remained elusive mainly because the Roche lobe formalism derived for circular binaries does not apply. Here we report the results of our smoothed particle hydrodynamic simulations of mass transfer in eccentric binaries using an alternate method in which we model only the outermost layers of the stars with appropriate boundary conditions. Using this technique, along with properly relaxed model stars, we characterize the mass transfer episodes of binaries with various orbital parameters. In particular, we show that these episodes can be described by Gaussians with an FWHM of {approx}0.12P{sub orb} and that the peak rates occur after periastron, at an orbital phase of {approx}0.58, independently of the eccentricity and mass of the stars. The accreted material is observed to form a rather sparse envelope around either or both stars. Although the fate of this envelope is not modeled in our simulations, we show that a constant fraction ({approx}5%) of the material transferred is ejected from the systems. We discuss this result in terms of the non-conservative mass transfer scenario. We suggest that our results could be incorporated in analytical and binary population synthesis studies to help better understand the evolution of eccentric binaries and the formation of exotic stellar populations.
Lai, Guanyu; Liu, Zhi; Zhang, Yun; Philip Chen, C L
2016-06-01
This paper is concentrated on the problem of adaptive fuzzy tracking control for an uncertain nonlinear system whose actuator is encountered by the asymmetric backlash behavior. First, we propose a new smooth inverse model which can approximate the asymmetric actuator backlash arbitrarily. By applying it, two adaptive fuzzy control scenarios, namely, the compensation-based control scheme and nonlinear decomposition-based control scheme, are then developed successively. It is worth noticing that the first fuzzy controller exhibits a better tracking control performance, although it recourses to a known slope ratio of backlash nonlinearity. The second one further removes the restriction, and also gets a desirable control performance. By the strict Lyapunov argument, both adaptive fuzzy controllers guarantee that the output tracking error is convergent to an adjustable region of zero asymptotically, while all the signals remain semiglobally uniformly ultimately bounded. Lastly, two comparative simulations are conducted to verify the effectiveness of the proposed fuzzy controllers. PMID:27187937
Adaptive Beam Smoothing with Plasma-Pinholes for Laser-Entrance-Hole Transmission Studies
NASA Astrophysics Data System (ADS)
Geissel, Matthias; Ruggles, Lawrence E.; Smith, Ian C.; Shores, Jonathn E.; Speas, C. Shane; Porter, John L.
2014-10-01
The concept of Magnetized Liner Inertial Fusion (MagLIF) requires the deposition of laser energy into a fuel-filled cylinder that is exposed to a magnetic field. To improve process, it is essential to optimize transmission through the foil covered laser entrance hole (LEH), which involves minimizing laser-plasma-instabilities (LPI). Laser beam smoothing is the most common approach to minimize LPI. It typically involves a Random-Phase-Plate (RPP) and smoothing by spectral dispersion (SSD). This approach can still cause LPI issues due to intensity ``hot-spots'' on a ps-time scale, and it inconveniently fixes the usable spot size. Changing laser spot sizes requires multiple dedicated RPPs. To study ideal spot sizes on a MagLIF LEH, the RPP/SSD approach gets cost prohibitive. As alternative, we use sacrificial thin foils (500 nm or less) at the laser focus, which instantly turn into a plasma-pinhole, acting as spatial filter. The smoothed laser spot size grows linearly with distance from best focus. We present experimental data for smoothing performance and resulting LEH transmission. Sandia is a multiprogram laboratory managed and operated by Sandia Corp., a wholly owned subsidiary of Lockheed Martin Corp., for the U.S. DOE's Nat'l Nucl. Sec. Admin. under Contract DE-AC04-94AL85000.
Galaxy mergers on a moving mesh: a comparison with smoothed particle hydrodynamics
NASA Astrophysics Data System (ADS)
Hayward, Christopher C.; Torrey, Paul; Springel, Volker; Hernquist, Lars; Vogelsberger, Mark
2014-08-01
Galaxy mergers have been investigated for decades using smoothed particle hydrodynamics (SPH), but recent work highlighting inaccuracies inherent in the traditional SPH technique calls into question the reliability of previous studies. We explore this issue by comparing a suite of GADGET-3 SPH simulations of idealized (i.e. non-cosmological) isolated discs and galaxy mergers with otherwise identical calculations performed using the moving-mesh code AREPO. When black hole (BH) accretion and active galactic nucleus (AGN) feedback are not included, the star formation histories (SFHs) obtained from the two codes agree well. When BHs are included, the code- and resolution-dependent variations in the SFHs are more significant, but the agreement is still good, and the stellar mass formed over the course of a simulation is robust to variations in the numerical method. During a merger, the gas morphology and phase structure are initially similar prior to the starburst phase. However, once a hot gaseous halo has formed from shock heating and AGN feedback (when included), the agreement is less good. In particular, during the post-starburst phase, the SPH simulations feature more prominent hot gaseous haloes and spurious clumps, whereas with AREPO, gas clumps and filaments are less apparent and the hot halo gas can cool more efficiently. We discuss the origin of these differences and explain why the SPH technique yields trustworthy results for some applications (such as the idealized isolated disc and galaxy merger simulations presented here) but not others (e.g. gas flows on to galaxies in cosmological hydrodynamical simulations).
NASA Astrophysics Data System (ADS)
Thomas, David M.; Wood, Matt A.
2015-04-01
Negative superhumps are believed to arise in cataclysmic variable systems when the accretion disk is tilted with respect to the orbital plane. Slow retrograde precession of the line-of-nodes results in a signal—the negative superhump—with a period slightly less than the orbital period. Previous studies have shown that tilted disks exhibit negative superhumps, but a consensus on how a disk initially tilts has not been reached. Analytical work by Lai (1999, ApJ, 524, 1030) suggests that a magnetic field on the primary can lead to a tilt instability in a disk when the dipole moment is offset in angle from the spin axis of the primary and when the primary’s spin axis is, itself, not aligned with the angular momentum axis of the binary orbit. However, Lai did not apply his work to the formation of negative superhumps. In this paper, we add Lai’s model to an existing smoothed particle hydrodynamics code. Using this code, we demonstrate the emergence of negative superhumps in the “light curve” for a range of magnetic dipole moments. We show that the period deficits calculated from these negative superhumps match those in simulations using manually tilted disks. When positive superhumps appear (q≲ 0.33), we show that the period excesses calculated from these signals are also consistent with previous results. Using examples, we show that the disks are tilted, though the tilt varies periodically, and that they precess in the retrograde direction. The magnetic fields found to lead to the emergence of negative superhumps lie in the kilogauss regime.
Smoothed Particle Hydrodynamics Modeling of Gravity Currents on a Dry Porous Medium
NASA Astrophysics Data System (ADS)
Daly, E.; Grimaldi, S.; Bui, H.
2014-12-01
Gravity currents flowing over porous media occur in many environmental processes and industrial applications, such as irrigation, benthic boundary layers, and oil spills. The coupling of the flow over the porous surface and the infiltration of the fluid in the porous media is complex and difficult to model. Of particular interest is the prediction of the position of the runoff front and the depth of the infiltration front. We present here a model for the flow of a finite volume of a highly viscous Newtonian fluid over a dry, homogenous porous medium. The Navier-Stokes equations describing the runoff flow are coupled to the Volume Averaged Navier-Stokes equations for the infiltration flow. The numerical solution of these equations is challenging because of the presence of two free surfaces (runoff and infiltration waves), the lack of fixed boundary conditions at the runoff front, and the difficulties in defining appropriate conditions at the surface of the porous medium. The first two challenges were addressed by using Smoothed Particle Hydrodynamics, which is a Lagrangian, mesh-free particle method particularly suitable for modelling free surface flows. Two different approaches were used to model the flow conditions at the surface of the porous medium. The Two Domain Approach (TDA) assumes that runoff and infiltration flows occur in two separate homogenous domains; here, we assume the continuity of velocity and stresses at the interface of the two domains. The One Domain Approach (ODA) models runoff and infiltration flows as occurring through a medium whose hydraulic properties vary continuously in space. The transition from the hydraulic properties of the atmosphere and the porous medium occur in a layer near the surface of the porous medium. Expressions listed in literature were used to compute the thickness of this transition layer and the spatial variation of porosity and permeability within it. Our results showed that ODA led to slower velocities of the runoff
Hydrodynamics in adaptive resolution particle simulations: Multiparticle collision dynamics
NASA Astrophysics Data System (ADS)
Alekseeva, Uliana; Winkler, Roland G.; Sutmann, Godehard
2016-06-01
A new adaptive resolution technique for particle-based multi-level simulations of fluids is presented. In the approach, the representation of fluid and solvent particles is changed on the fly between an atomistic and a coarse-grained description. The present approach is based on a hybrid coupling of the multiparticle collision dynamics (MPC) method and molecular dynamics (MD), thereby coupling stochastic and deterministic particle-based methods. Hydrodynamics is examined by calculating velocity and current correlation functions for various mixed and coupled systems. We demonstrate that hydrodynamic properties of the mixed fluid are conserved by a suitable coupling of the two particle methods, and that the simulation results agree well with theoretical expectations.
NASA Astrophysics Data System (ADS)
Seddik, Hassene
2014-12-01
Noise can occur during image capture, transmission, or processing phases. Image de-noising is a very important step in image processing, and many approaches are developed in order to achieve this goal such as the Gaussian filter which is efficient in noise removal. Its smoothing efficiency depends on the value of its standard deviation. The mask representing the filter presents generally static weights with invariant lobe. In this paper, an adaptive de-noising approach is proposed. The proposed approach uses a Gaussian kernel with variable width and direction called adaptive Gaussian kernel (AGK). In each processed window of the image, the smoothing strength changes according to the image content, noise kind, and intensity. In addition, the location of its lobe changes in eight different directions over the processed window. This directional variability avoids averaging details by the highest mask weights in order to preserve the edges and the borders. The recovered data is de-noised efficiently without introducing blur or losing details. A comparative study with the static Gaussian filter and other recent techniques is presented to prove the efficiency of the proposed approach.
NASA Astrophysics Data System (ADS)
Pan, Wenxiao; Bao, Jie; Tartakovsky, Alexandre
2013-11-01
A Continuous Boundary Force (CBF) method was developed for implementing Robin (Navier) boundary condition (BC) that can describe no-slip or slip conditions (slip length from zero to infinity) at the fluid-solid interface. In the CBF method the Robin BC is replaced by a homogeneous Neumann BC and an additional volumetric source term in the governing momentum equation. The formulation is derived based on an approximation of the sharp boundary with a diffuse interface of finite thickness, across which the BC is reformulated by means of a smoothed characteristic function. The CBF method is easy to be implemented in Lagrangian particle-based methods. We first implemented it in smoothed particle hydrodynamics (SPH) to solve numerically the Navier-Stokes equations subject to spatial-independent or dependent Robin BC in two and three dimensions. The numerical accuracy and convergence is examined through comparisons with the corresponding finite difference or finite element solutions. The CBF method is further implemented in smoothed dissipative particle dynamics (SDPD), a mesoscale scheme, for modeling slip flows commonly existent in micro/nano channels and microfluidic devices. The authors acknowledge the funding support by the ASCR Program of the Office of Science, U.S. Department of Energy.
NASA Astrophysics Data System (ADS)
Zhang, X.; Liang, S.; Wang, G.
2015-12-01
Incident solar radiation (ISR) over the Earth's surface plays an important role in determining the Earth's climate and environment. Generally, can be obtained from direct measurements, remotely sensed data, or reanalysis and general circulation models (GCMs) data. Each type of product has advantages and limitations: the surface direct measurements provide accurate but sparse spatial coverage, whereas other global products may have large uncertainties. Ground measurements have been normally used for validation and occasionally calibration, but transforming their "true values" spatially to improve the satellite products is still a new and challenging topic. In this study, an improved thin-plate smoothing spline approach is presented to locally "calibrate" the Global LAnd Surface Satellite (GLASS) ISR product using the reconstructed ISR data from surface meteorological measurements. The influences of surface elevation on ISR estimation was also considered in the proposed method. The point-based surface reconstructed ISR was used as the response variable, and the GLASS ISR product and the surface elevation data at the corresponding locations as explanatory variables to train the thin plate spline model. We evaluated the performance of the approach using the cross-validation method at both daily and monthly time scales over China. We also evaluated estimated ISR based on the thin-plate spline method using independent ground measurements at 10 sites from the Coordinated Enhanced Observation Network (CEON). These validation results indicated that the thin plate smoothing spline method can be effectively used for calibrating satellite derived ISR products using ground measurements to achieve better accuracy.
Model Adaptation for Prognostics in a Particle Filtering Framework
NASA Technical Reports Server (NTRS)
Saha, Bhaskar; Goebel, Kai Frank
2011-01-01
One of the key motivating factors for using particle filters for prognostics is the ability to include model parameters as part of the state vector to be estimated. This performs model adaptation in conjunction with state tracking, and thus, produces a tuned model that can used for long term predictions. This feature of particle filters works in most part due to the fact that they are not subject to the "curse of dimensionality", i.e. the exponential growth of computational complexity with state dimension. However, in practice, this property holds for "well-designed" particle filters only as dimensionality increases. This paper explores the notion of wellness of design in the context of predicting remaining useful life for individual discharge cycles of Li-ion batteries. Prognostic metrics are used to analyze the tradeoff between different model designs and prediction performance. Results demonstrate how sensitivity analysis may be used to arrive at a well-designed prognostic model that can take advantage of the model adaptation properties of a particle filter.
NASA Astrophysics Data System (ADS)
Rice, W. K. M.; Forgan, D. H.; Armitage, P. J.
2012-02-01
Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional smoothed particle hydrodynamics (SPH) code by Meru & Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time-scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time-scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.
Flight test results from a supercritical mission adaptive wing with smooth variable camber
NASA Technical Reports Server (NTRS)
Powers, Sheryll G.; Webb, Lannie D.; Friend, Edward L.; Lokos, William A.
1992-01-01
Results from the wing surface and boundary layer pressures, buffet studies and flight deflection measurement system for the advanced fighter technology integration F-111 mission adaptive wing program are presented. The different aerodynamic technologies studied on the aircraft, and their relationship with each other are described. The wingtip twist measurements provide an insight as to how dynamic pressures for positive normal accelerations affect the wingtip pressure profiles.
Wong Unhong; Wong Honcheng; Tang Zesheng
2010-05-21
The smoothed particle hydrodynamics (SPH), which is a class of meshfree particle methods (MPMs), has a wide range of applications from micro-scale to macro-scale as well as from discrete systems to continuum systems. Graphics hardware, originally designed for computer graphics, now provide unprecedented computational power for scientific computation. Particle system needs a huge amount of computations in physical simulation. In this paper, an efficient parallel implementation of a SPH method on graphics hardware using the Compute Unified Device Architecture is developed for fluid simulation. Comparing to the corresponding CPU implementation, our experimental results show that the new approach allows significant speedups of fluid simulation through handling huge amount of computations in parallel on graphics hardware.
Support Vector Machine Based on Adaptive Acceleration Particle Swarm Optimization
Abdulameer, Mohammed Hasan; Othman, Zulaiha Ali
2014-01-01
Existing face recognition methods utilize particle swarm optimizer (PSO) and opposition based particle swarm optimizer (OPSO) to optimize the parameters of SVM. However, the utilization of random values in the velocity calculation decreases the performance of these techniques; that is, during the velocity computation, we normally use random values for the acceleration coefficients and this creates randomness in the solution. To address this problem, an adaptive acceleration particle swarm optimization (AAPSO) technique is proposed. To evaluate our proposed method, we employ both face and iris recognition based on AAPSO with SVM (AAPSO-SVM). In the face and iris recognition systems, performance is evaluated using two human face databases, YALE and CASIA, and the UBiris dataset. In this method, we initially perform feature extraction and then recognition on the extracted features. In the recognition process, the extracted features are used for SVM training and testing. During the training and testing, the SVM parameters are optimized with the AAPSO technique, and in AAPSO, the acceleration coefficients are computed using the particle fitness values. The parameters in SVM, which are optimized by AAPSO, perform efficiently for both face and iris recognition. A comparative analysis between our proposed AAPSO-SVM and the PSO-SVM technique is presented. PMID:24790584
Support vector machine based on adaptive acceleration particle swarm optimization.
Abdulameer, Mohammed Hasan; Sheikh Abdullah, Siti Norul Huda; Othman, Zulaiha Ali
2014-01-01
Existing face recognition methods utilize particle swarm optimizer (PSO) and opposition based particle swarm optimizer (OPSO) to optimize the parameters of SVM. However, the utilization of random values in the velocity calculation decreases the performance of these techniques; that is, during the velocity computation, we normally use random values for the acceleration coefficients and this creates randomness in the solution. To address this problem, an adaptive acceleration particle swarm optimization (AAPSO) technique is proposed. To evaluate our proposed method, we employ both face and iris recognition based on AAPSO with SVM (AAPSO-SVM). In the face and iris recognition systems, performance is evaluated using two human face databases, YALE and CASIA, and the UBiris dataset. In this method, we initially perform feature extraction and then recognition on the extracted features. In the recognition process, the extracted features are used for SVM training and testing. During the training and testing, the SVM parameters are optimized with the AAPSO technique, and in AAPSO, the acceleration coefficients are computed using the particle fitness values. The parameters in SVM, which are optimized by AAPSO, perform efficiently for both face and iris recognition. A comparative analysis between our proposed AAPSO-SVM and the PSO-SVM technique is presented. PMID:24790584
NASA Astrophysics Data System (ADS)
Tong, Mingming; Browne, David J.
2012-01-01
Smoothed particle hydrodynamics is employed, for the first time, to develop a numerical model for the melting and fluid flow during laser welding process. In this meshlessLagrangian method the gas-melt two phase flow, heat transfer, surface tension, and melting of solid parent material are considered. This model was used to study the evolution of temperature field and fluid flow in the case study of laser spot welding in 2D. The simulation results show a strong influence of the melting process on the flow of liquid metal and a clear influence of the Marangoni flow on the heat transfer is also found.
NASA Astrophysics Data System (ADS)
Price, Daniel J.; Laibe, Guillaume
2015-07-01
We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves solving one extra diffusion equation for the dust fraction in addition to the usual equations of hydrodynamics. This `diffusion approximation for dust' is valid when the dust stopping time is smaller than the computational timestep. We present a numerical implementation using smoothed particle hydrodynamics that is conservative, accurate and fast. It does not require any implicit timestepping and can be straightforwardly ported into existing 3D codes.
Sinnott, M D; Cleary, P W; Arkwright, J W; Dinning, P G
2012-04-01
Complex relationships exist between gut contractility and the flow of digesta. We propose here a Smoothed Particle Hydrodynamics model coupling the flow of luminal content and wall flexure to help investigate these relationships. The model indicates that a zone of muscular relaxation preceding the contraction is an important element for transport. Low pressures in this zone generate positive thrust for low viscosity content. The viscosity of luminal content controls the localization of the flow and the magnitude of the radial pressure gradient and together with contraction amplitude they control the transport rate. For high viscosity content, high lumen occlusion is required for effective propulsion. PMID:22297431
Ansell, Thomas K; McFawn, Peter K; McLaughlin, Robert A; Sampson, David D; Eastwood, Peter R; Hillman, David R; Mitchell, Howard W; Noble, Peter B
2015-03-01
In isolated airway smooth muscle (ASM) strips, an increase or decrease in ASM length away from its current optimum length causes an immediate reduction in force production followed by a gradual time-dependent recovery in force, a phenomenon termed length adaptation. In situ, length adaptation may be initiated by a change in transmural pressure (Ptm), which is a primary physiological determinant of ASM length. The present study sought to determine the effect of sustained changes in Ptm and therefore, ASM perimeter, on airway function. We measured contractile responses in whole porcine bronchial segments in vitro before and after a sustained inflation from a baseline Ptm of 5 cmH2O to 25 cmH2O, or deflation to -5 cmH2O, for ∼50 min in each case. In one group of airways, lumen narrowing and stiffening in response to electrical field stimulation (EFS) were assessed from volume and pressure signals using a servo-controlled syringe pump with pressure feedback. In a second group of airways, lumen narrowing and the perimeter of the ASM in situ were determined by anatomical optical coherence tomography. In a third group of airways, active tension was determined under isovolumic conditions. Both inflation and deflation reduced the contractile response to EFS. Sustained Ptm change resulted in a further decrease in contractile response, which returned to baseline levels upon return to the baseline Ptm. These findings reaffirm the importance of Ptm in regulating airway narrowing. However, they do not support a role for ASM length adaptation in situ under physiological levels of ASM lengthening and shortening. PMID:25729015
NASA Astrophysics Data System (ADS)
Kumar, D.; Patra, A. K.; Pitman, E. B.; Chi, H.
2013-10-01
Smoothed Particle Hydrodynamics has been successfully used for various fluid-dynamics problems, such as breaking-waves, flooding etc., since it was originally proposed. While the Lagrangian approach is naturally suitable for free-surface flows, enforcing boundary conditions and poor approximations in the presence of discontinuities in the solution are major difficulties with the method. In this paper we present an enhanced conservative Godunov SPH based on the work of Inutsuka [S. Inutsuka, Reformulation of smoothed particle hydrodynamics with Riemann solver, Journal of Computational Physics 179 (2002) 238-267] that accurately resolves discontinuities without the need to use artificial viscosity, preserves partition of unity everywhere in the domain, correctly and flexibly enforces necessary essential and frictional slip boundary conditions to approximately solve free-surface granular flows. The development is motivated by the need to improve upon depth averaged grid based models of large scale debris flows and avalanches often characterized as granular flows. Simple validation of the results is obtained by comparison to table-top experiments.
NASA Astrophysics Data System (ADS)
Wang, Li; Sun, Xiaogang; Xing, Jian
2012-12-01
An inversion technique which combines the pattern search algorithm with the Tikhonov smoothing functional for retrieval of particle size distribution (PSD) by light extinction method is proposed. In the unparameterized shape-independent model, we first transform the PSD inversion problem into an optimization problem, with the Tikhonov smoothing functional employed to model the objective function. The optimization problem is then solved by the pattern search algorithm. To ensure good convergence rate and accuracy of the whole retrieval, a competitive strategy for determining the initial point of the pattern search algorithm is also designed. The accuracy and limitations of the proposed technique are tested by the inversion results of synthetic and real standard polystyrene particles immersed in water. In addition, the issues about the objective function and computation time are further discussed. Both simulation and experimental results show that the technique can be successfully applied to retrieve the PSD with high reliability and stability in the presence of random noise. Compared with the Phillips-Twomey method and genetic algorithm, the proposed technique has certain advantages in terms of reaching a more accurate and steady optimal solution with less computational effort, thus making this technique more suitable for quick and accurate measurement of PSD.
Adaptable Particle-in-Cell Algorithms for Graphical Processing Units
NASA Astrophysics Data System (ADS)
Decyk, Viktor; Singh, Tajendra
2010-11-01
Emerging computer architectures consist of an increasing number of shared memory computing cores in a chip, often with vector (SIMD) co-processors. Future exascale high performance systems will consist of a hierarchy of such nodes, which will require different algorithms at different levels. Since no one knows exactly how the future will evolve, we have begun development of an adaptable Particle-in-Cell (PIC) code, whose parameters can match different hardware configurations. The data structures reflect three levels of parallelism, contiguous vectors and non-contiguous blocks of vectors, which can share memory, and groups of blocks which do not. Particles are kept ordered at each time step, and the size of a sorting cell is an adjustable parameter. We have implemented a simple 2D electrostatic skeleton code whose inner loop (containing 6 subroutines) runs entirely on the NVIDIA Tesla C1060. We obtained speedups of about 16-25 compared to a 2.66 GHz Intel i7 (Nehalem), depending on the plasma temperature, with an asymptotic limit of 40 for a frozen plasma. We expect speedups of about 70 for an 2D electromagnetic code and about 100 for a 3D electromagnetic code, which have higher computational intensities (more flops/memory access).
Object tracking with adaptive HOG detector and adaptive Rao-Blackwellised particle filter
NASA Astrophysics Data System (ADS)
Rosa, Stefano; Paleari, Marco; Ariano, Paolo; Bona, Basilio
2012-01-01
Scenarios for a manned mission to the Moon or Mars call for astronaut teams to be accompanied by semiautonomous robots. A prerequisite for human-robot interaction is the capability of successfully tracking humans and objects in the environment. In this paper we present a system for real-time visual object tracking in 2D images for mobile robotic systems. The proposed algorithm is able to specialize to individual objects and to adapt to substantial changes in illumination and object appearance during tracking. The algorithm is composed by two main blocks: a detector based on Histogram of Oriented Gradient (HOG) descriptors and linear Support Vector Machines (SVM), and a tracker which is implemented by an adaptive Rao-Blackwellised particle filter (RBPF). The SVM is re-trained online on new samples taken from previous predicted positions. We use the effective sample size to decide when the classifier needs to be re-trained. Position hypotheses for the tracked object are the result of a clustering procedure applied on the set of particles. The algorithm has been tested on challenging video sequences presenting strong changes in object appearance, illumination, and occlusion. Experimental tests show that the presented method is able to achieve near real-time performances with a precision of about 7 pixels on standard video sequences of dimensions 320 × 240.
Smooth particle hydrodynamics: importance of correction terms in adaptive resolution algorithms
NASA Astrophysics Data System (ADS)
Alimi, J.-M.; Serna, A.; Pastor, C.; Bernabeu, G.
2003-11-01
We describe TREEASPH, a new code to evolve self-gravitating fluids, both with and without a collisionless component. In TREEASPH, gravitational forces are computed from a hierarchical tree algorithm (TREEcode), while hydrodynamic properties are computed by using a SPH method that includes the ∇h correction terms appearing when the spatial resolution h(t,r) is not a constant. Another important feature, which considerably increases the code efficiency on sequential and vectorial computers, is that time-stepping is performed from a PEC scheme (Predict-Evaluate-Correct) modified to allow for individual timesteps. Some authors have previously noted that the ∇h correction terms are needed to avoid the introduction on simulations of a non-physical entropy. By using TREEASPH we show here that, in cosmological simulations, this non-physical entropy has a negative sign. As a consequence, when the ∇h terms are neglected, the density peaks associated to shock fronts are overestimated. This in turn results in an overestimated efficiency of star-formation processes.
Pan, Wenxiao; Daily, Michael; Baker, Nathan A.
2015-05-07
Background: The calculation of diffusion-controlled ligand binding rates is important for understanding enzyme mechanisms as well as designing enzyme inhibitors. Methods: We demonstrate the accuracy and effectiveness of a Lagrangian particle-based method, smoothed particle hydrodynamics (SPH), to study diffusion in biomolecular systems by numerically solving the time-dependent Smoluchowski equation for continuum diffusion. Unlike previous studies, a reactive Robin boundary condition (BC), rather than the absolute absorbing (Dirichlet) BC, is considered on the reactive boundaries. This new BC treatment allows for the analysis of enzymes with “imperfect” reaction rates. Results: The numerical method is first verified in simple systems and thenmore » applied to the calculation of ligand binding to a mouse acetylcholinesterase (mAChE) monomer. Rates for inhibitor binding to mAChE are calculated at various ionic strengths and compared with experiment and other numerical methods. We find that imposition of the Robin BC improves agreement between calculated and experimental reaction rates. Conclusions: Although this initial application focuses on a single monomer system, our new method provides a framework to explore broader applications of SPH in larger-scale biomolecular complexes by taking advantage of its Lagrangian particle-based nature.« less
NASA Astrophysics Data System (ADS)
Avesani, Diego; Herrera, Paulo; Chiogna, Gabriele; Bellin, Alberto; Dumbser, Michael
2015-06-01
Most numerical schemes applied to solve the advection-diffusion equation are affected by numerical diffusion. Moreover, unphysical results, such as oscillations and negative concentrations, may emerge when an anisotropic dispersion tensor is used, which induces even more severe errors in the solution of multispecies reactive transport. To cope with this long standing problem we propose a modified version of the standard Smoothed Particle Hydrodynamics (SPH) method based on a Moving-Least-Squares-Weighted-Essentially-Non-Oscillatory (MLS-WENO) reconstruction of concentrations. This scheme formulation (called MWSPH) approximates the diffusive fluxes with a Rusanov-type Riemann solver based on high order WENO scheme. We compare the standard SPH with the MWSPH for different a few test cases, considering both homogeneous and heterogeneous flow fields and different anisotropic ratios of the dispersion tensor. We show that, MWSPH is stable and accurate and that it reduces the occurrence of negative concentrations compared to standard SPH. When negative concentrations are observed, their absolute values are several orders of magnitude smaller compared to standard SPH. In addition, MWSPH limits spurious oscillations in the numerical solution more effectively than classical SPH. Convergence analysis shows that MWSPH is computationally more demanding than SPH, but with the payoff a more accurate solution, which in addition is less sensitive to particles position. The latter property simplifies the time consuming and often user dependent procedure to define the initial dislocation of the particles.
Pan, Wenxiao; Daily, Michael D.; Baker, Nathan A.
2015-12-01
We demonstrate the accuracy and effectiveness of a Lagrangian particle-based method, smoothed particle hydrodynamics (SPH), to study diffusion in biomolecular systems by numerically solving the time-dependent Smoluchowski equation for continuum diffusion. The numerical method is first verified in simple systems and then applied to the calculation of ligand binding to an acetylcholinesterase monomer. Unlike previous studies, a reactive Robin boundary condition (BC), rather than the absolute absorbing (Dirichlet) boundary condition, is considered on the reactive boundaries. This new boundary condition treatment allows for the analysis of enzymes with "imperfect" reaction rates. Rates for inhibitor binding to mAChE are calculated at various ionic strengths and compared with experiment and other numerical methods. We find that imposition of the Robin BC improves agreement between calculated and experimental reaction rates. Although this initial application focuses on a single monomer system, our new method provides a framework to explore broader applications of SPH in larger-scale biomolecular complexes by taking advantage of its Lagrangian particle-based nature.
Pan, Wenxiao; Daily, Michael; Baker, Nathan A.
2015-05-07
Background: The calculation of diffusion-controlled ligand binding rates is important for understanding enzyme mechanisms as well as designing enzyme inhibitors. Methods: We demonstrate the accuracy and effectiveness of a Lagrangian particle-based method, smoothed particle hydrodynamics (SPH), to study diffusion in biomolecular systems by numerically solving the time-dependent Smoluchowski equation for continuum diffusion. Unlike previous studies, a reactive Robin boundary condition (BC), rather than the absolute absorbing (Dirichlet) BC, is considered on the reactive boundaries. This new BC treatment allows for the analysis of enzymes with “imperfect” reaction rates. Results: The numerical method is first verified in simple systems and then applied to the calculation of ligand binding to a mouse acetylcholinesterase (mAChE) monomer. Rates for inhibitor binding to mAChE are calculated at various ionic strengths and compared with experiment and other numerical methods. We find that imposition of the Robin BC improves agreement between calculated and experimental reaction rates. Conclusions: Although this initial application focuses on a single monomer system, our new method provides a framework to explore broader applications of SPH in larger-scale biomolecular complexes by taking advantage of its Lagrangian particle-based nature.
Gomme, Peter T; Prakash, Mahesh; Hunt, Ben; Stokes, Nick; Cleary, Paul; Tatford, Owen C; Bertolini, Joseph
2006-02-01
Using SPH (smoothed particle hydrodynamics), the motion of a lobe pump under load was simulated in order to predict the level of shear stress experienced by a protein solution. By varying the gap size between the lobes and pump housing, variations in pump efficiency and shear stress were determined. The simulations indicated that pump shear was dependent on gap size, with shear stress levels (0-40 Pa) correlating with those estimated in previous albumin-pumping studies. As gap size increased, the number of fluid particles experiencing low shear (<10 Pa) increased, whereas those experiencing high shear (>20 Pa) showed a decreasing trend. The pump efficiency, however, decreased with gap size, requiring more lobe revolutions to pass a unit volume. Taken together, these observations indicate that lobe pumps operated with increased gaps between the rotors and the housing result in larger number of particles within the fluid experiencing shear stresses. Moreover, the simulations indicate that it is best to use larger lobe pumps operated at high efficiency to transfer protein-containing solutions. PMID:16246177
NASA Astrophysics Data System (ADS)
S, Siva Prasad A. V.; Basu, Sumit
2015-10-01
Shock compaction of reactive powder mixtures to synthesize new materials is one of the oldest material processing techniques and has been studied extensively by several researchers over the past few decades. The quantitative connection between the shock energy imparted and the extent of reaction that can be completed in the small time window associated with the passage of the shock wave is complicated and depends on a large variety of parameters. In particular, our understanding of the complex interplay between the thermo-elasto-viscoplastic behaviour of the granular constituents and their temperature dependent, diffusion-limited reaction mechanism may be enriched through careful numerical simulations. A robust numerical model should be able to handle extremely large deformations coupled with diffusion mediated fast reaction kinetics. In this work, a meshfree discrete particle numerical method based on smoothed particle hydrodynamics (SPH) to simulate shock-induced chemical reactions (SICR) in reactive powder mixtures is proposed. We present a numerical strategy to carry out reactions between reactant powder particles and partition the obtained products between the particles in a manner that accounts for the requirement that the total mass of the entire system remains constant as the reactions occur. Instead of solving the reaction-diffusion problem, we propose a ‘pseudo-diffusion’ model in which a distance dependent reaction rate constant is defined to carry out chemical reaction kinetics. This approach mimics the actual reaction-diffusion process at short times. Our numerical model is demonstrated for the well-studied reaction system Nb + 2Si \\rightleftharpoons NbSi 2 . The predicted mass fractions of the product obtained from the simulations are in agreement with experimental observations. Finally, the effects of impact speed, particle arrangement and mixing ratio on the predicted product mass fractions are discussed.
NASA Astrophysics Data System (ADS)
Ganzenmüller, G. C.; Sauer, M.; May, M.; Hiermaier, S.
2016-03-01
We present a stabilization scheme for elastoplastic Smooth-Particle Hydrodynamics (SPH) which overcomes two major challenges: (i) the tensile instability inherent to the updated Lagrangian approach is suppressed and (ii) the rank-deficiency instability inherent to the nodal integration approach is cured. To achieve these goals, lessons learned from the Finite-Element Method are transferred to SPH. In particular, an analogue of hourglass control is derived for SPH, which locally linearizes the deformation field to obtain stable and accurate solutions, without the need to resort to stabilization via excessive artificial viscosity. The resulting SPH scheme combines the ability of updated Lagrangian SPH to model truly large deformations with the accuracy and stability needed to faithfully perform simulations. This claim is supported by the analysis of problematic cases and the simulation of an impact scenario.
NASA Astrophysics Data System (ADS)
Ganzenmüller, G. C.; Sauer, M.; May, M.; Hiermaier, S.
2016-04-01
We present a stabilization scheme for elastoplastic Smooth-Particle Hydrodynamics (SPH) which overcomes two major challenges: (i) the tensile instability inherent to the updated Lagrangian approach is suppressed and (ii) the rank-deficiency instability inherent to the nodal integration approach is cured. To achieve these goals, lessons learned from the Finite-Element Method are transferred to SPH. In particular, an analogue of hourglass control is derived for SPH, which locally linearizes the deformation field to obtain stable and accurate solutions, without the need to resort to stabilization via excessive artificial viscosity. The resulting SPH scheme combines the ability of updated Lagrangian SPH to model truly large deformations with the accuracy and stability needed to faithfully perform simulations. This claim is supported by the analysis of problematic cases and the simulation of an impact scenario.
NASA Astrophysics Data System (ADS)
Ganzenmüller, G. C.; Sauer, M.; May, M.; Hiermaier, S.
2016-05-01
We present a stabilization scheme for elastoplastic Smooth-Particle Hydrodynamics (SPH) which overcomes two major challenges: (i) the tensile instability inherent to the updated Lagrangian approach is suppressed and (ii) the rank-deficiency instability inherent to the nodal integration approach is cured. To achieve these goals, lessons learned from the Finite-Element Method are transferred to SPH. In particular, an analogue of hourglass control is derived for SPH, which locally linearizes the deformation field to obtain stable and accurate solutions, without the need to resort to stabilization via excessive artificial viscosity. The resulting SPH scheme combines the ability of updated Lagrangian SPH to model truly large deformations with the accuracy and stability needed to faithfully perform simulations. This claim is supported by the analysis of problematic cases and the simulation of an impact scenario.
NASA Astrophysics Data System (ADS)
Shi, Chuanqi; An, Yi; Wu, Qiang; Liu, Qingquan; Cao, Zhixian
2016-06-01
We simulate the generation of a landslide-induced impulse wave with a newly-developed soil-water coupling model in the smoothed particle hydrodynamics (SPH) framework. The model includes an elasto-plastic constitutive model for soil, a Navier-Stokes equation based model for water, and a bilateral coupling model at the interface. The model is tested with simulated waves induced by a slow and a fast landslide. Good agreement is obtained between simulation results and experimental data. The generated wave and the deformation of the landslide body can both be resolved satisfactorily. All parameters in our model have their physical meaning in soil mechanics and can be obtained from conventional soil mechanics experiments directly. The influence of the dilatancy angle of soil shows that the non-associated flow rule must be selected, and the value of the dilatancy angle should not be chosen arbitrarily, if it is not determined with relative experiments.
Margraf, J
2012-06-12
material flows through a still mesh. This is not typically done in an ALE3D analysis, especially if Lagrange elements exist. Deforming Lagrange elements would certainly tangle with a Eulerian mesh eventually. The best method in this case is to have an advecting mesh positioned as some relaxed version of the pre and post Lagrange step; this gives the best opportunity of modeling a high energy event with a combination of Lagrange and ALE elements. Dyne3D is another explicit dynamic analysis code, ParaDyn being the parallel version. ParaDyn is used for predicting the transient response of three dimensional structures using Lagrangian solid mechanics. Large deformation and mesh tangling is often resolved through the use of an element deletion scheme. This is useful to accommodate component failure, but if it is done purely as a means to preserve a useful mesh it can lead to problems because it does not maintain continuity of the material bulk response. Whatever medium exists between structural components is typically not modeled in ParaDyn. Instead, a structure either has a known loading profile applied or given initial conditions. The many included contact algorithms can calculate the loading response of materials if and when they collide. A recent implementation of an SPH module in which failed or deleted material nodes are converted to independent particles is currently being utilized for a variety of spall related problems and high velocity impact scenarios. Figure 4 shows an example of a projectile, given an initial velocity, and how it fails the first plate which generates SPH particles which then interact with and damage the second plate.
NASA Astrophysics Data System (ADS)
Suwa, T.; Imamura, F.; Sugawara, D.; Ogasawara, K.; Watanabe, M.; Hirahara, T.
2014-12-01
A tsunami simulator integrating a 3-D fluid simulation technology that runs on large-scale parallel computers using smoothed-particle hydrodynamics (SPH) method has been developed together with a 2-D tsunami propagation simulation technique using a nonlinear shallow water wave model. We use the 2-D simulation to calculate tsunami propagation of scale of about 1000km from epicenter to near shore. The 3-D SPH method can be used to calculate the water surface and hydraulic force that a tsunami can exert on a building, and to simulate flooding patterns at urban area of at most km scale. With our simulator we can also see three dimensional fluid feature such as complex changes a tsunami undergoes as it interacts with coastal topography or structures. As a result it is hoped that, e.g. , effect of the structures to dissipate waves energy passing over it can be elucidated. The authors utilize the simulator in the third of five fields of the Strategic Programs for Innovative Research, "Advanced Prediction Researches for Natural Disaster Prevention and Reduction," or the theme "Improvement of the tsunami forecasting system on the HPCI computer." The results of tsunami simulation using the K computer will be reported. We are going to apply it to a real problem of the disaster prevention in future.
NASA Astrophysics Data System (ADS)
Chang, Tsang-Jung; Chang, Kao-Hua; Kao, Hong-Ming
2014-11-01
A new approach to model weakly nonhydrostatic shallow water flows in open channels is proposed by using a Lagrangian meshless method, smoothed particle hydrodynamics (SPH). The Lagrangian form of the Boussinesq equations is solved through SPH to merge the local and convective derivatives as the material derivative. In the numerical SPH procedure, the present study uses a predictor-corrector method, in which the pure space derivative terms (the hydrostatic and source terms) are explicitly solved and the mixed space and time derivatives term (the material term of B1 and B2) is computed with an implicit scheme. It is thus a convenient tool in the processes of the space discretization compared to other Eulerian approaches. Four typical benchmark problems in weakly nonhydrostatic shallow water flows, including solitary wave propagation, nonlinear interaction of two solitary waves, dambreak flow propagation, and undular bore development, are selected to employ model validation under the closed and open boundary conditions. Numerical results are compared with the analytical solutions or published laboratory and numerical results. It is found that the proposed approach is capable of resolving weakly nonhydrostatic shallow water flows. Thus, the proposed SPH approach can supplement the lack of the SPH-Boussinesq researches in the literatures, and provide an alternative to model weakly nonhydrostatic shallow water flows in open channels.
NASA Astrophysics Data System (ADS)
Nordendale, Nikolas A.; Heard, William F.; Sherburn, Jesse A.; Basu, Prodyot K.
2016-03-01
The response of structural components of high-strength cementitious (HSC) materials to projectile impact is characterized by high-rate fragmentation resulting from strong compressive shock waves coupled with reflected tensile waves. Accurate modeling of armor panels of such brittle materials under high-velocity projectile impact is a complex problem requiring meticulous experimental characterization of material properties. In a recent paper by the authors, an approach to handle such problems based on a modified Advanced Fundamental Concrete (AFC) constitutive model was developed. In the HSC panels considered in this study, an analogous approach is applied, and the predictions are verified with ballistic impact test data. Traditional Lagrangian finite element analysis (FEA) of these problems tends to introduce errors and suffers from convergence issues resulting from large deformations at free surfaces. Also, FEA cannot properly account for the issues of secondary impact of spalled fragments when multiple armor panels are used. Smoothed particle hydrodynamics (SPH) is considered to be an attractive alternative to resolve these and other issues. However, SPH-based quantitative results have been found to be less accurate than the FEA-based ones when the deformations are not sufficiently large. This paper primarily focuses on a comparison of FEA and SPH models to predict high-velocity projectile impact on single and stacked HSC panels. Results are compared to recent ballistic experiments performed as a part of this research, and conclusions are drawn based on the findings.
Pan, Wenxiao; Bao, Jie; Tartakovsky, Alexandre M.
2014-02-15
Robin boundary condition for the Navier-Stokes equations is used to model slip conditions at the fluid-solid boundaries. A novel Continuous Boundary Force (CBF) method is proposed for solving the Navier-Stokes equations subject to Robin boundary condition. In the CBF method, the Robin boundary condition at boundary is replaced by the homogeneous Neumann boundary condition at the boundary and a volumetric force term added to the momentum conservation equation. Smoothed Particle Hydrodynamics (SPH) method is used to solve the resulting Navier-Stokes equations. We present solutions for two-dimensional and three-dimensional flows in domains bounded by flat and curved boundaries subject to various forms of the Robin boundary condition. The numerical accuracy and convergence are examined through comparison of the SPH-CBF results with the solutions of finite difference or finite element method. Taken the no-slip boundary condition as a special case of slip boundary condition, we demonstrate that the SPH-CBF method describes accurately both no-slip and slip conditions.
NASA Astrophysics Data System (ADS)
Marri, S.; White, S. D. M.
2003-10-01
We investigate a new implementation of the smoothed particle hydrodynamics technique designed to improve the realism with which galaxy formation can be simulated. In situations where cooling leads to the coexistence of phases of very different density and temperature, our method substantially reduces artificial overcooling near phase boundaries, prevents the exclusion of hot gas from the vicinity of cold `clouds' and allows relative motion of the two phases at each point. We demonstrate the numerical stability of our scheme in the presence of extremely steep density and temperature gradients, as well as in strong accretion shocks and cooling flows. In addition, we present new implementations of star formation and feedback which simulate the effect of energy injection into multiphase gas more successfully than previous schemes. Our feedback recipes deposit thermal energy separately in cold dense gas and hot diffuse gas, and can explicitly re-inject cold gas into the hot phase. They make it possible to dampen star formation effectively, to reheat cold gas, and to drive outflows into the galaxy halo and beyond. We show feedback effects to be strongest in small-mass objects where much of the gas can be expelled. After idealized tests, we carry out a first low-resolution study of galaxy formation in a Λ-cold dark matter universe. Feedback results in substantial and mass-dependent reductions in the total baryonic mass gathered on to the final object as well as in significant modulation of the star formation history.
Benchmark tests and spin adaptation for the particle-particle random phase approximation
Yang, Yang; Steinmann, Stephan N.; Peng, Degao; Aggelen, Helen van; Yang, Weitao
2013-11-07
The particle-particle random phase approximation (pp-RPA) provides an approximation to the correlation energy in density functional theory via the adiabatic connection [H. van Aggelen, Y. Yang, and W. Yang, Phys. Rev. A 88, 030501 (2013)]. It has virtually no delocalization error nor static correlation error for single-bond systems. However, with its formal O(N{sup 6}) scaling, the pp-RPA is computationally expensive. In this paper, we implement a spin-separated and spin-adapted pp-RPA algorithm, which reduces the computational cost by a substantial factor. We then perform benchmark tests on the G2/97 enthalpies of formation database, DBH24 reaction barrier database, and four test sets for non-bonded interactions (HB6/04, CT7/04, DI6/04, and WI9/04). For the G2/97 database, the pp-RPA gives a significantly smaller mean absolute error (8.3 kcal/mol) than the direct particle-hole RPA (ph-RPA) (22.7 kcal/mol). Furthermore, the error in the pp-RPA is nearly constant with the number of atoms in a molecule, while the error in the ph-RPA increases. For chemical reactions involving typical organic closed-shell molecules, pp- and ph-RPA both give accurate reaction energies. Similarly, both RPAs perform well for reaction barriers and nonbonded interactions. These results suggest that the pp-RPA gives reliable energies in chemical applications. The adiabatic connection formalism based on pairing matrix fluctuation is therefore expected to lead to widely applicable and accurate density functionals.
Moment preserving adaptive particle weights using octree velocity distributions for PIC simulations
Martin, Robert Scott; Cambier, Jean-Luc
2012-11-27
The ratio of computational to physical particles is of primary concern to statistical particle based simulations such as DSMC and PIC. An adaptive computational particle weight algorithm is presented that conserves mass, momentum, and energy. This algorithm is then enhanced with an octree adaptive mesh in velocity space to mitigate artificial thermalization. The new octree merge is compared to a merge that randomly selects merge partners for a bi-Maxwellian velocity distribution. Results for crossing beams in a fixed potential well along with an electrostatic PIC version with and without MCC collisions based ionizing breakdown show the advantages of the merge algorithm to both fixed particle weights and randomly selected merge partners.