Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Modeling and implementation of vertical excursion FFA in the Zgoubi ray-tracing code

Vertical Fixed Field Accelerators (vFFAs) feature complex and highly non-linear magnetic fields, which require simulation codes allowing step-wise particle tracking. Methods to model the 3D magnetic field of scaling vFFAs have been developed in the ray-tracing code Zgoubi. The field modeling and particle tracking methods include the field non-linearities, the fringe fields, and the field superposition of neighboring magnets. The procedure implements the vFFA analytical field expressions, allowing design studies and parameter optimizations using the Zgoubi built-in fit method. The vFFA procedure has been applied to a ten-fold symmetry ring with a triplet focusing structure designed to accelerate protons from 3 MeV to 12 MeV and studied under the ISIS-II proton driver prototype project. Results from particle tracking in externally generated 3D semi-analytical field maps and the developed vFFA analytical model are shown to be in excellent agreement.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Parallelized domain decomposition for multi-dimensional Lagrangian random walk mass-transfer particle tracking schemes

Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work. For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Sequential ensemble transform for Bayesian inverse problems

In this work, we present the Sequential Ensemble Transform (SET) method, an approach for generating approximate samples from a Bayesian posterior distribution. The method explores the posterior distribution by solving a sequence of discrete optimal transport problems to produce a series of transport plans which map prior samples to posterior samples. We prove that the sequence of Dirac mixture distributions produced by the SET method converges weakly to the true posterior as the sample size approaches infinity. Furthermore, our numerical results indicate that, when compared to standard Sequential Monte Carlo (SMC) methods, the SET approach is more robust to the choice of Markov mutation kernels and requires less computational efforts to reach a similar accuracy when used to explore complex posterior distributions. Finally, we describe adaptive schemes that allow to completely automate the use of the SET method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Compound pulse characteristics of a heterogeneous composite scintillator in a gamma-ray field

An important aspect of radiation detection in scintillators is the time dependence of emitted fluorescence, which results in characteristic measured pulse shape. In some materials, a relatively strong dependence of the scintillation pulse shape on stopping power exists, which can provide reliable identification of the incident particle. An alternative method to realize particle identification is by combining materials with different fluorescence lifetimes into heterogeneous structures. Such composite scintillators derive their properties from both the constituent material characteristics and the geometry of the composite structure. In composite scintillators with a high degree of heterogeneous loading, a superposition of fluorescence contributions originating in multiple materials, which is observed as single waveform with an intermediate pulse shape, may affect its ability to reject rays. We measure the time dependence of light output from a scintillator composed of 6 Li-containing glass shards and scintillating polyvinyl toluene and identify events that exhibit such compound behavior when exposed to rays and fast neutrons. We develop a modeling and simulation framework that reproduces the pulse shapes in heterogeneous scintillators and use it study the effect that the weight percentage of 6 Li glass has on rejection. The modeling framework is applied to the experimentally studied scintillator, finding a good agreement. Here, the developed modeling and simulation approach will help optimize the design of heterogeneous scintillators to meet the desired trade-o between neutron capture efficiency and -ray rejection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Electron Modulational Instability in the Strong Turbulent Regime for an Electron Beam Propagating in a Background Plasma

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. In conclusion, different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

Comparative studies of grain size separates of 60009

Five samples from 60009, the lower half of a double drive tube, were analyzed via grain-size methods, with particle types classified and counted in the coarser grain sizes. Studies were undertaken of particle types and distributions by petrographic methods, of magnetic fractions, of the size splits and magnetic splits as analyzed by ferromagnetic resonance (FMR) techniques, of maturity (based on agglutinate content, FMR index Is/FeO, mean size of sub-cm material, magnetic fraction), of possible reworking or mixing in situ, and of depositional history. Maturity indices are in substantial agreement for all of the five samples. Strong positive correlation of percent agglutinates and percent bedrock-derived lithic fragments, combined with negative correlation of those components with percent single crystal plagioclase, argue against in situ reworking of the same soil.

Mckay, D. S.↗
DOE OSTI2021

Multi-stage falling particle receivers

The present disclosure is directed to multi-stage falling particle receivers and methods of falling particle heating. As the particles fall through the receiver, the particles are periodically collected and released by flow retarding devices. The periodic catch-and-release of the particles falling through the receiver reduces particle flow dispersion, increases particle opacity and solar absorption, and reduces erosion and damage to surfaces caused by direct particle impingement.

Ho, Clifford K.↗
DOE OSTI2023

Multi-stage falling particle receivers

The present disclosure is directed to multi-stage falling particle receivers and methods of falling particle heating. As the particles fall through the receiver, the particles are periodically collected and released by flow retarding devices. The periodic catch-and-release of the particles falling through the receiver reduces particle flow dispersion, increases particle opacity and solar absorption, and reduces erosion and damage to surfaces caused by direct particle impingement.

Ho, Clifford K.↗
NASA NTRS1992

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie↗
DOE OSTI2024

Development and Experimental Optimization of High-Temperature Modeling Tools and Methods for Concentrated Solar Power Particle - Systems

A novel, open-source radiative modeling toolset was developed to extend the functionality of particle-based modeling software (e.g. discrete element method (DEM)) to environmental conditions relevant to concentrated solar power applications. This toolset was optimized for deployment on desktop workstations instead of high-performance computing systems, to render such tools more accessible to the research community. Both particle-based modeling and radiative exchange modeling are computationally expensive and often require specialized programming expertise, making these methods cumbersome to use. Recent developments in DEM software by DCS Computing have greatly reduced these challenges, providing a graphical-user-interface based platform and modeling optimization for desktop workstations, HPCs, and cloud computing. The University of Dayton leveraged the experience of DCS Computing in developing a user-friendly, open-source radiative heat transfer expansion for DEM modeling. The University of Dayton DEM+ radiative modeling toolset was developed using a combination of fundamental experimental measurements, modeling, and simplified flow experiments over a range of temperatures and flow conditions. The toolset provides researchers with access to multiple radiative models including an accelerated Monte-Carlo Ray Tracing (application agnostic, highly computationally expensive), an expanded database of distance-based approximations (application limited, computationally light), and a weighted blending of the two methods capable of achieving over 90% reduction in computation time with equivalent accuracy compared to Monte-Carlo Ray Tracing. Through a graphical user interface, users can customize the radiative models to match their desired accuracy and available computational resources, improving access to particle based modeling for the research community. Ceramic sintered bauxite proppants were used in modeling and experimentally as a baseline. Both the radiative heat transfer and flow properties for particulate systems were investigated at elevated temperatures up to 800 °C. The major accomplishments for this work include a verified, open-source radiative modeling toolset to be distributed amongst the research community and the fabrication of three small-scale test facilities to investigate particle behavior and tune DEM flow properties for operation up to 800 °C. The findings have been shared with the research community via conference modeling workshops, deployment of the tools in DCS Computing Aspherix®, and open-source access to the developed radiative modeling tool. The development of next-generation CSP facilities and thermal energy storage systems based on ceramic particles requires providing access to computationally efficient and accurate modeling tools. Particles will experience a wide range of environments (20-800 °C) and handling conditions (dilute curtains or dense packing), requiring specially designed and optimized equipment. Optimizing solid particle physics models and establishing best-practices for particle modeling in CSP environments will assist researchers with designing optimized equipment, accelerating the deployment of more economically-competitive CSP facilities.

14 SOLAR ENERGY↗
DOE OSTI2022

Characterizing particle-based thermal storage performance using optical methods for use in next generation concentrating solar power plants

Concentrating Solar Power (CSP) generation is an attractive option for low-emission power generation; however, the high costs of thermal storage associated with concentrating solar create a large barrier for their use and adaptation into modern life. Lowering their operation costs, while maintaining high thermal storage and transfer performance is essential. Solid particle-based heat exchange systems can reduce CSP cost but are often less efficient. Efforts to increase their performance have led to use of binary size particle mixes. Presented is an optical-based thermal analysis technique used to measure near-wall thermal conductivity of particle beds essential in determining their heat exchanger efficiency. Modulated Photothermal Radiometry is used to make dynamic temperature measurements, allowing for the extraction of the most relevant thermal properties like thermal conductivity, specific heat, and effusivity. The system uses a modulated laser source causing a damped periodic heat flux, resulting in a frequency and thermal property dependent surface temperature, of which is measured using radiometry. Lock-In techniques are used to extrapolate the amplitude of the signal. Plotting the amplitude against the root angular frequency allows for effusivity measurement by ratio to a known sample. Using specific heat measurements from literature and density measurements, the thermal conductivity of the particle mixes can be calculated. The simplicity of MPTR to probe through the depth of the bed is ideal for use in CSP for dynamic thermal performance monitoring.

Corona, Javier↗
NASA NTRS1975

A discharge-condensation method for the detection of traces of charged particles

A new discharge condensation method is suggested for detecting the traces of charged particles. The chamber designed on the basis of this method possesses a high time resolution (several microseconds) and a greater memory time (dozens of milliseconds). Photographed traces of cosmic particles in the discharge condensation chamber are given.

Mandzhavidze, Z. S.↗
NASA NTRS2012

Efficient Implementation of the Invariant Imbedding T-Matrix Method and the Separation of Variables Method Applied to Large Nonspherical Inhomogeneous Particles

Three terms, ''Waterman's T-matrix method'', ''extended boundary condition method (EBCM)'', and ''null field method'', have been interchangeable in the literature to indicate a method based on surface integral equations to calculate the T-matrix. Unlike the previous method, the invariant imbedding method (IIM) calculates the T-matrix by the use of a volume integral equation. In addition, the standard separation of variables method (SOV) can be applied to compute the T-matrix of a sphere centered at the origin of the coordinate system and having a maximal radius such that the sphere remains inscribed within a nonspherical particle. This study explores the feasibility of a numerical combination of the IIM and the SOV, hereafter referred to as the IIMþSOV method, for computing the single-scattering properties of nonspherical dielectric particles, which are, in general, inhomogeneous. The IIMþSOV method is shown to be capable of solving light-scattering problems for large nonspherical particles where the standard EBCM fails to converge. The IIMþSOV method is flexible and applicable to inhomogeneous particles and aggregated nonspherical particles (overlapped circumscribed spheres) representing a challenge to the standard superposition T-matrix method. The IIMþSOV computational program, developed in this study, is validated against EBCM simulated spheroid and cylinder cases with excellent numerical agreement (up to four decimal places). In addition, solutions for cylinders with large aspect ratios, inhomogeneous particles, and two-particle systems are compared with results from discrete dipole approximation (DDA) computations, and comparisons with the improved geometric-optics method (IGOM) are found to be quite encouraging.

Invariant imbedding method↗
NASA NTRS1979

Computational methods for analyzing the transmission characteristics of a beta particle magnetic analysis system

Computational methods were developed to study the trajectories of beta particles (positrons) through a magnetic analysis system as a function of the spatial distribution of the radionuclides in the beta source, size and shape of the source collimator, and the strength of the analyzer magnetic field. On the basis of these methods, the particle flux, their energy spectrum, and source-to-target transit times have been calculated for Na-22 positrons as a function of the analyzer magnetic field and the size and location of the target. These data are in studies requiring parallel beams of positrons of uniform energy such as measurement of the moisture distribution in composite materials. Computer programs for obtaining various trajectories are included.

Singh, J. J.↗
NASA NTRS1991

Light scattering by nonspherical particles - A refinement to the coupled-dipole method

In the coupled-dipole method, an arbitrary particle is modeled as an array of N polarizable subunits, each of which gives rise to only electric dipole radiation. The Clausius-Mosotti relation is widely used to calculate the polarizability of the subunits that correspond to the dielectric function of the particle that the array represents. The Clausius-Mosotti relation is replaced here with an exact expression for the electric dipole polarizability, and improvement is found in extinction calculations for spheres as compared with Mie theory. Near a Froehlich frequency the coupled-dipole method yields extinction cross sections for spheres and spheroids that compare favorably with the continuous distribution of ellipsoids method and measured values.

Dungey, Clifton E.↗
DOE OSTI2021

Enabling particle applications for exascale computing platforms

The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four “sub-motifs”: short-range particle–particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle–particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or “apps”) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable “lessons learned” that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA’s ECP application partner codes are also addressed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Energy storage material and method of producing the same

The instant invention includes a spherical porous secondary silicon-based particle and methods for producing the same. The spherical porous secondary silicon-based particle is comprised of agglomerated primary silicon-based nanoparticles. The secondary particle comprises a carbon coating that reduces the effective exposed surface area of the primary particles to the electrolyte, thus improving first cycle efficiency. The secondary particle further comprises porous regions that enable the silicon nanoparticles to expand during lithiation. Advantages include ease of castability with micron-sized spherical particles, ease of mixing spherical particles, ease of flow for spherical particles in various processing steps, and ease with obtaining higher loading, which translates to higher areal capacity and overall energy density of the cell. A readily scalable process for producing the particles using low-cost materials and low-cost processing methods is disclosed.

Holt, Christopher T.↗