Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle Methods”

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 793 records · Page 44

Spaceflight Holography Investigation in a Virtual Apparatus

The movement of a particle in a fluid is one of the most fundamental processes in physics and plays an important role in materials science. For example, the physics of crystal growth from a nucleation site in solution can be dominated by movement of the solution over the crystal surface. Gravity usually dominates the equations of motion, but in microgravity other terms can dominate, making the equation much more complex. Until recently, the equations were solved only by numerical methods and/or by neglecting terms. During this study, we discovered an exact solution to the equations, which shows that the usually neglected terms become extremely important in microgravity. We also developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the virtual apparatus . We will record holograms of particle distributions in motion in microgravity and bring them back to earth for analysis, allowing the study of the full three-dimensional motion of sets of particles, allowing us to test the new analytical solutions. The experiment will also provide accurate measurements of the quasi-steady acceleration of the space platform and other interesting microgravity and g-jitter effects. This program will produce the flight definition for an experiment in the microgravity environment of space to validate the theoretical model. We will design an experiment with the help of the theoretical model that is optimized for testing the model, measuring g, g-jitter, and other microgravity phenomena.

Source record↗

Swept wing ice accretion modeling

An effort to develop a three-dimensional modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flow fields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane normal to the wing leading edge and compared to 2-D results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that the flow field over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3-D calculation.

Potapczuk, Mark G.↗

A Particle Batch Smoother Approach to Snow Water Equivalent Estimation

This paper presents a newly proposed data assimilation method for historical snow water equivalent SWE estimation using remotely sensed fractional snow-covered area fSCA. The newly proposed approach consists of a particle batch smoother (PBS), which is compared to a previously applied Kalman-based ensemble batch smoother (EnBS) approach. The methods were applied over the 27-yr Landsat 5 record at snow pillow and snow course in situ verification sites in the American River basin in the Sierra Nevada (United States). This basin is more densely vegetated and thus more challenging for SWE estimation than the previous applications of the EnBS. Both data assimilation methods provided significant improvement over the prior (modeling only) estimates, with both able to significantly reduce prior SWE biases. The prior RMSE values at the snow pillow and snow course sites were reduced by 68%-82% and 60%-68%, respectively, when applying the data assimilation methods. This result is encouraging for a basin like the American where the moderate to high forest cover will necessarily obscure more of the snow-covered ground surface than in previously examined, less-vegetated basins. The PBS generally outperformed the EnBS: for snow pillows the PBSRMSE was approx.54%of that seen in the EnBS, while for snow courses the PBSRMSE was approx.79%of the EnBS. Sensitivity tests show relative insensitivity for both the PBS and EnBS results to ensemble size and fSCA measurement error, but a higher sensitivity for the EnBS to the mean prior precipitation input, especially in the case where significant prior biases exist.

EnBS↗

Does Orbital Angular Momentum Have Effect on Laser’s Scattering by Molecular Atmosphere?

Lasers with orbital angular momentum (OAM) have potential applications in communication technology, manipulation of particles, and remote sensing. Because of its unusual light-scattering properties, the OAM laser’s interaction with a molecular atmosphere must be studied to ensure that it is not lossy for communication or remote-sensing applications that involve its transmission through an atmospheric environment. In this study, the finite-difference time-domain (FDTD) method [21] is applied to calculate the light scattering of the purely azimuthal (the radial mode number is assumed to be zero) Laguerre-Gaussian (LG) beams with OAM by very small dielectric particles. Not like Lorentz-Mie solutions, the FDTD method can calculate for particles off the central axis of the LG beam. It is found that when the particles are very small, and the topological charge number of the OAM of a laser is not extremely large, the laser’s OAM has little effect on the scattering phase function. This suggests that Rayleigh theory can be applied directly to calculate the light scattering by atmospheric molecules. The transmission of a laser beam with OAM in a molecular atmosphere is not different from that of a regular Gaussian beam.

Molecular atmosphere↗

Light scattering by size-shape distributions of randomly oriented axially symmetric particles of a size comparable to a wavelength

Rigorously light scattering by size-shape distributions of randomly oriented axially symmetric particles are calculated by the T-matrix method, as extended to randomly oriented scatterers. The computational scheme is described along with a newly developed convergence procedure that makes it possible to substantially reduce computer time and storage requirements. The elements of the Stokes scattering matrix for a power-law size distribution of randomly oriented moderately aspherical spheroids are shown to be much smoother than and differ substantially from those of equivalent monodisperse spheroids; averaging over orientations does not eliminate the necessity of averaging over particle sizes. The angular-scattering behavior of the ensembles of nonspherical particles is found to be significantly different from that of the equivalent polydisperse spheres.

Mishchenko, Michael I.↗

Fuzzy logic particle tracking velocimetry

Fuzzy logic has proven to be a simple and robust method for process control. Instead of requiring a complex model of the system, a user defined rule base is used to control the process. In this paper the principles of fuzzy logic control are applied to Particle Tracking Velocimetry (PTV). Two frames of digitally recorded, single exposure particle imagery are used as input. The fuzzy processor uses the local particle displacement information to determine the correct particle tracks. Fuzzy PTV is an improvement over traditional PTV techniques which typically require a sequence (greater than 2) of image frames for accurately tracking particles. The fuzzy processor executes in software on a PC without the use of specialized array or fuzzy logic processors. A pair of sample input images with roughly 300 particle images each, results in more than 200 velocity vectors in under 8 seconds of processing time.

Wernet, Mark P.↗

A Generalized Dual-Frequency Ratio (DFR) Approach for Rain Retrievals

The dual-frequency ratio of radar reflectivity factors (DFR) has been shown to be a useful quantity as it is independent of the number concentration of the particle size distribution and primarily a function of the mass-weighted particle diameter, Dm. A drawback of DFR-related methods for rain estimation, however, is the non-unique relationship between Dm and DFR. At Ku- and Ka-band frequencies, two solutions for Dm exist when DFR is less than zero. This ambiguity generates multiple solutions for the range profiles of the particle size parameters. We investigate characteristics of these solutions for both the initial-value (forward) and final-value (backward) forms of the equations. To choose one of among many possible range profiles of Dm, number concentration and rain rate, R, independently measured path attenuations are used. For the backward approach the possibility exists of dispensing with externally measured path attenuations by achieving consistency between the input and output path attenuations. The methods are tested by means of a simulation based on disdrometer-measured raindrop size distributions and the results are compared with a simplified version of the operational R-Dm method.

Robert Meneghini↗

Particle displacement tracking applied to air flows

Electronic Particle Image Velocimeter (PIV) techniques offer many advantages over conventional photographic PIV methods such as fast turn around times and simplified data reduction. A new all electronic PIV technique was developed which can measure high speed gas velocities. The Particle Displacement Tracking (PDT) technique employs a single cw laser, small seed particles (1 micron), and a single intensified, gated CCD array frame camera to provide a simple and fast method of obtaining two-dimensional velocity vector maps with unambiguous direction determination. Use of a single CCD camera eliminates registration difficulties encountered when multiple cameras are used to obtain velocity magnitude and direction information. An 80386 PC equipped with a large memory buffer frame-grabber board provides all of the data acquisition and data reduction operations. No array processors of other numerical processing hardware are required. Full video resolution (640x480 pixel) is maintained in the acquired images, providing high resolution video frames of the recorded particle images. The time between data acquisition to display of the velocity vector map is less than 40 sec. The new electronic PDT technique is demonstrated on an air nozzle flow with velocities less than 150 m/s.

Wernet, Mark P.↗

Particle displacement tracking applied to air flows

Electronic Particle Image Velocimetric (PIV) techniques offer many advantages over conventional photographic PIV methods such as fast turn around times and simplified data reduction. A new all electronic PIV technique was developed which can measure high speed gas velocities. The Particle Displacement Tracking (PDT) technique employs a single CW laser, small seed particles (1 micron), and a single intensified, gated CCD array frame camera to provide a simple and fast method of obtaining two-dimensional velocity vector maps with unambiguous direction determination. Use of a single CCD camera eliminates registration difficulties encountered when multiple cameras are used to obtain velocity magnitude and direction information. An 80386 PC equipped with a large memory buffer frame-grabber board provides all of the data acquisition and data reduction operations. No array processors of other numerical processing hardware are required. Full video resolution (640 x 480 pixel) is maintained in the acquired images, providing high resolution video frames of the recorded particle images. The time between data acquisition to display of the velocity vector map is less than 40 sec. The new electronic PDT technique is demonstrated on an air nozzle flow with velocities less than 150 m/s.

Wernet, Mark P.↗

An Efficient and Accurate Algorithm for Computing Grid-Averaged Solar Fluxes for Horizontally Inhomogeneous Clouds

A computationally efficient method is presented to account for the horizontal cloud inhomogeneity by using a radiatively equivalent plane parallel homogeneous (PPH) cloud. The algorithm can accurately match the calculations of the reference (rPPH) independent column approximation (ICA) results, but use only the same computational time required for a single plane parallel computation. The effective optical depth of this synthetic sPPH cloud is derived by exactly matching the direct transmission to that of the inhomogeneous ICA cloud. The ffective9 scattering asymmetry factor is found from a pre-calculated albedo inverse look-up-table that is allowed to vary over the range from -1.0 to 1.0. In the special cases of conservative scattering and total absorption, the synthetic method is exactly equivalent to the ICA, with only a small bias (about 0.2% in flux) relative to ICA due to imperfect interpolation in using the look-up tables. In principle, the ICA albedo can be approximated accurately regardless of cloud inhomogeneity. For a more complete comparison, the broadband shortwave albedo and transmission calculated from the synthetic sPPH cloud and averaged over all incident directions, have the RMS biases of 0.26% and 0.76%, respectively, for inhomogeneous clouds over a wide variation of particle size. The advantages of the synthetic PPH method are that (1) it is not required that all the cloud subcolumns have uniform microphysical characteristic, (2) it is applicable to any 1D radiative transfer scheme, and (3) it can handle arbitrary cloud optical depth distributions and an arbitrary number of cloud subcolumns with uniform computational efficiency.

cloud inhomogeneity↗

Vibration-translation energy transfer in anharmonic diatomic molecules. 1: A critical evaluation of the semiclassical approximation

The semiclassical approximation is applied to anharmonic diatomic oscillators in excited initial states. Multistate numerical solutions giving the vibrational transition probabilities for collinear collisions with an inert atom are compared with equivalent, exact quantum-mechanical calculations. Several symmetrization methods are shown to correlate accurately the predictions of both theories for all initial states, transitions, and molecular types tested, but only if coupling of the oscillator motion and the classical trajectory of the incident particle is considered. In anharmonic heteronuclear molecules, the customary semiclassical method of computing the classical trajectory independently leads to transition probabilities with anomalous low-energy resonances. Proper accounting of the effects of oscillator compression and recoil on the incident particle trajectory removes the anomalies and restores the applicability of the semiclassical approximation.

Mckenzie, R. L.↗

Vibration-translation energy transfer in anharmonic diatomic molecules. I - A comparative evaluation of the semiclassical approximation

The semiclassical approximation (quantum oscillator, classical path) is applied to anharmonic diatomic oscillators in excited initial states. Multistate numerical solutions giving the vibrational transition probabilities for collinear collisions with an inert atom are compared with equivalent, exact quantum-mechanical calculations. Several symmetrization methods are shown to correlate accurately the predictions of both theories for all initial states, transitions, and molecular types tested, but only if coupling of the oscillator motion and the classical trajectory of the incident particle is considered. In anharmonic heteronuclear molecules, the customary semiclassical method of computing the classical trajectory independently leads to transition probabilities with anomalous low-energy resonances. Proper accounting of the effects of oscillator compression and recoil on the incident particle trajectory removes the anomalies and restores the applicability of the semiclassical approximation.

Mckenzie, R. L.↗

Double Power-law Formation by Sequential Particle Acceleration

Spectral double power laws are common in solar high-energy phenomena such as flares and interplanetary energetic-electron events. However, the physical mechanism that produces the changes in power-law index within a single spectrum is unclear. We developed a fully analytical method of forming single power-law spectra from sequential acceleration of particles orbiting inside and hopping between simulated large-scale magnetic islands formed by flare reconnection. Here, we extend the analytical method to the formation of double power-law spectra by assuming sequential acceleration in two successive regions with different acceleration and particle-transport rates. The resulting spectral distribution is continuous and smooth, with a flattening at low energies, two power-law regions at mid-energies, and a steep rollover at high energies. The model provides analytical expressions for the spectral indices, all energy breaks, and normalization constants as functions of just three physical parameters of each acceleration region: (1) the energy gain in each accelerator, (2) the percentage of particles transferred between accelerators, and (3) the number of accelerators visited. One of the most salient predictions of our work is that the spectral index at high (low) energies is determined by the parameters of the first "seed" (second) acceleration region. By constructing the spectral distribution through an iterative analytical process, the evolution toward a double power law is easily characterized and explained. Our analytical model provides tools to interpret space- and ground-based observations from RHESSI, FOXSI, NuSTAR, Solar Orbiter/STIX, EOVSA, and future high-energy missions.

Solar flares↗

Size of metallic and polyethylene debris particles in failed cemented total hip replacements

Reports of differing failure rates of total hip prostheses made of various metals prompted us to measure the size of metallic and polyethylene particulate debris around failed cemented arthroplasties. We used an isolation method, in which metallic debris was extracted from the tissues, and a non-isolation method of routine preparation for light and electron microscopy. Specimens were taken from 30 cases in which the femoral component was of titanium alloy (10), cobalt-chrome alloy (10), or stainless steel (10). The mean size of metallic particles with the isolation method was 0.8 to 1.0 microns by 1.5 to 1.8 microns. The non-isolation method gave a significantly smaller mean size of 0.3 to 0.4 microns by 0.6 to 0.7 microns. For each technique the particle sizes of the three metals were similar. The mean size of polyethylene particles was 2 to 4 microns by 8 to 13 microns. They were larger in tissue retrieved from failed titanium-alloy implants than from cobalt-chrome and stainless-steel implants. Our results suggest that factors other than the size of the metal particles, such as the constituents of the alloy, and the amount and speed of generation of debris, may be more important in the failure of hip replacements.

Non-NASA Center↗

The reliability of macroscopic identification of lunar coarse fines particles and the petrogenesis of 2-4 mm particles in Apennine Front sample 15243

Particles from the 15243 Apennine Front sample (Station 6) were studied using macroscopic, neutron activation, and petrographic methods. It was found that macroscopic identifications allow for good identification of the particles. It was shown that the Station 6 regolith must consist of a fine-grained, comparatively KREEP-poor material and a coarser component consisting of Apollo 15 KREEP basalts and KREEP-like impact melt fragments.

Ryder, Graham↗

Equilibrium properties of a one-dimensional kinetic system.

One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.

Williams, J. H.↗

The application of nonintrusive optical methods for physical measurements in combustion

Laser optical diagnostic techniques have proven to be effective for making physical measurements in hostile combustion environments. Two major and complementary approaches have been applied; namely, single event imaging and multi-signal detection methods. Single event imaging methods (e.g., pulsed laser holography) can provide the essentially instantaneous observation of physical events, e.g., fuel atomization, droplet breakup, particle and surface combustion, condensed phase formation, and flow field visualization. Multi-signal detection methods, which involve the processing of signals produced by the light scattered from the interaction of particles and/or droplets with an incident coherent ray, can provide information on flow velocity, turbulence, fluid shear stresses, and particle/droplet size and velocity. Recent results using both approaches are presented.

Wuerer, J. E.↗

Multi-Scale Modelling of the Bound Metal Deposition Manufacturing of Ti6Al4V

Nonlinear shrinkage of the metal part during manufacturing by bound metal deposition, both on the ground and under microgravity, is considered. A multi-scale physics-based approach is developed to address the problem. It spans timescales from atomistic dynamics on the order of nanoseconds to full-part shrinkage on the order of hours. This approach enables estimation of the key parameters of the problem, including the widths of grain boundaries, the coefficient of surface diffusion, the initial redistribution of particles during the debinding stage, the evolution of the microstructure from round particles to densely-packed grains, the corresponding changes in the total and chemical free energies, and the sintering stress. The method has been used to predict shrinkage at the levels of two particles, of the filament cross-section, of the sub-model, and of the whole green, brown, and metal parts.

Nonlinear shrinkage↗