SEARCH · Engineering Papers
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.
Compositions, methods, and articles of manufacture relating to a transparent/translucent insulating material formed by a hybrid silica gel composition made of porous silica particles joined with each other by silane crosslinking through a three dimensional silica network. In one embodiment, the hybrid silica gel includes small porous silica particles, and is made by a process that avoids energy intensive techniques.
Abstract. A key constraint of particle-based methods for modeling cloud microphysics is the conservation of total particle number, which is required for computational tractability. The process of collisional breakup poses a particular challenge to this framework, as breakup events often produce many droplet fragments of varying sizes, which would require creating new particles in the system. This work introduces a representation of collisional breakup in the so-called “superdroplet” method which conserves the total number of superdroplets in the system. This representation extends an existing stochastic collisional-coalescence scheme and samples from a fragment size distribution in an additional Monte Carlo step. This method is demonstrated in a set of idealized box model and single-column warm-rain simulations. We further discuss the effects of the breakup dynamic and fragment size distribution on the particle size distribution, hydrometeor population, and microphysical process rates. Box model experiments serve to characterize the impacts of properties such as coalescence efficiency and fragmentation function on the relative roles of collisional breakup and coalescence. The results demonstrate that this representation of collisional breakup can produce a stationary particle size distribution, in which breakup and coalescence rates are approximately equal, and that it recovers expected behavior such as a reduction in precipitate-sized particles in the column model. The breakup algorithm presented here contributes to an open-source pythonic implementation of the superdroplet method, PySDM, which will facilitate future research using particle-based microphysics.
A method for evaluating a target, the target having a surface, includes pulsing a defined, energetic particle beam through the surface and into the target such that particle energy deposition from the particle beam is concentrated in a subsurface target volume within a target medium of the target. The deposited particle energy induces a thermoelastic expansion of the target medium in the target volume that generates a corresponding acoustic wave. The method further includes detecting the acoustic wave from the target medium.
Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.
Recent research showed that fracture of sand particles plays a significant role in determining the plastic bulk volumetric changes of granular materials under different loading conditions. One of the major tools used to better understand the influence of particle fracture on the behavior of granular materials is discrete-element modeling (DEM). This paper employed the bonded block model (BBM) to simulate the fracture behavior of sand. Each sand particle is modeled as an agglomerate of rigid blocks bonded at their contacts using the linear-parallel contact model, which can transmit both moment and force. DEM simulated particles closely matched the actual three-dimensional (3D) shape of sand particles acquired using high-resolution 3D synchrotron microcomputed tomography (SMT). Results from unconfined one-dimensional (1D) compression of a single synthetic silica cube were used to calibrate the model parameters. Particle fracture was investigated for specimens composed of three sand particles that were loaded under confined 1D compression. Breakage energy measured from DEM models matched well with that measured experimentally. The paper studied the effects of contact loading condition and particle interaction on the fracture mode of particles using BBM that can closely capture the 3D shape of real sand particles.
In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.
A parallel, relativistic, three-dimensional particle-in-cell code SPACE has been developed for the simulation of electromagnetic fields, relativistic particle beams, and plasmas. In addition to the standard second-order Particle-in-Cell (PIC) algorithm, SPACE includes efficient novel algorithms to resolve atomic physics processes such as multi-level ionization of plasma atoms, recombination, and electron attachment to dopants in dense neutral gases. SPACE also contains a highly adaptive particle-based method, called Adaptive Particle-in-Cloud (AP-Cloud), for solving the Vlasov-Poisson problems. It eliminates the traditional Cartesian mesh of PIC and replaces it with an adaptive octree data structure. The code's algorithms, structure, capabilities, parallelization strategy, and performance have been discussed. Additionally, typical examples of SPACE applications to accelerator science and engineering problems are described.
The identification of prompt and isolated muons, as well asmuons from heavy-flavour hadron decays, is an important task. Wedeveloped two multivariate techniques to provide highly efficientidentification for muons with transverse momentum greater than10 GeV. One provides a continuous variable as an alternative to acut-based identification selection and offers a betterdiscrimination power against misidentified muons. The other oneselects prompt and isolated muons by using isolation requirements toreduce the contamination from nonprompt muons arising inheavy-flavour hadron decays. Both algorithms are developed using59.7 fb$^{-1}$ of proton-proton collisions data at a centre-of-massenergy of √(s)=13 TeV collected in 2018 with the CMSexperiment at the CERN LHC.
Neutron invariant-mass spectroscopy (NIMS) is a key technique for studying unbound and weakly bound nuclei at the limits of stability. At the Facility for Rare Isotope Beams (FRIB), such measurements are performed using the Sweeper spectrometer, a large-gap, high-rigidity dipole system coupled to the MoNA-LISA neutron detector arrays. To meet the demands imposed by higher beam energies (>130 MeV/u) and the broad cocktail-beam selection available at FRIB, the spectrometer has recently been upgraded to improve particle-identification and detection performance. Upstream of the reaction target, a plastic scintillator with Silicon photomultiplier (SiPM) readout provides the global trigger and time reference, two parallel plate avalanche counters (PPACs) track the trajectories of incoming beam particles, and a silicon PIN detector measures the energy loss, ΔE, for charge (Z) identification. After the Sweeper magnet, the trajectories of the reaction products are tracked by two micro-pattern drift chambers (MPDCs), their charge (Z) is identified by a Frisch-grid ionization chamber (FG-IC), and their mass-to-charge ratio (A/Q) is deduced by time-of-flight measurement using a fast plastic scintillator read out by an array of photomultiplier tubes (PMTs). The detection system also incorporates the Modular Neutron Array (MoNA) for neutron detection and the CAESium-iodide scintillator ARray (CAESAR) for high-efficiency γ-ray measurements to enable full kinematic reconstruction. Performance was evaluated using a cocktail beam around 37 Al accelerated at E ≈ 130 MeV/u during the first FRIB campaign, demonstrating the readiness of the upgraded system for future studies of nuclei at and beyond the neutron drip line.
The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.
Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.
A method of reducing electrophoretic mobility of macromolecules, particles, cells, and other substances is provided which comprises interacting in a conventional electrophoretic separating procedure, the substances with a polymer-linked affinity compound comprised of a hydrophilic neutral polymer such as polyethylene glycol bound to a second component such as a hydrophobic compound, an immunocompound such as an antibody or antibody active fragment, or a ligand such as a hormone, drug, antigen, or a hapten. The reduction of electrophoretic mobility achieved is directly proportional to the concentration of the polymer-linked affinity compound employed, and such reduction can comprise up to 100 percent for particular particles and cells. The present invention is advantageous in that electrophoretic separation can now be achieved for substances whose native surface charge structure had prevented them from being separated by normal electrophoretic means. Depending on the affinity component utilized, separation can be achieved on the basis of the specific/irreversible, specific/reversible, semi-specific/reversible, relatively nonspecific/reversible, or relatively nonspecific/irreversible ligand-substance interactions.
Studies of the electrostatic interaction between micrometer-sized particles and polymer surfaces are of great interest to NASA's planetary exploration program. The unmanned landing missions to Mars planned for this decade as well as the possible manned missions that might take place during the second decade of this century require a better understanding of the electrostatic response of the materials used in landing crafts and equipment when exposed to wind-blown dust or to surface dust and sand particles. We report on preliminary experiments designed to measure the electrostatic charge developed on five polymer surfaces as they are impacted simultaneously by Mars simulant particles less than 5 micrometers in diameter moving at 20 m/s. Experiments were performed in a CO2 atmosphere at 10 mbars of pressure using a particle delivery method that propels the particles with contact. Experiments were also performed in dry air at atmospheric pressures using a pressurized particle delivery system. The five polymer surfaces, commonly used in space applications, were chosen so that they span the triboelectric series.
Methods for measuring electrostatic potentials on and near dielectric surfaces charged to several kilovolts are studied. Secondary emission from those charged dielectrics is measured. Candidates for potential measurement include the induced charge, from which potential is calculated; the trajectory endpoints of either high or low energy particles traversing the region near the surface; trajectory impact on the surface; and creating ions at points of interest near the surface. Some of the methods require computer simulations and iterative calculation if potential maps are to be generated. Several approaches are described and compared. A method using a half-cylinder as a test chamber and low-energy probing beams is adapted for the measurement of seconary emission.
Many problems in materials science and biology involve particles interacting with strong, short-ranged bonds that can break and form on experimental timescales. Treating such bonds as constraints can significantly speed up sampling their equilibrium distribution, and there are several methods to sample probability distributions subject to fixed constraints. We introduce a Monte Carlo method to handle the case when constraints can break and form. More generally, the method samples a probability distribution on a stratification: a collection of manifolds of different dimensions, where the lower-dimensional manifolds lie on the boundaries of the higher-dimensional manifolds. We show several applications of the method in polymer physics, self-assembly of colloids, and volume calculation in high dimensions.
An apparatus and method for destructively removing particles from a flowing gas containing the particles is described. In the specific embodiments disclosed the apparatus is adapted to remove carbon particles from diesel engine exhaust products. The exhaust products are directed to a predetermined location where they are rapidly vaporized and combine with oxygen in the exhaust products to form carbon dioxide. Vaporization in one embodiment is effected by a discharge grid located within an exhaust conduit, the grid being chosen so that alternate conductors defining the grid are spaced apart a distance approximately 125 times the mean diameter of the particles to be removed. A voltage differential of approximately 690 volts is applied across adjacent conductors.
Provided is a particle that includes a first porous region and a second porous region that differs from the first porous region. Also provided is a particle that has a wet etched porous region and that does have a nucleation layer associated with wet etching. Methods of making porous particles are also provided.
Methods and apparatus for measuring gravitational and inertial forces, magnetic fields, or wave or radiant energy acting on an object or fluid in space provide an electric tunneling current through a gap between an electrode and that object or fluid in space and vary that gap with any selected one of such forces, magnetic fields, or wave or radiant energy acting on that object or fluid. These methods and apparatus sense a corresponding variation in an electric property of that gap and determine the latter force, magnetic fields, or wave or radiant energy in response to that corresponding variation, and thereby sense or measure such parameters as acceleration, position, particle mass, velocity, magnetic field strength, presence or direction, or wave or radiant energy intensity, presence or direction.