Science.gov

Sample records for positive-energy particle models

  1. DEM Particle Fracture Model

    SciTech Connect

    Zhang, Boning; Herbold, Eric B.; Homel, Michael A.; Regueiro, Richard A.

    2015-12-01

    An adaptive particle fracture model in poly-ellipsoidal Discrete Element Method is developed. The poly-ellipsoidal particle will break into several sub-poly-ellipsoids by Hoek-Brown fracture criterion based on continuum stress and the maximum tensile stress in contacts. Also Weibull theory is introduced to consider the statistics and size effects on particle strength. Finally, high strain-rate split Hopkinson pressure bar experiment of silica sand is simulated using this newly developed model. Comparisons with experiments show that our particle fracture model can capture the mechanical behavior of this experiment very well, both in stress-strain response and particle size redistribution. The effects of density and packings o the samples are also studied in numerical examples.

  2. Modeling atmospheric particle deposition

    NASA Astrophysics Data System (ADS)

    Jackson, Msafiri M.

    Experimentally determined dry deposition velocities for atmospheric particles in the size range of 5-80 μm in diameter have been shown to be greater than predictions made with the current state-of-the-art (Sehmel-Hodgson) model which is based on wind tunnel experiment, particularly at higher wind speed. In this research, a model to predict the atmospheric dry deposition velocities of particles has been developed that is similar to a model developed for particle deposition in vertical pipes. The model uses a sigmoid curve to correlate nondimensional inertial deposition velocity (Vdi+) with dimensionless particle relaxation time (/tau+) and flow Reynolds number (Re). Vdi+ obtained from data collected in the atmosphere with particle size classifier system and a flat greased plate, Re, and /tau+ for particles between 1 and 100 μm diameter were fit with a sigmoid curve using the least square procedure to obtain coefficients for the sigmoid curve. Deposition velocities data for particles between 0.06 and 4 μm diameter developed by Sehmel-Hodgson model were used to introduce a Schmidt number (Sc) term to take care of Brownian diffusion. The atmospheric plate deposition velocity model is a function of Vst (Stokes settling velocity), V* (friction velocity), /tau+, Re, and Sc. Model application to 62 atmospheric data set revealed that: generated flux predictions agreed well with atmospheric measurements, and its performance is better than Sehmel-Hodgson model. By comparing the sigmoid curve coefficients developed for vertical pipe data with the coefficients developed for atmospheric data it is concluded that, the two types of deposition are similar when the effects of Re and /tau+ are properly considered. Sensitivity analysis for the model has revealed three distinct regions based on particle size. Of the three physical parameters (/tau+, Re, Sc) in the model, not more than two controls the deposition in any of the identified regions. The plate deposition model which is

  3. Particle bed reactor modeling

    NASA Technical Reports Server (NTRS)

    Sapyta, Joe; Reid, Hank; Walton, Lew

    1993-01-01

    The topics are presented in viewgraph form and include the following: particle bed reactor (PBR) core cross section; PBR bleed cycle; fuel and moderator flow paths; PBR modeling requirements; characteristics of PBR and nuclear thermal propulsion (NTP) modeling; challenges for PBR and NTP modeling; thermal hydraulic computer codes; capabilities for PBR/reactor application; thermal/hydralic codes; limitations; physical correlations; comparison of predicted friction factor and experimental data; frit pressure drop testing; cold frit mask factor; decay heat flow rate; startup transient simulation; and philosophy of systems modeling.

  4. Gyrokinetic particle simulation model

    SciTech Connect

    Lee, W.W.

    1986-07-01

    A new type of particle simulation model based on the gyrophase-averaged Vlasov and Poisson equations is presented. The reduced system, in which particle gyrations are removed from the equations of motion while the finite Larmor radius effects are still preserved, is most suitable for studying low frequency microinstabilities in magnetized plasmas. It is feasible to simulate an elongated system (L/sub parallel/ >> L/sub perpendicular/) with a three-dimensional grid using the present model without resorting to the usual mode expansion technique, since there is essentially no restriction on the size of ..delta..x/sub parallel/ in a gyrokinetic plasma. The new approach also enables us to further separate the time and spatial scales of the simulation from those associated with global transport through the use of multiple spatial scale expansion. Thus, the model can be a very efficient tool for studying anomalous transport problems related to steady-state drift-wave turbulence in magnetic confinement devices. It can also be applied to other areas of plasma physics.

  5. MODELING DEPOSITION OF INHALED PARTICLES

    EPA Science Inventory

    Modeling Deposition of Inhaled Particles: ABSTRACT

    The mathematical modeling of the deposition and distribution of inhaled aerosols within human lungs is an invaluable tool in predicting both the health risks associated with inhaled environmental aerosols and the therapeut...

  6. Patchy Particle Model for Vitrimers

    NASA Astrophysics Data System (ADS)

    Smallenburg, Frank; Leibler, Ludwik; Sciortino, Francesco

    2013-11-01

    Vitrimers—a recently invented new class of polymers—consist of covalent networks that can rearrange their topology via a bond shuffling mechanism, preserving the total number of network links. We introduce a patchy particle model whose dynamics directly mimic the bond exchange mechanism and reproduce the observed glass-forming ability. We calculate the free energy of this model in the limit of strong (chemical) bonds between the particles, both via the Wertheim thermodynamic perturbation theory and using computer simulations. The system exhibits an entropy-driven phase separation between a network phase and a dilute cluster gas, bringing new insight into the swelling behavior of vitrimers in solvents.

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

  8. Modeling particle loss in ventilation ducts

    SciTech Connect

    Sippola, Mark R.; Nazaroff, William W.

    2003-04-01

    Empirical equations were developed and applied to predict losses of 0.01-100 {micro}m airborne particles making a single pass through 120 different ventilation duct runs typical of those found in mid-sized office buildings. For all duct runs, losses were negligible for submicron particles and nearly complete for particles larger than 50 {micro}m. The 50th percentile cut-point diameters were 15 {micro}m in supply runs and 25 {micro}m in return runs. Losses in supply duct runs were higher than in return duct runs, mostly because internal insulation was present in portions of supply duct runs, but absent from return duct runs. Single-pass equations for particle loss in duct runs were combined with models for predicting ventilation system filtration efficiency and particle deposition to indoor surfaces to evaluate the fates of particles of indoor and outdoor origin in an archetypal mechanically ventilated building. Results suggest that duct losses are a minor influence for determining indoor concentrations for most particle sizes. Losses in ducts were of a comparable magnitude to indoor surface losses for most particle sizes. For outdoor air drawn into an unfiltered ventilation system, most particles smaller than 1 {micro}m are exhausted from the building. Large particles deposit within the building, mostly in supply ducts or on indoor surfaces. When filters are present, most particles are either filtered or exhausted. The fates of particles generated indoors follow similar trends as outdoor particles drawn into the building.

  9. Discrete Element Modeling of Triboelectrically Charged Particles

    NASA Technical Reports Server (NTRS)

    Hogue, Michael D.; Calle, Carlos I.; Weitzman, Peter S.; Curry, David R.

    2008-01-01

    Tribocharging of particles is common in many processes including fine powder handling and mixing, printer toner transport and dust extraction. In a lunar environment with its high vacuum and lack of water, electrostatic forces are an important factor to consider when designing and operating equipment. Dust mitigation and management is critical to safe and predictable performance of people and equipment. The extreme nature of lunar conditions makes it difficult and costly to carry out experiments on earth which are necessary to better understand how particles gather and transfer charge between each other and with equipment surfaces. DEM (Discrete Element Modeling) provides an excellent virtual laboratory for studying tribocharging of particles as well as for design of devices for dust mitigation and for other purposes related to handling and processing of lunar regolith. Theoretical and experimental work has been performed pursuant to incorporating screened Coulombic electrostatic forces into EDEM, a commercial DEM software package. The DEM software is used to model the trajectories of large numbers of particles for industrial particulate handling and processing applications and can be coupled with other solvers and numerical models to calculate particle interaction with surrounding media and force fields. While simple Coulombic force between two particles is well understood, its operation in an ensemble of particles is more complex. When the tribocharging of particles and surfaces due to frictional contact is also considered, it is necessary to consider longer range of interaction of particles in response to electrostatic charging. The standard DEM algorithm accounts for particle mechanical properties and inertia as a function of particle shape and mass. If fluid drag is neglected, then particle dynamics are governed by contact between particles, between particles and equipment surfaces and gravity forces. Consideration of particle charge and any tribocharging and

  10. Exploring the Standard Model of Particles

    ERIC Educational Resources Information Center

    Johansson, K. E.; Watkins, P. M.

    2013-01-01

    With the recent discovery of a new particle at the CERN Large Hadron Collider (LHC) the Higgs boson could be about to be discovered. This paper provides a brief summary of the standard model of particle physics and the importance of the Higgs boson and field in that model for non-specialists. The role of Feynman diagrams in making predictions for…

  11. Modeling and simulation of bubbles and particles

    NASA Astrophysics Data System (ADS)

    Dorgan, Andrew James

    The interaction of particles, drops, and bubbles with a fluid (gas or liquid) is important in a number of engineering problems. The present works seeks to extend the understanding of these interactions through numerical simulation. To model many of these relevant flows, it is often important to consider finite Reynolds number effects on drag, lift, torque and history force. Thus, the present work develops an equation of motion for spherical particles with a no-slip surface based on theoretical analysis, experimental data and surface-resolved simulations which is appropriate for dispersed multiphase flows. The equation of motion is then extended to account for finite particle size. This extension is critical for particles which will have a size significantly larger than the grid cell size, particularly important for bubbles and low-density particles. The extension to finite particle size is accomplished through spatial-averaging (both volume-based and surface-based) of the continuous flow properties. This averaging is consistent with the Faxen limit for solid spheres at small Reynolds numbers and added mass and fluid stress forces at inviscid limits. Further work is needed for more quantitative assessment of the truncation terms in complex flows. The new equation of motion is then used to assess the relative importance of each force in the context of two low-density particles (an air bubble and a sand particle) in a boundary layer of water. This relative importance is measured by considering effects on particle concentration, visualization of particle-fluid interactions, diffusion rates, and Lagrangian statistics collected along the particle trajectory. Strong added mass and stress gradient effects are observed for the bubble but these were found to have little effect on a sand particle of equal diameter. Lift was shown to be important for both conditions provided the terminal velocity was aligned with the flow direction. The influence of lift was found to be

  12. Lagrangian particle model for multiphase flows

    SciTech Connect

    Tartakovsky, Alexandre M.; Ferris, Kim F.; Meakin, Paul

    2009-10-01

    A Lagrangian particle model for multiphase multicomponent fluid flow, based on smoothed particle hydrodynamics (SPH), was developed and used to simulate the flow of an emulsion consisting of bubbles of a non-wetting liquid surrounded by a wetting liquid. In SPH simulations, fluids are represented by sets of particles that are used as discretization points to solve the Navier-Stokes fluid dynamics equations. In the multiphase multicomponent SPH model, a modified van der Waals equation of state is used to close the system of flow equations. The combination of the momentum conservation equation with the van der Waals equation of state results in a particle equation of motion in which the total force acting on each particle consists of many-body repulsive and viscous forces, two-body (particle-particle) attractive forces, and body forces such as gravitational forces. Similarly to molecular dynamics, for a given fluid component the combination of repulsive and attractive forces causes a phase separation. The surface tension at liquid-liquid interfaces is imposed through component dependent attractive forces. The wetting behavior of the fluids is controlled by phase dependent attractive interactions between the fluid particles and stationary particles that represent the solid phase. The dynamics of fluids away from interface is governed by purely hydrodynamic forces. Comparison with analytical solutions for static conditions and relatively simple flows demonstrates the accuracy of the SPH model.

  13. Observations and Modeling of Geospace Energetic Particles

    NASA Astrophysics Data System (ADS)

    Li, Xinlin

    2016-07-01

    Comprehensive measurements of energetic particles and electric and magnetic fields from state-of-art instruments onboard Van Allen Probes, in a geo-transfer-like orbit, revealed new features of the energetic particles and the fields in the inner magnetosphere and impose new challenges to any quantitative modeling of the physical processes responsible for these observations. Concurrent measurements of energetic particles by satellites in highly inclined low Earth orbits and plasma and fields by satellites in farther distances in the magnetospheres and in the up stream solar wind are the critically needed information for quantitative modeling and for leading to eventual accurate forecast of the variations of the energetic particles in the magnetosphere. In this presentation, emphasis will be on the most recent advance in our understanding of the energetic particles in the magnetosphere and the missing links for significantly advance in our modeling and forecasting capabilities.

  14. Polarizable water model for Dissipative Particle Dynamics

    NASA Astrophysics Data System (ADS)

    Pivkin, Igor; Peter, Emanuel

    2015-11-01

    Dissipative Particle Dynamics (DPD) is an efficient particle-based method for modeling mesoscopic behavior of fluid systems. DPD forces conserve the momentum resulting in a correct description of hydrodynamic interactions. Polarizability has been introduced into some coarse-grained particle-based simulation methods; however it has not been done with DPD before. We developed a new polarizable coarse-grained water model for DPD, which employs long-range electrostatics and Drude oscillators. In this talk, we will present the model and its applications in simulations of membrane systems, where polarization effects play an essential role.

  15. Modelling of pyrolysis of large wood particles.

    PubMed

    Sadhukhan, Anup Kumar; Gupta, Parthapratim; Saha, Ranajit Kumar

    2009-06-01

    A fully transient mathematical model has been developed to describe the pyrolysis of large biomass particles. The kinetic model consists of both primary and secondary reactions. The heat transfer model includes conductive and internal convection within the particle and convective and radiative heat transfer between the external surface and the bulk. An implicit Finite Volume Method (FVM) with Tridiagonal Matrix Algorithm (TDMA) is employed to solve the energy conservation equation. Experimental investigations are carried out for wood fines and large wood cylinder and sphere in an electrically heated furnace under inert atmosphere. The model predictions for temperature and mass loss histories are in excellent agreement with experimental results. The effect of internal convection and particle shrinkage on pyrolysis behaviour is investigated and found to be significant. Finally, simulation studies are carried out to analyze the effect of bulk temperature and particle size on total pyrolysis time and the final yield of char. PMID:19231172

  16. Sugar cane. Positive energy source for alcohol

    SciTech Connect

    Polack, J.A.; Birkett, H.S.; West, M.D.

    1981-06-01

    Sugar cane stands out as a renewable resource for fuel alcohol production, thanks to its unique, highly positive energy balance. It supplies its own processing fuel, bagasse. Net liquid fuel usage is only that consumed on the farm, amounting to a maximum of 0.3 volume per volume of ethanol produced. In some locations, the net liquid fuel consumption of the farm is as low as 0.12 volume/volume produced. This small debit may be offset by generating electric power and by foreseeable processing improvements. In view of the very favorable fuel balance for sugar cane, a decision to employ it as a renewable source of ethanol depends wholly on economic and political factors, which in turn are highly location-dependent.

  17. Modeling of particle agglomeration in nanofluids

    NASA Astrophysics Data System (ADS)

    Krishna, K. Hari; Neti, S.; Oztekin, A.; Mohapatra, S.

    2015-03-01

    Agglomeration strongly influences the stability or shelf life of nanofluid. The present computational and experimental study investigates the rate of agglomeration quantitatively. Agglomeration in nanofluids is attributed to the net effect of various inter-particle interaction forces. For the nanofluid considered here, a net inter-particle force depends on the particle size, volume fraction, pH, and electrolyte concentration. A solution of the discretized and coupled population balance equations can yield particle sizes as a function of time. Nanofluid prepared here consists of alumina nanoparticles with the average particle size of 150 nm dispersed in de-ionized water. As the pH of the colloid was moved towards the isoelectric point of alumina nanofluids, the rate of increase of average particle size increased with time due to lower net positive charge on particles. The rate at which the average particle size is increased is predicted and measured for different electrolyte concentration and volume fraction. The higher rate of agglomeration is attributed to the decrease in the electrostatic double layer repulsion forces. The rate of agglomeration decreases due to increase in the size of nano-particle clusters thus approaching zero rate of agglomeration when all the clusters are nearly uniform in size. Predicted rates of agglomeration agree adequate enough with the measured values; validating the mathematical model and numerical approach is employed.

  18. Modeling of particle agglomeration in nanofluids

    SciTech Connect

    Krishna, K. Hari; Neti, S.; Oztekin, A.; Mohapatra, S.

    2015-03-07

    Agglomeration strongly influences the stability or shelf life of nanofluid. The present computational and experimental study investigates the rate of agglomeration quantitatively. Agglomeration in nanofluids is attributed to the net effect of various inter-particle interaction forces. For the nanofluid considered here, a net inter-particle force depends on the particle size, volume fraction, pH, and electrolyte concentration. A solution of the discretized and coupled population balance equations can yield particle sizes as a function of time. Nanofluid prepared here consists of alumina nanoparticles with the average particle size of 150 nm dispersed in de-ionized water. As the pH of the colloid was moved towards the isoelectric point of alumina nanofluids, the rate of increase of average particle size increased with time due to lower net positive charge on particles. The rate at which the average particle size is increased is predicted and measured for different electrolyte concentration and volume fraction. The higher rate of agglomeration is attributed to the decrease in the electrostatic double layer repulsion forces. The rate of agglomeration decreases due to increase in the size of nano-particle clusters thus approaching zero rate of agglomeration when all the clusters are nearly uniform in size. Predicted rates of agglomeration agree adequate enough with the measured values; validating the mathematical model and numerical approach is employed.

  19. a Conduction Model Describing Particle-Particle Interaction in the Case of Surface Conducting Particles

    NASA Astrophysics Data System (ADS)

    Gonon, P.; Foulc, J.-N.; Atten, P.

    We propose an analytical conduction model describing particle-particle interactions for the case of electrorheological fluids based on surface conducting particles. The system consisting of two contacting spheres immersed in a dielectric liquid is modeled by a distributed impedances network, from which we derive analytical expressions for the potential at the spheres surface, for the electric field in the liquid phase, and finally for the interaction force. The theoretical interaction force is compared with experimental results obtained on insulating spheres coated with a thin conducting polyaniline film. A good agreement is found between the theory and experiment.

  20. Modeling Deposition of Inhaled Particles

    EPA Science Inventory

    The mathematical modeling of the deposition and distribution of inhaled aerosols within human lungs is an invaluable tool in predicting both the health risks associated with inhaled environmental aerosols and the therapeutic dose delivered by inhaled pharmacological drugs. Howeve...

  1. Probabilistic Solar Energetic Particle Models

    NASA Technical Reports Server (NTRS)

    Adams, James H., Jr.; Dietrich, William F.; Xapsos, Michael A.

    2011-01-01

    To plan and design safe and reliable space missions, it is necessary to take into account the effects of the space radiation environment. This is done by setting the goal of achieving safety and reliability with some desired level of confidence. To achieve this goal, a worst-case space radiation environment at the required confidence level must be obtained. Planning and designing then proceeds, taking into account the effects of this worst-case environment. The result will be a mission that is reliable against the effects of the space radiation environment at the desired confidence level. In this paper we will describe progress toward developing a model that provides worst-case space radiation environments at user-specified confidence levels. We will present a model for worst-case event-integrated solar proton environments that provide the worst-case differential proton spectrum. This model is based on data from IMP-8 and GOES spacecraft that provide a data base extending from 1974 to the present. We will discuss extending this work to create worst-case models for peak flux and mission-integrated fluence for protons. We will also describe plans for similar models for helium and heavier ions.

  2. Polarizable protein model for Dissipative Particle Dynamics

    NASA Astrophysics Data System (ADS)

    Peter, Emanuel; Lykov, Kirill; Pivkin, Igor

    2015-11-01

    In this talk, we present a novel polarizable protein model for the Dissipative Particle Dynamics (DPD) simulation technique, a coarse-grained particle-based method widely used in modeling of fluid systems at the mesoscale. We employ long-range electrostatics and Drude oscillators in combination with a newly developed polarizable water model. The protein in our model is resembled by a polarizable backbone and a simplified representation of the sidechains. We define the model parameters using the experimental structures of 2 proteins: TrpZip2 and TrpCage. We validate the model on folding of five other proteins and demonstrate that it successfully predicts folding of these proteins into their native conformations. As a perspective of this model, we will give a short outlook on simulations of protein aggregation in the bulk and near a model membrane, a relevant process in several Amyloid diseases, e.g. Alzheimer's and Diabetes II.

  3. Model of resonance scattering of composite particles

    SciTech Connect

    Kuperin, Yu.A.; Makarov, K.A.; Pavlov, B.S.

    1987-04-01

    A model of binary reactions in a system of particles having a nontrivial internal structure is constructed by the theory of extensions of a nonrelativistic Hamiltonian with the addition of a space of internal degrees of freedom. The model is used to describe hadron-hadron scattering at low and intermediate energies.

  4. Particle-Based Multidimensional Multispecies Biofilm Model

    PubMed Central

    Picioreanu, Cristian; Kreft, Jan-Ulrich; van Loosdrecht, Mark C. M.

    2004-01-01

    In this paper we describe a spatially multidimensional (two-dimensional [2-D] and three-dimensional [3-D]) particle-based approach for modeling the dynamics of multispecies biofilms growing on multiple substrates. The model is based on diffusion-reaction mass balances for chemical species coupled with microbial growth and spreading of biomass represented by hard spherical particles. Effectively, this is a scaled-up version of a previously proposed individual-based biofilm model. Predictions of this new particle-based model were quantitatively compared with those obtained with an established one-dimensional (1-D) multispecies model for equivalent problems. A nitrifying biofilm containing aerobic ammonium and nitrite oxidizers, anaerobic ammonium oxidizers, and inert biomass was chosen as an example. The 2-D and 3-D models generally gave the same results. If only the average flux of nutrients needs to be known, 2-D and 1-D models are very similar. However, the behavior of intermediates, which are produced and consumed in different locations within the biofilm, is better described in 2-D and 3-D models because of the multidirectional concentration gradients. The predictions of 2-D or 3-D models are also different from those of 1-D models for slowly growing or minority species in the biofilm. This aspect is related to the mechanism of biomass spreading or advection implemented in the models and should receive more attention in future experimental studies. PMID:15128564

  5. Modeling pollutant transport using a meshless-lagrangian particle model

    SciTech Connect

    Carrington, D. B.; Pepper, D. W.

    2002-01-01

    A combined meshless-Lagrangian particle transport model is used to predict pollutant transport over irregular terrain. The numerical model for initializing the velocity field is based on a meshless approach utilizing multiquadrics established by Kansa. The Lagrangian particle transport technique uses a random walk procedure to depict the advection and dispersion of pollutants over any type of surface, including street and city canyons

  6. Interspecies modeling of inhaled particle deposition patterns

    SciTech Connect

    Martonen, T.B.; Zhang, Z.; Yang, Y.

    1992-01-01

    To evaluate the potential toxic effects of ambient contaminants or therapeutic effects of airborne drugs, inhalation exposure experiments can be performed with surrogate laboratory animals. Herein, an interspecies particle deposition theory is presented for physiologically based pharmacokinetic modeling. It is derived to improve animal testing protocols. The computer code describes the behavior and fate of particles in the lungs of human subjects and a selected surrogate, the laboratory rat. In the simulations CO2 is integrated with exposure chamber atmospheres, and its concentrations regulated to produce rat breathing profiles corresponding to selected levels of human physical activity. The dosimetric model is used to calculate total, compartmental (i.e., tracheobronchial and pulmonary), and localized distribution patterns of inhaled particles in rats and humans for comparable ventilatory conditions. It is demonstrated that the model can be used to predetermine the exposure conditions necessary to produce deposition patterns in rats that are equivalent to those in humans at prescribed physical activities.

  7. Model of Image Artifacts from Dust Particles

    NASA Technical Reports Server (NTRS)

    Willson, Reg

    2008-01-01

    A mathematical model of image artifacts produced by dust particles on lenses has been derived. Machine-vision systems often have to work with camera lenses that become dusty during use. Dust particles on the front surface of a lens produce image artifacts that can potentially affect the performance of a machine-vision algorithm. The present model satisfies a need for a means of synthesizing dust image artifacts for testing machine-vision algorithms for robustness (or the lack thereof) in the presence of dust on lenses. A dust particle can absorb light or scatter light out of some pixels, thereby giving rise to a dark dust artifact. It can also scatter light into other pixels, thereby giving rise to a bright dust artifact. For the sake of simplicity, this model deals only with dark dust artifacts. The model effectively represents dark dust artifacts as an attenuation image consisting of an array of diffuse darkened spots centered at image locations corresponding to the locations of dust particles. The dust artifacts are computationally incorporated into a given test image by simply multiplying the brightness value of each pixel by a transmission factor that incorporates the factor of attenuation, by dust particles, of the light incident on that pixel. With respect to computation of the attenuation and transmission factors, the model is based on a first-order geometric (ray)-optics treatment of the shadows cast by dust particles on the image detector. In this model, the light collected by a pixel is deemed to be confined to a pair of cones defined by the location of the pixel s image in object space, the entrance pupil of the lens, and the location of the pixel in the image plane (see Figure 1). For simplicity, it is assumed that the size of a dust particle is somewhat less than the diameter, at the front surface of the lens, of any collection cone containing all or part of that dust particle. Under this assumption, the shape of any individual dust particle artifact

  8. Computer Models Simulate Fine Particle Dispersion

    NASA Technical Reports Server (NTRS)

    2010-01-01

    Through a NASA Seed Fund partnership with DEM Solutions Inc., of Lebanon, New Hampshire, scientists at Kennedy Space Center refined existing software to study the electrostatic phenomena of granular and bulk materials as they apply to planetary surfaces. The software, EDEM, allows users to import particles and obtain accurate representations of their shapes for modeling purposes, such as simulating bulk solids behavior, and was enhanced to be able to more accurately model fine, abrasive, cohesive particles. These new EDEM capabilities can be applied in many industries unrelated to space exploration and have been adopted by several prominent U.S. companies, including John Deere, Pfizer, and Procter & Gamble.

  9. Particle production beyond the thermal model

    NASA Astrophysics Data System (ADS)

    Wolschin, Georg

    2016-07-01

    The sources of particle production in relativistic heavy-ion collisions are investigated from RHIC to LHC energies. Whereas charged-hadron production in the fragmentation sources follows a ln(sNN/s0) law, particle production in the mid-rapidity low-x gluon-gluon source exhibits a much stronger dependence ∝ ln3(sNN/s0), and becomes dominant between RHIC and LHC energies. The equilibration of the three sources is investigated in a relativistic diffusion model (RDM). It agrees with the thermal model only for t → ∞.

  10. Lagrangian Trajectory Modeling of Lunar Dust Particles

    NASA Technical Reports Server (NTRS)

    Lane, John E.; Metzger, Philip T.; Immer, Christopher D.

    2008-01-01

    Apollo landing videos shot from inside the right LEM window, provide a quantitative measure of the characteristics and dynamics of the ejecta spray of lunar regolith particles beneath the Lander during the final 10 [m] or so of descent. Photogrammetry analysis gives an estimate of the thickness of the dust layer and angle of trajectory. In addition, Apollo landing video analysis divulges valuable information on the regolith ejecta interactions with lunar surface topography. For example, dense dust streaks are seen to originate at the outer rims of craters within a critical radius of the Lander during descent. The primary intent of this work was to develop a mathematical model and software implementation for the trajectory simulation of lunar dust particles acted on by gas jets originating from the nozzle of a lunar Lander, where the particle sizes typically range from 10 micron to 500 micron. The high temperature, supersonic jet of gas that is exhausted from a rocket engine can propel dust, soil, gravel, as well as small rocks to high velocities. The lunar vacuum allows ejected particles to travel great distances unimpeded, and in the case of smaller particles, escape velocities may be reached. The particle size distributions and kinetic energies of ejected particles can lead to damage to the landing spacecraft or to other hardware that has previously been deployed in the vicinity. Thus the primary motivation behind this work is to seek a better understanding for the purpose of modeling and predicting the behavior of regolith dust particle trajectories during powered rocket descent and ascent.

  11. Modeling light scattering from diesel soot particles

    SciTech Connect

    Hull, Patricia; Shepherd, Ian; Hunt, Arlon

    2002-07-16

    The Mie model is widely used to analyze light scattering from particulate aerosols. The Diesel Particle Scatterometer (DPS), for example, determines the size and optical properties of diesel exhaust particles that are characterized by measuring three angle-dependent elements of the Mueller scattering matrix. These elements are then fitted using Mie calculations with a Levenburg-Marquardt optimization program. This approach has achieved good fits for most experimental data. However, in many cases, the predicted real and imaginary parts of the index of refraction were less than that for solid carbon. To understand this result and explain the experimental data, we present an assessment of the Mie model by use of a light scattering model based on the coupled dipole approximation. The results indicate that the Mie calculation can be used to determine the largest dimension of irregularly shaped particles at sizes characteristic of Diesel soot and, for particles of known refractive index, tables can be constructed to determine the average porosity of the particles from the predicted index of refraction.

  12. Bonded-cell model for particle fracture.

    PubMed

    Nguyen, Duc-Hanh; Azéma, Emilien; Sornay, Philippe; Radjai, Farhang

    2015-02-01

    Particle degradation and fracture play an important role in natural granular flows and in many applications of granular materials. We analyze the fracture properties of two-dimensional disklike particles modeled as aggregates of rigid cells bonded along their sides by a cohesive Mohr-Coulomb law and simulated by the contact dynamics method. We show that the compressive strength scales with tensile strength between cells but depends also on the friction coefficient and a parameter describing cell shape distribution. The statistical scatter of compressive strength is well described by the Weibull distribution function with a shape parameter varying from 6 to 10 depending on cell shape distribution. We show that this distribution may be understood in terms of percolating critical intercellular contacts. We propose a random-walk model of critical contacts that leads to particle size dependence of the compressive strength in good agreement with our simulation data. PMID:25768494

  13. Bonded-cell model for particle fracture

    NASA Astrophysics Data System (ADS)

    Nguyen, Duc-Hanh; Azéma, Emilien; Sornay, Philippe; Radjai, Farhang

    2015-02-01

    Particle degradation and fracture play an important role in natural granular flows and in many applications of granular materials. We analyze the fracture properties of two-dimensional disklike particles modeled as aggregates of rigid cells bonded along their sides by a cohesive Mohr-Coulomb law and simulated by the contact dynamics method. We show that the compressive strength scales with tensile strength between cells but depends also on the friction coefficient and a parameter describing cell shape distribution. The statistical scatter of compressive strength is well described by the Weibull distribution function with a shape parameter varying from 6 to 10 depending on cell shape distribution. We show that this distribution may be understood in terms of percolating critical intercellular contacts. We propose a random-walk model of critical contacts that leads to particle size dependence of the compressive strength in good agreement with our simulation data.

  14. Modeling of Particle Agglomeration in Nanofluids

    NASA Astrophysics Data System (ADS)

    Kanagala, Hari Krishna

    Nanofluids are colloidal dispersions of nano sized particles (<100nm in diameter) in dispersion mediums. They are of great interest in industrial applications as heat transfer fluids owing to their enhanced thermal conductivities. Stability of nanofluids is a major problem hindering their industrial application. Agglomeration and then sedimentation are some reasons, which drastically decrease the shelf life of these nanofluids. Current research addresses the agglomeration effect and how it can affect the shelf life of a nanofluid. The reasons for agglomeration in nanofluids are attributable to the interparticle interactions which are quantified by the various theories. By altering the governing properties like volume fraction, pH and electrolyte concentration different nanofluids with instant agglomeration, slow agglomeration and no agglomeration can be produced. A numerical model is created based on the discretized population balance equations which analyses the particle size distribution at different times. Agglomeration effects have been analyzed for alumina nanoparticles with average particle size of 150nm dispersed in de-ionized water. As the pH was moved towards the isoelectric point of alumina nanofluids, the particle size distribution became broader and moved to bigger sizes rapidly with time. Particle size distributions became broader and moved to bigger sizes more quickly with time with increase in the electrolyte concentration. The two effects together can be used to create different temporal trends in the particle size distributions. Faster agglomeration is attributed to the decrease in the electrostatic double layer repulsion forces which is due to decrease in the induced charge and the double layer thickness around the particle. Bigger particle clusters show lesser agglomeration due to reaching the equilibrium size. The procedures and processes described in this work can be used to generate more stable nanofluids.

  15. Impact modeling with Smooth Particle Hydrodynamics

    SciTech Connect

    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.

  16. Validation of a Lagrangian particle model

    NASA Astrophysics Data System (ADS)

    Brzozowska, Lucyna

    2013-05-01

    In this paper a custom-developed model of dispersion of pollutants is presented. The proposed approach is based on both a Lagrangian particle model and an urban-scale diagnostic model of the air velocity field. Both models constitute a part of an operational air quality assessment system. The proposed model is validated by comparing its computed results with the results of measurements obtained in a wind tunnel reflecting conditions of the Mock Urban Setting Test (MUST) experiment. Commonly used measures of errors and model concordance are employed and the results obtained are additionally compared with those obtained by other authors for CFD and non-CFD class models. The obtained results indicate that the validity of the model presented in this paper is acceptable.

  17. Beyond the standard model of particle physics.

    PubMed

    Virdee, T S

    2016-08-28

    The Large Hadron Collider (LHC) at CERN and its experiments were conceived to tackle open questions in particle physics. The mechanism of the generation of mass of fundamental particles has been elucidated with the discovery of the Higgs boson. It is clear that the standard model is not the final theory. The open questions still awaiting clues or answers, from the LHC and other experiments, include: What is the composition of dark matter and of dark energy? Why is there more matter than anti-matter? Are there more space dimensions than the familiar three? What is the path to the unification of all the fundamental forces? This talk will discuss the status of, and prospects for, the search for new particles, symmetries and forces in order to address the open questions.This article is part of the themed issue 'Unifying physics and technology in light of Maxwell's equations'. PMID:27458261

  18. Particle Image Velocimetry in Lung Bifurcation Models

    NASA Astrophysics Data System (ADS)

    Theunissen, Raf; Riethmuller, Michel L.

    To better understand the human pulmonary system and governing aerosol deposition mechanisms within the lung, an accurate description of the airflow in the conductive and respiratory pulmonary airways is considered to be essential. In-vivo measurements are deemed impossible due to the small scales in the lung structure. Though numerical simulations can improve the insight, validation of the results is mostly lacking, especially in the extraction of particle trajectories. Numerical and experimental studies have been performed at the von Karman Institute for Fluid Dynamics on single and multiple bifurcation models representing simplifications of the lung system. Both steady and oscillating flows have been studied for the upper lung airways, while steady flow conditions were imposed in the modeling of the alveolar zones. Further experiments consisted in the extraction of particle trajectories in the models of the respiratory airways. This chapter presents an overview of the conducted measurement campaigns complemented with the main observations and results.

  19. Probabilistic Models for Solar Particle Events

    NASA Technical Reports Server (NTRS)

    Adams, James H., Jr.; Dietrich, W. F.; Xapsos, M. A.; Welton, A. M.

    2009-01-01

    Probabilistic Models of Solar Particle Events (SPEs) are used in space mission design studies to provide a description of the worst-case radiation environment that the mission must be designed to tolerate.The models determine the worst-case environment using a description of the mission and a user-specified confidence level that the provided environment will not be exceeded. This poster will focus on completing the existing suite of models by developing models for peak flux and event-integrated fluence elemental spectra for the Z>2 elements. It will also discuss methods to take into account uncertainties in the data base and the uncertainties resulting from the limited number of solar particle events in the database. These new probabilistic models are based on an extensive survey of SPE measurements of peak and event-integrated elemental differential energy spectra. Attempts are made to fit the measured spectra with eight different published models. The model giving the best fit to each spectrum is chosen and used to represent that spectrum for any energy in the energy range covered by the measurements. The set of all such spectral representations for each element is then used to determine the worst case spectrum as a function of confidence level. The spectral representation that best fits these worst case spectra is found and its dependence on confidence level is parameterized. This procedure creates probabilistic models for the peak and event-integrated spectra.

  20. Modelling the motion of particles around choanoflagellates

    NASA Astrophysics Data System (ADS)

    Orme, Belinda; Pettitt, Michala; Otto, Steve; Blake, John

    2001-11-01

    The three-dimensional particle paths due to a helical beat pattern of the flagellum of a sessile choanoflagellate, Salpingoeca Amphoridium (SA), are modelled and compared against experimental observations. The organism's main components are a flagellum and cell body which are situated above a substrate such that the interaction between these entities is crucial in determining the fluid flow around the choanoflagellate. The flow of fluid in the organism's environment can be characterised as Stokes flow and a flow field analogous to one created by the flagellum is generated by a distribution of stokeslets and dipoles along a helical curve. The model describing the flow considers interactions between a slender flagellum, an infinite flat plane (modelling the substrate) and a sphere (modelling the cell body). The use of image systems appropriate to Green's functions for a sphere and plane boundary are described. The computations depict particle paths representing passive tracers from experiments and their motion illustrates overall flow patterns. Figures are presented comparing recorded experimental data with numerically generated results for a number of particle paths. The principal results show good agreement between the experiments and theory.

  1. Online maintaining appearance model using particle filter

    NASA Astrophysics Data System (ADS)

    Chen, Siying; Lan, Tian; Wang, Jianyu; Ni, Guoqiang

    2008-03-01

    Tracking by foreground matching heavily depends on the appearance model to establish object correspondences among frames and essentially, the appearance model should encode both the difference part between the object and background to guarantee the robustness and the stable part to ensure tracking consistency. This paper provides a solution for online maintaining appearance models by adjusting features in the model. Object appearance is co-modeled by a subset of Haar features selected from the over-complete feature dictionary which encodes the discriminative part of object appearance and the color histogram which describes the stable appearance. During the particle filtering process, feature values both from background patches and object observations are sampled efficiently by the aid of "foreground" and "background" particles respectively. Based on these sampled values, top-ranked discriminative features are added and invalid features are removed out to ensure the object being distinguishable from current background according to the evolving appearance model. The tracker based on this online appearance model maintaining technique has been tested on people and car tracking tasks and promising experimental results are obtained.

  2. Modeling Deep Burn TRISO particle nuclear fuel

    NASA Astrophysics Data System (ADS)

    Besmann, T. M.; Stoller, R. E.; Samolyuk, G.; Schuck, P. C.; Golubov, S. I.; Rudin, S. P.; Wills, J. M.; Coe, J. D.; Wirth, B. D.; Kim, S.; Morgan, D. D.; Szlufarska, I.

    2012-11-01

    Under the DOE Deep Burn program TRISO fuel is being investigated as a fuel form for consuming plutonium and minor actinides, and for greater efficiency in uranium utilization. The result will thus be to drive TRISO particulate fuel to very high burn-ups. In the current effort the various phenomena in the TRISO particle are being modeled using a variety of techniques. The chemical behavior is being treated utilizing thermochemical analysis to identify phase formation/transformation and chemical activities in the particle, including kernel migration. Density functional theory is being used to understand fission product diffusion within the plutonia oxide kernel, the fission product's attack on the SiC coating layer, as well as fission product diffusion through an alternative coating layer, ZrC. Finally, a multiscale approach is being used to understand thermal transport, including the effect of radiation damage induced defects, in a model SiC material.

  3. Potentiometric detection of model bioaerosol particles.

    PubMed

    Sarantaridis, Dimitris; Caruana, Daren J

    2010-09-15

    A new technique for the detection of bioaerosols is presented, utilizing particle combustion/ionization in a premixed hydrogen/oxygen/nitrogen flame plasma, followed by gas phase electrochemical detection. Bermuda grass pollen (Cynodon dactylon, one of the most common causes of pollen allergy) and black walnut pollen (Juglans nigra) were used as model bioaerosol particles. We demonstrate that single particle detection can be comfortably achieved by zero current potential measurements between two platinum electrodes, giving potential signals of over 800 mV and unique fragmentation features which may be used for differentiating between species. The high sensitivity is due to the inherent amplification through flame fragmentation, gasification and ionization; a single pollen grain of 25 μm diameter can give a plume of combustion products measuring 4 mm in diameter. The physical basis of the potential difference is a mixed interfacial potential with an additive diffusion/junction potential due to the increase in ionization from the pollen combustion. The results suggest this methodology may be applied to the detection of particulates composed of ionizable species (organic or inorganic) in gaseous environments, such as bacteria, viruses, pollen grains, and dust. Its effectiveness will depend on the propensity of the target particle to combust and generate voltages under specific flame and electrode conditions. PMID:20738107

  4. Modeling positronium beyond the single particle approximation

    NASA Astrophysics Data System (ADS)

    Zubiaga, A.; Ervasti, M. M.; Makkonen, I.; Harju, A.; Tuomisto, F.; Puska, M. J.

    2016-03-01

    Understanding the properties of the positronium atom in matter is of interest for the interpretation of positron annihilation experiments. This technique has a unique capability for the investigation of nanometer sized voids and pores in soft molecular materials (polymers, liquids or biostructures) and porous materials. However, detailed interpretations of the experimental data rely on modeling of the annihilation properties of positronium in the host material. New applications of the technique are being developed but the computational models still are based on single particle approaches and there is no way to address the influence of the electronic properties of the host material. In this work we discuss new directions of research.

  5. A Classical Model for Virtual Particle Exchange

    NASA Astrophysics Data System (ADS)

    Bullock, Daniel; Bixler, David

    2008-10-01

    As raindrops splash on the surface of a parking lot, a bubble of air may form briefly in a puddle. Some bubbles are very small and others can be rather large. They also vary in how long they last before releasing the air trapped inside in a manner reminiscent of high-energy particle collisions. When a bubble is formed, it essentially draws energy from the ``vacuum'' or surrounding medium, and the energy must be deposited back into the medium within a predictable time. The lifetime of a bubble may follow an uncertainty principle that determines the size or energy of the bubble. This research project attempts to formulate the uncertainty in the energy and lifetime of these bubbles and model the four fundamental forces based on the range of interaction. Just as spandex has been used and disproved as a model for the ``fabric of space-time,'' this bubble model may provide some insight into how elementary particles make up the fundamental forces of nature. This model will also be tested against a change in interaction medium to better correlate the data with known uncertainties.

  6. Abnormal behaviors detection using particle motion model

    NASA Astrophysics Data System (ADS)

    Chen, Yutao; Zhang, Hong; Cheng, Feiyang; Yuan, Ding; You, Yuhu

    2015-03-01

    Human abnormal behaviors detection is one of the most challenging tasks in the video surveillance for the public security control. Interaction Energy Potential model is an effective and competitive method published recently to detect abnormal behaviors, but their model of abnormal behaviors is not accurate enough, so it has some limitations. In order to solve this problem, we propose a novel Particle Motion model. Firstly, we extract the foreground to improve the accuracy of interest points detection since the complex background usually degrade the effectiveness of interest points detection largely. Secondly, we detect the interest points using the graphics features. Here, the movement of each human target can be represented by the movements of detected interest points of the target. Then, we track these interest points in videos to record their positions and velocities. In this way, the velocity angles, position angles and distance between each two points can be calculated. Finally, we proposed a Particle Motion model to calculate the eigenvalue of each frame. An adaptive threshold method is proposed to detect abnormal behaviors. Experimental results on the BEHAVE dataset and online videos show that our method could detect fight and robbery events effectively and has a promising performance.

  7. Particle-based model for skiing traffic

    NASA Astrophysics Data System (ADS)

    Holleczek, Thomas; Tröster, Gerhard

    2012-05-01

    We develop and investigate a particle-based model for ski slope traffic. Skiers are modeled as particles with a mass that are exposed to social and physical forces, which define the riding behavior of skiers during their descents on ski slopes. We also report position and speed data of 21 skiers recorded with GPS-equipped cell phones on two ski slopes. A comparison of these data with the trajectories resulting from computer simulations of our model shows a good correspondence. A study of the relationship among the density, speed, and flow of skiers reveals that congestion does not occur even with arrival rates of skiers exceeding the maximum ski lift capacity. In a sensitivity analysis, we identify the kinetic friction coefficient of skis on snow, the skier mass, the range of repelling social forces, and the arrival rate of skiers as the crucial parameters influencing the simulation results. Our model allows for the prediction of speed zones and skier densities on ski slopes, which is important in the prevention of skiing accidents.

  8. Langangian Particle Model of Friction Stir Welding

    SciTech Connect

    Tartakovsky, Alexandre M.

    2006-12-13

    Since its invention fifteen years ago, Friction Stir Welding (FSW) has found commercial application in the marine, aerospace, rail, and now automotive industries. Development of the FSW process for each new application, however, has remained largely empirical. Few detailed numerical modeling techniques have been developed that can explain and predict important features of the process physics. This is particularly true in the areas of material flow, mixing mechanisms, and void prediction. In this paper we present a novel modeling approach to simulate FSW processes that may have significant advantages over current traditional finite element or finite difference based methods. The proposed model is based on the Smoothed Particle Hydrodynamics (SPH) method. Unlike traditional grid-based methods, Lagrangian particle methods such as SPH can simulate the dynamics of interfaces, large material deformations, void formations and the material's strain and temperature history without employing complex tracking schemes. Two- and three-dimensional FSW simulations for different tool designs are presented. Preliminary numerical results are in qualitative agreement with experimental observations. Detailed comparisons between experimental measurements and larger scale FSW simulations are required to further validate and calibrate the SPH based FSW model.

  9. Particle-Surface Interaction Model and Method of Determining Particle-Surface Interactions

    NASA Technical Reports Server (NTRS)

    Hughes, David W. (Inventor)

    2012-01-01

    A method and model of predicting particle-surface interactions with a surface, such as the surface of a spacecraft. The method includes the steps of: determining a trajectory path of a plurality of moving particles; predicting whether any of the moving particles will intersect a surface; predicting whether any of the particles will be captured by the surface and/or; predicting a reflected trajectory and velocity of particles reflected from the surface.

  10. Trajectory Model of Lunar Dust Particles

    NASA Technical Reports Server (NTRS)

    2008-01-01

    The goal of this work was to predict the trajectories of blowing lunar regolith (soil) particles when a spacecraft lands on or launches from the Moon. The blown regolith is known to travel at very high velocity and to damage any hardware located nearby on the Moon. It is important to understand the trajectories so we can develop technologies to mitigate the blast effects for the launch and landing zones at a lunar outpost. A mathematical model was implemented in software to predict the trajectory of a single spherical mass acted on by the gas jet from the nozzle of a lunar lander.

  11. Model-independent particle accelerator tuning

    SciTech Connect

    Scheinker, Alexander; Pang, Xiaoying; Rybarcyk, Larry

    2013-10-21

    We present a new model-independent dynamic feedback technique, rotation rate tuning, for automatically and simultaneously tuning coupled components of uncertain, complex systems. The main advantages of the method are: 1) It has the ability to handle unknown, time-varying systems, 2) It gives known bounds on parameter update rates, 3) We give an analytic proof of its convergence and its stability, and 4) It has a simple digital implementation through a control system such as the Experimental Physics and Industrial Control System (EPICS). Because this technique is model independent it may be useful as a real-time, in-hardware, feedback-based optimization scheme for uncertain and time-varying systems. In particular, it is robust enough to handle uncertainty due to coupling, thermal cycling, misalignments, and manufacturing imperfections. As a result, it may be used as a fine-tuning supplement for existing accelerator tuning/control schemes. We present multi-particle simulation results demonstrating the scheme’s ability to simultaneously adaptively adjust the set points of twenty two quadrupole magnets and two RF buncher cavities in the Los Alamos Neutron Science Center Linear Accelerator’s transport region, while the beam properties and RF phase shift are continuously varying. The tuning is based only on beam current readings, without knowledge of particle dynamics. We also present an outline of how to implement this general scheme in software for optimization, and in hardware for feedback-based control/tuning, for a wide range of systems.

  12. Geometrical Scaling of Hall Thruster Particle Model

    SciTech Connect

    Taccogna, Francesco; Longo, Savino; Capitelli, Mario; Schneider, Ralf

    2005-05-16

    Non-Maxwellian behaviour and plasma-wall interaction are key processes in the physics of Hall thrusters. For this purpose, a 2D{l_brace}r,z{r_brace}-3V axisymmetric fully kinetic Particle-in-Cell/Monte Carlo Collision (PIC-MCC) model of the acceleration channel including the process of secondary electron emission (SEE) from the dielectric walls has been developed. In order to make the simulation possible with regard to the computational time, a reduction of the thruster dimension was done. This was derived from a new physics-based scaling law. This model has demonstrated its outstanding capability in improving the physics insight into the processes in Stationary Plasma Thruster (SPT) and in reproducing accurately well experimental data.

  13. Toward a descriptive model of solar particles in the heliosphere

    NASA Technical Reports Server (NTRS)

    Shea, M. A.; Smart, D. F.; Adams, James H., Jr.; Chenette, D.; Feynman, Joan; Hamilton, Douglas C.; Heckman, G. R.; Konradi, A.; Lee, Martin A.; Nachtwey, D. S.

    1988-01-01

    During a workshop on the interplanetary charged particle environment held in 1987, a descriptive model of solar particles in the heliosphere was assembled. This model includes the fluence, composition, energy spectra, and spatial and temporal variations of solar particles both within and beyong 1 AU. The ability to predict solar particle fluences was also discussed. Suggestions for specific studies designed to improve the basic model were also made.

  14. Multiscale modeling with smoothed dissipative particle dynamics.

    PubMed

    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

  15. Microscale Simulations of Shock Interaction with Large Assembly of Particles for Developing Point-Particle Models

    NASA Astrophysics Data System (ADS)

    Thakur, Siddharth; Neal, Chris; Mehta, Yash; Sridharan, Prashanth; Jackson, Tom; Balachandar, S.; University of Florida Team

    2015-06-01

    Micrsoscale simulations are being conducted for developing point-particle models that are needed for macroscale simulations of explosive dispersal of particles. These particle models are required to compute instantaneous force and heat transfer between particles and surroundings. A strategy for a sequence of microscale simulations has been devised for systematic development of hybrid surrogate models that are applicable at conditions representative of explosive dispersal. The microscale simulations examine particle force dependence on: Mach number, Reynolds number, and volume fraction (particle arrangements such as cubic, face-centered cubic, body-centered cubic and random). Future plans include investigation of sequences of fully-resolved microscale simulations consisting of an array of particles subjected to more realistic time-dependent flows that progressively better approximate the problem of explosive dispersal. Additionally, effects of particle shape, size, and number as well as the transient particle deformation dependence on parameters including: (a) particle material, (b) medium material, (c) multiple particles, (d) incoming shock pressure and speed, (e) medium to particle impedance ratio, (f) particle shape and orientation to shock, etc. are being investigated.

  16. Model-independent particle accelerator tuning

    DOE PAGESBeta

    Scheinker, Alexander; Pang, Xiaoying; Rybarcyk, Larry

    2013-10-21

    We present a new model-independent dynamic feedback technique, rotation rate tuning, for automatically and simultaneously tuning coupled components of uncertain, complex systems. The main advantages of the method are: 1) It has the ability to handle unknown, time-varying systems, 2) It gives known bounds on parameter update rates, 3) We give an analytic proof of its convergence and its stability, and 4) It has a simple digital implementation through a control system such as the Experimental Physics and Industrial Control System (EPICS). Because this technique is model independent it may be useful as a real-time, in-hardware, feedback-based optimization scheme formore » uncertain and time-varying systems. In particular, it is robust enough to handle uncertainty due to coupling, thermal cycling, misalignments, and manufacturing imperfections. As a result, it may be used as a fine-tuning supplement for existing accelerator tuning/control schemes. We present multi-particle simulation results demonstrating the scheme’s ability to simultaneously adaptively adjust the set points of twenty two quadrupole magnets and two RF buncher cavities in the Los Alamos Neutron Science Center Linear Accelerator’s transport region, while the beam properties and RF phase shift are continuously varying. The tuning is based only on beam current readings, without knowledge of particle dynamics. We also present an outline of how to implement this general scheme in software for optimization, and in hardware for feedback-based control/tuning, for a wide range of systems.« less

  17. Multiscale modelling of nucleosome core particle aggregation

    NASA Astrophysics Data System (ADS)

    Lyubartsev, Alexander P.; Korolev, Nikolay; Fan, Yanping; Nordenskiöld, Lars

    2015-02-01

    The nucleosome core particle (NCP) is the basic building block of chromatin. Under the influence of multivalent cations, isolated mononucleosomes exhibit a rich phase behaviour forming various columnar phases with characteristic NCP-NCP stacking. NCP stacking is also a regular element of chromatin structure in vivo. Understanding the mechanism of nucleosome stacking and the conditions leading to self-assembly of NCPs is still incomplete. Due to the complexity of the system and the need to describe electrostatics properly by including the explicit mobile ions, novel modelling approaches based on coarse-grained (CG) methods at the multiscale level becomes a necessity. In this work we present a multiscale CG computer simulation approach to modelling interactions and self-assembly of solutions of NCPs induced by the presence of multivalent cations. Starting from continuum simulations including explicit three-valent cobalt(III)hexammine (CoHex3+) counterions and 20 NCPs, based on a previously developed advanced CG NCP model with one bead per amino acid and five beads per two DNA base pair unit (Fan et al 2013 PLoS One 8 e54228), we use the inverse Monte Carlo method to calculate effective interaction potentials for a ‘super-CG’ NCP model consisting of seven beads for each NCP. These interaction potentials are used in large-scale simulations of up to 5000 NCPs, modelling self-assembly induced by CoHex3+. The systems of ‘super-CG’ NCPs form a single large cluster of stacked NCPs without long-range order in agreement with experimental data for NCPs precipitated by the three-valent polyamine, spermidine3+.

  18. Shock interaction with a deformable particle: Direct numerical simulation and point-particle modeling

    NASA Astrophysics Data System (ADS)

    Ling, Y.; Haselbacher, A.; Balachandar, S.; Najjar, F. M.; Stewart, D. S.

    2013-01-01

    The interaction of shock waves with deformable particles is an important fundamental problem. In some applications, e.g., the detonation of explosives loaded with metal particles, the pressure behind the shock wave can be significantly larger than the yield strength of the particle material. This means that particles can deform severely during their interaction with the shock wave. The experimental and theoretical studies of shock interaction with deformable particles (SIDP) are extremely challenging because of its highly transient nature. As a result, no accurate model exists yet that can be used in simulations. The objective of this paper is to develop a simple point-particle model that accurately captures the unsteady force and heat-transfer in SIDP. In the development of this model, we build on earlier models by Ling et al. (Int. J. Multiphase Flow 37, 1026-1044 (2011)) for the unsteady force and heat-transfer contributions for rigid particles. Insights gained from direct numerical simulations (DNS) guide the extension of these models to deforming particles. Results obtained with the extended model for the interaction of a deforming particle with a shock wave and a Chapman-Jouguet detonation wave compare well with DNS results and therefore offer significant improvements over standard models.

  19. Modeling Particle Acceleration and Transport at CIRs

    NASA Astrophysics Data System (ADS)

    Li, G.; Zhao, L.; Ebert, R. W.; Desai, M. I.; Dayeh, M. A.; Mason, G. M.; Chen, Y.; Wu, Z.

    2014-12-01

    CIRs are a major site for particle acceleration during solar minimum. Earlier Ulysses observations have found that particles can be accelerated at both the forward and the reverse shocks that often form at a few AUs. The accelerated particles then propagate back to the Earth along Parker's field line. Theoretical calculations predicted a modulation of the spectrum at low energies, qualitatively agreed with obsevations at 1 AU. However, this picture was recently challenged by STEREO observations, where local accelerations near 1 AU were inferred in many events. In this work, we perform a detailed numerical calculation to study particle acceleration and transport in one CIR event which was observed by both ACE and STEREO spacecraft. We obtain particle currents at different heliocentric distances and different longitudes, as well as particle anisotropy. These values are compared with observations and the implication on the acceleration site and the interplanetary turbulence spectrum is discussed.

  20. Modeling Particle Exposure in US Trucking Terminals

    PubMed Central

    Davis, ME; Smith, TJ; Laden, F; Hart, JE; Ryan, LM; Garshick, E

    2007-01-01

    Multi-tiered sampling approaches are common in environmental and occupational exposure assessment, where exposures for a given individual are often modeled based on simultaneous measurements taken at multiple indoor and outdoor sites. The monitoring data from such studies is hierarchical by design, imposing a complex covariance structure that must be accounted for in order to obtain unbiased estimates of exposure. Statistical methods such as structural equation modeling (SEM) represent a useful alternative to simple linear regression in these cases, providing simultaneous and unbiased predictions of each level of exposure based on a set of covariates specific to the exposure setting. We test the SEM approach using data from a large exposure assessment of diesel and combustion particles in the US trucking industry. The exposure assessment includes data from 36 different trucking terminals across the United States sampled between 2001 and 2005, measuring PM2.5 and its elemental carbon (EC), organic carbon (OC) components, by personal monitoring, and sampling at two indoor work locations and an outdoor “background” location. Using the SEM method, we predict: 1) personal exposures as a function of work related exposure and smoking status; 2) work related exposure as a function of terminal characteristics, indoor ventilation, job location, and background exposure conditions; and 3) background exposure conditions as a function of weather, nearby source pollution, and other regional differences across terminal sites. The primary advantage of SEMs in this setting is the ability to simultaneously predict exposures at each of the sampling locations, while accounting for the complex covariance structure among the measurements and descriptive variables. The statistically significant results and high R2 values observed from the trucking industry application supports the broader use of this approach in exposure assessment modeling. PMID:16856739

  1. Multiscale Modeling of Metallic Materials Containing Embedded Particles

    NASA Technical Reports Server (NTRS)

    Phillips, Dawn R.; Iesulauro, Erin; Glaessgen, Edward H.

    2004-01-01

    Multiscale modeling at small length scales (10(exp -9) to 10(exp -3) m) is discussed for aluminum matrices with embedded particles. A configuration containing one particle surrounded by about 50 grains and subjected to uniform tension and lateral constraint is considered. The analyses are performed to better understand the effects of material configuration on the initiation and progression of debonding of the particles from the surrounding aluminum matrix. Configurational parameters considered include particle aspect ratio and orientation within the surrounding matrix. Both configurational parameters are shown to have a significant effect on the behavior of the materials as a whole. For elliptical particles with the major axis perpendicular to the direction of loading, a particle with a 1:1 aspect ratio completely debonds from the surrounding matrix at higher loads than particles with higher aspect ratios. As the particle major axis is aligned with the direction of the applied load, increasing amounts of load are required to completely debond the particles.

  2. Cluster kinetics model of particle separation in vibrated granular media

    NASA Astrophysics Data System (ADS)

    McCoy, Benjamin J.; Madras, Giridhar

    2006-01-01

    We model the Brazil-nut effect (BNE) by hypothesizing that granules form clusters that fragment and aggregate. This provides a heterogeneous medium in which the immersed intruder particle rises (BNE) or sinks (reverse BNE) according to relative convection currents and buoyant and drag forces. A simple relationship proposed for viscous drag in terms of the vibrational intensity and the particle to grain density ratio allows simulation of published experimental data for rise and sink times as functions of particle radius, initial depth of the particle, and particle-grain density ratio. The proposed model correctly describes the experimentally observed maximum in risetime.

  3. Cluster kinetics model of particle separation in vibrated granular media.

    PubMed

    McCoy, Benjamin J; Madras, Giridhar

    2006-01-01

    We model the Brazil-nut effect (BNE) by hypothesizing that granules form clusters that fragment and aggregate. This provides a heterogeneous medium in which the immersed intruder particle rises (BNE) or sinks (reverse BNE) according to relative convection currents and buoyant and drag forces. A simple relationship proposed for viscous drag in terms of the vibrational intensity and the particle to grain density ratio allows simulation of published experimental data for rise and sink times as functions of particle radius, initial depth of the particle, and particle-grain density ratio. The proposed model correctly describes the experimentally observed maximum in risetime. PMID:16486131

  4. Continuum modeling of hydrodynamic particle-particle interactions in microfluidic high-concentration suspensions.

    PubMed

    Ley, Mikkel W H; Bruus, Henrik

    2016-03-23

    A continuum model is established for numerical studies of hydrodynamic particle-particle interactions in microfluidic high-concentration suspensions. A suspension of microparticles placed in a microfluidic channel and influenced by an external force, is described by a continuous particle-concentration field coupled to the continuity and Navier-Stokes equation for the solution. The hydrodynamic interactions are accounted for through the concentration dependence of the suspension viscosity, of the single-particle mobility, and of the momentum transfer from the particles to the suspension. The model is applied on a magnetophoretic and an acoustophoretic system, respectively, and based on the results, we illustrate three main points: (1) for relative particle-to-fluid volume fractions greater than 0.01, the hydrodynamic interaction effects become important through a decreased particle mobility and an increased suspension viscosity. (2) At these high particle concentrations, particle-induced flow rolls occur, which can lead to significant deviations of the advective particle transport relative to that of dilute suspensions. (3) Which interaction mechanism that dominates, depends on the specific flow geometry and the specific external force acting on the particles. PMID:26948344

  5. Experimental and Modeling Studies of Particles in Plasmas

    NASA Astrophysics Data System (ADS)

    Daugherty, John Edward

    Particles that are generated during plasma processing are an important source of contamination in microelectronic device fabrication. In this work we investigate the transport of plasma-generated particles in a radiofrequency (rf) diode plasma reactor, and we characterize this transport in terms of the forces that act on particles in the plasma environment. In this way we determine the mechanisms by which particles transport to critical processing surfaces. Aluminum and copper particles are formed in a radiofrequency argon sputtering system. The spatial distribution of these particles is observed with elastic laser light scattering. The distribution of particles in the discharge is found to depend most strongly on discharge power and particle size. Large particles and particles in high power discharges tend to segregate near the plasma sheath boundary. Small particles and particles in very low power discharges tend to accumulate in the center of the discharge. The spatial distributions that are observed experimentally are consistent with a model of the forces on particles in a plasma environment. The model includes the electrostatic force, the thermophoretic force, gravity, and momentum transfer from drifting plasma ions and from drifting neutral gas. In addition, the charging of small particles is investigated. We find that the particle charge can usually be predicted using simplified theory rather than the full plasma kinetic theory. Also, the electrostatic potential distribution in the vicinity of a charged particle is found to resemble a Debye-Huckel potential profile. Particles are observed to affect the discharge structure. Spatially- and temporally-resolved optical emission intensities take on characteristics that are associated with discharges in electronegative gases. We take this as evidence that the particles act as heavy, multiply charged negative ions. Plasmas with particles have a smaller self bias, and they are more resistive than discharges in pure

  6. Particle hopping vs. fluid-dynamical models for traffic flow

    SciTech Connect

    Nagel, K.

    1995-12-31

    Although particle hopping models have been introduced into traffic science in the 19509, their systematic use has only started recently. Two reasons for this are, that they are advantageous on modem computers, and that recent theoretical developments allow analytical understanding of their properties and therefore more confidence for their use. In principle, particle hopping models fit between microscopic models for driving and fluiddynamical models for traffic flow. In this sense, they also help closing the conceptual gap between these two. This paper shows connections between particle hopping models and traffic flow theory. It shows that the hydrodynamical limits of certain particle hopping models correspond to the Lighthill-Whitham theory for traffic flow, and that only slightly more complex particle hopping models produce already the correct traffic jam dynamics, consistent with recent fluid-dynamical models for traffic flow. By doing so, this paper establishes that, on the macroscopic level, particle hopping models are at least as good as fluid-dynamical models. Yet, particle hopping models have at least two advantages over fluid-dynamical models: they straightforwardly allow microscopic simulations, and they include stochasticity.

  7. Saturn's rings - Particle size distributions for thin layer model

    NASA Technical Reports Server (NTRS)

    Zebker, H. A.; Marouf, E. A.; Tyler, G. L.

    1985-01-01

    A model incorporating limited interaction between the incident energy and particles in the ring is considered which appears to be consistent with the multiple scattering process in Saturn's rings. The model allows for the small physical thickness of the rings and can be used to relate Voyager 1 observations of 3.6- and 13-cm wavelength microwave scatter from the rings to the ring particle size distribution function for particles with radii ranging from 0.001 to 20 m. This limited-scatter model yields solutions for particle size distribution functions for eight regions in the rings, which exhibit approximately inverse-cubic power-law behavior.

  8. INTERSPECIES MODELING OF INHALED PARTICLE DEPOSITION

    EPA Science Inventory

    To evaluate the potential toxic effects of ambient contaminants or therapeutic effects of airborne drugs, inhalation exposure experiments can be performed with surrogate laboratory animals. erein, an interspecies particle deposition theory is presented for physiologically based p...

  9. Experimental validation of different modeling approaches for solid particle receivers.

    SciTech Connect

    Khalsa, Siri Sahib S.; Amsbeck, Lars , Spain and Stuttgart, Germany); Roger, Marc , Spain and Stuttgart, Germany); Siegel, Nathan Phillip; Kolb, Gregory J.; Buck, Reiner , Spain and Stuttgart, Germany); Ho, Clifford Kuofei

    2009-07-01

    Solid particle receivers have the potential to provide high-temperature heat for advanced power cycles, thermochemical processes, and thermal storage via direct particle absorption of concentrated solar energy. This paper presents two different models to evaluate the performance of these systems. One model is a detailed computational fluid dynamics model using FLUENT that includes irradiation from the concentrated solar flux, two-band re-radiation and emission within the cavity, discrete-phase particle transport and heat transfer, gas-phase convection, wall conduction, and radiative and convective heat losses. The second model is an easy-to-use and fast simulation code using Matlab that includes solar and thermal radiation exchange between the particle curtain, cavity walls, and aperture, but neglects convection. Both models were compared to unheated particle flow tests and to on-sun heating tests. Comparisons between measured and simulated particle velocities, opacity, particle volume fractions, particle temperatures, and thermal efficiencies were found to be in good agreement. Sensitivity studies were also performed with the models to identify parameters and modifications to improve the performance of the solid particle receiver.

  10. White dwarfs constraints on dark sector models with light particles

    SciTech Connect

    Ubaldi, Lorenzo

    2014-06-24

    The white dwarf luminosity function is well understood in terms of standard model physics and leaves little room for exotic cooling mechanisms related to the possible existence of new weakly interacting light particles. This puts significant constraints on the parameter space of models that contain a massive dark photon and light dark sector particles.

  11. An incompressible two-dimensional multiphase particle-in-cell model for dense particle flows

    SciTech Connect

    Snider, D.M.; O`Rourke, P.J.; Andrews, M.J.

    1997-06-01

    A two-dimensional, incompressible, multiphase particle-in-cell (MP-PIC) method is presented for dense particle flows. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Difficulties associated with calculating interparticle interactions for dense particle flows with volume fractions above 5% have been eliminated by mapping particle properties to a Eulerian grid and then mapping back computed stress tensors to particle positions. This approach utilizes the best of Eulerian/Eulerian continuum models and Eulerian/Lagrangian discrete models. The solution scheme allows for distributions of types, sizes, and density of particles, with no numerical diffusion from the Lagrangian particle calculations. The computational method is implicit with respect to pressure, velocity, and volume fraction in the continuum solution thus avoiding courant limits on computational time advancement. MP-PIC simulations are compared with one-dimensional problems that have analytical solutions and with two-dimensional problems for which there are experimental data.

  12. Dirac equation for electrodynamic model particles

    NASA Astrophysics Data System (ADS)

    Zheng-Johansson, J. X.

    2008-08-01

    We set up the Maxwell's equations and subsequently the classical wave equations for the electromagnetic waves which together with their generating source, an oscillatory charge of zero rest mass in general travelling, make up a particle travelling similarly as the source at velocity ν in the field of an external scalar and vector potentials. The direct solutions in constant external field are Doppler-displaced plane waves propagating at the velocity of light c; at the de Broglie wavelength scale and expressed in terms of the dynamically equivalent and appropriate geometric mean wave variables, these render as functions identical to the space-time functions of a corresponding Dirac spinor, and in turn identical to de Broglie phase waves previously obtained from explicit superposition. For two spin-half particles of a common set of space-time functions constrained with antisymmetric spin functions as follows the Pauli principle for same charges and as separately indirectly induced based on experiment for opposite charges, the complete wave functions are identical to the Dirac spinor. The back-substitution of the so explicitly determined complete wave functions in the corresponding classical wave equations of the two particles, subjected further to reductions appropriate for the stationary-state particle motion and to rotation invariance when in three dimensions, give a Dirac equation set; the procedure and conclusion are directly extendible to arbitrarily varying potentials by use of the Furious theorem and to particle motions in three dimensions by virtue of the characteristics of de Broglie particle motion. Through the derivation of the Dirac equation, the study hopes to lend insight into the connections between the Dirac wave functions and the electrodynamic components of simple particles under the government by the well established basic laws of electrodynamics.

  13. Particle Physics Primer: Explaining the Standard Model of Matter.

    ERIC Educational Resources Information Center

    Vondracek, Mark

    2002-01-01

    Describes the Standard Model, a basic model of the universe that describes electromagnetic force, weak nuclear force radioactivity, and the strong nuclear force responsible for holding particles within the nucleus together. (YDS)

  14. Modeling of Oxidation of Molybdenum Particles during Plasma Spray Deposition

    SciTech Connect

    Fincke, James Russell; Wan, Y. P.; Jiang, X. Y.; Sampath, S.; Prasad, V.; Herman, H.

    2001-06-01

    An oxidation model for molybdenum particles during the plasma spray deposition process is presented. Based on a well-verified model for plasma chemistry and the heating and phase change of particles in a plasma plume, this model accounts for the oxidant diffusion around the surface of particles or splats, oxidation on the surface, as well as oxygen diffusion in molten molybdenum. Calculations are performed for a single molybdenum particle sprayed under Metco-9MB spraying conditions. The oxidation features of particles during the flight are compared with those during the deposition. The result shows the dominance of oxidation of a molybdenum particle during the flight, as well as during deposition when the substrate temperature is high (above 400 °C).

  15. Predictive modeling of particle-laden turbulent flows. Final report

    SciTech Connect

    Shaffer, F.; Bolio, E.J.; Hrenya, C.M.

    1993-12-31

    Earlier work of Sinclair and Jackson which treats the laminar flow of gas-solid suspensions is extended to model dilute turbulent flow. The random particle motion, often exceeding the turbulent fluctuations in the gas, is obtained using a model based on kinetic theory of granular materials. A two-equation low Reynolds number turbulence model is, modified to account for the presence of the dilute particle phase. Comparisons of the model predictions with available experimental data for the mean and fluctuating velocity profiles for both phases indicate that the resulting theory captures many of the flow features observed in the pneumatic transport of large particles. The model predictions did not manifest an extreme sensitivity to the degree of inelasticity in the particle-particle collisions for the range of solid loading ratios investigated.

  16. Investigation of particles size effects in Dissipative Particle Dynamics (DPD) modelling of colloidal suspensions

    NASA Astrophysics Data System (ADS)

    Mai-Duy, N.; Phan-Thien, N.; Khoo, B. C.

    2015-04-01

    In the Dissipative Particle Dynamics (DPD) simulation of suspension, the fluid (solvent) and colloidal particles are replaced by a set of DPD particles and therefore their relative sizes (as measured by their exclusion zones) can affect the maximal packing fraction of the colloidal particles. In this study, we investigate roles of the conservative, dissipative and random forces in this relative size ratio (colloidal/solvent). We propose a mechanism of adjusting the DPD parameters to properly model the solvent phase (the solvent here is supposed to have the same isothermal compressibility to that of water).

  17. Particle Modelling of Fluid Phenomena in Three -

    NASA Astrophysics Data System (ADS)

    Mahmoudi, Mohsen

    A new, numerical approach is developed to simulate fluid phenomena by means of molecular type behavior. First, consider the large number of molecules to be approximated by a smaller number of aggregates called particles. Then, let the particles interact with each other according to a classical molecular type force vec F whose magnitude F is given (Hirchfelder, Curtiss and Bird (1954)) by: F = rm -{Gover r^{p}} + {Hover r^{q}}, in which G, H, p, q, are positive constants and r is the distance between two particles. The acceleration of each particle is related to the force by the Newtonian dynamical equations vec F = m vec a. Displacement, velocity, and acceleration of each particle are then approximated by the "Leap Frog" formulas. The CRAY X-MP/24 is used to solve numerically the resulting large system of nonlinear, ordinary differential equations. We then study fluid phenomena in the following order. Part 1. Generation of particle fluids in a cylindrical region. Part 2. Verification of basic fluid properties. Part 3. Simulation of surface motion. In this part, we simulate three phenomena, which can be observed physically, by dropping a small object into a container filled with liquid. First, there is a backdrop. Then, a wave will be generated and going outward from the point of entry of the object into the container. Last, a reaction which can be recorded only with a high speed camera (Trefethen (1972)) is that very small drops of the container fluid may actually pinch off from the backdrop. Part 4. Simulation of surface tension. This phenomena can be observed by performing the following experiment. A small needle placed gently upon a water surface will not be sunk but will be supported by the molecular forces in the liquid surface. The molecules in the surface are depressed slightly in the process.

  18. Hydrodynamic model for particle size segregation in granular media

    NASA Astrophysics Data System (ADS)

    Trujillo, Leonardo; Herrmann, Hans J.

    2003-12-01

    We present a hydrodynamic theoretical model for “Brazil nut” size segregation in granular materials. We give analytical solutions for the rise velocity of a large intruder particle immersed in a medium of monodisperse fluidized small particles. We propose a new mechanism for this particle size-segregation due to buoyant forces caused by density variations which come from differences in the local “granular temperature”. The mobility of the particles is modified by the energy dissipation due to inelastic collisions and this leads to a different behavior from what one would expect for an elastic system. Using our model we can explain the size ratio dependence of the upward velocity.

  19. Merging for Particle-Mesh Complex Particle Kinetic Modeling of the Multiple Plasma Beams

    NASA Technical Reports Server (NTRS)

    Lipatov, Alexander S.

    2011-01-01

    We suggest a merging procedure for the Particle-Mesh Complex Particle Kinetic (PMCPK) method in case of inter-penetrating flow (multiple plasma beams). We examine the standard particle-in-cell (PIC) and the PMCPK methods in the case of particle acceleration by shock surfing for a wide range of the control numerical parameters. The plasma dynamics is described by a hybrid (particle-ion-fluid-electron) model. Note that one may need a mesh if modeling with the computation of an electromagnetic field. Our calculations use specified, time-independent electromagnetic fields for the shock, rather than self-consistently generated fields. While a particle-mesh method is a well-verified approach, the CPK method seems to be a good approach for multiscale modeling that includes multiple regions with various particle/fluid plasma behavior. However, the CPK method is still in need of a verification for studying the basic plasma phenomena: particle heating and acceleration by collisionless shocks, magnetic field reconnection, beam dynamics, etc.

  20. Confined turbulent fluid-particle flow modeling using multiple-realization particle trajectory schemes

    NASA Astrophysics Data System (ADS)

    Adeniji-Fashola, A. A.

    1988-07-01

    A multiple-realization particle trajectory scheme has been developed and applied to the numerical prediction of confined turbulent fluid-particle flows. The example flows investigated include the vertical pipe upflow experimental data of Tsuji et al. and the experimental data of Leavitt for a coaxial jet flow, comprising a particle-laden central jet and a clean annular jet, into a large recirculation chamber. The results obtained from the numerical scheme agree well with the experimental data, lending confidence to the modeling approach. The multiple-realization particle trajectory turbulent flow modeling scheme is believed to be a more elegant and accurate approach to the extension of single-particle hydrodynamics to dilute multi-particle systems than the more commonly employed two-fluid modeling approach. It is also better able to incorporate additional force items such as lift, virtual mass and Bassett history terms directly into the particle equation of motion as appropriate. This makes it a suitable candidate for particle migration studies and an extension to situations involving liquid particulate phases with possible propulsion applications, such as in spray combustion, follows naturally.

  1. Confined turbulent fluid-particle flow modeling using multiple-realization particle trajectory schemes

    NASA Technical Reports Server (NTRS)

    Adeniji-Fashola, A. A.

    1988-01-01

    A multiple-realization particle trajectory scheme has been developed and applied to the numerical prediction of confined turbulent fluid-particle flows. The example flows investigated include the vertical pipe upflow experimental data of Tsuji et al. and the experimental data of Leavitt for a coaxial jet flow, comprising a particle-laden central jet and a clean annular jet, into a large recirculation chamber. The results obtained from the numerical scheme agree well with the experimental data, lending confidence to the modeling approach. The multiple-realization particle trajectory turbulent flow modeling scheme is believed to be a more elegant and accurate approach to the extension of single-particle hydrodynamics to dilute multi-particle systems than the more commonly employed two-fluid modeling approach. It is also better able to incorporate additional force items such as lift, virtual mass and Bassett history terms directly into the particle equation of motion as appropriate. This makes it a suitable candidate for particle migration studies and an extension to situations involving liquid particulate phases with possible propulsion applications, such as in spray combustion, follows naturally.

  2. Application of The Stochastic Particle Tracking Model To Evaluate Particle Movement Uncertainty in Extreme Flows

    NASA Astrophysics Data System (ADS)

    Tsai, C.; Lin, E.

    2014-12-01

    In this study, modeling of suspended sediment particle movement in extreme flows is proposed by stochastic particle tracking modeling approaches. The proposed stochastic model is governed by a stochastic differential equation (SDE) composed of two random processes (a Wiener process and a Poisson process), and a random variable (i.e., flow magnitude) simulated by the extreme value Type I distribution. An extreme flow is defined as a hydrologic flow event (such as a flash flood) or a large flow perturbation with a low probability of occurrence and a high impact on its ambient flow environment. In the proposed particle tracking model, a random term mainly caused by fluid eddy motions is modeled as a Wiener process, while the random occurrences of a sequence of extreme flows can be modeled as a Poisson process. Following previous work by Oh and Tsai (2010)[1] and Tsai et al. (2014)[2], this study is intended to modify the jump term, which models the abrupt changes of particle position in the extreme flow environments. It is proposed that the probabilistic magnitude of extreme events can be simulated by the extreme value type I (EV I) distribution. The ensemble mean and variance of particle trajectory can be obtained from the proposed stochastic models via simulations. Our findings suggest that the ability to consider the probabilistic magnitude of extreme events can provide a more comprehensive and realistic estimate of the uncertainty of particle movement when extreme flow events occur. It is also found that the variance of particle position may be attributed to both the random magnitudes and occurrences of particle jumps in the presence of extreme flow events. It is demonstrated from this study that the proposed model can more explicitly quantify the uncertainty of particle movement by taking into considerations both the random arrival process of extreme flows and the variability of the extreme flow magnitude. [1] Oh, J. S., and Tsai, C.W.(2010). "A stochastic jump

  3. A Particle-based model for water simulation

    SciTech Connect

    Max, N., LLNL

    1998-01-01

    The Smooth-Particle Applied Mechanics (SPAM) model is a relatively recent physical modeling technique It can model both fluids and solids using free-moving particles An implemented SPAM model is described that solved the compressible Navier-Stokes equations to produce animations of splashing and pooling water Because the particle positions are known explicitly each timestep, the SPAM technique produces data amenable to visualization A ray-tracing renderer is also described It samples the underwater light-field distribution and stores tbe information into a Light Accumulation Lattice which is used for scattered light calculations and caustics.

  4. Viscous FRW Models with Particle Creation in Early Universe

    NASA Astrophysics Data System (ADS)

    Singh, C. P.

    2012-05-01

    We discuss the dynamical effects of bulk viscosity and particle creation on the early evolution of the Friedmann-Robertson-Walker model in the framework of open thermodynamical systems. We consider bulk viscosity and particle creation as separate irreversible processes. Exact solutions of the Einstein field equations are obtained by using the "gamma-law" equation of state p = (γ-1)ρ, where the adiabatic parameter γ varies with scale factor of the metric. We consider the cosmological model to study the early phases of the evolution of the universe as it goes from an inflationary phase to a radiation-dominated era in the presence of bulk viscosity and particle creation. Analytical solutions are obtained for particle number density and entropy for all models. It is found that, by choosing appropriate functions for particle creation rate and bulk viscous coefficient, the models exhibit singular and non-singular beginnings.

  5. Dirac Equation for Electrodynamic Model Particle

    NASA Astrophysics Data System (ADS)

    Zheng-Johansson, J. X.

    2008-03-01

    We set up the Maxwell's equations and subsequently the classical wave equations for the electromagnetic waves which together with their generating source, an oscillatory charge of zero rest mass, make up a particle travelling at velocity v as with the charge in the fields of an external scalar and vector potentials. The direct solutions in constant external field are Doppler-displaced plane waves propagating at the velocity of light c; at the de Broglie wavelength scale and expressed in terms of the dynamically equivalent and appropriate geometric mean wave variables, these render as functons identical to the space-time functions of the Dirac spinor, and these are identical to the de Broglie phase waves given previously from explicit superposition. For two spin-half particles of a common set of space-time functions constrained with antisymmetric spin functions as follows the Pauli principle for same charges and as separately indirectly induced based on experiment for opposite charges, the complete wave functions are identical to a Dirac spinor. The back-substitution of the so explicitly determined complete wave functions in the corresponding classical wave equations of the two particles, subjected further to reductions appropriate for the stationary- state particle motion and to rotation invariance when in three dimensions, give a Dirac equation set; the procedure and conclusion are directly extendible to arbitrarily varying potentials by use of the Furious theorem and to three dimensions (full paper: QTS5).

  6. Fractal Rings and Composite Elementary Particles (FRACEP): A Picture of Composite Standard Model Fundamental Particles

    NASA Astrophysics Data System (ADS)

    Giannini, Judith

    2016-03-01

    The object of this work was to study the feasibility of identifying a minimum set of fundamental particles that could be used to build up composite fermions and bosons that exhibit the same properties and behavior as the Standard Model (SM) fundamental particles. The spontaneous decay of most of the SM fermions suggests the possibility that they are composite in nature. The results of this arithmetically-based conceptual model identify a minimum set of only two fundamental particles (with equal and opposite mass) that combine in fractal-like configurations to form Intermediate Building Blocks (IMB). The IMBs then combine to form all of the SM fundamental particles and their anti-counterparts. These composite (bright universe) particles agree with the SM particles in mass, spin, electric charge, decay products and maximum classical radius (indicated by the scattering cross-section). Further, FRACEP identifies an equal set of dark universe particles, based primarily on its negative fundamental particle, which could represent the dark matter and energy understood to be the cause of the expansion of our (bright) universe.

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

  8. Hysteresis modeling and measurement for two-dimensional particle assemblies

    NASA Astrophysics Data System (ADS)

    Hauser, H.; Fulmek, P. L.; Grössinger, R.

    2002-04-01

    The increasing accuracy of circuit data storage simulations demands reliable models for the magnetic behaviour of the magnetic storage material. This paper introduces and compares the results of measurements and the results of model calculations by applying the Jiles-Atherton model, and the energetic model of ferromagnetic hysteresis by Hauser. The results show good agreement for uniaxial particle assemblies.

  9. A Simplified Model for the Acceleration of Cosmic Ray Particles

    ERIC Educational Resources Information Center

    Gron, Oyvind

    2010-01-01

    Two important questions concerning cosmic rays are: Why are electrons in the cosmic rays less efficiently accelerated than nuclei? How are particles accelerated to great energies in ultra-high energy cosmic rays? In order to answer these questions we construct a simple model of the acceleration of a charged particle in the cosmic ray. It is not…

  10. Optical properties of soot particles: measurement - model comparison

    NASA Astrophysics Data System (ADS)

    Forestieri, S.; Lambe, A. T.; Lack, D.; Massoli, P.; Cross, E. S.; Dubey, M.; Mazzoleni, C.; Olfert, J.; Freedman, A.; Davidovits, P.; Onasch, T. B.; Cappa, C. D.

    2013-12-01

    Soot, a product of incomplete combustion, plays an important role in the earth's climate system through the absorption and scattering of solar radiation. In order to accurately model the direct radiative impact of black carbon (BC), the refractive index and shape dependent scattering and absorption characteristics must be known. At present, the assumed shape remains highly uncertain because BC particles are fractal-like, being agglomerates of smaller (20-40 nm) spherules, yet traditional optical models such as Mie theory typically assume a spherical particle morphology. To investigate the ability of various optical models to reproduce observed BC optical properties, we measured light absorption and extinction coefficients of methane and ethylene flame soot particles. Optical properties were measured by multiple instruments: absorption by a dual cavity ringdown photoacoustic spectrometer (CRD-PAS), absorption and scattering by a 3-wavelength photoacoustic/nephelometer spectrometer (PASS-3) and extinction and scattering by a cavity attenuated phase shift spectrometer (CAPS). Soot particle mass was quantified using a centrifugal particle mass analyzer (CPMA) and mobility size was measured with a scanning mobility particle sizer (SMPS). Measurements were made for nascent soot particles and for collapsed soot particles following coating with dioctyl sebacate or sulfuric acid and thermal denuding to remove the coating. Wavelength-dependent refractive indices for the sampled particles were derived by fitting the observed absorption and extinction cross-sections to spherical particle Mie theory and Rayleigh-Debye-Gans theory. The Rayleigh-Debye-Gans approximation assumes that the absorption properties of soot are dictated by the individual spherules and neglects interaction between them. In general, Mie theory reproduces the observed absorption and extinction cross-sections for particles with volume equivalent diameters (VED) < ~160 nm, but systematically predicts lower

  11. Calibration of TSI model 3025 ultrafine condensation particle counter

    SciTech Connect

    Kesten, J.; Reineking, A.; Porstendoerfer, J. )

    1991-01-01

    The registration efficiency of the TSI model 3025 ultrafine condensation particle counter for Ag and NaCl particles of between 2 and 20 nm in diameter was determined. Taking into account the different shapes of the input aerosol size distributions entering the differential mobility analyzer (DMA) and the transfer function of the DMA, the counting efficiencies of condensation nucleus counters (CNC) for monodisperse Ag and NaCl particles were estimated. In addition, the dependence of the CNC registration efficiency on the particle concentration was investigated.

  12. Lagrangian modeling of large volcanic particles: Application to Vulcanian explosions

    NASA Astrophysics Data System (ADS)

    de'Michieli Vitturi, M.; Neri, A.; Esposti Ongaro, T.; Lo Savio, S.; Boschi, E.

    2010-08-01

    A new 2D/3D Lagrangian particle model (named LPAC) for the dynamics of clasts ejected during explosive eruptions is presented. The novelty of the model lies in the one-way coupling of the carrier flow field, given by a Eulerian multiphase flow code, and the particles. The model is based on a simplification of the Basset-Boussinesq-Oseen equation, expressing the Lagrangian equation of a particle as the sum of the forces exerted on it along its trajectory. It is assumed that particles are non-interacting and do not affect the background carrier flow and that the drag coefficient is constant. The model was applied to large clasts produced by Vulcanian explosions, in particular those occurring in August 1997 at Soufrière Hills Volcano, Montserrat (West Indies, UK). Simulation results allowed parametric studies as well as semi-quantitative comparisons between modeling results and field evidence. Major results include (1) the carrier flow was found to play a fundamental role even for meter-sized particles—a 1 m diameter block is predicted to reach a distance that is about 70% greater than that predicted without the effect of the carrier flow (assuming the same initial velocity), (2) assumption of the initial velocity of the particle was dropped thanks to the description of both the acceleration and deceleration phases along the particle trajectory, (3) by adopting experimentally based drag coefficients, large particles were able to reach greater distances with respect to smaller particles consistently with field observations and (4) the initial depth of the particle in the conduit was found to mainly influence the ejection velocity while the initial radial position with respect to the conduit axis was found to play a major role on the distance reached by the particle.

  13. Quasilinear Model for Energetic Particles Interacting with TAE Modes

    NASA Astrophysics Data System (ADS)

    Ghantous, Katy; Gorelenkov, Nikolai; Berk, Herbert

    2010-11-01

    TAE instabilities are thought to be a major source of Energetic Particle transport which could set limits on operational scenarios, especially for burning plasmas, and causes damage to the first wall. The quasilinear model proposed by Berk et al.ootnotetextH. L. Berk et al, Nucl. Fusion, 35:1661, 1995. relies on diffusion mechanisms for particle dynamics to captures the evolution of the energetic particle distribution function and the associated mode amplitude. Using the bump-on-tail as a paradigm, we analyze the dynamics near the resonances for accurate diffusion coefficient representation. We verify the model to get the predicted single mode saturation levels and benchmark the case of multimode overlap against particle codes. Using the TAE mode structures computed by the ideal MHD code NOVA, we generalize this method to relax energetic particles' profiles in the full 3D phase space.

  14. Derivation of manifestly covariant quantum models for spinning relativistic particles

    NASA Astrophysics Data System (ADS)

    Marnelius, Robert; Mårtensson, Ulf

    1990-05-01

    A method to construct manifestly covariant models for relativistic spinning particles is given. The models involve manifestly covariant internal variables leading to discrete state spaces apart from the coordinate and momentum variables. The precise form of the models are extracted from the Bargmann-Wigner conditions on the Pauli-Lubanski operator, together with some consistency conditions. Several simple models are derived and analysed, some of which are new. Also manifestly conformally invariant models for particles with arbitrary spins are derived using a condition of Bracken and Jessup.

  15. Modelling of externally mixed particles in the atmosphere

    NASA Astrophysics Data System (ADS)

    ZHU, Shupeng; Sartelet, Karine; Seigneur, Christian

    2014-05-01

    Particles present in the atmosphere have significant impacts on climate as well as on human health. Thus, it is important to accurately simulate and forecast their concentrations. Most commonly used air quality models assume that particles are internally mixed, largely for computational reasons. However, this assumption is disproved by measurements, especially close to sources. In fact, the externally-mixed properties of particles are important for aerosol source identification, radiative effects and particle evolution. In this study, a new size-composition resolved aerosol model is developed. It can solve the aerosol dynamic evolution for external mixtures taking into account the processes of coagulation, condensation and nucleation. Both the size of particles and the mass fraction of each chemical compound are discretized. For a given particle size, particles of different chemical composition may co-exist. Aerosol dynamics is solved in each grid cell by splitting coagulation and condensation/evaporation-nucleation processes. For the condensation/evaporation, surface equilibrium between gas and aerosol is calculated based on ISORROPIA and the newly developed H2O (Hydrophilic/Hydrophobic Organic) Model. Because size and chemical composition sections evolve during condensation/evaporation, concentrations need to be redistributed on fixed sections after condensation/evaporation to be able to use the model in 3 dimensions. This is done based on the numerical scheme HEMEN, which was initially developed for size redistribution. Chemical components can be grouped into several aggregates to reduce computational cost. The 0D model is validated by comparison to results obtained for internally mixed particles and the effect of mixing is investigated for up to 31 species and 4 aggregates. The model will be integrated into the air quality modeling platform POLYPHEMUS to investigate its performance in modeling air quality by comparing with observations during the MEGAPOLI

  16. Model independence of constraints on particle dark matter

    SciTech Connect

    Griest, K.; Sadoulet, B.

    1989-03-01

    The connection between the annihilation, elastic, and production cross sections is reviewed, showing how a general lower limit on the interaction rate in a detector is obtained from the requirement that a particle be the dark matter. High energy production experiments further constrain models, making very light dark matter particles unlikely. Special attention is paid to the uncertainties, loopholes and model dependencies that go into the arguments and several examples are given. 12 refs., 6 figs.

  17. The ESyS_Particle: A New 3-D Discrete Element Model with Single Particle Rotation

    NASA Astrophysics Data System (ADS)

    Wang, Yucang; Mora, Peter

    In this paper, the Discrete Element Model (DEM) is reviewed, and the ESyS_Particle, our new version of DEM, is introduced. We particularly highlight some of the major physical concerns about DEMs and major differences between our model and most current DEMs. In the new model, single particle rotation is introduced and represented by a unit quaternion. For each 3-D particle, six degrees of freedom are employed: three for translational motion, and three for orientation. Six kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial, two shearing forces, twisting and two bending torques. The relative rotation between two particles is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive rigid bodies can be uniquely determined. This algorithm can give more accurate results because physical principles are obeyed. A theoretical analysis about how to choose the model parameters is presented. Several numerical tests have been carried out, the results indicate that most laboratory tests can be well reproduced using our model.

  18. A hybrid mathematical model for controlling particle size, particle size distribution, and color properties of toner particles

    NASA Astrophysics Data System (ADS)

    Ataeefard, Maryam; Shadman, Alireza; Saeb, Mohammad Reza; Mohammadi, Yousef

    2016-08-01

    A mathematical modeling approach was proposed combining the capabilities of response surface methodology (RSM) and desirability function (DF) and implemented successfully in production of printing toner particles. Toner powders were systematically synthesized through suspension copolymerization process. Applying RSM, a series of experiments were designed and toner particles were prepared and the effects of monomer ratio, colorant and surfactant content on the particle size (PS), particle size distribution (PSD), thermal and colorimetric properties (∆ E) of the resulting toner were monitored and discussed. The second-order models corresponding to each target characteristic, i.e., PS, PSD, and ∆ E of different types of toner powders, were obtained by individual optimization to express variation of each property in terms of polymerization parameters. Applying statistical calculations, the best reduced models were identified to be fed in the second step of optimization. Since toners with appropriate PS, PSD, and CP were needed, we applied multi-objective optimization based on DF approach. The results show that exact tuning of toner properties is closely possible with the aid of hybrid mathematical model developed in this work. Noticeably, desirabilities are very close to 100 %.

  19. Relativistic models in nuclear and particle physics

    SciTech Connect

    Coester, F.

    1988-01-01

    A comparative overview is presented of different approaches to the construction of phenomenological dynamical models that respect basic principles of quantum theory and relativity. Wave functions defined as matrix elements of products of field operators on one hand and wave functions that are defined as representatives of state vectors in model Hilbert spaces are related differently to observables and dynamical models for these wave functions have each distinct advantages and disadvantages 34 refs.

  20. Constant Energy of Time Involute Particles of Biharmonic Particles in Bianchi Type-I Cosmological Model Spacetime

    NASA Astrophysics Data System (ADS)

    Körpinar, Talat; Turhan, Essin

    2015-05-01

    In this paper, we study energy of time involute particles of biharmonic particles in Bianchi type-I (B-I) cosmological model spacetime. We give a geometrical description of energy for a Frenet vector fields of timelike biharmonic particle. Finally, using the Frenet frame of the given particle, we obtain different cases for this particles and give important characterizations about them in Bianchi type-I (B-I) cosmological model spacetime.

  1. Surface tension models for particle laden thin films

    NASA Astrophysics Data System (ADS)

    Wong, Jeffrey; Wang, Li; Bertozzi, Andrea

    We study viscous slurries on an incline, for which particles migrate in a fluid due to a combination of gravity-induced settling and shear-induced migration. The lubrication model for the bulk of the fluid is a hyperbolic system of conservation laws for the film height and particle concentration which exhibits in interesting behavior, including singular shock solutions corresponding to accumulation of particles at the front. The addition of surface tension to the model produces a a capillary ridge that is affected by the particle accumulation and in two dimensions leads to fingering instabilities. We compare this model to experimental results. This work is supported by NSF Grants DMS-1312543 and DMS-1045536.

  2. Statistical modeling for particle impact noise detection testing

    SciTech Connect

    Prairie, R.R. ); Zimmer, W.J. )

    1990-01-01

    Particle Impact Noise Detection (PIND) testing is widely used to test electronic devices for the presence of conductive particles which can cause catastrophic failure. This paper develops a statistical model based on the rate of particles contaminating the part, the rate of particles induced by the test vibration, the escape rate, and the false alarm rate. Based on data from a large number of PIND tests for a canned transistor, the model is shown to fit the observed results closely. Knowledge of the parameters for which this fit is made is important in evaluating the effectiveness of the PIND test procedure and for developing background judgment about the performance of the PIND test. Furthermore, by varying the input parameters to the model, the resulting yield, failure rate and percent fallout can be examined and used to plan and implement PIND test programs.

  3. A particle-grid air quality modeling approach

    SciTech Connect

    Chock, D.P.; Winkler, S.L.

    1996-12-31

    A particle-grid air quality modeling approach that can incorporate chemistry is proposed as an alternative to the conventional PDF-grid air quality modeling. The particle trajectory model can accurately describe advection of air pollutants without introducing artificial diffusion, generating negative concentrations or distorting the concentration distributions. It also accurately describes the dispersion of emissions from point sources and is capable of retaining subgrid-scale information. Inhomogeneous turbulence necessitates use of a small timestep, say, 10 s to describe vertical dispersion of particles in convective conditions. A timestep as large as 200 s can be used to simulate horizontal dispersion. A time-splitting scheme can be used to couple the horizontal and vertical dispersion in a 3D simulation, and about 2000-3000 particles per cell of size 5 km x 5 km X 50 m is sufficient to yield a highly accurate simulation of 3D dispersion. Use of an hourly-averaged concentration further reduces the demand of particle per cell to 500. The particle-grid method is applied to a system of ten reacting chemical species in a two-dimensional rotating flow field with and without diffusion. A chemistry grid within which reactions are assumed to take place can be decoupled from the grid describing the flow field. Two types of chemistry grids are used to describe the chemical reactions: a fixed coarse grid and a moving (the advection case) or stationary (the advection plus diffusion case) fine grid. Two particle-number densities are also used: 256 and 576 particles per fixed coarse grid cell. The species mass redistributed back to the particle after each reaction step is assumed to be proportional to the species mass in the particle before the reaction. The simulation results are very accurate, especially in the advection-chemistry case. Accuracy improves with the use of a fine grid.

  4. Mesoscopic dynamics of colloids simulated with dissipative particle dynamics and fluid particle model.

    PubMed

    Dzwinel, Witold; Yuen, David A; Boryczko, Krzysztof

    2002-01-01

    We report results of numerical simulations of complex fluids, using a combination of discrete-particle methods. Our molecular modeling repertoire comprises three simulation techniques: molecular dynamics (MD), dissipative particle dynamics (DPD), and the fluid particle model (FPM). This type of model can depict multi-resolution molecular structures found in complex fluids ranging from single micelle, colloidal crystals, large-scale colloidal aggregates up to the mesoscale processes of hydrodynamical instabilities in the bulk of colloidal suspensions. We can simulate different colloidal structures in which the colloidal beds are of comparable size to the solvent particles. This undertaking is accomplished with a two-level discrete particle model consisting of the MD paradigm with a Lennard-Jones (L-J) type potential for defining the colloidal particle system and DPD or FPM for modeling the solvent. We observe the spontaneous emergence of spherical or rod-like micelles and their crystallization in stable hexagonal or worm-like structures, respectively. The ordered arrays obtained by using the particle model are similar to the 2D colloidal crystals observed in laboratory experiments. The micelle shape and its hydrophobic or hydrophilic character depend on the ratio between the scaling factors of the interactions between colloid-colloid to colloid-solvent. Unlike the miscellar arrays, the colloidal aggregates involve the colloid-solvent interactions prescribed by the DPD forces. Different from the assumption of equilibrium growth, the two-level particle model can display much more realistic molecular physics, which allows for the simulation of aggregation for various types of colloids and solvent liquids over a very broad range of conditions. We discuss the potential prospects of combining MD, DPD, and FPM techniques in a single three-level model. Finally, we present results from large-scale simulation of the Rayleigh-Taylor instability and dispersion of colloidal slab

  5. Particle-Turbulence Interaction Model for Aluminum Combustion

    NASA Astrophysics Data System (ADS)

    Sinha, Neeraj; Calhoon, William; Tomes, Jeremy

    2011-06-01

    Particle-turbulence interactions will have a substantial impact on the performance of thermobaric explosives that rely on the particle combustion for secondary heat release. Modeling these interactions from a fundamental perspective is very difficult and intractable for large-scale problems of practical interest. Alternatively, these interactions may be modeled from a macroscopic perspective that seeks to account for the probability distribution function (PDF) of variables within the modeled laminar burning rate for the particulates. Such a formulation would account for the first order effect of turbulent fluctuations on the burning rate within a computationally affordable model. This paper will describe the development of such a model for aluminum particle combustion in both the diffusion and kinetic burning regimes. This formulation is based on an assumed PDF method that may be parameterized into a database that may be deployed within a flow solver. As a result, the formulation is computational efficient and affordable for large-scale simulations.

  6. A micromechanical modelling approach for predicting particle dislodgement

    NASA Astrophysics Data System (ADS)

    Shih, W.; Diplas, P.

    2015-03-01

    In recent years, the impulse and energy concepts have been proposed as a means of improving the prediction of sediment particle entrainment. Because they account for temporal aspects of the applied hydrodynamic forces, impulse-based or energy-based approaches provide a better representation of a highly fluctuating phenomenon compared to the more traditional time-averaged Shields parameter or constant force magnitude-based model. In the present paper, we elaborate further on the simulation of particle migration by applying a theoretical model based on hydrodynamic pressure data. Three common movement scenarios and the dynamic critical drag force, allowed to vary during particle motion, are qualitatively discussed here. Also, the possible inaccuracy in entrainment prediction of three different models is addressed with reference to our simulation results. It is demonstrated that the multi-parcel impulse model has the best performance when compared to the single-parcel impulse and force magnitude-based models.

  7. One-dimensional particle models for heat transfer analysis

    NASA Astrophysics Data System (ADS)

    Bufferand, H.; Ciraolo, G.; Ghendrih, Ph; Tamain, P.; Bagnoli, F.; Lepri, S.; Livi, R.

    2010-11-01

    For a better understanding of Spitzer-Härm closure restrictions and for estimating the relevancy of this expression when collisionnality decreases, an effort is done in developing simple models that aim at catching the physics of the transition from conductive to free-streaming heat flux. In that perspective, one-dimensional particle models are developed to study heat transfer properties in the direction parallel to the magnetic field in tokamaks. These models are based on particles that carry energy at a specific velocity and that can interact with each other or with heat sources. By adjusting the particle dynamics and particle interaction properties, it is possible to generate a broad range of models of growing complexity. The simplest models can be solved analytically and are used to link particle behavior to general macroscopic heat transfer properties. In particular, some configurations recover Fourier's law and make possible to investigate the dependance of thermal conductivity on temperature. Besides, some configurations where local balance is lost require defining non local expression for heat flux. These different classes of models could then be linked to different plasma configurations and used to study transition from collisional to non-collisional plasma.

  8. ANFIS modeling for prediction of particle motions in fluid flows

    NASA Astrophysics Data System (ADS)

    Safdari, Arman; Kim, Kyung Chun

    2015-11-01

    Accurate dynamic analysis of parcel of solid particles driven in fluid flow system is of interest for many natural and industrial applications such as sedimentation process, study of cloud particles in atmosphere, etc. In this paper, numerical modeling of solid particles in incompressible flow using Eulerian-Lagrangian approach is carried out to investigate the dynamic behavior of particles in different flow conditions; channel and cavity flow. Although modern computers have been well developed, the high computational time and costs for this kind of problems are still demanded. The Lattice Boltzmann Method (LBM) is used to simulate fluid flows and combined with the Lagrangian approach to predict the motion of particles in the range of masses. Some particles are selected, and subjected to Adaptive-network-based fuzzy inference system (ANFIS) to predict the trajectory of moving solid particles. Using a hybrid learning procedure from computational particle movement, the ANFIS can construct an input-output mapping based on fuzzy if-then rules and stipulated computational fluid dynamics prediction pairs. The obtained results from ANFIS algorithm is validated and compared with the set of benchmark data provided based on point-like approach coupled with the LBM method.

  9. The Modeling of Pickup Ion or Energetic Particle Mediated Plasmas

    NASA Astrophysics Data System (ADS)

    Zank, G. P.; Mostafavi, P.; Hunana, P.

    2016-05-01

    Suprathermal energetic particles, such as solar energetic particles (SEPs) in the inner heliosphere and pickup ions (PUIs) in the outer heliosphere and the very local interstellar medium, often form a thermodynamically dominant component in their various environments. In the supersonic solar wind beyond > 10 AU, in the inner heliosheath (IHS), and in the very local interstellar medium (VLISM), PUIs do not equilibrate collisionally with the background plasma. Similarly, SEPs do not equilibrate collisionally with the background solar wind in the inner heliosphere. In the absence of equilibration between plasma components, a separate coupled plasma description for the energetic particles is necessary. Using a collisionless Chapman-Enskog expansion, we derive a closed system of multi-component equations for a plasma comprised of thermal protons and electrons, and suprathermal particles (SEPs, PUIs). The energetic particles contribute an isotropic scalar pressure to leading order, a collisionless heat flux at the next order, and a collisionless stress tensor at the second-order. The collisionless heat conduction and viscosity in the multi-fluid description results from a nonisotropic energetic particle distribution. A simpler single-fluid MHD-like system of equations with distinct equations of state for both the background plasma and the suprathermal particles is derived. We note briefly potential pitfalls that can emerge in the numerical modeling of collisionless plasma flows that contain a dynamically important energetic particle component.

  10. Dissipative particle dynamics model for colloid transport in porous media

    SciTech Connect

    Pan, W.; Tartakovsky, A. M.

    2013-08-01

    We present that the transport of colloidal particles in porous media can be effectively modeled with a new formulation of dissipative particle dynamics, which augments standard DPD with non-central dissipative shear forces between particles while preserving angular momentum. Our previous studies have demonstrated that the new formulation is able to capture accurately the drag forces as well as the drag torques on colloidal particles that result from the hydrodynamic retardation effect. In the present work, we use the new formulation to study the contact efficiency in colloid filtration in saturated porous media. Note that the present model include all transport mechanisms simultaneously, including gravitational sedimentation, interception and Brownian diffusion. Our results of contact efficiency show a good agreement with the predictions of the correlation equation proposed by Tufenkji and EliMelech, which also incorporate all transport mechanisms simultaneously without the additivity assumption.

  11. Multi-particle FEM modeling on microscopic behavior of 2D particle compaction

    NASA Astrophysics Data System (ADS)

    Zhang, Y. X.; An, X. Z.; Zhang, Y. L.

    2015-03-01

    In this paper, the discrete random packing and various ordered packings such as tetragonal and hexagonal close packed structures generated by discrete element method and honeycomb, which is manually generated were input as the initial packing structures into the multi-particle finite element model (FEM) to study their densification during compaction, where each particle is discretized as a FEM mesh. The macro-property such as relative density and micro-properties such as local morphology, stress, coordination number and densification mechanism obtained from various initial packings are characterized and analyzed. The results show that the coupling of discrete feature in particle scale with the continuous FEM in macro-scale can effectively conquer the difficulties in traditional FEM modeling, which provides a reasonable way to reproduce the compaction process and identify the densification mechanism more accurately and realistically.

  12. Space Charge Models for Particle Tracking on Long Time Scales

    SciTech Connect

    Holmes, Jeffrey A; Cousineau, Sarah M; Shishlo, Andrei P; Potts III, Robert E

    2013-01-01

    In order to efficiently track charged particles over long times, most tracking codes use either analytic charge distributions or particle-in-cell (PIC) methods based on fast Fourier transforms (FFTs). While useful for theoretical studies, analytic distribution models do not allow accurate simulation of real machines. PIC calculations can utilize realistic space charge distributions, but these methods suffer from the presence of discretization errors. We examine the situation for particle tracking with space charge over long times, and consider possible ideas to improve the accuracy of such calculations.

  13. Particle production within the quark meson coupling model

    SciTech Connect

    Panda, P. K.; Menezes, D. P.; Providencia, C.

    2009-07-15

    Quark-meson coupling (QMC) models can be successfully applied to the description of compact star properties in nuclear astrophysics as well as to nuclear matter. In the regime of hot hadronic matter very few calculations exist using the QMC model, in particular when applied to particle yields in heavy ion collisions. In the present work, we identify the free energy of the bag with the effective mass of the baryons and we calculate the particle production yields on a Au+Au collision at the BNL Relativistic Heavy Ion Collider (RHIC) with the QMC model and compare them with results obtained previously with other relativistic models. A smaller temperature for the fireball, T=132 MeV, is obtained because of the smaller effective baryon masses predicted by QMC. QMC was also applied to the description of particle yields at the CERN Super Proton Synchrotron (SPS) in Pb+Pb collisions.

  14. Isotopic overabundances and the energetic particle model of solar flares

    NASA Technical Reports Server (NTRS)

    Perez-Enriquez, R.; Bravo, S.

    1985-01-01

    According to the energetic particle model of solar flares particles are efficiently accelerated in the magnetic field loop of an active region (AR) by hydromagnetic turbulence. It is demonstrated that the isotopic overabundances observed in some flares are not a consequence of the flare itself but are characteristic of the plasma in the AR. Only when a flare releases the plasma into the interplanetary space it is possible to observe this anomalous composition at spacecraft locations.

  15. Smoothed Particle Hydrodynamics Model for Reactive Transport and Mineral Precipitation

    SciTech Connect

    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.

  16. THEORETICAL MODEL OF SOILING OF SURFACES BY AIRBORNE PARTICLES

    EPA Science Inventory

    A model is developed which can be used to predict the change in reflectance from a surface as a function of time. Reflectance change is a measure of soiling caused by the deposition of particles on a surface. The major inputs to the model are the parameters to a bimodal distribut...

  17. Modeling of Fine-Particle Formation in Turbulent Flames

    NASA Astrophysics Data System (ADS)

    Raman, Venkat; Fox, Rodney O.

    2016-01-01

    The generation of nanostructured particles in high-temperature flames is important both for the control of emissions from combustion devices and for the synthesis of high-value chemicals for a variety of applications. The physiochemical processes that lead to the production of fine particles in turbulent flames are highly sensitive to the flow physics and, in particular, the history of thermochemical compositions and turbulent features they encounter. Consequently, it is possible to change the characteristic size, structure, composition, and yield of the fine particles by altering the flow configuration. This review describes the complex multiscale interactions among turbulent fluid flow, gas-phase chemical reactions, and solid-phase particle evolution. The focus is on modeling the generation of soot particles, an unwanted pollutant from automobile and aircraft engines, as well as metal oxides, a class of high-value chemicals sought for specialized applications, including emissions control. Issues arising due to the numerical methods used to approximate the particle number density function, the modeling of turbulence-chemistry interactions, and model validation are also discussed.

  18. Effect of Particle Shape on Mechanical Behaviors of Rocks: A Numerical Study Using Clumped Particle Model

    PubMed Central

    Rong, Guan; Liu, Guang; Zhou, Chuang-bing

    2013-01-01

    Since rocks are aggregates of mineral particles, the effect of mineral microstructure on macroscopic mechanical behaviors of rocks is inneglectable. Rock samples of four different particle shapes are established in this study based on clumped particle model, and a sphericity index is used to quantify particle shape. Model parameters for simulation in PFC are obtained by triaxial compression test of quartz sandstone, and simulation of triaxial compression test is then conducted on four rock samples with different particle shapes. It is seen from the results that stress thresholds of rock samples such as crack initiation stress, crack damage stress, and peak stress decrease with the increasing of the sphericity index. The increase of sphericity leads to a drop of elastic modulus and a rise in Poisson ratio, while the decreasing sphericity usually results in the increase of cohesion and internal friction angle. Based on volume change of rock samples during simulation of triaxial compression test, variation of dilation angle with plastic strain is also studied. PMID:23997677

  19. Sticky Particles: Modeling Rigid Aggregates in Dense Planetary Rings

    NASA Astrophysics Data System (ADS)

    Perrine, Randall P.; Richardson, D. C.; Scheeres, D. J.

    2008-09-01

    We present progress on our study of planetary ring dynamics. We use local N-body simulations to examine small patches of dense rings in which self-gravity and mutual collisions dominate the dynamics of the ring material. We use the numerical code pkdgrav to model the motions of 105-7 ring particles, using a sliding patch model with modified periodic boundary conditions. The exact nature of planetary ring particles is not well understood. If covered in a frost-like layer, such irregular surfaces may allow for weak cohesion between colliding particles. Thus we have recently added new functionality to our model, allowing "sticky particles” to lock into rigid aggregates while in a rotating reference frame. This capability allows particles to adhere to one another, forming irregularly shaped aggregates that move as rigid bodies. (The bonds between particles can subsequently break, given sufficient stress.) These aggregates have greater strength than gravitationally bound "rubble piles,” and are thus able to grow larger and survive longer under similar stresses. This new functionality allows us to explore planetary ring properties and dynamics in a new way, by self-consistently forming (and destroying) non-spherical aggregates and moonlets via cohesive forces, while in a rotating frame, subjected to planetary tides. (We are not aware of any similar implementations in other existing models.) These improvements allow us to study the many effects that particle aggregation may have on the rings, such as overall ring structure; wake formation; equilibrium properties of non-spherical particles, like pitch angle, orientation, shape, size distribution, and spin; and the surface properties of the ring material. We present test cases and the latest results from this new model. This work is supported by a NASA Earth and Space Science Fellowship.

  20. Modelling new particle formation events in the South African savannah

    SciTech Connect

    Gierens, Rosa; Laakso, Lauri; Mogensen, Ditte; Vakkari, Ville; Buekes, Johan P.; Van Zyl, Pieter; Hakola, H.; Guenther, Alex B.; Pienaar, J. J.; Boy, Michael

    2014-05-28

    Africa is one of the less studied continents with respect to atmospheric aerosols. Savannahs are complex dynamic systems sensitive to climate and land-use changes, but the interaction of these systems with the atmosphere is not well understood. Atmospheric particles, called aerosols, affect the climate on regional and global scales, and are an important factor in air quality. In this study, measurements from a relatively clean savannah environment in South Africa were used to model new particle formation and growth. There already are some combined long-term measurements of trace gas concentrations together with aerosol and meteorological variables available, but to our knowledge this is the first detailed simulation that includes all the main processes relevant to particle formation. The results show that both of the particle formation mechanisms investigated overestimated the dependency of the formation rates on sulphuric acid. From the two particle formation mechanisms tested in this work, the approach that included low volatile organic compounds to the particle formation process was more accurate in describing the nucleation events than the approach that did not. To obtain a reliable estimate of aerosol concentration in simulations for larger scales, nucleation mechanisms would need to include organic compounds, at least in southern Africa. This work is the first step in developing a more comprehensive new particle formation model applicable to the unique environment in southern Africa. Such a model will assist in better understanding and predicting new particle formation – knowledge which could ultimately be used to mitigate impacts of climate change and air quality.

  1. A contoured continuum surface force model for particle methods

    NASA Astrophysics Data System (ADS)

    Duan, Guangtao; Koshizuka, Seiichi; Chen, Bin

    2015-10-01

    A surface tension model is essential to simulate multiphase flows with deformed interfaces. This study develops a contoured continuum surface force (CCSF) model for particle methods. A color function that varies sharply across the interface to mark different fluid phases is smoothed in the transition region, where the local contour curvature can be regarded as the interface curvature. The local contour passing through each reference particle in the transition region is extracted from the local profile of the smoothed color function. The local contour curvature is calculated based on the Taylor series expansion of the smoothed color function, whose derivatives are calculated accurately according to the definition of the smoothed color function. Two schemes are proposed to specify the smooth radius: fixed scheme, where 2 ×re (re = particle interaction radius) is assigned to all particles in the transition region; and varied scheme, where re and 2 ×re are assigned to the central and edged particles in the transition region respectively. Numerical examples, including curvature calculation for static circle and ellipse interfaces, deformation of square droplet to a circle (2D and 3D), droplet deformation in shear flow, and droplet coalescence, are simulated to verify the CCSF model and compare its performance with those of other methods. The CCSF model with the fixed scheme is proven to produce the most accurate curvature and lowest parasitic currents among the tested methods.

  2. Modeling nanoscale hydrodynamics by smoothed dissipative particle dynamics

    SciTech Connect

    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.

  3. Hybrid Modeling Method for a DEP Based Particle Manipulation

    PubMed Central

    Miled, Mohamed Amine; Gagne, Antoine; Sawan, Mohamad

    2013-01-01

    In this paper, a new modeling approach for Dielectrophoresis (DEP) based particle manipulation is presented. The proposed method fulfills missing links in finite element modeling between the multiphysic simulation and the biological behavior. This technique is amongst the first steps to develop a more complex platform covering several types of manipulations such as magnetophoresis and optics. The modeling approach is based on a hybrid interface using both ANSYS and MATLAB to link the propagation of the electrical field in the micro-channel to the particle motion. ANSYS is used to simulate the electrical propagation while MATLAB interprets the results to calculate cell displacement and send the new information to ANSYS for another turn. The beta version of the proposed technique takes into account particle shape, weight and its electrical properties. First obtained results are coherent with experimental results. PMID:23364197

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

  5. Modeling reactive transport with particle tracking and kernel estimators

    NASA Astrophysics Data System (ADS)

    Rahbaralam, Maryam; Fernandez-Garcia, Daniel; Sanchez-Vila, Xavier

    2015-04-01

    Groundwater reactive transport models are useful to assess and quantify the fate and transport of contaminants in subsurface media and are an essential tool for the analysis of coupled physical, chemical, and biological processes in Earth Systems. Particle Tracking Method (PTM) provides a computationally efficient and adaptable approach to solve the solute transport partial differential equation. On a molecular level, chemical reactions are the result of collisions, combinations, and/or decay of different species. For a well-mixed system, the chem- ical reactions are controlled by the classical thermodynamic rate coefficient. Each of these actions occurs with some probability that is a function of solute concentrations. PTM is based on considering that each particle actually represents a group of molecules. To properly simulate this system, an infinite number of particles is required, which is computationally unfeasible. On the other hand, a finite number of particles lead to a poor-mixed system which is limited by diffusion. Recent works have used this effect to actually model incomplete mix- ing in naturally occurring porous media. In this work, we demonstrate that this effect in most cases should be attributed to a defficient estimation of the concentrations and not to the occurrence of true incomplete mixing processes in porous media. To illustrate this, we show that a Kernel Density Estimation (KDE) of the concentrations can approach the well-mixed solution with a limited number of particles. KDEs provide weighting functions of each particle mass that expands its region of influence, hence providing a wider region for chemical reactions with time. Simulation results show that KDEs are powerful tools to improve state-of-the-art simulations of chemical reactions and indicates that incomplete mixing in diluted systems should be modeled based on alternative conceptual models and not on a limited number of particles.

  6. Modeling oceanic multiphase flow by using Lagrangian particle tracking

    NASA Astrophysics Data System (ADS)

    Matsumura, Y.

    2014-12-01

    While the density of seawater is basically determined by its temperature, salinity and pressure, the effective density becomes higher when the water mass contains suspended sediment. On the other hands, effective density declines when water mass contains fine scale materials of lower density such as bubbles and ice crystals. Such density anomaly induced by small scale materials suspended in water masses sometimes plays important roles in the sub-mesoscale ocean physics. To simulate these small scale oceanic multiphase flow, a new modeling framework using an online Lagrangian particle tracking method is developed. A Lagrangian particle tracking method has substantial advantages such as an explicit treatment of buoyancy force acting on each individual particle, no numerical diffusion and dissipation, high dynamic range and an ability to track the history and each individual particle. However, its numerical cost causes difficulty when we try to simulate a large number of particles. In the present study we implement a numerically efficient particle tracking scheme using linked-list data structure, which is coupled with a nonhydrostatic dynamical core. This newly developed model successfully reproduces characteristics of some interesting small scale multiphase processes, for example hyperpycnal flow (a sediment-rich river water plume trapped at ocean floor) and grease ice cover (a slurry mixture of frazil ice crystals and seawater).

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

  8. Simple model for fine particle (dust) clouds in plasmas

    NASA Astrophysics Data System (ADS)

    Totsuji, Hiroo

    2016-04-01

    In the cloud of fine particles (dusts) in plasmas, the charge neutrality can be much enhanced due to large charge numbers of fine particles. The required condition is not difficult to satisfy even when their charge density is substantially smaller than electrons or ions. Based on this fact, a simple model of fine particle clouds is proposed and the cloud radius is related to the half-width, the radius where the density of surrounding plasmas drops by half, in cylindrically and spherically symmetric cases under microgravity. When fine particles are gradually introduced with parameters of surrounding plasma especially the half-width being fixed, the size of clouds first increases and then saturates at the value determined by the plasma half-width, giving a possibility to control the size and density of clouds independently.

  9. Investigation of self-oscillation using particle balance model

    SciTech Connect

    Bae, Inshik; Na, Byungkeun Chang, Hongyoung

    2015-08-15

    Self-oscillation obtained using a DC-only power supply under specific anode voltage conditions is investigated in a cylindrical system with thermal electrons using tungsten filaments. Analysis of the obtained oscillation profiles reveals that the experimental data are consistent with a model derived from the particle balance model. The self-oscillation period characteristics with respect to the pressure and gas species are also analyzed. As the physics and particle motion of self-oscillation near the plasma transition region are analyzed from different perspectives, this paper may advance the study of this phenomenon.

  10. Investigation of Self-Oscillation using Particle Balance Model

    NASA Astrophysics Data System (ADS)

    Bae, Inshik; Na, Byungkeun; Chang, Hongyoung

    2015-09-01

    Self-oscillation, which is obtained by using a DC-only power supply with specific anode voltage conditions, is investigated in a cylindrical system with thermal electrons using tungsten filaments. From analysis of the obtained oscillation profiles, the experimental data is consistent with the model derived from the particle balance model. The self-oscillation period characteristics with respect to the pressure and gas species are also analyzed. As the physics and particle motion of self-oscillation near the electron avalanche is analyzed in different perspective, this study may advance the understanding of this phenomenon. This research was supported by the Ministry of Knowledge Economy (MKE) of Korea (Grant No. 10041681).

  11. Kinetic model of particle-inhibited grain growth

    NASA Astrophysics Data System (ADS)

    Thompson, Gary Scott

    The effects of second phase particles on matrix grain growth kinetics were investigated using Al2O3-SiC as a model system. In particular, the validity of the conclusion drawn from a previous kinetic analysis that the kinetics of particle-inhibited grain growth in Al2 O3-SiC samples with an intermediate volume fraction of second phase could be well quantified by a modified-Zener model was investigated. A critical analysis of assumptions made during the previous kinetic analysis revealed oversimplifications which affect the validity of the conclusion. Specifically, the degree of interaction between particles and grain boundaries was assumed to be independent of the mean second phase particle size and size distribution. In contrast, current measurements indicate that the degree of interaction in Al2O3-SiC is dependent on these parameters. An improved kinetic model for particle-inhibited grain growth in Al 2O3-SiC was developed using a modified-Zener approach. The comparison of model predictions with experimental grain growth data indicated that significant discrepancies (as much as 4--5 orders of magnitude) existed. Based on this, it was concluded that particles had a much more significant effect on grain growth kinetics than that caused by a simple reduction of the boundary driving force due to the removal of boundary area. Consequently, it was also concluded that the conclusion drawn from the earlier kinetic analysis regarding the validity of a modified-Zener model was incorrect. Discrepancies between model and experiment were found to be the result of a significant decrease in experimental growth rate constant not predicted by the model. Possible physical mechanisms for such a decrease were investigated. The investigation of a small amount of SiO2 on grain growth in Al2O3 indicated that the decrease was not the result of a decrease in grain boundary mobility due to impurity contamination by particles. By process of elimination and based on previous observations

  12. NMR relaxation induced by iron oxide particles: testing theoretical models.

    PubMed

    Gossuin, Y; Orlando, T; Basini, M; Henrard, D; Lascialfari, A; Mattea, C; Stapf, S; Vuong, Q L

    2016-04-15

    Superparamagnetic iron oxide particles find their main application as contrast agents for cellular and molecular magnetic resonance imaging. The contrast they bring is due to the shortening of the transverse relaxation time T 2 of water protons. In order to understand their influence on proton relaxation, different theoretical relaxation models have been developed, each of them presenting a certain validity domain, which depends on the particle characteristics and proton dynamics. The validation of these models is crucial since they allow for predicting the ideal particle characteristics for obtaining the best contrast but also because the fitting of T 1 experimental data by the theory constitutes an interesting tool for the characterization of the nanoparticles. In this work, T 2 of suspensions of iron oxide particles in different solvents and at different temperatures, corresponding to different proton diffusion properties, were measured and were compared to the three main theoretical models (the motional averaging regime, the static dephasing regime, and the partial refocusing model) with good qualitative agreement. However, a real quantitative agreement was not observed, probably because of the complexity of these nanoparticulate systems. The Roch theory, developed in the motional averaging regime (MAR), was also successfully used to fit T 1 nuclear magnetic relaxation dispersion (NMRD) profiles, even outside the MAR validity range, and provided a good estimate of the particle size. On the other hand, the simultaneous fitting of T 1 and T 2 NMRD profiles by the theory was impossible, and this occurrence constitutes a clear limitation of the Roch model. Finally, the theory was shown to satisfactorily fit the deuterium T 1 NMRD profile of superparamagnetic particle suspensions in heavy water. PMID:26933908

  13. NMR relaxation induced by iron oxide particles: testing theoretical models

    NASA Astrophysics Data System (ADS)

    Gossuin, Y.; Orlando, T.; Basini, M.; Henrard, D.; Lascialfari, A.; Mattea, C.; Stapf, S.; Vuong, Q. L.

    2016-04-01

    Superparamagnetic iron oxide particles find their main application as contrast agents for cellular and molecular magnetic resonance imaging. The contrast they bring is due to the shortening of the transverse relaxation time T 2 of water protons. In order to understand their influence on proton relaxation, different theoretical relaxation models have been developed, each of them presenting a certain validity domain, which depends on the particle characteristics and proton dynamics. The validation of these models is crucial since they allow for predicting the ideal particle characteristics for obtaining the best contrast but also because the fitting of T 1 experimental data by the theory constitutes an interesting tool for the characterization of the nanoparticles. In this work, T 2 of suspensions of iron oxide particles in different solvents and at different temperatures, corresponding to different proton diffusion properties, were measured and were compared to the three main theoretical models (the motional averaging regime, the static dephasing regime, and the partial refocusing model) with good qualitative agreement. However, a real quantitative agreement was not observed, probably because of the complexity of these nanoparticulate systems. The Roch theory, developed in the motional averaging regime (MAR), was also successfully used to fit T 1 nuclear magnetic relaxation dispersion (NMRD) profiles, even outside the MAR validity range, and provided a good estimate of the particle size. On the other hand, the simultaneous fitting of T 1 and T 2 NMRD profiles by the theory was impossible, and this occurrence constitutes a clear limitation of the Roch model. Finally, the theory was shown to satisfactorily fit the deuterium T 1 NMRD profile of superparamagnetic particle suspensions in heavy water.

  14. Charged Particle Environment Definition for NGST: Model Development

    NASA Technical Reports Server (NTRS)

    Blackwell, William C.; Minow, Joseph I.; Evans, Steven W.; Hardage, Donna M.; Suggs, Robert M.

    2000-01-01

    NGST will operate in a halo orbit about the L2 point, 1.5 million km from the Earth, where the spacecraft will periodically travel through the magnetotail region. There are a number of tools available to calculate the high energy, ionizing radiation particle environment from galactic cosmic rays and from solar disturbances. However, space environment tools are not generally available to provide assessments of charged particle environment and its variations in the solar wind, magnetosheath, and magnetotail at L2 distances. An engineering-level phenomenology code (LRAD) was therefore developed to facilitate the definition of charged particle environments in the vicinity of the L2 point in support of the NGST program. LRAD contains models tied to satellite measurement data of the solar wind and magnetotail regions. The model provides particle flux and fluence calculations necessary to predict spacecraft charging conditions and the degradation of materials used in the construction of NGST. This paper describes the LRAD environment models for the deep magnetotail (XGSE < -100 Re) and solar wind, and presents predictions of the charged particle environment for NGST.

  15. Generalized slave-particle method for extended Hubbard models

    NASA Astrophysics Data System (ADS)

    Georgescu, Alexandru B.; Ismail-Beigi, Sohrab

    2015-12-01

    We introduce a set of generalized slave-particle models for extended Hubbard models that treat localized electronic correlations using slave-boson decompositions. Our models automatically include two slave-particle methods of recent interest, the slave-rotor and slave-spin methods, as well as a ladder of new intermediate models where one can choose which of the electronic degrees of freedom (e.g., spin or orbital labels) are treated as correlated degrees of freedom by the slave bosons. In addition, our method removes the aberrant behavior of the slave-rotor model, where it systematically overestimates the importance of electronic correlation effects for weak interaction strength, by removing the contribution of unphysical states from the bosonic Hilbert space. The flexibility of our formalism permits one to separate and isolate the effect of correlations on the key degrees of freedom.

  16. Modeling Correlation Effects in Nickelates with Slave Particles

    NASA Astrophysics Data System (ADS)

    Georgescu, Alexandru Bogdan; Ismail-Beigi, Sohrab

    Nickelate interfaces display interesting electronic properties including orbital ordering similar to that of cuprate superconductors and thickness dependent metal-insulator transitions. One-particle band theory calculations do not include dynamic localized correlation effects on the nickel sites and thus often incorrectly predict metallic systems or incorrect ARPES spectra. Building on two previous successful slave-particle treatments of local correlations, we present a generalized slave-particle method that includes prior models and allows us to produce new intermediate models. The computational efficiency of these slave-boson methods means that one can readily study correlation effects in complex heterostructures. We show some predictions of these methods for the electronic structure of bulk and thin film nickelates. Work supported by NSF Grant MRSEC DMR-1119826.

  17. Explicit particle-dynamics model for granular materials

    SciTech Connect

    Walton, O.R.

    1982-05-01

    Discrete-particle simulation of granular-material motion is developing into a viable method for studying how various interparticulate forces affect the bulk behavior of granular solids. A two-dimensional, polygonal-particle computer model, developed from the ideas of Cundall (1976), and incorporating other techniques from molecular dynamics, is being used in a study of the flow behavior of rubblized oil shale. Direct comparison with physical tests involving multiblock systems have verified the model's ability to predict the motion of real materials. Computer generated movies and high-speed motion pictures of physical tests involving gravity flow of 2-dimensional polygonal particles show formation of temporary arches followed by dynamic rupture and reformation of new arches. Direct shear tests on oil-shale rubble involving very large displacements indicate significant circulatory motion in the rubble. Computer simulation of the direct shear tests show similar behavior.

  18. Quasilinear Line Broadened Model for Energetic Particle Transport

    NASA Astrophysics Data System (ADS)

    Ghantous, Katy; Gorelenkov, Nikolai; Berk, Herbert

    2011-10-01

    We present a self-consistent quasi-linear model that describes wave-particle interaction in toroidal geometry and computes fast ion transport during TAE mode evolution. The model bridges the gap between single mode resonances, where it predicts the analytically expected saturation levels, and the case of multiple modes overlapping, where particles diffuse across phase space. Results are presented in the large aspect ratio limit where analytic expressions are used for Fourier harmonics of the power exchange between waves and particles, . Implemention of a more realistic mode structure calculated by NOVAK code are also presented. This work is funded by DOE contract DE-AC02-09CH11466.

  19. A description of rotations for DEM models of particle systems

    NASA Astrophysics Data System (ADS)

    Campello, Eduardo M. B.

    2015-06-01

    In this work, we show how a vector parameterization of rotations can be adopted to describe the rotational motion of particles within the framework of the discrete element method (DEM). It is based on the use of a special rotation vector, called Rodrigues rotation vector, and accounts for finite rotations in a fully exact manner. The use of fictitious entities such as quaternions or complicated structures such as Euler angles is thereby circumvented. As an additional advantage, stick-slip friction models with inter-particle rolling motion are made possible in a consistent and elegant way. A few examples are provided to illustrate the applicability of the scheme. We believe that simple vector descriptions of rotations are very useful for DEM models of particle systems.

  20. Modeling of long range frequency sweeping for energetic particle modes

    SciTech Connect

    Nyqvist, R. M.; Breizman, B. N.

    2013-04-15

    Long range frequency sweeping events are simulated numerically within a one-dimensional, electrostatic bump-on-tail model with fast particle sources and collisions. The numerical solution accounts for fast particle trapping and detrapping in an evolving wave field with a fixed wavelength, and it includes three distinct collisions operators: Drag (dynamical friction on the background electrons), Krook-type collisions, and velocity space diffusion. The effects of particle trapping and diffusion on the evolution of holes and clumps are investigated, and the occurrence of non-monotonic (hooked) frequency sweeping and asymptotically steady holes is discussed. The presented solution constitutes a step towards predictive modeling of frequency sweeping events in more realistic geometries.

  1. Advances in Bayesian Model Based Clustering Using Particle Learning

    SciTech Connect

    Merl, D M

    2009-11-19

    Recent work by Carvalho, Johannes, Lopes and Polson and Carvalho, Lopes, Polson and Taddy introduced a sequential Monte Carlo (SMC) alternative to traditional iterative Monte Carlo strategies (e.g. MCMC and EM) for Bayesian inference for a large class of dynamic models. The basis of SMC techniques involves representing the underlying inference problem as one of state space estimation, thus giving way to inference via particle filtering. The key insight of Carvalho et al was to construct the sequence of filtering distributions so as to make use of the posterior predictive distribution of the observable, a distribution usually only accessible in certain Bayesian settings. Access to this distribution allows a reversal of the usual propagate and resample steps characteristic of many SMC methods, thereby alleviating to a large extent many problems associated with particle degeneration. Furthermore, Carvalho et al point out that for many conjugate models the posterior distribution of the static variables can be parametrized in terms of [recursively defined] sufficient statistics of the previously observed data. For models where such sufficient statistics exist, particle learning as it is being called, is especially well suited for the analysis of streaming data do to the relative invariance of its algorithmic complexity with the number of data observations. Through a particle learning approach, a statistical model can be fit to data as the data is arriving, allowing at any instant during the observation process direct quantification of uncertainty surrounding underlying model parameters. Here we describe the use of a particle learning approach for fitting a standard Bayesian semiparametric mixture model as described in Carvalho, Lopes, Polson and Taddy. In Section 2 we briefly review the previously presented particle learning algorithm for the case of a Dirichlet process mixture of multivariate normals. In Section 3 we describe several novel extensions to the original

  2. Applying dispersive changes to Lagrangian particles in groundwater transport models

    USGS Publications Warehouse

    Konikow, Leonard F.

    2010-01-01

    Method-of-characteristics groundwater transport models require that changes in concentrations computed within an Eulerian framework to account for dispersion be transferred to moving particles used to simulate advective transport. A new algorithm was developed to accomplish this transfer between nodal values and advecting particles more precisely and realistically compared to currently used methods. The new method scales the changes and adjustments of particle concentrations relative to limiting bounds of concentration values determined from the population of adjacent nodal values. The method precludes unrealistic undershoot or overshoot for concentrations of individual particles. In the new method, if dispersion causes cell concentrations to decrease during a time step, those particles in the cell having the highest concentration will decrease the most, and those with the lowest concentration will decrease the least. The converse is true if dispersion is causing concentrations to increase. Furthermore, if the initial concentration on a particle is outside the range of the adjacent nodal values, it will automatically be adjusted in the direction of the acceptable range of values. The new method is inherently mass conservative.

  3. Applying Dispersive Changes to Lagrangian Particles in Groundwater Transport Models

    USGS Publications Warehouse

    Konikow, L.F.

    2010-01-01

    Method-of-characteristics groundwater transport models require that changes in concentrations computed within an Eulerian framework to account for dispersion be transferred to moving particles used to simulate advective transport. A new algorithm was developed to accomplish this transfer between nodal values and advecting particles more precisely and realistically compared to currently used methods. The new method scales the changes and adjustments of particle concentrations relative to limiting bounds of concentration values determined from the population of adjacent nodal values. The method precludes unrealistic undershoot or overshoot for concentrations of individual particles. In the new method, if dispersion causes cell concentrations to decrease during a time step, those particles in the cell having the highest concentration will decrease the most, and those with the lowest concentration will decrease the least. The converse is true if dispersion is causing concentrations to increase. Furthermore, if the initial concentration on a particle is outside the range of the adjacent nodal values, it will automatically be adjusted in the direction of the acceptable range of values. The new method is inherently mass conservative. ?? US Government 2010.

  4. MODELING OF PARTICLE FORMATION AND DYNAMICS IN A FLAME INCINERATOR

    EPA Science Inventory

    A model has been developed to predict the formation and growth of metallic particles in a flame incinerator system. Flow fields and temperature profiles in a cylindrical laminar jet flame have been used to determine the position and physical conditions of the species along the fl...

  5. A note on the multispecies model for identical particles

    SciTech Connect

    Yu, M. Y.; Luo Huaqiang

    2008-02-15

    It is pointed out that the multispecies model for identical particles is in general not suitable for considering quasistationary nonlinear plasma structures. For phenomena on long time scales, the a priori partitioning of electrons into hot and cold species is often unrealistic unless the latter are actually separated in the phase space by external or self-consistent fields.

  6. Thermodynamic model for bouncing charged particles inside a capacitor

    NASA Astrophysics Data System (ADS)

    Rezaeizadeh, Amin; Mameghani, Pooya

    2013-08-01

    We introduce an equation of state for a conducting particle inside a charged parallel-plate capacitor and show that it is similar to the equation of state for an ideal gas undergoing an adiabatic process. We describe a simple experiment that shows reasonable agreement with the theoretical model.

  7. Neural Networks for Modeling and Control of Particle Accelerators

    DOE PAGESBeta

    Edelen, A. L.; Biedron, S. G.; Chase, B. E.; Edstrom, D.; Milton, S. V.; Stabile, P.

    2016-04-01

    Myriad nonlinear and complex physical phenomena are host to particle accelerators. They often involve a multitude of interacting systems, are subject to tight performance demands, and should be able to run for extended periods of time with minimal interruptions. Often times, traditional control techniques cannot fully meet these requirements. One promising avenue is to introduce machine learning and sophisticated control techniques inspired by artificial intelligence, particularly in light of recent theoretical and practical advances in these fields. Within machine learning and artificial intelligence, neural networks are particularly well-suited to modeling, control, and diagnostic analysis of complex, nonlinear, and time-varying systems,more » as well as systems with large parameter spaces. Consequently, the use of neural network-based modeling and control techniques could be of significant benefit to particle accelerators. For the same reasons, particle accelerators are also ideal test-beds for these techniques. Moreover, many early attempts to apply neural networks to particle accelerators yielded mixed results due to the relative immaturity of the technology for such tasks. For the purpose of this paper is to re-introduce neural networks to the particle accelerator community and report on some work in neural network control that is being conducted as part of a dedicated collaboration between Fermilab and Colorado State University (CSU). We also describe some of the challenges of particle accelerator control, highlight recent advances in neural network techniques, discuss some promising avenues for incorporating neural networks into particle accelerator control systems, and describe a neural network-based control system that is being developed for resonance control of an RF electron gun at the Fermilab Accelerator Science and Technology (FAST) facility, including initial experimental results from a benchmark controller.« less

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

  9. Micromechanical Modeling of Storage Particles in Lithium Ion Batteries

    NASA Astrophysics Data System (ADS)

    Purkayastha, Rajlakshmi Tarun

    The effect of stress on storage particles within a lithium ion battery, while acknowledged, is not understood very well. In this work three non-dimensional parameters were identified which govern the stress response within a spherical storage particle. These parameters are developed using material properties such as the diffusion coefficient, particle radius, partial molar volume and Young's modulus. Stress maps are then generated for various values of these parameters for fixed rates of insertion, applying boundary conditions similar to those found in a battery. Stress and concentration profiles for various values of these parameters show the coupling between stress and concentration is magnified depending on the values of the parameters. These maps can be used for different materials, depending on the value of the dimensionless parameters. The value of maximum stress generated is calculated for extraction as well as insertion of lithium into the particle. The model was then used to study to ellipsoidal particles in order to ascertain the effect of geometry on the maximum stress within the particle. By performing a parameter study, we can identify those materials for which particular aspect ratios of ellipsoids are more beneficial, in terms of reducing stress. We find that the stress peaks at certain aspect ratios, mostly at 2 and 1/ 2 . A parameter study was also performed on cubic particle. The values of maximum stresses for both insertion and extraction of lithium were plotted as contour plots. It was seen that the material parameters influenced the location of the maximum stress, with the maximum stress occurring either at the center of the edge between two faces or the point at the center of a face. Newer materials such as silicon are being touted as new lithium storage materials for batteries due to their higher capacity. Their tendency to rapidly loose capacity in a short period of time has led to a variety designs such are the use of carbon nanotubes or

  10. Constants of motion of the four-particle Calogero model

    SciTech Connect

    Saghatelian, A.

    2012-10-15

    We present the explicit expressions of the complete set of constants of motion of four-particle Calogero model with excluded center of mass, i.e. of the A{sub 3} rational Calogero model. Then we find the constants of motion of its spherical part, defining two-dimensional 12-center spherical oscillator, with the force centers located at the vertexes of cuboctahedron.

  11. Modelling of critical phenomena for ignition of metal particles

    NASA Astrophysics Data System (ADS)

    Shchepakina, E.; Sobolev, V.

    2008-11-01

    The singularly perturbed system of differential equations describing metal particle ignition is analyzed. Two general types of combustion reactions are recognized, namely subcritical and supercritical. It is shown that there exists a so-called critical regime, which separates domains of subcritical and supercritical regimes. The critical regime is modelled by the special phase trajectory of the system, which includes the repulsive part with slow variation of the temperature of a metal particle and the thickness of the oxide film. The asymptotic formulae for the calculation of critical values of the modified Semenov number are obtained in the cases of parabolic and cubic laws of oxidation kinetics.

  12. Discrete element method modeling of the triboelectric charging of polyethylene particles: Can particle size distribution and segregation reduce the charging?

    NASA Astrophysics Data System (ADS)

    Konopka, Ladislav; Kosek, Juraj

    2015-10-01

    Polyethylene particles of various sizes are present in industrial gas-dispersion reactors and downstream processing units. The contact of the particles with a device wall as well as the mutual particle collisions cause electrons on the particle surface to redistribute in the system. The undesirable triboelectric charging results in several operational problems and safety risks in industrial systems, for example in the fluidized-bed polymerization reactor. We studied the charging of polyethylene particles caused by the particle-particle interactions in gas. Our model employs the Discrete Element Method (DEM) describing the particle dynamics and incorporates the ‘Trapped Electron Approach’ as the physical basis for the considered charging mechanism. The model predicts the particle charge distribution for systems with various particle size distributions and various level of segregation. Simulation results are in a qualitative agreement with experimental observations of similar particulate systems specifically in two aspects: 1) Big particles tend to gain positive charge and small particles the negative one. 2) The wider the particle size distribution is, the more pronounced is the charging process. Our results suggest that not only the size distribution, but also the effect of the spatial segregation of the polyethylene particles significantly influence the resulting charge distribution ‘generated’ in the system. The level of particle segregation as well as the particle size distribution of polyethylene particles can be in practice adjusted by the choice of supported catalysts, by the conditions in the fluidized-bed polymerization reactor and by the fluid dynamics. We also attempt to predict how the reactor temperature affects the triboelectric charging of particles.

  13. The 5th Generation model of Particle Physics

    NASA Astrophysics Data System (ADS)

    Lach, Theodore

    2009-05-01

    The Standard model of Particle Physics is able to account for all known HEP phenomenon, yet it is not able to predict the masses of the quarks or leptons nor can it explain why they have their respective values. The Checker Board Model (CBM) predicts that there are 5 generation of quarks and leptons and shows a pattern to those masses, namely each three quarks or leptons (within adjacent generations or within a generation) are related to each other by a geometric mean relationship. A 2D structure of the nucleus can be imaged as 2D plate spinning on its axis, it would for all practical circumstances appear to be a 3D object. The masses of the hypothesized ``up'' and ``dn'' quarks determined by the CBM are 237.31 MeV and 42.392 MeV respectively. These new quarks in addition to a lepton of 7.4 MeV make up one of the missing generations. The details of this new particle physics model can be found at the web site: checkerboard.dnsalias.net. The only areas were this theory conflicts with existing dogma is in the value of the mass of the Top quark. The particle found at Fermi Lab must be some sort of composite particle containing Top quarks.

  14. Cirrus cloud model parameterizations: Incorporating realistic ice particle generation

    NASA Technical Reports Server (NTRS)

    Sassen, Kenneth; Dodd, G. C.; Starr, David OC.

    1990-01-01

    Recent cirrus cloud modeling studies have involved the application of a time-dependent, two dimensional Eulerian model, with generalized cloud microphysical parameterizations drawn from experimental findings. For computing the ice versus vapor phase changes, the ice mass content is linked to the maintenance of a relative humidity with respect to ice (RHI) of 105 percent; ice growth occurs both with regard to the introduction of new particles and the growth of existing particles. In a simplified cloud model designed to investigate the basic role of various physical processes in the growth and maintenance of cirrus clouds, these parametric relations are justifiable. In comparison, the one dimensional cloud microphysical model recently applied to evaluating the nucleation and growth of ice crystals in cirrus clouds explicitly treated populations of haze and cloud droplets, and ice crystals. Although these two modeling approaches are clearly incompatible, the goal of the present numerical study is to develop a parametric treatment of new ice particle generation, on the basis of detailed microphysical model findings, for incorporation into improved cirrus growth models. For example, the relation between temperature and the relative humidity required to generate ice crystals from ammonium sulfate haze droplets, whose probability of freezing through the homogeneous nucleation mode are a combined function of time and droplet molality, volume, and temperature. As an example of this approach, the results of cloud microphysical simulations are presented showing the rather narrow domain in the temperature/humidity field where new ice crystals can be generated. The microphysical simulations point out the need for detailed CCN studies at cirrus altitudes and haze droplet measurements within cirrus clouds, but also suggest that a relatively simple treatment of ice particle generation, which includes cloud chemistry, can be incorporated into cirrus cloud growth.

  15. Determining Trajectory of Triboelectrically Charged Particles, Using Discrete Element Modeling

    NASA Technical Reports Server (NTRS)

    2008-01-01

    The Kennedy Space Center (KSC) Electrostatics and Surface Physics Laboratory is participating in an Innovative Partnership Program (IPP) project with an industry partner to modify a commercial off-the-shelf simulation software product to treat the electrodynamics of particulate systems. Discrete element modeling (DEM) is a numerical technique that can track the dynamics of particle systems. This technique, which was introduced in 1979 for analysis of rock mechanics, was recently refined to include the contact force interaction of particles with arbitrary surfaces and moving machinery. In our work, we endeavor to incorporate electrostatic forces into the DEM calculations to enhance the fidelity of the software and its applicability to (1) particle processes, such as electrophotography, that are greatly affected by electrostatic forces, (2) grain and dust transport, and (3) the study of lunar and Martian regoliths.

  16. Particle systems for adaptive, isotropic meshing of CAD models

    PubMed Central

    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

  17. Noise, Bifurcations, and Modeling of Interacting Particle Systems

    PubMed Central

    Mier-y-Teran-Romero, Luis; Forgoston, Eric; Schwartz, Ira B.

    2011-01-01

    We consider the stochastic patterns of a system of communicating, or coupled, self-propelled particles in the presence of noise and communication time delay. For sufficiently large environmental noise, there exists a transition between a translating state and a rotating state with stationary center of mass. Time delayed communication creates a bifurcation pattern dependent on the coupling amplitude between particles. Using a mean field model in the large number limit, we show how the complete bifurcation unfolds in the presence of communication delay and coupling amplitude. Relative to the center of mass, the patterns can then be described as transitions between translation, rotation about a stationary point, or a rotating swarm, where the center of mass undergoes a Hopf bifurcation from steady state to a limit cycle. Examples of some of the stochastic patterns will be given for large numbers of particles. PMID:22124204

  18. Simulation of Cell Adhesion using a Particle Transport Model

    NASA Astrophysics Data System (ADS)

    Chesnutt, Jennifer

    2005-11-01

    An efficient computational method for simulation of cell adhesion through protein binding forces is discussed. In this method, the cells are represented by deformable elastic particles, and the protein binding is represented by a rate equation. The method is first developed for collision and adhesion of two similar cells impacting on each other from opposite directions. The computational method is then applied in a particle-transport model for a cloud of interacting and colliding cells, each of which are represented by particles of finite size. One application might include red blood cells adhering together to form rouleaux, which are chains of red blood cells that are found in different parts of the circulatory system. Other potential applications include adhesion of platelets to a blood vessel wall or mechanical heart valve, which is a precursor of thrombosis formation, or adhesion of cancer cells to organ walls in the lymphatic, circulatory, digestive or pulmonary systems.

  19. An asymptotic model of particle deposition at an airway bifurcation

    PubMed Central

    Zierenberg, Jennifer R.; Halpern, David; Filoche, Marcel; Sapoval, Bernard; Grotberg, James B.

    2013-01-01

    Particle transport and deposition associated with flow over a wedge is investigated as a model for particle transport and flow at the carina of an airway bifurcation during inspiration. Using matched asymptotics, a uniformly valid solution is obtained to represent the high Reynolds number flow over a wedge that considers the viscous boundary layer near the wedge and the outer inviscid region and is then used to solve the particle transport equations. Sometimes particle impaction on the wedge is prevented due to the boundary layer. We call this boundary layer shielding (BLS). This effect can be broken down into different types: rejection, trapping and deflection that are described by what happens to the particle’s initial negative velocity normal to the wall either changing sign, reaching zero, or remaining negative in the boundary layer region. The deposition efficiency depends on the critical Stokes number but exhibits a weak dependence on Reynolds number. Deposition efficiency for Sc in the range 0 < Sc < 0.4 yields the following relationship De ≈ (1.867 Sc1.78− 0.016) sin(βπ/2) at large Reynolds numbers, where βπ is the wedge angle. For a specific deposition efficiency, Sc decreases as βπ increases. The distribution of impacted particles was also computed and revealed that particles primarily impact within one airway diameter of the carina, consistent with computational fluid dynamics approaches. This work provides a new insight that the BLS inherent to the wedge component of the structure is the dominant reason for the particle distribution. This finding is important in linking aerosol deposition to the location of airway disease as well as target sites for therapeutic deposition. PMID:22378463

  20. Physics of Hall-effect thruster by particle model

    SciTech Connect

    Taccogna, Francesco; Minelli, Pierpaolo; Capitelli, Mario; Longo, Savino

    2012-11-27

    A realistic three-dimensional fully kinetic particle simulation of a Hall-effect thruster discharge has been developed. The model consists of a Particle-in-Cell methodology tracking electrons, Xe{sup +} and Xe{sup ++} ions in their selfconsistent electric field. A detailed secondary electron emission from lateral walls is also implemented in addition with electron-atom and electron-ion volume collisions. The model is able to capture the most relevant features of axial, radial and azimuthal behaviors of the start-up transient and steady state phases detecting inverted sheaths and azimuthal modulation in the acceleration region. The model has the potentiality to investigate the driving mechanisms at the origin of the electron anomalous cross-field transport.

  1. Particle Swarm Based Collective Searching Model for Adaptive Environment

    SciTech Connect

    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.

  2. Particle Swarm Based Collective Searching Model for Adaptive Environment

    SciTech Connect

    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.

  3. Nonlinear EEG Decoding Based on a Particle Filter Model

    PubMed Central

    Hong, Jun

    2014-01-01

    While the world is stepping into the aging society, rehabilitation robots play a more and more important role in terms of both rehabilitation treatment and nursing of the patients with neurological diseases. Benefiting from the abundant contents of movement information, electroencephalography (EEG) has become a promising information source for rehabilitation robots control. Although the multiple linear regression model was used as the decoding model of EEG signals in some researches, it has been considered that it cannot reflect the nonlinear components of EEG signals. In order to overcome this shortcoming, we propose a nonlinear decoding model, the particle filter model. Two- and three-dimensional decoding experiments were performed to test the validity of this model. In decoding accuracy, the results are comparable to those of the multiple linear regression model and previous EEG studies. In addition, the particle filter model uses less training data and more frequency information than the multiple linear regression model, which shows the potential of nonlinear decoding models. Overall, the findings hold promise for the furtherance of EEG-based rehabilitation robots. PMID:24949420

  4. Particle-hole duality, integrability, and Russian doll BCS model

    NASA Astrophysics Data System (ADS)

    Bork, L. V.; Pogosov, W. V.

    2015-08-01

    We address a generalized Richardson model (Russian doll BCS model), which is characterized by the breaking of time-reversal symmetry. This model is known to be exactly solvable and integrable. We point out that the Russian doll BCS model, on the level of Hamiltonian, is also particle-hole symmetric. This implies that the same state can be expressed both in the particle and hole representations with two different sets of Bethe roots. We then derive exact relations between Bethe roots in the two representations, which can hardly be obtained staying on the level of Bethe equations. In a quasi-classical limit, similar identities for usual Richardson model, known from literature, are recovered from our results. We also show that these relations for Richardson roots take a remarkably simple form at half-filling and for a symmetric with respect to the middle of the interaction band distribution of one-body energy levels, since, in this special case, the rapidities in the particle and hole representations up to the translation satisfy the same system of equations.

  5. Theory and modeling of particles with DNA-mediated interactions

    NASA Astrophysics Data System (ADS)

    Licata, Nicholas A.

    2008-05-01

    In recent years significant attention has been attracted to proposals which utilize DNA for nanotechnological applications. Potential applications of these ideas range from the programmable self-assembly of colloidal crystals, to biosensors and nanoparticle based drug delivery platforms. In Chapter I we introduce the system, which generically consists of colloidal particles functionalized with specially designed DNA markers. The sequence of bases on the DNA markers determines the particle type. Due to the hybridization between complementary single-stranded DNA, specific, type-dependent interactions can be introduced between particles by choosing the appropriate DNA marker sequences. In Chapter II we develop a statistical mechanical description of the aggregation and melting behavior of particles with DNA-mediated interactions. In Chapter III a model is proposed to describe the dynamical departure and diffusion of particles which form reversible key-lock connections. In Chapter IV we propose a method to self-assemble nanoparticle clusters using DNA scaffolds. A natural extension is discussed in Chapter V, the programmable self-assembly of nanoparticle clusters where the desired cluster geometry is encoded using DNA-mediated interactions. In Chapter VI we consider a nanoparticle based drug delivery platform for targeted, cell specific chemotherapy. In Chapter VII we present prospects for future research: the connection between DNA-mediated colloidal crystallization and jamming, and the inverse problem in self-assembly.

  6. Modelling New Particle Formation Events in the South African Savannah

    SciTech Connect

    Gierens, Rosa; Laakso, Lauri; Mogensen, Ditte; Vakkari, Ville; Beukes, J. P.; Van Zyl, Pieter; Hakola, H.; Guenther, Alex B.; Pienaar, J. J.; Boy, Michael

    2014-01-01

    Africa is one of the less studied continents with respect to atmospheric aerosols. Savannahs are complex dynamic systems sensitive to climate and land-use changes, but the interaction with the atmosphere is not well understood. Atmospheric particles, aka aerosols, affect the climate on regional and global scale, and are an important factor in air quality. In this study measurements from a relatively clean savannah environment in South Africa were used to model new particle formation and growth. There are already some combined long-term measurements of trace gas concentrations together with aerosol and meteorological variables available, but to our knowledge this is the first time detailed simulations, that include all the main processes relevant to particle formation, were done. The results show that both investigated particle formation mechanisms overestimated the formation rates dependency on sulphuric acid. The approach including low volatile organic compounds to the particle formation process was more accurate in describing the nucleation events. To get reliable estimation of aerosol concentration in simulations for larger scales, nucleation mechanisms would need to include organic compounds, at least in southern Africa.

  7. Particle-based ablation model for faint meteors

    NASA Astrophysics Data System (ADS)

    Stokan, E.; Campbell-Brown, M.

    2014-07-01

    Modeling the ablation of meteoroids as they enter the atmosphere is a way of determining their physical structure and elemental composition. This can provide insight into the structure of parent bodies when combined with an orbit computed from observations. The Canadian Automated Meteor Observatory (CAMO) is a source of new, high-resolution observations of faint meteors [1]. These faint objects tend to have pre-atmospheric masses around 10^{-5} kg, corresponding to a radius of 1 mm. A wide-field camera with a 28° field of view provides guidance to a high-resolution camera that tracks meteors in flight with 1.5° field of view. Meteors are recorded with a scale of 4 m per pixel at a range of 135 km, at 110 frames per second, allowing us to investigate detailed meteor morphology. This serves as an important new constraint for ablation models, in addition to meteor brightness (lightcurves) and meteoroid deceleration. High-resolution observations of faint meteors have revealed that contemporary ablation models are not able to predict meteor morphology, even while matching the observed lightcurve and meteoroid deceleration [2]. This implies that other physical processes, in addition to fragmentation, must be considered for faint meteor ablation. We present a new, particle-based approach to modeling the ablation of small meteoroids. In this model, we simulate the collisions between atmospheric particles and the meteoroid to determine the rate of evaporation and deceleration. Subsequent collisions simulated between evaporated meteoroid particles and ambient atmospheric particles then produce light that would be observed by high-resolution cameras. Preliminary results show simultaneous agreement with meteor morphology, lightcurves, and decelerations recorded with CAMO. A sample comparison of simulated and observed meteor morphology is given in the attached figure. Several meteoroids are well-represented as solid, stony bodies, but some require modeling as a dustball [3

  8. A theoretical model for the Lorentz force particle analyzer

    NASA Astrophysics Data System (ADS)

    Moreau, René; Tao, Zhen; Wang, Xiaodong

    2016-07-01

    In a previous paper [X. Wang et al., J. Appl. Phys. 120, 014903 (2016)], several experimental devices have been presented, which demonstrate the efficiency of electromagnetic techniques for detecting and sizing electrically insulating particles entrained in the flow of a molten metal. In each case, a non-uniform magnetic field is applied across the flow of the electrically conducting liquid, thereby generating a braking Lorentz force on this moving medium and a reaction force on the magnet, which tends to be entrained in the flow direction. The purpose of this letter is to derive scaling laws for this Lorentz force from an elementary theoretical model. For simplicity, as in the experiments, the flowing liquid is modeled as a solid body moving with a uniform velocity U. The eddy currents in the moving domain are derived from the classic induction equation and Ohm's law, and expressions for the Lorentz force density j ×B and for its integral over the entire moving domain follow. The insulating particles that are eventually present and entrained with this body are then treated as small disturbances in a classic perturbation analysis, thereby leading to scaling laws for the pulses they generate in the Lorentz force. The purpose of this letter is both to illustrate the eddy currents without and with insulating particles in the electrically conducting liquid and to derive a key relation between the pulses in the Lorentz force and the main parameters (particle volume and dimensions of the region subjected to the magnetic field).

  9. Tabletop Traffic Jams: Modeling Traffic Jams using Self Propelled Particles

    NASA Astrophysics Data System (ADS)

    Yadav, Vikrant; Kudrolli, Arshad

    2015-03-01

    We model behavior of traffic using Self Propelled Particles (SPPs). Granular rods with asymmetric mass distribution confined to move in a circular channel on a vibrated substrate and interact with each other through inelastic collision serve as our model vehicle. Motion of a single vehicle is observed to be composed of 2 parts, a linear velocity in the direction of lighter end of particle and a non-Gaussian random velocity. We find that the collective mean speed of the SPPs is constant over a wide range of line densities before decreasing rapidly as the maximum packing is approached indicating the spontaneous formation of Phantom jams. This decrease in speed is observed to be far greater than any small differences in the mean drift speed of individual SPPs , and occurs as the collision frequency between SPPs increase exponentially with line density. However the random velocity component of SPPs remain super-diffusive over entire range of line densities. While the collective motion at low densities is characterized by caravan following behind the slowest particle leading to clustering, at higher densities we see formation of jamming waves travelling in direction opposite to that of motion of particles.

  10. Dissipative-particle-dynamics model of biofilm growth

    SciTech Connect

    Xu, Zhijie; Meakin, Paul; Tartakovsky, Alexandre M.; Scheibe, Timothy D.

    2011-06-13

    A dissipative particle dynamics (DPD) model for the quantitative simulation of biofilm growth controlled by substrate (nutrient) consumption, advective and diffusive substrate transport, and hydrodynamic interactions with fluid flow (including fragmentation and reattachment) is described. The model was used to simulate biomass growth, decay, and spreading. It predicts how the biofilm morphology depends on flow conditions, biofilm growth kinetics, the rheomechanical properties of the biofilm and adhesion to solid surfaces. The morphology of the model biofilm depends strongly on its rigidity and the magnitude of the body force that drives the fluid over the biofilm.

  11. Modeling convection onset in colloidal suspensions of particles

    NASA Astrophysics Data System (ADS)

    Hadji, Layachi

    2011-09-01

    A particulate medium model is used to investigate the onset of Rayleigh-Bénard convection in a colloidal suspension of inert solid particles. The model accounts for the effects of thermophoresis, sedimentation and Brownian diffusion. Depending on the size of the particles, the problem has up to four time scales. These are due to thermal diffusion, particle diffusion, particle migration due to thermophoresis and particle sedimentation. The ratios of these time scales lead to the emergence of three parameters, one of which is the Lewis number τ. The smallness of the latter makes the differential eigenvalue system governing convection onset singular. The other two are the density number Γ and the dimensionless migration velocity β. For a given experimental set-up, β has a quadratic functional dependence on the particles radius. A combination of asymptotics and numerical computation is used to capture the effect of the resulting thin particle concentration boundary layers on the leading order instability thresholds. Results, which are depicted as function of Γ and β, reveal a non-monotonic dependence of the critical Rayleigh number ℛ c on β. The curve ℛ c (β) is bimodal and it exhibits a maximum , the value of which increases very sharply with |Γ| while the critical wavelength decreases. Experimental conditions can be fixed so that occurs at β values that correspond to nanoparticles. Therefore, experimental parameters can be controlled so that the mixing of a small amount of nano-size particles has a substantial stabilizing effect. An investigation of the flow patterns in the range of parameters associated with this sharp increase in ℛ c reveals a decrease in the effective

  12. Particle dispersion models and drag coefficients for particles in turbulent flows

    NASA Technical Reports Server (NTRS)

    Crowe, C. T.; Chung, J. N.; Troutt, T. R.

    1988-01-01

    Some of the concepts underlying particle dispersion due to turbulence are reviewed. The traditional approaches to particle dispersion in homogeneous, stationary turbulent fields are addressed, and recent work on particle dispersion in large scale turbulent structures is reviewed. The state of knowledge of particle drag coefficients in turbulent gas-particle flows is also reviewed.

  13. 3D flare particle model for ShipIR/NTCS

    NASA Astrophysics Data System (ADS)

    Ramaswamy, Srinivasan; Vaitekunas, David A.

    2016-05-01

    A key component in any soft-kill response to an incoming guided missile is the flare /chaff decoy used to distract or seduce the seeker homing system away from the naval platform. This paper describes a new 3D flare particle model in the naval threat countermeasure simulator (NTCS) of the NATO-standard ship signature model (ShipIR), which provides independent control over the size and radial distribution of its signature. The 3D particles of each flare sub-munition are modelled stochastically and rendered using OpenGL z-buffering, 2D projection, and alpha-blending to produce a unique and time varying signature. A sensitivity analysis on each input parameter provides the data and methods needed to synthesize a model from an IR measurement of a decoy. The new model also eliminated artifacts and deficiencies in our previous model which prevented reliable tracks from the adaptive track gate algorithm already presented by Ramaswamy and Vaitekunas (2015). A sequence of scenarios are used to test and demonstrate the new flare model during a missile engagement.

  14. Modelling the dispersion of particle numbers in five European cities

    NASA Astrophysics Data System (ADS)

    Kukkonen, J.; Karl, M.; Keuken, M. P.; Denier van der Gon, H. A. C.; Denby, B. R.; Singh, V.; Douros, J.; Manders, A.; Samaras, Z.; Moussiopoulos, N.; Jonkers, S.; Aarnio, M.; Karppinen, A.; Kangas, L.; Lützenkirchen, S.; Petäjä, T.; Vouitsis, I.; Sokhi, R. S.

    2016-02-01

    We present an overview of the modelling of particle number concentrations (PNCs) in five major European cities, namely Helsinki, Oslo, London, Rotterdam, and Athens, in 2008. Novel emission inventories of particle numbers have been compiled both on urban and European scales. We used atmospheric dispersion modelling for PNCs in the five target cities and on a European scale, and evaluated the predicted results against available measured concentrations. In all the target cities, the concentrations of particle numbers (PNs) were mostly influenced by the emissions originating from local vehicular traffic. The influence of shipping and harbours was also significant for Helsinki, Oslo, Rotterdam, and Athens, but not for London. The influence of the aviation emissions in Athens was also notable. The regional background concentrations were clearly lower than the contributions originating from urban sources in Helsinki, Oslo, and Athens. The regional background was also lower than urban contributions in traffic environments in London, but higher or approximately equal to urban contributions in Rotterdam. It was numerically evaluated that the influence of coagulation and dry deposition on the predicted PNCs was substantial for the urban background in Oslo. The predicted and measured annual average PNCs in four cities agreed within approximately ≤ 26 % (measured as fractional biases), except for one traffic station in London. This study indicates that it is feasible to model PNCs in major cities within a reasonable accuracy, although major challenges remain in the evaluation of both the emissions and atmospheric transformation of PNCs.

  15. Modeling of Solar Energetic Particles in Interplanetary Space

    NASA Astrophysics Data System (ADS)

    Vainio, Rami; Agueda, Neus; Aran, Angels; Lario, David

    Solar energetic particles (SEPs) in the interplanetary (IP) medium are transported under the influence of electromagnetic fields of the solar wind. These fields consists of the smooth background fields, which can be modeled by the MHD equations governing the expansion of the solar wind, and of the small-scale fluctuations (waves or turbulence) that scatter the particles in pitch angle and act as agents enabling their acceleration at IP shock waves. We review theoretical models of SEP transport and acceleration in the IP medium. We start from the simple analytical approaches (diffusion models), which assume quasi-isotropic particle distributions, and then continue to the more accurate numerical approaches based on the focused transport equation, not making this simplifying assumption. A careful analysis of two SEP events, an impulsive and a gradual one, is presented and the spatial scaling of their peak intensities, differential fluences and time-integrated net fluxes is discussed. We conclude that rather simple scaling laws for these quantities can be obtained for impulsive events but no simple scaling laws can be expected to govern the gradual SEP events

  16. Multi-point Observations and Modeling of Particle Injections

    NASA Astrophysics Data System (ADS)

    Henderson, M. G.; Morley, S.; Reeves, G. D.; Larsen, B.; Skoug, R. M.; Funsten, H. O.; Spence, H. E.

    2014-12-01

    Dispersionless and dispersed particle injections associated with substorms have been studied for many years based on observations acquired primarily at geosynchronous orbit. A general picture that has emerged is that particles are energized and rapidly transported/organized behind an "injection boundary" that penetrates closer to Earth in some magnetic local time sector (e.g. the so-called double-spiral injection boundary model). While this picture provides a very good description of injections at geosynchronous orbit, with the recent launchof the Van Allen Probes mission, we are now able to explore the evolution of injection signatures well inside of geosynchronous orbit at multiple locations as well. Here we report initial results of injection boundary modeling based on simultaneous multi-point measurements at both geosynchronous orbit (from the LANL/GEO spacecraft) and inside (from Van Allen Probes spacecraft). It is shown that many of the complex dispersion features observed in Van Allen Probes particle data are reproduced by the injection boundary model.

  17. Object Oriented Design and the Standard Model of particle physics

    NASA Astrophysics Data System (ADS)

    Lipovaca, Samir

    2007-04-01

    Inspired by the computer as both tool and metaphor, a new path emerges toward understanding life, physics, and existence. The path leads throughout all of nature, from the interior of cells to inside black holes. This view of science is based on the idea that information is the ultimate ``substance'' from which all things are made. Exploring this view, we will focus on Object - Oriented (OO) design as one of the most important designs in software development. The OO design views the world as composed of objects with well defined properties. The dynamics is pictured as interactions among objects. Interactions are mediated by messages that objects exchange with each other. This description closely resembles the view of the elementary particles world created by the Standard Model of particle physics. The object model (OM) provides a theoretical foundation upon which the OO design is built. The OM is based on the principles of abstraction, encapsulation, modularity and hierarchy. We will show that the Standard Model of particle physics follows the OM principles.

  18. Modelling of strongly coupled particle growth and aggregation

    NASA Astrophysics Data System (ADS)

    Gruy, F.; Touboul, E.

    2013-02-01

    The mathematical modelling of the dynamics of particle suspension is based on the population balance equation (PBE). PBE is an integro-differential equation for the population density that is a function of time t, space coordinates and internal parameters. Usually, the particle is characterized by a unique parameter, e.g. the matter volume v. PBE consists of several terms: for instance, the growth rate and the aggregation rate. So, the growth rate is a function of v and t. In classical modelling, the growth and the aggregation are independently considered, i.e. they are not coupled. However, current applications occur where the growth and the aggregation are coupled, i.e. the change of the particle volume with time is depending on its initial value v0, that in turn is related to an aggregation event. As a consequence, the dynamics of the suspension does not obey the classical Von Smoluchowski equation. This paper revisits this problem by proposing a new modelling by using a bivariate PBE (with two internal variables: v and v0) and by solving the PBE by means of a numerical method and Monte Carlo simulations. This is applied to a physicochemical system with a simple growth law and a constant aggregation kernel.

  19. Receptor modeling application framework for particle source apportionment.

    PubMed

    Watson, John G; Zhu, Tan; Chow, Judith C; Engelbrecht, Johann; Fujita, Eric M; Wilson, William E

    2002-12-01

    Receptor models infer contributions from particulate matter (PM) source types using multivariate measurements of particle chemical and physical properties. Receptor models complement source models that estimate concentrations from emissions inventories and transport meteorology. Enrichment factor, chemical mass balance, multiple linear regression, eigenvector. edge detection, neural network, aerosol evolution, and aerosol equilibrium models have all been used to solve particulate air quality problems, and more than 500 citations of their theory and application document these uses. While elements, ions, and carbons were often used to apportion TSP, PM10, and PM2.5 among many source types, many of these components have been reduced in source emissions such that more complex measurements of carbon fractions, specific organic compounds, single particle characteristics, and isotopic abundances now need to be measured in source and receptor samples. Compliance monitoring networks are not usually designed to obtain data for the observables, locations, and time periods that allow receptor models to be applied. Measurements from existing networks can be used to form conceptual models that allow the needed monitoring network to be optimized. The framework for using receptor models to solve air quality problems consists of: (1) formulating a conceptual model; (2) identifying potential sources; (3) characterizing source emissions; (4) obtaining and analyzing ambient PM samples for major components and source markers; (5) confirming source types with multivariate receptor models; (6) quantifying source contributions with the chemical mass balance; (7) estimating profile changes and the limiting precursor gases for secondary aerosols; and (8) reconciling receptor modeling results with source models, emissions inventories, and receptor data analyses. PMID:12492167

  20. Particle Tracking Model and Abstraction of Transport Processes

    SciTech Connect

    B. Robinson

    2000-04-07

    The purpose of the transport methodology and component analysis is to provide the numerical methods for simulating radionuclide transport and model setup for transport in the unsaturated zone (UZ) site-scale model. The particle-tracking method of simulating radionuclide transport is incorporated into the FEHM computer code and the resulting changes in the FEHM code are to be submitted to the software configuration management system. This Analysis and Model Report (AMR) outlines the assumptions, design, and testing of a model for calculating radionuclide transport in the unsaturated zone at Yucca Mountain. In addition, methods for determining colloid-facilitated transport parameters are outlined for use in the Total System Performance Assessment (TSPA) analyses. Concurrently, process-level flow model calculations are being carrier out in a PMR for the unsaturated zone. The computer code TOUGH2 is being used to generate three-dimensional, dual-permeability flow fields, that are supplied to the Performance Assessment group for subsequent transport simulations. These flow fields are converted to input files compatible with the FEHM code, which for this application simulates radionuclide transport using the particle-tracking algorithm outlined in this AMR. Therefore, this AMR establishes the numerical method and demonstrates the use of the model, but the specific breakthrough curves presented do not necessarily represent the behavior of the Yucca Mountain unsaturated zone.

  1. Current models of the intensely ionizing particle environment in space

    NASA Technical Reports Server (NTRS)

    Adams, James H., Jr.

    1988-01-01

    The Cosmic Ray Effects on MicroElectronics (CREME) model that is currently in use to estimate single event effect rates in spacecraft is described. The CREME model provides a description of the radiation environment in interplanetary space near the orbit of the earth that contains no major deficiencies. The accuracy of the galactic cosmic ray model is limited by the uncertainties in solar modulation. The model for solar energetic particles could be improved by making use of all the data that has been collected on solar energetic particle events. There remain major uncertainties about the environment within the earth's magnetosphere, because of the uncertainties over the charge states of the heavy ions in the anomalous component and solar flares, and because of trapped heavy ions. The present CREME model is valid only at 1 AU, but it could be extended to other parts of the heliosphere. There is considerable data on the radiation environment from 0.2 to 35 AU in the ecliptic plane. This data could be used to extend the CREME model.

  2. Generation of Random Particle Packings for Discrete Element Models

    NASA Astrophysics Data System (ADS)

    Abe, S.; Weatherley, D.; Ayton, T.

    2012-04-01

    An important step in the setup process of Discrete Element Model (DEM) simulations is the generation of a suitable particle packing. There are quite a number of properties such a granular material specimen should ideally have, such as high coordination number, isotropy, the ability to fill arbitrary bounding volumes and the absence of locked-in stresses. An algorithm which is able to produce specimens fulfilling these requirements is the insertion based sphere packing algorithm originally proposed by Place and Mora, 2001 [2] and extended in this work. The algorithm works in two stages. First a number of "seed" spheres are inserted into the bounding volume. In the second stage the gaps between the "seed" spheres are filled by inserting new spheres in a way so they have D+1 (i.e. 3 in 2D, 4 in 3D) touching contacts with either other spheres or the boundaries of the enclosing volume. Here we present an implementation of the algorithm and a systematic statistical analysis of the generated sphere packings. The analysis of the particle radius distribution shows that they follow a power-law with an exponent ≈ D (i.e. ≈3 for a 3D packing and ≈2 for 2D). Although the algorithm intrinsically guarantees coordination numbers of at least 4 in 3D and 3 in 2D, the coordination numbers realized in the generated packings can be significantly higher, reaching beyond 50 if the range of particle radii is sufficiently large. Even for relatively small ranges of particle sizes (e.g. Rmin = 0.5Rmax) the maximum coordination number may exceed 10. The degree of isotropy of the generated sphere packing is also analysed in both 2D and 3D, by measuring the distribution of orientations of vectors joining the centres of adjacent particles. If the range of particle sizes is small, the packing algorithm yields moderate anisotropy approaching that expected for a face-centred cubic packing of equal-sized particles. However, once Rmin < 0.3Rmax a very high degree of isotropy is demonstrated in

  3. Effect of particle volume fraction on the settling velocity of volcanic ash particles: implications for ash dispersion models

    NASA Astrophysics Data System (ADS)

    Del Bello, E.; Taddeucci, J.; De'Michieli Vitturi, M.; Scarlato, P.; Andronico, D.; Scollo, S.; Kueppers, U.

    2015-12-01

    We present the first report of experimental measurements of the enhanced settling velocity of volcanic particles as function of particle volume fraction. In order to investigate the differences in the aerodynamic behavior of ash particles when settling individually or in mass, we performed systematic large-scale ash settling experiments using natural basaltic and phonolitic ash. By releasing ash particles at different, controlled volumetric flow rates, in an unconstrained open space and at minimal air movement, we measured their terminal velocity, size, and particle volume fraction with a high-speed camera at 2000 fps. Enhanced settling velocities of individual particles increase with increasing particle volume fraction. This suggests that particle clustering during fallout may be one reason explaining larger than theoretical depletion rates of fine particles from volcanic ash clouds. We provide a quantitative empirical model that allows to calculate, from a given particle size and density, the enhanced velocity resulting from a given particle volume fraction. The proposed model has the potential to serve as a simple tool for the prediction of the terminal velocity of ash of an hypothetical distribution of ash of known particle size and volume fraction. This is of particular importance for advection-diffusion transport model of ash where generally a one-way coupling is adopted, considering only the flow effects on particles. To better quantify the importance of the enhanced settling velocity in ash dispersal, we finally introduced the new formulation in a Lagrangian model calculating for realistic eruptive conditions the resulting ash concentration in the atmosphere and on the ground.

  4. Force models for particle-dynamics simulations of granular materials

    SciTech Connect

    Walton, O.R.

    1994-12-01

    Engineering-mechanics contact models are utilized to describe the inelastic, frictional interparticle forces acting in dry granular systems. Simple analyses based on one-dimensional chains are utilized to illustrate wave propagation phenomena in dense and dilute discrete particulates. The variation of restitution coefficient with impact velocity is illustrated for a variety of viscous and hysteretic normal force models. The effects of interparticle friction on material strength in discrete-particle simulations are much closer to measured values than are theories that do not allow article rotations.

  5. Modelling the internal structure of nascent soot particles

    SciTech Connect

    Totton, Tim S.; Sander, Markus; Kraft, Markus; Chakrabarti, Dwaipayan; Wales, David J.; Misquitta, Alston J.

    2010-05-15

    In this paper we present studies of clusters assembled from polycyclic aromatic hydrocarbon (PAH) molecules similar in size to small soot particles. The clusters studied were comprised of coronene (C{sub 24}H{sub 12}) or pyrene (C{sub 16}H{sub 10}) molecules and represent the types of soot precursor molecule typically found in flame environments. A stochastic 'basin-hopping' global optimisation scheme was used to locate low-lying local minima on the potential energy surface of the molecular clusters. TEM-style projections of the resulting geometries show similarities with those observed experimentally in TEM images of soot particles. The mass densities of these clusters have also been calculated and are lower than bulk values of the pure crystalline PAH structures. They are also significantly lower than the standard value of 1.8 g/cm{sup 3} used in our soot models. Consequently we have varied the mass density between 1.0 g/cm{sup 3} and 1.8 g/cm{sup 3} to examine the effects of varying soot density on our soot model and observed how the shape of the particle size distribution changes. Based on similarities between nascent soot particles and PAH clusters a more accurate soot density is likely to be significantly lower than 1.8 g/cm{sup 3}. As such, for modelling purposes, we recommend that the density of nascent soot should be taken to be the value obtained for our coronene cluster of 1.12 g/cm{sup 3}. (author)

  6. Theory and modeling of particles with DNA-mediated interactions

    NASA Astrophysics Data System (ADS)

    Licata, Nicholas A.

    In recent years significant attention has been attracted to proposals which utilize DNA for nanotechnological applications. Potential applications of these ideas range from the programmable self-assembly of colloidal crystals, to biosensors and nanoparticle based drug delivery platforms. In Chapter I we introduce the system, which generically consists of colloidal particles functionalized with specially designed DNA markers. The sequence of bases on the DNA markers determines the particle type. Due to the hybridization between complementary single-stranded DNA, specific, type-dependent interactions can be introduced between particles by choosing the appropriate DNA marker sequences. In Chapter II we develop a statistical mechanical description of the aggregation and melting behavior of particles with DNA-mediated interactions. A quantitative comparison between the theory and experiments is made by calculating the experimentally observed melting profile. In Chapter III a model is proposed to describe the dynamical departure and diffusion of particles which form reversible key-lock connections. The model predicts a crossover from localized to diffusive behavior. The random walk statistics for the particles' in plane diffusion is discussed. The lateral motion is analogous to dispersive transport in disordered semiconductors, ranging from standard diffusion with a renormalized diffusion coefficient to anomalous, subdiffusive behavior. In Chapter IV we propose a method to self-assemble nanoparticle clusters using DNA scaffolds. An optimal concentration ratio is determined for the experimental implementation of our self-assembly proposal. A natural extension is discussed in Chapter V, the programmable self-assembly of nanoparticle clusters where the desired cluster geometry is encoded using DNA-mediated interactions. We determine the probability that the system self-assembles the desired cluster geometry, and discuss the connections to jamming in granular and colloidal

  7. Precision Modeling of Solar Energetic Particle Intensity and Anisotropy Profiles

    NASA Astrophysics Data System (ADS)

    Ruffolo, D.; Sáiz, A.; Bieber, J. W.; Evenson, P.; Pyle, R.; Rujiwarodom, M.; Tooprakai, P.; Wechakama, M.; Khumlumlert, T.

    2006-12-01

    A focused transport equation for solar energetic particles is sufficiently complex that simple analytic approximations are generally inadequate, but the physics is sufficiently well established to permit precise numerical modeling of high energy particle observations at various distances from the Sun. We demonstrate how observed profiles of intensity and anisotropy vs. time can be quantitatively fit to determine an optimal injection profile at the Sun, scattering mean free path λ, and magnetic configuration. For several ground level enhancements (GLE) of solar energetic particles at energies ~ 1 GeV, the start time of injection has been determined to 1 or 2 minutes. In each case this start time coincides, within that precision, to the soft X-ray peak time, when the flare's primary energy release has ended. This is not inconsistent with acceleration at a coronal mass ejection (CME)-driven shock, though the rapid timescale is challenging to understand. For the GLE of 2005 January 20, λ decreases substantially over ~ 10 minutes, which is consistent with concepts of proton-amplified waves. The GLE of 2000 July 14 is properly fit only when a magnetic bottleneck beyond Earth is taken into account, a feature later confirmed by NEAR observations. The long-standing puzzle of the 1989 October 22 event can now be explained by simultaneous injection of relativistic solar particles along both legs of a closed interplanetary magnetic loop, while other reasonable explanations fail the test of quantitative fitting. The unusually long λ (confirming many previous reports) and a low turbulent spectral index hint at unusual properties of turbulence in the loop. While the early GLE peak on 2003 October 28 remains a mystery, the main peak's strong anisotropy is inconsistent with a suggestion of injection along the far leg of a magnetic loop; quantitative fitting fails because of reverse focusing during Sunward motion. With these modeling capabilities, one is poised to take full

  8. Modelling of aircrew radiation exposure during solar particle events

    NASA Astrophysics Data System (ADS)

    Al Anid, Hani Khaled

    show a very different response during anisotropic events, leading to variations in aircrew radiation doses that may be significant for dose assessment. To estimate the additional exposure due to solar flares, a model was developed using a Monte-Carlo radiation transport code, MCNPX. The model transports an extrapolated particle spectrum based on satellite measurements through the atmosphere using the MCNPX analysis. This code produces the estimated flux at a specific altitude where radiation dose conversion coefficients are applied to convert the particle flux into effective and ambient dose-equivalent rates. A cut-off rigidity model accounts for the shielding effects of the Earth's magnetic field. Comparisons were made between the model predictions and actual flight measurements taken with various types of instruments used to measure the mixed radiation field during Ground Level Enhancements 60 and 65. An anisotropy analysis that uses neutron monitor responses and the pitch angle distribution of energetic solar particles was used to identify particle anisotropy for a solar event in December 2006. In anticipation of future commercial use, a computer code has been developed to implement the radiation dose assessment model for routine analysis. Keywords: Radiation Dosimetry, Radiation Protection, Space Physics.

  9. Entrainment of coarse grains using a discrete particle model

    SciTech Connect

    Valyrakis, Manousos; Arnold, Roger B. Jr.

    2014-10-06

    Conventional bedload transport models and incipient motion theories relying on a time-averaged boundary shear stress are incapable of accounting for the effects of fluctuating near-bed velocity in turbulent flow and are therefore prone to significant errors. Impulse, the product of an instantaneous force magnitude and its duration, has been recently proposed as an appropriate criterion for quantifying the effects of flow turbulence in removing coarse grains from the bed surface. Here, a discrete particle model (DPM) is used to examine the effects of impulse, representing a single idealized turbulent event, on particle entrainment. The results are classified according to the degree of grain movement into the following categories: motion prior to entrainment, initial dislodgement, and energetic displacement. The results indicate that in all three cases the degree of particle motion depends on both the force magnitude and the duration of its application and suggest that the effects of turbulence must be adequately accounted for in order to develop a more accurate method of determining incipient motion. DPM is capable of simulating the dynamics of grain entrainment and is an appropriate tool for further study of the fundamental mechanisms of sediment transport.

  10. Comprehensive computer model for magnetron sputtering. II. Charged particle transport

    SciTech Connect

    Jimenez, Francisco J. Dew, Steven K.; Field, David J.

    2014-11-01

    Discharges for magnetron sputter thin film deposition systems involve complex plasmas that are sensitively dependent on magnetic field configuration and strength, working gas species and pressure, chamber geometry, and discharge power. The authors present a numerical formulation for the general solution of these plasmas as a component of a comprehensive simulation capability for planar magnetron sputtering. This is an extensible, fully three-dimensional model supporting realistic magnetic fields and is self-consistently solvable on a desktop computer. The plasma model features a hybrid approach involving a Monte Carlo treatment of energetic electrons and ions, along with a coupled fluid model for thermalized particles. Validation against a well-known one-dimensional system is presented. Various strategies for improving numerical stability are investigated as is the sensitivity of the solution to various model and process parameters. In particular, the effect of magnetic field, argon gas pressure, and discharge power are studied.

  11. A Refined Model for Solid Particle Rock Erosion

    NASA Astrophysics Data System (ADS)

    Momber, A. W.

    2016-02-01

    A procedure for the estimation of distribution parameters of a Weibull distribution model K 1 = f( K Ic 12/4 / σ C 23/4 ) for solid particle erosion, as recently suggested in Rock Mech Rock Eng, doi: 10.1007/s00603-014-0658-x, 2014, is derived. The procedure is based on examinations of elastic-plastically responding rocks (rhyolite, granite) and plastically responding rocks (limestone, schist). The types of response are quantified through SEM inspections of eroded surfaces. Quantitative numbers for the distribution parameter K 1 are calculated for 30 rock materials, which cover a wide range of mechanical properties. The ranking according to the parameter K 1 is related to qualitative rock classification schemes. A modified proposal for the erosion of schist due to solid particle impingement at normal incidence is introduced.

  12. A Generalized Brownian Motion Model for Turbulent Relative Particle Dispersion

    NASA Astrophysics Data System (ADS)

    Shivamoggi, Bhimsen

    2015-11-01

    A generalized Brownian motion model has been applied to the turbulent relative particle dispersion problem (Shivamoggi). The fluctuating pressure forces acting on a fluid particle are taken to follow an Uhlenbeck-Ornstein process while it appears plausible to take their correlation time to have a power-law dependence on the flow Reynolds number Re. This ansatz provides an insight into the result that the Richardson-Obukhov scaling holds only in the infinite-Re limit and disappears otherwise. It provides a determination of the Richardson-Obukhov constant g as a function of Re, with an asymptotic constant value in the infinite-Re limit. This ansatz is further shown to be in quantitative agreement, in the small-Re limit, with the Batchelor-Townsend ansatz for the rate of change of the mean square interparticle separation in 3D FDT. My thanks to The Netherlands Organization for Scientific Research for Support.

  13. From flow and particle transport modeling to vibration isolation

    NASA Astrophysics Data System (ADS)

    Ellison, Joseph Fabian

    2001-08-01

    This thesis is composed of two parts. Part I is devoted to the analysis of particle transport and deposition. In this part, flow and particle transport and deposition in a furnace and under microgravity conditions are analyzed. In the first study of Part I, fluid flow, combustion, heat and mass transfer involved in a methane/air furnace were studied. The purpose of this study was to investigate variations in the flow field and thermal conditions in the furnace and to develop methods for improving its efficiency. The analysis of the combustor model was performed using an unstructured grid model developed with the Gambit grid generator of FLUENT version 5. In the second study of Part I, particle dispersion in a liquid filled box under orbital g-jitter excitation is analyzed. The study investigated particle motion experiments that were performed aboard the orbiting shuttle. The experiments have provided confusing data as to the nature of particle dispersion in the orbital environment. To obtain a better understanding of the dynamics involved, a series of numerical simulations are performed to study the dispersion of suspended particles subject to g-jitter excitations. Part II of the thesis is devoted to the analysis of vibration and vibration isolation problems. Following along the lines of vibration effecting system performance, a study of vibration isolation used to protect avionics equipment from adverse aircraft vibration environments was conducted. Passive isolation is the simplest means to achieve this goal. The system used here consisted of a circular steel ring with a lump mass on top and exposed to base excitation. Sinusoidal and filtered zero-mean Gaussian white noise were used to excite the structure and the acceleration response spectra at the top of the ring were computed. An experiment was performed to identify the natural frequencies and modal damping of the circular ring. The polished homogeneity measurement of large optics mounted in a vertical ring

  14. A particle based simulation model for glacier dynamics

    NASA Astrophysics Data System (ADS)

    Åström, J. A.; Riikilä, T. I.; Tallinen, T.; Zwinger, T.; Benn, D.; Moore, J. C.; Timonen, J.

    2013-10-01

    A particle-based computer simulation model was developed for investigating the dynamics of glaciers. In the model, large ice bodies are made of discrete elastic particles which are bound together by massless elastic beams. These beams can break, which induces brittle behaviour. At loads below fracture, beams may also break and reform with small probabilities to incorporate slowly deforming viscous behaviour in the model. This model has the advantage that it can simulate important physical processes such as ice calving and fracturing in a more realistic way than traditional continuum models. For benchmarking purposes the deformation of an ice block on a slip-free surface was compared to that of a similar block simulated with a Finite Element full-Stokes continuum model. Two simulations were performed: (1) calving of an ice block partially supported in water, similar to a grounded marine glacier terminus, and (2) fracturing of an ice block on an inclined plane of varying basal friction, which could represent transition to fast flow or surging. Despite several approximations, including restriction to two-dimensions and simplified water-ice interaction, the model was able to reproduce the size distributions of the debris observed in calving, which may be approximated by universal scaling laws. On a moderate slope, a large ice block was stable and quiescent as long as there was enough of friction against the substrate. For a critical length of frictional contact, global sliding began, and the model block disintegrated in a manner suggestive of a surging glacier. In this case the fragment size distribution produced was typical of a grinding process.

  15. Reduced quasilinear models for energetic particles interaction with Alfvenic eigenmodes

    NASA Astrophysics Data System (ADS)

    Ghantous, Katy

    The Line Broadened Quasilinear (LBQ) and the 1.5D reduced models are able to predict the effect of Alfvenic eigenmodes' interaction with energetic particles in burning plasmas. This interaction can result in energetic-particle losses that can damage the first wall, deteriorate the plasma performance, and even prevent ignition. The 1.5D model assumes a broad spectrum of overlapping modes and, based on analytic expressions for the growth and damping rates, calculates the pressure profiles that the energetic particles relax to upon interacting with the modes. 1.5D is validated with DIII-D experiments and predicted neutron losses consistent with observation. The model is employed to predict alpha-particle fusion-product losses in a large-scale operational parameter-space for burning plasmas. The LBQ model captures the interaction both in the regime of isolated modes as well as in the conventional regime of overlapping modes. Rules were established that allow quasilinear equations to replicate the expected steady-state saturation levels of isolated modes. The fitting formula is improved and the model is benchmarked with a Vlasov code, BOT. The saturation levels are accurately predicted and the mode evolution is well-replicated in the case of steady-state evolution where the collisions are high enough that coherent structures do not form. When the collisionality is low, oscillatory behavior can occur. LBQ can also exhibit non-steady behavior, but the onset of oscillations occurs for much higher collisional rates in BOT than in LBQ. For certain parameters of low collisionality, hole-clump creation and frequency chirping can occur which are not captured by the LBQ model. Also, there are cases of non-steady evolution without chirping which is possible for LBQ to study. However the results are inconclusive since the periods and amplitudes of the oscillations in the mode evolution are not well-replicated. If multiple modes exist, they can grow to the point of overlap which

  16. Fluctuating Nonlinear Spring Model of Mechanical Deformation of Biological Particles

    PubMed Central

    Kononova, Olga; Snijder, Joost; Kholodov, Yaroslav; Marx, Kenneth A.; Wuite, Gijs J. L.; Roos, Wouter H.; Barsegov, Valeri

    2016-01-01

    The mechanical properties of virus capsids correlate with local conformational dynamics in the capsid structure. They also reflect the required stability needed to withstand high internal pressures generated upon genome loading and contribute to the success of important events in viral infectivity, such as capsid maturation, genome uncoating and receptor binding. The mechanical properties of biological nanoparticles are often determined from monitoring their dynamic deformations in Atomic Force Microscopy nanoindentation experiments; but a comprehensive theory describing the full range of observed deformation behaviors has not previously been described. We present a new theory for modeling dynamic deformations of biological nanoparticles, which considers the non-linear Hertzian deformation, resulting from an indenter-particle physical contact, and the bending of curved elements (beams) modeling the particle structure. The beams’ deformation beyond the critical point triggers a dynamic transition of the particle to the collapsed state. This extreme event is accompanied by a catastrophic force drop as observed in the experimental or simulated force (F)-deformation (X) spectra. The theory interprets fine features of the spectra, including the nonlinear components of the FX-curves, in terms of the Young’s moduli for Hertzian and bending deformations, and the structural damage dependent beams’ survival probability, in terms of the maximum strength and the cooperativity parameter. The theory is exemplified by successfully describing the deformation dynamics of natural nanoparticles through comparing theoretical curves with experimental force-deformation spectra for several virus particles. This approach provides a comprehensive description of the dynamic structural transitions in biological and artificial nanoparticles, which is essential for their optimal use in nanotechnology and nanomedicine applications. PMID:26821264

  17. Fluctuating Nonlinear Spring Model of Mechanical Deformation of Biological Particles.

    PubMed

    Kononova, Olga; Snijder, Joost; Kholodov, Yaroslav; Marx, Kenneth A; Wuite, Gijs J L; Roos, Wouter H; Barsegov, Valeri

    2016-01-01

    The mechanical properties of virus capsids correlate with local conformational dynamics in the capsid structure. They also reflect the required stability needed to withstand high internal pressures generated upon genome loading and contribute to the success of important events in viral infectivity, such as capsid maturation, genome uncoating and receptor binding. The mechanical properties of biological nanoparticles are often determined from monitoring their dynamic deformations in Atomic Force Microscopy nanoindentation experiments; but a comprehensive theory describing the full range of observed deformation behaviors has not previously been described. We present a new theory for modeling dynamic deformations of biological nanoparticles, which considers the non-linear Hertzian deformation, resulting from an indenter-particle physical contact, and the bending of curved elements (beams) modeling the particle structure. The beams' deformation beyond the critical point triggers a dynamic transition of the particle to the collapsed state. This extreme event is accompanied by a catastrophic force drop as observed in the experimental or simulated force (F)-deformation (X) spectra. The theory interprets fine features of the spectra, including the nonlinear components of the FX-curves, in terms of the Young's moduli for Hertzian and bending deformations, and the structural damage dependent beams' survival probability, in terms of the maximum strength and the cooperativity parameter. The theory is exemplified by successfully describing the deformation dynamics of natural nanoparticles through comparing theoretical curves with experimental force-deformation spectra for several virus particles. This approach provides a comprehensive description of the dynamic structural transitions in biological and artificial nanoparticles, which is essential for their optimal use in nanotechnology and nanomedicine applications. PMID:26821264

  18. Closure Models for Turbulent Particle-laden Flows from Particle-resolved Direct Numerical Simulation

    NASA Astrophysics Data System (ADS)

    Subramaniam, Shankar; Tenneti, Sudheer; Mehrabadi, Mohammad; Garg, Rahul

    2012-11-01

    Gas-phase velocity fluctuations in fixed particle beds and freely evolving suspensions are quantified using a particle-resolved direct numerical simulation (PR-DNS). The flow regime corresponds to gas-solid systems typically encountered in fluidized bed risers, with high solid to gas density ratio and particle diameter being greater than the dissipative length scales. The kinetic energy associated with gas-phase velocity fluctuations in homogeneous monodisperse fixed beds is characterized as a function of solid volume fraction φ and the Reynolds number based on the mean slip velocity Re. A simple scaling analysis is used to explain the dependence of k on ɛ and Re. The steady value of k results from the balance between the source of k due to interphase transfer of kinetic energy, and the dissipation rate (ɛ) of k in the gas-phase. It is found that the dissipation rate of k in gas-solid flows can be modeled using a length scale that is analogous to the Taylor microscale used in single-phase turbulence. Using the PR-DNS data for k and ɛ we also infer an eddy viscosity for gas-solid flow. For the parameter values considered here, the level of gas-phase velocity fluctuations in freely evolving suspensions differs by only 10% from the value for the corresponding fixed beds. Funded in part by the US Department of Energy's National Energy Technology Laboratory Grant DE-FC26-07NT43098 (Advanced Research) and the National Science Foundation's grant CBET 1134500.

  19. Dynamics of magnetic particles near a surface: model and experiments on field-induced disaggregation.

    PubMed

    van Reenen, A; Gao, Y; de Jong, A M; Hulsen, M A; den Toonder, J M J; Prins, M W J

    2014-04-01

    Magnetic particles are widely used in biological research and bioanalytical applications. As the corresponding tools are progressively being miniaturized and integrated, the understanding of particle dynamics and the control of particles down to the level of single particles become important. Here, we describe a numerical model to simulate the dynamic behavior of ensembles of magnetic particles, taking account of magnetic interparticle interactions, interactions with the liquid medium and solid surfaces, as well as thermal diffusive motion of the particles. The model is verified using experimental data of magnetic field-induced disaggregation of magnetic particle clusters near a physical surface, wherein the magnetic field properties, particle size, cluster size, and cluster geometry were varied. Furthermore, the model clarifies how the cluster configuration, cluster alignment, magnitude of the field gradient, and the field repetition rate play a role in the particle disaggregation process. The simulation model will be very useful for further in silico studies on magnetic particle dynamics in biotechnological tools. PMID:24827250

  20. Statistical thermal model analysis of particle production at LHC

    NASA Astrophysics Data System (ADS)

    Karasu Uysal, A.; Vardar, N.

    2016-04-01

    A successful description of the particle ratios measured in heavy-ion collisions has been achieved in the framework of thermal models. In such a way, a large number of observables can be reproduced with a small number of parameters, namely the temperature, baryo-chemical potential and a factor measuring the degree of strangeness saturation. The comparison of experimental data at and the model estimations has made possible to define the thermodynamic parameters of strongly interacting matter at chemical freeze-out temperature. The detailed study of hadron and meson production including resonances using the statistical-thermal model is discussed. Their ratios are compared with the existing experimental data and predictions are made for pp and heavy-ion collisions at RHIC and LHC energies.

  1. Density-independent model of self-propelled particles.

    PubMed

    Schubring, Daniel; Ohmann, Paul R

    2013-09-01

    We examine a density-independent modification of the Vicsek model in which a particle interacts with neighbors defined by Delaunay triangulation. To feasibly simulate the model, an algorithm for repairing the triangulation over time was developed. This algorithm may also be applied to any time varying two-dimensional Delaunay triangulation. This model exhibits a continuous phase transition with noise, and a distinct set of critical exponents were measured which satisfy a hyperscaling relationship. The critical exponents are found to vary between a low and high velocity regime, but they are robust under the inclusion of a repulsive interaction. We present evidence that the correlation length approximately scales with the size of the system even in the ordered phase. PMID:24125215

  2. Modeling crowd turbulence by many-particle simulations

    NASA Astrophysics Data System (ADS)

    Yu, Wenjian; Johansson, Anders

    2007-10-01

    A recent study [D. Helbing, A. Johansson, and H. Z. Al-Abideen, Phys. Rev. E 75, 046109 (2007)] has revealed a “turbulent” state of pedestrian flows, which is characterized by sudden displacements and causes the falling and trampling of people. However, turbulent crowd motion is not reproduced well by current many-particle models due to their insufficient representation of the local interactions in areas of extreme densities. In this contribution, we extend the repulsive force term of the social force model to reproduce crowd turbulence. We perform numerical simulations of pedestrians moving through a bottleneck area with this model. The transitions from laminar to stop-and-go and turbulent flows are observed. The empirical features characterizing crowd turbulence, such as the structure function and the probability density function of velocity increments, are reproduced well; i.e., they are well compatible with an analysis of video data during the annual Muslim pilgrimage.

  3. Emergent smectic order in simple active particle models

    NASA Astrophysics Data System (ADS)

    Romanczuk, Pawel; Chaté, Hugues; Chen, Leiming; Ngo, Sandrine; Toner, John

    2016-06-01

    Novel ‘smectic-P’ behavior, in which self-propelled particles form rows and move on average along them, occurs generically within the orientationally ordered phase of simple models that we simulate. Both apolar (head–tail symmetric) and polar (head–tail asymmetric) models with aligning and repulsive interactions exhibit slow algebraic decay of smectic order with system size up to some finite length scale, after which faster decay occurs. In the apolar case, this scale is that of an undulation instability of the rows. In the polar case, this instability is absent, but traveling fluctuations disrupt the rows in large systems and motion and smectic order may spontaneously globally rotate. These observations agree with a new hydrodynamic theory which we present here. Variants of our models also exhibit active smectic ‘A’ and ‘C’ order, with motion orthogonal and oblique to the layers respectively.

  4. Unsteady Free-Wake Vortex Particle Model for HAWT

    NASA Astrophysics Data System (ADS)

    Bogateanu, R.; Frunzulicǎ, F.; Cardos, V.

    2010-09-01

    In the design of horizontal axis wind turbines (HAWT) one problem is to determine the aeroelastic behaviour of the rotor blades for the various wind inflow conditions. A step in this process is to predict with accuracy the aerodynamic loads on the blades. The Vortex Lattice Method (VLM) provides a transparent investigation concerning the role of various physical parameters which influence the aerodynamic problem. In this paper we present a method for the calculation of the non-uniform induced downwash of a HAWT rotor using the vortex ring model for the lifting surface coupled with an unsteady free-wake vortex particle model. Comparative studies between results obtained with different models of wake for a generic HAWT are presented.

  5. Modeling of mesoscopic electrokinetic phenomena using charged dissipative particle dynamics

    NASA Astrophysics Data System (ADS)

    Deng, Mingge; Li, Zhen; Karniadakis, George

    2015-11-01

    In this work, we propose a charged dissipative particle dynamics (cDPD) model for investigation of mesoscopic electrokinetic phenomena. In particular, this particle-based method was designed to simulate micro- or nano- flows which governing by Poisson-Nernst-Planck (PNP) equation coupled with Navier-Stokes (NS) equation. For cDPD simulations of wall-bounded fluid systems, a methodology for imposing correct Dirichlet and Neumann boundary conditions for both PNP and NS equations is developed. To validate the present cDPD model and the corresponding boundary method, we perform cDPD simulations of electrostatic double layer (EDL) in the vicinity of a charged wall, and the results show good agreement with the mean-field theoretical solutions. The capacity density of a parallel plate capacitor in salt solution is also investigated with different salt concentration. Moreover, we utilize the proposed methodology to study the electroosmotic and electroosmotic/pressure-driven flow in a micro-channel. In the last, we simulate the dilute polyelectrolyte solution both in bulk and micro-channel, which show the flexibility and capability of this method in studying complex fluids. This work was sponsored by the Collaboratory on Mathematics for Mesoscopic Modeling of Materials (CM4) supported by DOE.

  6. Particle-gamma and particle-particle correlations in nuclear reactions using Monte Carlo Hauser-Feshback model

    SciTech Connect

    Kawano, Toshihiko; Talou, Patrick; Watanabe, Takehito; Chadwick, Mark

    2010-01-01

    Monte Carlo simulations for particle and {gamma}-ray emissions from an excited nucleus based on the Hauser-Feshbach statistical theory are performed to obtain correlated information between emitted particles and {gamma}-rays. We calculate neutron induced reactions on {sup 51}V to demonstrate unique advantages of the Monte Carlo method. which are the correlated {gamma}-rays in the neutron radiative capture reaction, the neutron and {gamma}-ray correlation, and the particle-particle correlations at higher energies. It is shown that properties in nuclear reactions that are difficult to study with a deterministic method can be obtained with the Monte Carlo simulations.

  7. Efficiency of a statistical transport model for turbulent particle dispersion

    NASA Astrophysics Data System (ADS)

    Litchford, Ron J.; Jeng, San-Mou

    1992-05-01

    In developing its theory for turbulent dispersion transport, the Litchford and Jeng (1991) statistical transport model for turbulent particle dispersion took a generalized approach in which the perturbing influence of each turbulent eddy on consequent interactions was transported through all subsequent eddies. Nevertheless, examinations of this transport relation shows it to be able to decay rapidly: this implies that additional computational efficiency may be obtained via truncation of unneccessary transport terms. Attention is here given to the criterion for truncation, as well as to expected efficiency gains.

  8. Efficiency of a statistical transport model for turbulent particle dispersion

    SciTech Connect

    Litchford, R.J.; Jeng, San-Mou )

    1992-05-01

    In developing its theory for turbulent dispersion transport, the Litchford and Jeng (1991) statistical transport model for turbulent particle dispersion took a generalized approach in which the perturbing influence of each turbulent eddy on consequent interactions was transported through all subsequent eddies. Nevertheless, examinations of this transport relation shows it to be able to decay rapidly: this implies that additional computational efficiency may be obtained via truncation of unneccessary transport terms. Attention is here given to the criterion for truncation, as well as to expected efficiency gains. 2 refs.

  9. Efficiency of a statistical transport model for turbulent particle dispersion

    NASA Technical Reports Server (NTRS)

    Litchford, Ron J.; Jeng, San-Mou

    1992-01-01

    In developing its theory for turbulent dispersion transport, the Litchford and Jeng (1991) statistical transport model for turbulent particle dispersion took a generalized approach in which the perturbing influence of each turbulent eddy on consequent interactions was transported through all subsequent eddies. Nevertheless, examinations of this transport relation shows it to be able to decay rapidly: this implies that additional computational efficiency may be obtained via truncation of unneccessary transport terms. Attention is here given to the criterion for truncation, as well as to expected efficiency gains.

  10. How to model the interaction of charged Janus particles

    NASA Astrophysics Data System (ADS)

    Hieronimus, Reint; Raschke, Simon; Heuer, Andreas

    2016-08-01

    We analyze the interaction of charged Janus particles including screening effects. The explicit interaction is mapped via a least square method on a variable number n of systematically generated tensors that reflect the angular dependence of the potential. For n = 2 we show that the interaction is equivalent to a model previously described by Erdmann, Kröger, and Hess (EKH). Interestingly, this mapping is for n = 2 not able to capture the subtleties of the interaction for small screening lengths. Rather, a larger number of tensors has to be used. We find that the characteristics of the Janus type interaction plays an important role for the aggregation behavior. We obtained cluster structures up to the size of 13 particles for n = 2 and 36 and screening lengths κ-1 = 0.1 and 1.0 via Monte Carlo simulations. The influence of the screening length is analyzed and the structures are compared to results for an electrostatic-type potential and for the multipole-expanded Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. We find that a dipole-like potential (EKH or dipole DLVO approximation) is not able to sufficiently reproduce the anisotropy effects of the potential. Instead, a higher order expansion has to be used to obtain cluster structures that are compatible with experimental observations. The resulting minimum-energy clusters are compared to those of sticky hard sphere systems. Janus particles with a short-range screened interaction resemble sticky hard sphere clusters for all considered particle numbers, whereas for long-range screening even very small clusters are structurally different.

  11. Event-based total suspended sediment particle size distribution model

    NASA Astrophysics Data System (ADS)

    Thompson, Jennifer; Sattar, Ahmed M. A.; Gharabaghi, Bahram; Warner, Richard C.

    2016-05-01

    One of the most challenging modelling tasks in hydrology is prediction of the total suspended sediment particle size distribution (TSS-PSD) in stormwater runoff generated from exposed soil surfaces at active construction sites and surface mining operations. The main objective of this study is to employ gene expression programming (GEP) and artificial neural networks (ANN) to develop a new model with the ability to more accurately predict the TSS-PSD by taking advantage of both event-specific and site-specific factors in the model. To compile the data for this study, laboratory scale experiments using rainfall simulators were conducted on fourteen different soils to obtain TSS-PSD. This data is supplemented with field data from three construction sites in Ontario over a period of two years to capture the effect of transport and deposition within the site. The combined data sets provide a wide range of key overlooked site-specific and storm event-specific factors. Both parent soil and TSS-PSD in runoff are quantified by fitting each to a lognormal distribution. Compared to existing regression models, the developed model more accurately predicted the TSS-PSD using a more comprehensive list of key model input parameters. Employment of the new model will increase the efficiency of deployment of required best management practices, designed based on TSS-PSD, to minimize potential adverse effects of construction site runoff on aquatic life in the receiving watercourses.

  12. Particle Tracking Model and Abstraction of Transport Processes

    SciTech Connect

    B. Robinson

    2004-10-21

    The purpose of this report is to document the abstraction model being used in total system performance assessment (TSPA) model calculations for radionuclide transport in the unsaturated zone (UZ). The UZ transport abstraction model uses the particle-tracking method that is incorporated into the finite element heat and mass model (FEHM) computer code (Zyvoloski et al. 1997 [DIRS 100615]) to simulate radionuclide transport in the UZ. This report outlines the assumptions, design, and testing of a model for calculating radionuclide transport in the UZ at Yucca Mountain. In addition, methods for determining and inputting transport parameters are outlined for use in the TSPA for license application (LA) analyses. Process-level transport model calculations are documented in another report for the UZ (BSC 2004 [DIRS 164500]). Three-dimensional, dual-permeability flow fields generated to characterize UZ flow (documented by BSC 2004 [DIRS 169861]; DTN: LB03023DSSCP9I.001 [DIRS 163044]) are converted to make them compatible with the FEHM code for use in this abstraction model. This report establishes the numerical method and demonstrates the use of the model that is intended to represent UZ transport in the TSPA-LA. Capability of the UZ barrier for retarding the transport is demonstrated in this report, and by the underlying process model (BSC 2004 [DIRS 164500]). The technical scope, content, and management of this report are described in the planning document ''Technical Work Plan for: Unsaturated Zone Transport Model Report Integration'' (BSC 2004 [DIRS 171282]). Deviations from the technical work plan (TWP) are noted within the text of this report, as appropriate. The latest version of this document is being prepared principally to correct parameter values found to be in error due to transcription errors, changes in source data that were not captured in the report, calculation errors, and errors in interpretation of source data.

  13. Dissipative particle dynamics modeling of blood flow in arterial bifurcations

    NASA Astrophysics Data System (ADS)

    Li, Xuejin; Lykov, Kirill; Pivkin, Igor V.; Karniadakis, George Em

    2013-11-01

    The motion of a suspension of red blood cells (RBCs) flowing in bifurcations is investigated using both low-dimensional RBC (LD-RBC) and multiscale RBC (MS-RBC) models based on dissipative particle dynamics (DPD). The blood flow is first simulated in a symmetric geometry between the diverging and converging channels to satisfy the periodic flow assumption along the flow direction. The results show that the flowrate ratio of the daughter channels and the feed hematocrit level has considerable influence on blood-plasma separation. We also propose a new method to model the inflow and outflow boundaries for the blood flow simulations: the inflow at the inlet is duplicated from a fully developed flow generated by DPD fluid with periodic boundary conditions; the outflow in two adjacent regions near the outlet is controlled by adaptive forces to keep the flowrate and velocity gradient equal, while the particles leaving the microfluidic channel at the outlet at each time step are removed from the system. The simulation results of the developing flow match analytical solutions from continuum theory. Plasma skimming and the all-or-nothing phenomenon of RBCs in bifurcation have been investigated in the simulations. The simulation results are consistent with previous experimental results and theoretical predictions. This work is supported by the NIH Grant R01HL094270.

  14. Sterile Neutrinos in Non-Standard Cosmologies and Particle Models

    NASA Astrophysics Data System (ADS)

    Osoba, Efunwande

    2010-12-01

    The discovery of neutrino masses suggests that the Standard Model should be supplemented with new gauge-singlet fermions, often called sterile neutrinos. These sterile neutrinos can shed new light on open questions in cosmology. I will highlight some interesting contributions that sterile neutrinos bring to the understanding of cosmology. In this dissertation, I will show a novel way in which sterile neutrinos could be a dark matter candidate in the form of "Inert-Sterile" neutrinos. In usual particle models, sterile neutrinos can account for the dark matter of the Universe only if they have niasses in the keV range and are warm dark matter. Stringent cosmological and astrophysical bounds, in particular imposed by X-ray observations, apply to them. I will point out that in a particular variation of the Inert Doublet Model, sterile neutrinos can account for the dark matter in the Universe and may be either cold or warm dark matter candidates, even for masses much larger than the keV range. These "Inert-Sterile" neutrinos, produced non-thermally in the early Universe, would be stable and have very small couplings to Standard Model particles, rendering very difficult their detection in either direct or indirect dark matter searches. They could be, in principle, revealed in colliders by discovering other particles in the model. I also show how the existence of the sterile neutrino may force us to rethink the standard cosmology. It is commonly assumed that the cosmological and astrophysical bounds on the mixings of sterile with active neutrinos are much more stringent than those obtained from laboratory measurements. In this dissertation, I show that in scenarios with a very low reheating temperature at the end of (the last episode of) inflation or entropy creation, the abundance of heavy (> 1 MeV) sterile neutrinos becomes largely suppressed with respect to that obtained within the standard framework. Thus, in this case cosmological bounds become much less stringent

  15. Particle Acceleration in the Low Corona Over Broad Longitudes: Coupling Between 3D Magnetohydrodynamic and Energetic Particle Models

    NASA Astrophysics Data System (ADS)

    Gorby, M.; Schwadron, N.; Lee, M. A.; Booth, A. C.; Spence, H.; Torok, T.; Downs, C.; Lionello, R.; Linker, J.; Titov, V. S.; Mikic, Z.; Riley, P.; Desai, M. I.; Dayeh, M. A.; Kozarev, K. A.

    2013-12-01

    Recent work on the coupling between 3D energetic particle models (e.g., the Energetic Particle Radiation Environment Model, EPREM) and magnetohydrodynamic (MHD) models of Coronal Mass Ejections (CMEs, e.g., the PSI MAS model) has demonstrated the efficacy of compression regions in front of fast CMEs for particle acceleration from remarkably low in the corona (3-6 solar radii). Typically particle acceleration becomes rapid beyond 3Rs and often in regions where compression regions have not yet formed active shocks. The challenge for forming large SEP events in such compression-acceleration scenarios is to have enhanced scattering within the acceleration region while also allowing for efficient escape of accelerated particles downstream (away from the Sun) from the compression region. Simulations show rapid particle acceleration in the range of 3-8 Rs over a broad longitudinal region (80°) resulting from the pile-up of magnetic flux in the compression and the subsequent expansion of these fields. These results have important implications for multi-instrument observations that will allow Solar Probe Plus and Solar Orbiter to test the developing paradigm of SEP acceleration and transport from coronal compressions. We present here recent coupled simulations for SEP acceleration and transport, including energetic particle and CME plasma profiles. The broadness of the longitudinal profile from such events may be a key observational test of compression acceleration in the low corona.

  16. Fractality à la carte: a general particle aggregation model.

    PubMed

    Nicolás-Carlock, J R; Carrillo-Estrada, J L; Dossetti, V

    2016-01-01

    In nature, fractal structures emerge in a wide variety of systems as a local optimization of entropic and energetic distributions. The fractality of these systems determines many of their physical, chemical and/or biological properties. Thus, to comprehend the mechanisms that originate and control the fractality is highly relevant in many areas of science and technology. In studying clusters grown by aggregation phenomena, simple models have contributed to unveil some of the basic elements that give origin to fractality, however, the specific contribution from each of these elements to fractality has remained hidden in the complex dynamics. Here, we propose a simple and versatile model of particle aggregation that is, on the one hand, able to reveal the specific entropic and energetic contributions to the clusters' fractality and morphology, and, on the other, capable to generate an ample assortment of rich natural-looking aggregates with any prescribed fractal dimension. PMID:26781204

  17. A computer model of particle balance in ECR ion sources

    NASA Astrophysics Data System (ADS)

    Shirkov, G. D.; Philippov, A. V.

    2008-12-01

    The investigation of the widespread model of particle balance and energy transport [1-5] for calculating the ion charge state distribution (CSD) in an electron cyclotron resonance (ECR) ion source [6] is given. The modification of this model that makes it possible to describe the confinement and accumulation processes of highly charged ions in ECR plasma for the case of gas mixing is more precisely discussed. The discussion of the new technique for calculating the time confinement of ions and electrons, which is based on the Pastukhov theory [7, 8], is given: calculation of confinement times during two step minimization of special type functionals. The preliminary results obtained with this approach have been compared with available experimental data.

  18. Fractality à la carte: a general particle aggregation model

    NASA Astrophysics Data System (ADS)

    Nicolás-Carlock, J. R.; Carrillo-Estrada, J. L.; Dossetti, V.

    2016-01-01

    In nature, fractal structures emerge in a wide variety of systems as a local optimization of entropic and energetic distributions. The fractality of these systems determines many of their physical, chemical and/or biological properties. Thus, to comprehend the mechanisms that originate and control the fractality is highly relevant in many areas of science and technology. In studying clusters grown by aggregation phenomena, simple models have contributed to unveil some of the basic elements that give origin to fractality, however, the specific contribution from each of these elements to fractality has remained hidden in the complex dynamics. Here, we propose a simple and versatile model of particle aggregation that is, on the one hand, able to reveal the specific entropic and energetic contributions to the clusters’ fractality and morphology, and, on the other, capable to generate an ample assortment of rich natural-looking aggregates with any prescribed fractal dimension.

  19. Fractality à la carte: a general particle aggregation model

    PubMed Central

    Nicolás-Carlock, J. R.; Carrillo-Estrada, J. L.; Dossetti, V.

    2016-01-01

    In nature, fractal structures emerge in a wide variety of systems as a local optimization of entropic and energetic distributions. The fractality of these systems determines many of their physical, chemical and/or biological properties. Thus, to comprehend the mechanisms that originate and control the fractality is highly relevant in many areas of science and technology. In studying clusters grown by aggregation phenomena, simple models have contributed to unveil some of the basic elements that give origin to fractality, however, the specific contribution from each of these elements to fractality has remained hidden in the complex dynamics. Here, we propose a simple and versatile model of particle aggregation that is, on the one hand, able to reveal the specific entropic and energetic contributions to the clusters’ fractality and morphology, and, on the other, capable to generate an ample assortment of rich natural-looking aggregates with any prescribed fractal dimension. PMID:26781204

  20. Model of cosmology and particle physics at an intermediate scale

    SciTech Connect

    Bastero-Gil, M.; Di Clemente, V.; King, S. F.

    2005-05-15

    We propose a model of cosmology and particle physics in which all relevant scales arise in a natural way from an intermediate string scale. We are led to assign the string scale to the intermediate scale M{sub *}{approx}10{sup 13} GeV by four independent pieces of physics: electroweak symmetry breaking; the {mu} parameter; the axion scale; and the neutrino mass scale. The model involves hybrid inflation with the waterfall field N being responsible for generating the {mu} term, the right-handed neutrino mass scale, and the Peccei-Quinn symmetry breaking scale. The large scale structure of the Universe is generated by the lightest right-handed sneutrino playing the role of a coupled curvaton. We show that the correct curvature perturbations may be successfully generated providing the lightest right-handed neutrino is weakly coupled in the seesaw mechanism, consistent with sequential dominance.

  1. Particle-based simulation of ellipse-shaped particle aggregation as a model for vascular network formation

    NASA Astrophysics Data System (ADS)

    Palachanis, Dimitrios; Szabó, András; Merks, Roeland M. H.

    2015-12-01

    Computational modeling is helpful for elucidating the cellular mechanisms driving biological morphogenesis. Previous simulation studies of blood vessel growth based on the cellular Potts model proposed that elongated, adhesive or mutually attractive endothelial cells suffice for the formation of blood vessel sprouts and vascular networks. Because each mathematical representation of a model introduces potential artifacts, it is important that model results are reproduced using alternative modeling paradigms. Here, we present a lattice-free, particle-based simulation of the cell elongation model of vasculogenesis. The new, particle-based simulations confirm the results obtained from the previous cellular Potts simulations. Furthermore, our current findings suggest that the emergence of order is possible with the application of a high enough attractive force or, alternatively, a longer attraction radius. The methodology will be applicable to a range of problems in morphogenesis and noisy particle aggregation in which cell shape is a key determining factor.

  2. A new coincidence model for single particle counters, Part I: Theory and experimental verification.

    PubMed

    Knapp, J Z; Abramson, L R

    1994-01-01

    The prerequisites for estimating the effect of signal coincidence on both particle undercounting and the injection of false counts in the implementation of U.S.P. 788 contaminating particle assays by light extinction particle counters are defined. These include a particle concentration measure that varies with particle size and a new model of the counting process. Both prerequisites have been verified empirically: a single normalized equation describes the coincidence effect in all single particle counters. The single parameter of the normalized equation is the number of effective detector volumes per milliliter. A maximum undercount limit of 5% is proposed based on adequately suspended particles. Using the SVP U.S.P. XXII acceptance limits of 10,000 particles per container or the PMA propose 6,000 particles per container maximum for particles > 10 microns in U.S.P. XXIII, undercount errors are estimated for the smallest container sizes. The large concentration of particles below the controlled 10 microns particle size, that has been documented in injectable solutions, can pose an additional 788 measurement hazard. A Poisson model is used to estimate and control the injection of false particle counts into the mandated measurement through particle coincidence. Acceptable counting accuracy limits with present particle counting systems can be achieved by understanding the capabilities of the particle counter measurement system and using a dilution technique when appropriate. The new model of the counting process and the new particle concentration measures can result in standard, conservative, instrument specifications for use in Pharmacopeial contamination testing and in GLP user evaluation tests. Part I of this paper includes the theory of the coincidence effect on particle counting and the particle size distribution measured. A summary of the experimental verification employed to determine coincidence count loss as a function of particle concentration for single

  3. Predictive modeling of particle-laden, turbulent flows

    SciTech Connect

    Sinclair, J.L.

    1992-01-01

    The successful prediction of particle-laden, turbulent flows relies heavily on the representation of turbulence in the gas phase. Several types of turbulence models for single-phase gas flow have been developed which compare reasonably well with experimental data. In the present work, a low-Reynolds'' k-[epsilon], closure model is chosen to describe the Reynolds stresses associated with gas-phase turbulence. This closure scheme, which involves transport equations for the turbulent kinetic energy and its dissipation rate, is valid in the turbulent core as well as the viscous sublayer. Several versions of the low-Reynolds k-[epsilon] closure are documented in the literature. However, even those models which are similar in theory often differ considerably in their quantitative and qualitative predictions, making the selection of such a model a difficult task. The purpose of this progress report is to document our findings on the performance of ten different versions of the low-Reynolds k-[epsilon] model on predicting fully developed pipe flow. The predictions are compared with the experimental data of Schildknecht, et al. (1979). With the exception of the model put forth by Hoffman (1975), the predictions of all the closures show reasonable agreement for the mean velocity profile. However, important quantitative differences exist for the turbulent kinetic energy profile. In addition, the predicted eddy viscosity profile and the wall-region profile of the turbulent kinetic energy dissipation rate exhibit both quantitative and qualitative differences. An effort to extend the present comparisons to include experimental measurements of other researchers is recommended in order to further evaluate the performance of the models.

  4. Model for boiling and dryout in particle beds. [LMFBR

    SciTech Connect

    Lipinski, R. J.

    1982-06-01

    Over the last ten years experiments and modeling of dryout in particle beds have produced over fifty papers. Considering only volume-heated beds, over 250 dryout measurements have been made, and are listed in this work. In addition, fifteen models to predict dryout have been produced and are discussed. A model is developed in this report for one-dimensional boiling and dryout in a porous medium. It is based on conservation laws for mass, momentum, and energy. The initial coupled differential equations are reduced to a single first-order differential equation with an algebraic equation for the upper boundary condition. The model includes the effects of both laminar and turbulent flow, two-phase friction, and capillary force. The boundary condition at the bed bottom includes the possibility of inflowing liquid and either an adiabatic or a bottom-cooled support structure. The top of the bed may be either channeled or subcooled. In the first case the channel length and the saturation at the base of the channels are predicted. In the latter case, a criterion for penetration of the subcooled zone by channels is obtained.

  5. Independent-particle models for light negative atomic ions

    NASA Technical Reports Server (NTRS)

    Ganas, P. S.; Talman, J. D.; Green, A. E. S.

    1980-01-01

    For the purposes of astrophysical, aeronomical, and laboratory application, a precise independent-particle model for electrons in negative atomic ions of the second and third period is discussed. The optimum-potential model (OPM) of Talman et al. (1979) is first used to generate numerical potentials for eight of these ions. Results for total energies and electron affinities are found to be very close to Hartree-Fock solutions. However, the OPM and HF electron affinities both depart significantly from experimental affinities. For this reason, two analytic potentials are developed whose inner energy levels are very close to the OPM and HF levels but whose last electron eigenvalues are adjusted precisely with the magnitudes of experimental affinities. These models are: (1) a four-parameter analytic characterization of the OPM potential and (2) a two-parameter potential model of the Green, Sellin, Zachor type. The system O(-) or e-O, which is important in upper atmospheric physics is examined in some detail.

  6. Modeling Proton Dissociation and Transfer Using Dissipative Particle Dynamics Simulation.

    PubMed

    Lee, Ming-Tsung; Vishnyakov, Aleksey; Neimark, Alexander V

    2015-09-01

    We suggest a coarse-grained model for dissipative particle dynamics (DPD) simulations of solutions with dissociated protons. The model uses standard short-range soft repulsion and smeared charge electrostatic potentials between the beads, representing solution components. The proton is introduced as a separate charged bead that forms dissociable bonds with proton receptor base beads, such as water or deprotonated acid anions. The proton-base bonds are described by Morse potentials. When the proton establishes the Morse bonds with two bases, they form an intermediate complex, and the proton is able to "hop" between the bases artificially mimicking the Grotthuss diffusion mechanism. By adjusting the Morse potential parameters, one can regulate the potential barrier associated with intermediate complex formation and breakup and control the hopping frequency. This makes the proposed model applicable to simulations of proton mobility and reaction equilibria between protonated and deprotonated acid forms in aqueous solutions. The proposed model provides quantitative agreement with experiments for the proton self-diffusion coefficient and hopping frequency, as well as for the degree of dissociation of benzenesulfonic acid. PMID:26575931

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

  8. Validation of chemistry models employed in a particle simulation method

    NASA Technical Reports Server (NTRS)

    Haas, Brian L.; Mcdonald, Jeffrey D.

    1991-01-01

    The chemistry models employed in a statistical particle simulation method, as implemented in the Intel iPSC/860 multiprocessor computer, are validated and applied. Chemical relaxation of five-species air in these reservoirs involves 34 simultaneous dissociation, recombination, and atomic-exchange reactions. The reaction rates employed in the analytic solutions are obtained from Arrhenius experimental correlations as functions of temperature for adiabatic gas reservoirs in thermal equilibrium. Favorable agreement with the analytic solutions validates the simulation when applied to relaxation of O2 toward equilibrium in reservoirs dominated by dissociation and recombination, respectively, and when applied to relaxation of air in the temperature range 5000 to 30,000 K. A flow of O2 over a circular cylinder at high Mach number is simulated to demonstrate application of the method to multidimensional reactive flows.

  9. Model-independent analyses of dark-matter particle interactions

    DOE PAGESBeta

    Anand, Nikhil; Fitzpatrick, A. Liam; Haxton, W. C.

    2015-03-24

    A model-independent treatment of dark-matter particle elastic scattering has been developed, yielding the most general interaction for WIMP-nucleon low-energy scattering, and the resulting amplitude has been embedded into the nucleus, taking into account the selection rules imposed by parity and time-reversal. One finds that, in contrast to the usual spin-independent/spin-dependent (SI/SD) formulation, the resulting cross section contains six independent nuclear response functions, three of which are associated with possible velocity-dependent interactions. We find that current experiments are four orders of magnitude more sensitive to derivative couplings than is apparent in the standard SI/SD treatment, which necessarily associated such interactions withmore » cross sections proportional to v2T ~ 10⁻⁶, where vT is the WIMP velocity relative to the center of mass of the nuclear target.« less

  10. Model-independent analyses of dark-matter particle interactions

    SciTech Connect

    Anand, Nikhil; Fitzpatrick, A. Liam; Haxton, W. C.

    2015-03-24

    A model-independent treatment of dark-matter particle elastic scattering has been developed, yielding the most general interaction for WIMP-nucleon low-energy scattering, and the resulting amplitude has been embedded into the nucleus, taking into account the selection rules imposed by parity and time-reversal. One finds that, in contrast to the usual spin-independent/spin-dependent (SI/SD) formulation, the resulting cross section contains six independent nuclear response functions, three of which are associated with possible velocity-dependent interactions. We find that current experiments are four orders of magnitude more sensitive to derivative couplings than is apparent in the standard SI/SD treatment, which necessarily associated such interactions with cross sections proportional to v2T ~ 10⁻⁶, where vT is the WIMP velocity relative to the center of mass of the nuclear target.

  11. A Simple Mathematical Model for Standard Model of Elementary Particles and Extension Thereof

    NASA Astrophysics Data System (ADS)

    Sinha, Ashok

    2016-03-01

    An algebraically (and geometrically) simple model representing the masses of the elementary particles in terms of the interaction (strong, weak, electromagnetic) constants is developed, including the Higgs bosons. The predicted Higgs boson mass is identical to that discovered by LHC experimental programs; while possibility of additional Higgs bosons (and their masses) is indicated. The model can be analyzed to explain and resolve many puzzles of particle physics and cosmology including the neutrino masses and mixing; origin of the proton mass and the mass-difference between the proton and the neutron; the big bang and cosmological Inflation; the Hubble expansion; etc. A novel interpretation of the model in terms of quaternion and rotation in the six-dimensional space of the elementary particle interaction-space - or, equivalently, in six-dimensional spacetime - is presented. Interrelations among particle masses are derived theoretically. A new approach for defining the interaction parameters leading to an elegant and symmetrical diagram is delineated. Generalization of the model to include supersymmetry is illustrated without recourse to complex mathematical formulation and free from any ambiguity. This Abstract represents some results of the Author's Independent Theoretical Research in Particle Physics, with possible connection to the Superstring Theory. However, only very elementary mathematics and physics is used in my presentation.

  12. Particle dispersion in homogeneous turbulence using the one-dimensional turbulence model

    SciTech Connect

    Sun, Guangyuan; Lignell, David O.; Hewson, John C.; Gin, Craig R.

    2014-10-09

    Lagrangian particle dispersion is studied using the one-dimensional turbulence (ODT) model in homogeneous decaying turbulence configurations. The ODT model has been widely and successfully applied to a number of reacting and nonreacting flow configurations, but only limited application has been made to multiphase flows. We present a version of the particle implementation and interaction with the stochastic and instantaneous ODT eddy events. The model is characterized by comparison to experimental data of particle dispersion for a range of intrinsic particle time scales and body forces. Particle dispersion, velocity, and integral time scale results are presented. Moreover, the particle implementation introduces a single model parameter β p , and sensitivity to this parameter and behavior of the model are discussed. Good agreement is found with experimental data and the ODT model is able to capture the particle inertial and trajectory crossing effects. Our results serve as a validation case of the multiphase implementations of ODT for extensions to other flow configurations.

  13. Particle dispersion in homogeneous turbulence using the one-dimensional turbulence model

    SciTech Connect

    Sun, Guangyuan Lignell, David O.; Hewson, John C.; Gin, Craig R.

    2014-10-15

    Lagrangian particle dispersion is studied using the one-dimensional turbulence (ODT) model in homogeneous decaying turbulence configurations. The ODT model has been widely and successfully applied to a number of reacting and nonreacting flow configurations, but only limited application has been made to multiphase flows. Here, we present a version of the particle implementation and interaction with the stochastic and instantaneous ODT eddy events. The model is characterized by comparison to experimental data of particle dispersion for a range of intrinsic particle time scales and body forces. Particle dispersion, velocity, and integral time scale results are presented. The particle implementation introduces a single model parameter β{sub p}, and sensitivity to this parameter and behavior of the model are discussed. Good agreement is found with experimental data and the ODT model is able to capture the particle inertial and trajectory crossing effects. These results serve as a validation case of the multiphase implementations of ODT for extensions to other flow configurations.

  14. A generalized Brownian motion model for turbulent relative particle dispersion

    NASA Astrophysics Data System (ADS)

    Shivamoggi, B. K.

    2016-08-01

    There is speculation that the difficulty in obtaining an extended range with Richardson-Obukhov scaling in both laboratory experiments and numerical simulations is due to the finiteness of the flow Reynolds number Re in these situations. In this paper, a generalized Brownian motion model has been applied to describe the relative particle dispersion problem in more realistic turbulent flows and to shed some light on this issue. The fluctuating pressure forces acting on a fluid particle are taken to be a colored noise and follow a stationary process and are described by the Uhlenbeck-Ornstein model while it appears plausible to take their correlation time to have a power-law dependence on Re, thus introducing a bridge between the Lagrangian quantities and the Eulerian parameters for this problem. This ansatz is in qualitative agreement with the possibility of a connection speculated earlier by Corrsin [26] between the white-noise representation for the fluctuating pressure forces and the large-Re assumption in the Kolmogorov [4] theory for the 3D fully developed turbulence (FDT) as well as a similar argument of Monin and Yaglom [23] and a similar result of Sawford [13] and Borgas and Sawford [24]. It also provides an insight into the result that the Richardson-Obukhov scaling holds only in the infinite-Re limit and disappears otherwise. This ansatz further provides a determination of the Richardson-Obukhov constant g as a function of Re, with an asymptotic constant value in the infinite-Re limit. It is shown to lead to full agreement, in the small-Re limit as well, with the Batchelor-Townsend [27] scaling for the rate of change of the mean square interparticle separation in 3D FDT, hence validating its soundness further.

  15. GRAVITATIONAL LENS MODELING WITH GENETIC ALGORITHMS AND PARTICLE SWARM OPTIMIZERS

    SciTech Connect

    Rogers, Adam; Fiege, Jason D.

    2011-02-01

    Strong gravitational lensing of an extended object is described by a mapping from source to image coordinates that is nonlinear and cannot generally be inverted analytically. Determining the structure of the source intensity distribution also requires a description of the blurring effect due to a point-spread function. This initial study uses an iterative gravitational lens modeling scheme based on the semilinear method to determine the linear parameters (source intensity profile) of a strongly lensed system. Our 'matrix-free' approach avoids construction of the lens and blurring operators while retaining the least-squares formulation of the problem. The parameters of an analytical lens model are found through nonlinear optimization by an advanced genetic algorithm (GA) and particle swarm optimizer (PSO). These global optimization routines are designed to explore the parameter space thoroughly, mapping model degeneracies in detail. We develop a novel method that determines the L-curve for each solution automatically, which represents the trade-off between the image {chi}{sup 2} and regularization effects, and allows an estimate of the optimally regularized solution for each lens parameter set. In the final step of the optimization procedure, the lens model with the lowest {chi}{sup 2} is used while the global optimizer solves for the source intensity distribution directly. This allows us to accurately determine the number of degrees of freedom in the problem to facilitate comparison between lens models and enforce positivity on the source profile. In practice, we find that the GA conducts a more thorough search of the parameter space than the PSO.

  16. Aspects of Particle Physics Beyond the Standard Model

    NASA Astrophysics Data System (ADS)

    Lu, Xiaochuan

    This dissertation describes a few aspects of particles beyond the Standard Model, with a focus on the remaining questions after the discovery of a Standard Model-like Higgs boson. In specific, three topics are discussed in sequence: neutrino mass and baryon asymmetry, naturalness problem of Higgs mass, and placing constraints on theoretical models from precision measurements. First, the consequence of the neutrino mass anarchy on cosmology is studied. Attentions are paid in particular to the total mass of neutrinos and baryon asymmetry through leptogenesis. With the assumption of independence among mass matrix entries in addition to the basis independence, Gaussian measure is the only choice. On top of Gaussian measure, a simple approximate U(1) flavor symmetry makes leptogenesis highly successful. Correlations between the baryon asymmetry and the light-neutrino quantities are investigated. Also discussed are possible implications of recently suggested large total mass of neutrinos by the SDSS/BOSS data. Second, the Higgs mass implies fine-tuning for minimal theories of weak-scale supersymmetry (SUSY). Non-decoupling effects can boost the Higgs mass when new states interact with the Higgs, but new sources of SUSY breaking that accompany such extensions threaten naturalness. I will show that two singlets with a Dirac mass can increase the Higgs mass while maintaining naturalness in the presence of large SUSY breaking in the singlet sector. The modified Higgs phenomenology of this scenario, termed "Dirac NMSSM", is also studied. Finally, the sensitivities of future precision measurements in probing physics beyond the Standard Model are studied. A practical three-step procedure is presented for using the Standard Model effective field theory (SM EFT) to connect ultraviolet (UV) models of new physics with weak scale precision observables. With this procedure, one can interpret precision measurements as constraints on the UV model concerned. A detailed explanation is

  17. Modeling of hydrogen production methods: Single particle model and kinetics assessment

    SciTech Connect

    Miller, R.S.; Bellan, J.

    1996-10-01

    The investigation carried out by the Jet Propulsion Laboratory (JPL) is devoted to the modeling of biomass pyrolysis reactors producing an oil vapor (tar) which is a precursor to hydrogen. This is an informal collaboration with NREL whereby JPL uses the experimentally-generated NREL data both as initial and boundary conditions for the calculations, and as a benchmark for model validation. The goal of this investigation is to find drivers of biomass fast-pyrolysis in the low temperature regime. The rationale is that experimental observations produce sparse discrete conditions for model validation, and that numerical simulations produced with a validated model are an economic way to find control parameters and an optimal operation regime, thereby circumventing costly changes in hardware and tests. During this first year of the investigation, a detailed mathematical model has been formulated for the temporal and spatial accurate modeling of solid-fluid reactions in biomass particles. These are porous particles for which volumetric reaction rate data is known a priori and both the porosity and the permeability of the particle are large enough to allow for continuous gas phase flow. The methodology has been applied to the pyrolysis of spherically symmetric biomass particles by considering previously published kinetics schemes for both cellulose and wood. The results show that models which neglect the thermal and species boundary layers exterior to the particle will generally over predict both the pyrolysis rates and experimentally obtainable tar yields. An evaluation of the simulation results through comparisons with experimental data indicates that while the cellulose kinetics is reasonably accurate, the wood pyrolysis kinetics is not accurate; particularly at high reactor temperatures. Current effort in collaboration with NREL is aimed at finding accurate wood kinetics.

  18. Supports and modified nano-particles for designing model catalysts.

    PubMed

    O'Brien, C P; Dostert, K-H; Hollerer, M; Stiehler, C; Calaza, F; Schauermann, S; Shaikhutdinov, S; Sterrer, M; Freund, H-J

    2016-07-01

    In order to design catalytic materials, we need to understand the essential causes for material properties resulting from its composite nature. In this paper we discuss two, at first sight, diverse aspects: (a) the effect of the oxide-metal interface on metal nanoparticle properties and (b) the consequences of metal particle modification after activation on the selectivity of hydrogenation reactions. However, these two aspects are intimately linked. The metal nanoparticle's electronic structure changes at the interface as a catalyst is brought to different reaction temperatures due to morphological modifications in the metal and, as we will discuss, these changes in the chemistry lead to changes in the reaction path. As the morphology of the particle varies, facets of different orientations and sizes are exposed, which may lead to a change in the surface chemistry as well. We use two specific reactions to address these issues in some detail. To the best of our knowledge, the present paper reports the first observations of this kind for well-defined model systems. The changes in the electronic structure of Au nanoparticles due to their size and interaction with a supporting oxide are revealed as a function of temperature using CO2 activation as a probe. The presence of spectator species (oxopropyl), formed during an activation step of acrolein hydrogenation, strongly controls the selectivity of the reaction towards hydrogenation of the unsaturated C[double bond, length as m-dash]O bond vs. the C[double bond, length as m-dash]C bond on Pd(111) when compared with oxide-supported Pd nanoparticles. PMID:27064816

  19. A cyclic model for particle motion in the pulmonary acinus

    NASA Astrophysics Data System (ADS)

    Haber, S.; Tsuda, A.

    2006-11-01

    A simplified model for the pulmonary alveolus that imitates the rhythmical expansion of the alveolus and the periodic shear flow in the adjacent airway is explored. The model consists of two eccentric cylinders and incompressible fluid that occupies the gap between them. The two cylinders undergo a simultaneous rhythmical expansion and contraction (mimicking the alveolus expansion) while the inner cylinder performs a periodic rotation about its axis (inducing shear flow mimicking airway ductal flow). An analytical solution is obtained for the creeping flow induced by the simultaneously expanding cylinders. It is shown that above a certain critical value of rotation to expansion velocity ratio, the flow exhibits characteristic features such as a saddle point and closed streamlines about a centre, similar to those existing inside a single alveolus during inhalation and exhalation. Poincaré maps of the trajectories of fluid particles demonstrate that, under various flow conditions, chaotic trajectories may exist, provided that expansion and rotation are slightly out of phase. This is similar to normal breathing conditions where the periodic expansion of the alveolus and the tidal flow (i.e. shear flow above the mouth of the alveolus) may be slightly out of phase. A novel definition of overall convective mixing efficiency is also suggested that inherently discounts reversible processes that do not contribute to mixing. It is demonstrated that two different convective mechanisms, related to the irreversibility of exhalation and inhalation and the onset of chaos, govern mixing efficiency in lung alveoli.

  20. Treating asphericity in fuel particle pressure vessel modeling

    NASA Astrophysics Data System (ADS)

    Miller, Gregory K.; Wadsworth, Derek C.

    1994-07-01

    The prototypical nuclear fuel of the New Production Modular High Temperature Gas-Cooled Reactor (NP-MHTGR) consists of spherical TRISO-coated particles suspended in graphite cylinders. The coating layers surrounding the fuel kernels in these particles consist of pyrolytic carbon layers and a silicon carbide layer. These coating layers act as a pressure vessel which retains fission product gases. In the operating conditions of the NP-MHTGR, a small percentage of these particles (pressure vessels) are expected to fail due to the pressure loading. The fuel particles of the NP-MHTGR deviate to some degree from a true spherical shape, which may have some effect on the failure percentages. A method is presented that treats the asphericity of the particles in predicting failure probabilities for particle samples. It utilizes a combination of finite element analysis and Monte Carlo sampling and is based on the Weibull statistical theory. The method is used here to assess the effects of asphericity in particles of two common geometric shapes, i.e. faceted particles and ellipsoidal particles. The method presented could be used to treat particle anomalies other than asphericity.

  1. Implicit frictional-contact model for soft particle systems

    NASA Astrophysics Data System (ADS)

    Nezamabadi, Saeid; Radjai, Farhang; Averseng, Julien; Delenne, Jean-Yves

    2015-10-01

    We introduce a novel numerical approach for the simulation of soft particles interacting via frictional contacts. This approach is based on an implicit formulation of the Material Point Method, allowing for large particle deformations, combined with the Contact Dynamics method for the treatment of unilateral frictional contacts between particles. This approach is both precise due to the treatment of contacts with no regularization and artificial damping parameters, and robust due to implicit time integration of both bulk degrees of freedom and relative contact velocities at the nodes representing the contact points. By construction, our algorithm is capable of handling arbitrary particle shapes and deformations. We illustrate this approach by two simple 2D examples: a Hertz contact and a rolling particle on an inclined plane. We also investigate the compaction of a packing of circular particles up to a solid fraction well above the jamming limit of hard particles. We find that, for the same level of deformation, the solid fraction in a packing of frictional particles is above that of a packing of frictionless particles as a result of larger particle shape change.

  2. Effect of Microstructural Geometry for Computing Closure Models in Multiscale Modeling of Shocked Particle Laden Flow

    NASA Astrophysics Data System (ADS)

    Sen, Oishik; Udaykumar, H. S.; Jacobs, Gustaaf

    Interaction of a shock wave with dust particles is a complex physical phenomenon. A computational model for studying this two-phase system is the Particle-Source in Cell (PSIC) approach. In this method, the dust particles are tracked as point particles in a Lagrangian frame of reference immersed in a compressible fluid. Two-way interaction between the carrier and the dispersed phases is ensured by coupling the momentum and energy transfer between the two phases as source terms in the respective governing equations. These source terms (e.g. drag force on particles) may be computed from resolved numerical simulations by treating each macroscopic point particle as an ensemble of cylinders immersed in a compressed fluid. However the drag so computed must be independent of the geometry of the mesoscale. In this work, the effect of the stochasticity of the microstructural geometry in construction of drag laws from resolved mesoscale computations is studied. Several different arrangement of cylinders are considered and the mean drag law as a function of Mach Number and Volume Fraction for each arrangement is computed using the Dynamic Kriging Method. The uncertainty in the drag forces arising because of the arrangement of the cylinders for a given volume fraction is quantified as 90% credible sets and the effect of the uncertainty on PSIC computations is studied.

  3. Particle models for discrete element modeling of bulk grain properties of wheat kernels

    Technology Transfer Automated Retrieval System (TEKTRAN)

    Recent research has shown the potential of discrete element method (DEM) in simulating grain flow in bulk handling systems. Research has also revealed that simulation of grain flow with DEM requires establishment of appropriate particle models for each grain type. This research completes the three-p...

  4. Coupled Particle Transport and Pattern Formation in a Nonlinear Leaky-Box Model

    NASA Technical Reports Server (NTRS)

    Barghouty, A. F.; El-Nemr, K. W.; Baird, J. K.

    2009-01-01

    Effects of particle-particle coupling on particle characteristics in nonlinear leaky-box type descriptions of the acceleration and transport of energetic particles in space plasmas are examined in the framework of a simple two-particle model based on the Fokker-Planck equation in momentum space. In this model, the two particles are assumed coupled via a common nonlinear source term. In analogy with a prototypical mathematical system of diffusion-driven instability, this work demonstrates that steady-state patterns with strong dependence on the magnetic turbulence but a rather weak one on the coupled particles attributes can emerge in solutions of a nonlinearly coupled leaky-box model. The insight gained from this simple model may be of wider use and significance to nonlinearly coupled leaky-box type descriptions in general.

  5. Critical Examination of the Colloidal Particle Model of Globular Proteins

    PubMed Central

    Sarangapani, Prasad S.; Hudson, Steven D.; Jones, Ronald L.; Douglas, Jack F.; Pathak, Jai A.

    2015-01-01

    Recent studies of globular protein solutions have uniformly adopted a colloidal view of proteins as particles, a perspective that neglects the polymeric primary structure of these biological macromolecules, their intrinsic flexibility, and their ability to sample a large configurational space. While the colloidal perspective often serves as a useful idealization in many cases, the macromolecular identity of proteins must reveal itself under thermodynamic conditions in which the native state is no longer stable, such as denaturing solvents and high protein concentrations where macromolecules tend to have screened excluded volume, charge, and hydrodynamic interactions. Under extreme pH conditions, charge repulsion interactions within the protein chain can overcome the attractive hydrogen-bonding interactions, holding it in its native globular state. Conformational changes can therefore be expected to have great significance on the shear viscosity and other rheological properties of protein solutions. These changes are not envisioned in conventional colloidal protein models and we have initiated an investigation of the scattering and rheological properties of model proteins. We initiate this effort by considering bovine serum albumin because it is a globular protein whose solution properties have also been extensively investigated as a function of pH, temperature, ionic strength, and concentration. As we anticipated, near-ultraviolet circular dichroism measurements and intrinsic viscosity measurements clearly indicate that the bovine serum albumin tertiary structure changes as protein concentration and pH are varied. Our findings point to limited validity of the colloidal protein model and to the need for further consideration and quantification of the effects of conformational changes on protein solution viscosity, protein association, and the phase behavior. Small-angle Neutron Scattering measurements have allowed us to assess how these conformational changes

  6. Modelling the surface deposition of meteoric smoke particles

    NASA Astrophysics Data System (ADS)

    Brooke, James S. A.; Feng, Wuhu; Mann, Graham W.; Dhomse, Sandip S.; Bardeen, Charles G.; Plane, John M. C.

    2016-04-01

    The flux of meteoric smoke particles (MSPs) in Greenland and Antarctica has been measured using Ir and Pt observations in ice cores, by Gabrielli et al. [1,2]. They obtained MSP deposition fluxes of 1.5 ± 0.45 × 10‑4 g m‑2 yr‑1 (209 ± 63 t d‑1) in Greenland and 3.9 ± 1.4 × 10‑5 g m‑2 yr‑1 (55 ± 19 t d‑1) in Antarctica, where the values in parentheses are total atmospheric inputs, assuming a uniform global deposition rate. These results show reasonable agreement with those of Lanci et al. [3], who used ice core magnetisation measurements, resulting in MSP fluxes of 1.7 ± 0.23 × 10‑4 g m‑2 yr‑1 (236 ± 50 t d‑1) (Greenland) and 2.0 ± 0.52 × 10‑5 g m‑2 yr‑1 (29 ± 5.0 t d‑1) (Antarctica). Atmospheric modelling studies have been performed to assess the transport and deposition of MSPs, using WACCM (Whole Atmosphere Community Climate Model), and the CARMA (Community Aerosol and Radiation Model) aerosol microphysics package. An MSP input function totalling 44 t d‑1 was added between about 80 and 105 km. Several model runs have been performed in which the aerosol scavenging by precipitation was varied. Wet deposition is expected (and calculated here) to be the main deposition process; however, rain and snow aerosol scavenging coefficients have uncertainties spanning up to two and three orders of magnitude, respectively [4]. The model experiments that we have carried out include simple adjustments of the scavenging coefficients, full inclusion of a parametrisation reported by Wang et al. [4], and a scheme based on aerosol removal where relative humidity > 100 %. The MSP fluxes obtained vary between 1.4 × 10‑5 and 2.6 × 10‑5 g m‑2 yr‑1 for Greenland, and 5.1 × 10‑6 and 1.7 × 10‑5 g m‑2 yr‑1 for Antarctica. These values are about an order of magnitude lower than the Greenland observations, but show reasonable agreement for Antarctica. The UM (Unified Model), UKCA (United Kingdom Chemistry and Aerosols Model

  7. Particle methods to solve modelling problems in wound healing and tumor growth

    NASA Astrophysics Data System (ADS)

    Vermolen, F. J.

    2015-12-01

    The paper deals with a compilation of several of our modelling studies on particle methods used for simulation of wound healing and tumor growth processes. The paper serves as an introduction of our modelling approaches to researchers with interest in biological cell-based models that use particle-based modelling approaches. The particles that we consider in the present models mimic either cells or points on cell boundaries that are allowed to migrate as a result of several chemical and mechanical factors. A distinct feature of our modelling frameworks with respect to conventional particle models, is that cells, mimicked by particles, are allowed to divide, differentiate and to die as a result of apoptosis or any causes for cell death. The paper is merely descriptive, rather than written in full mathematical rigor, and will show some of the potentials of the applied modelling.

  8. Microbial interactions lead to rapid micro-scale successions on model marine particles.

    PubMed

    Datta, Manoshi S; Sliwerska, Elzbieta; Gore, Jeff; Polz, Martin F; Cordero, Otto X

    2016-01-01

    In the ocean, organic particles harbour diverse bacterial communities, which collectively digest and recycle essential nutrients. Traits like motility and exo-enzyme production allow individual taxa to colonize and exploit particle resources, but it remains unclear how community dynamics emerge from these individual traits. Here we track the taxon and trait dynamics of bacteria attached to model marine particles and demonstrate that particle-attached communities undergo rapid, reproducible successions driven by ecological interactions. Motile, particle-degrading taxa are selected for during early successional stages. However, this selective pressure is later relaxed when secondary consumers invade, which are unable to use the particle resource but, instead, rely on carbon from primary degraders. This creates a trophic chain that shifts community metabolism away from the particle substrate. These results suggest that primary successions may shape particle-attached bacterial communities in the ocean and that rapid community-wide metabolic shifts could limit rates of marine particle degradation. PMID:27311813

  9. Simulation and Modeling of Explosive Dispersal of Compressed Gas-Particle Suspensions

    NASA Astrophysics Data System (ADS)

    Ling, Yue; Haselbacher, Andreas; Balachandar, S.

    2008-11-01

    Some effects of a detonating particle-laden explosive can be modeled through the explosion of a compressed gas-particle suspension. The suspension can be considered to be located in a spherical container and released through the instantaneous removal of the container similar to a shock tube. The transient flow generated by the removal of the container is well understood in the case without particles. Comparatively little is known about the dispersal of the particles and the unsteady interactions between the various waves and the particles. The objective of this work is to use Eulerian-Lagrangian simulations to improve our understanding of the dispersal of particles. Particular attention is focused on the interactions of the particles with the expansion fan and their effect on the flow behavior at large times. Parametric studies are carried out to systematically investigate the effects of particle size, particle density, and mass fraction.

  10. Microbial interactions lead to rapid micro-scale successions on model marine particles

    PubMed Central

    Datta, Manoshi S.; Sliwerska, Elzbieta; Gore, Jeff; Polz, Martin F.; Cordero, Otto X.

    2016-01-01

    In the ocean, organic particles harbour diverse bacterial communities, which collectively digest and recycle essential nutrients. Traits like motility and exo-enzyme production allow individual taxa to colonize and exploit particle resources, but it remains unclear how community dynamics emerge from these individual traits. Here we track the taxon and trait dynamics of bacteria attached to model marine particles and demonstrate that particle-attached communities undergo rapid, reproducible successions driven by ecological interactions. Motile, particle-degrading taxa are selected for during early successional stages. However, this selective pressure is later relaxed when secondary consumers invade, which are unable to use the particle resource but, instead, rely on carbon from primary degraders. This creates a trophic chain that shifts community metabolism away from the particle substrate. These results suggest that primary successions may shape particle-attached bacterial communities in the ocean and that rapid community-wide metabolic shifts could limit rates of marine particle degradation. PMID:27311813