A Hybrid Finite Volume-Smoothed Particle Hydrodynamics Approach for Shock Capturing
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Particle fuel forms generally consist of a dispersion of fuel, such as tristructural isotropic (TRISO) particles, within a refractory matrix (e.g., graphite or silicon carbide). The density of matrix materials for particle fuel forms is of interest for modeling fuel form strength and thermal properties and may be specified as a quality control parameter, depending on reactor design. Some of the uncertainty associated with traditional, manual approaches can be eliminated by performing x-ray computed tomography (XCT) on the fuel forms and applying image processing methods to generate a precise count of the number of particles. This also removes the need to include determination of particle count within each individual fuel form during fabrication. Unfortunately, reconstruction artifacts from high-Z uranium-bearing kernels prevent accurate measurement of individual particle volumes using this approach, so the use of mean particle mass and volume are still necessary for computation of average fuel form matrix density. This method of using XCT to count particles in individual fuel form for determination of average matrix density was applied to three archived compacts from the Advanced Gas Reactor Fuel Development and Qualification (AGR)-1 campaign, four archived compacts with uranium carbide/uranium oxide (UCO) TRISO from the AGR-2 campaign, and three archived UO 2 -TRISO compacts from the AGR-2 campaign. The resulting density values were compared with those previously reported, showing slight changes due to uncertainties in the previously used number of particles in each of these cylindrical, graphite matrix compacts.
Particle-laden flows in solar photovoltaic (PV) systems are inevitable, where wind-swept debris in open environments are carried by high winds and turbulence, coating panel surfaces or damaging structures. Particle deposition, or soiling, is a well-known issue for large-scale plants which rely on uninhibited solar rays for optimal production. But understanding the mechanisms leading to soiling requires a physical and fluid dynamics-centered focus, since turbulence dominates PV panel wakes and is also known to alter particle concentration and trajectories. This study presents an experimental campaign toward consequences of particle-laden flow between two model PV panels using time-resolved particle image velocimetry. The model array was subjected to varied particle volume fractions, including a tracer particle case and a water droplet case. Characterization of mean velocity, turbulence statistics, and mean kinetic energy within the single phase and, separately, particle phase flows showed modified features due to particle inertia. Images captured at a frequency of 1 kHz in the near wake of the upstream panel allow for a first experimental look at vorticity and convective velocity of vortex structures for single-phase and particle-phase flows which are crucial to debris transport and soiling in PV environments.
Two problems of similar mathematical structure are presented: the thermocapillary motion of bubbles and the electrophoresis of colloidal particles. For both problems, it is shown that in a cloud of n particles surrounded by an infinite expanse of fluid, the velocity of each sphere under creeping flow conditions is equal to the velocity of an isolated particle, unchanged by interactions between the particles. However, when the cloud fills a container, conservation of mass shows that this result cannot continue to hold, and the average translational velocity must be calculated subject to a constraint on the mass flux. It is concluded that the average thermocapillary or electrophoretic translational velocity of a particle in the cloud is related to the effective conductivity of the cloud over the whole range of particle volume fractions, provided that the particles are identical, non-conducting and, for the thermocapillary problem, inviscid.
Particle dispersion and the influence that particle momentum exchange has on the properties of a turbulent carrier flow in micro-gravity environments challenge present understanding and predictive schemes. The objective of this effort has been to develop and assess high-fidelity simulation tools for predicting particle transport within micro-gravity environments suspended in turbulent flows. The computational technique is based on Direct Numerical Simulation (DNS) of the incompressible Navier-Stokes equations. The particular focus of the present work is on the class of dilute flows in which particle volume fractions and inter-particle collisions are negligible. Particle motion is assumed to be governed by drag with particle relaxation times ranging from the Kolmogorov scale to the Eulerian timescale of the turbulence and particle mass loadings up to one. The velocity field was made statistically stationary by forcing the low wavenumbers of the flow. The calculations were performed using 96(exp 3) collocation points and the Taylor-scale Reynolds number for the stationary flow was 62. The effect of particles on the turbulence was included in the Navier-Stokes equations using the point-force approximation in which 96(exp 3) particles were used in the calculations. DNS results show that particles increasingly dissipate fluid kinetic energy with increased loading, with the reduction in kinetic energy being relatively independent of the particle relaxation time. Viscous dissipation in the fluid decreases with increased loading and is larger for particles with smaller relaxation times. Fluid energy spectra show that there is a non-uniform distortion of the turbulence with a relative increase in small-scale energy. The non-uniform distortion significantly affects the transport of the dissipation rate, with the production and destruction of dissipation exhibiting completely different behaviors. The spectrum of the fluid-particle energy exchange rate shows that the fluid drags particles at low wavenumbers while the converse is true at high wavenumbers for small particles. A spectral analysis shows that the increase of the high wavenumber portion of the fluid energy spectrum can be attributed to transfer of the fluid-particle covariance by the fluid turbulence. This in turn explains the relative increase of small-scale energy caused by small particles observed in the present simulations as well as those of others.
The suspension balance model (SBM) for viscous Stokes flow has been well studied for the case of identical, monodisperse spherical particles in a channel, however, more work remains to be done to explore the bidisperse and polydisperse SBM. Here, we present a simple extension of the SBM that allows for modeling of suspensions of particles of bidisperse size. The comparison with avaiable experiments and direct simulations is found to be good. Additionally, we present a range of simulations to justify the assumptions made in the polydisperse suspension balance model for non-Brownian, neutrally buoyant particles that vary in size with a bidisperse distribution. The aim is to study the effects of moderate size variation on the rheology of the suspension and the distribution of the particles across the flow, especially in the near-wall regions. It is shown that, at the size ratios considered, bidispersity in the particle size does not significantly affect rheological parameters including the particle volume fraction and relative viscosity in the bulk of a Couette flow. The particle phase stresses are distributed between small and large particles in proportion to the relative volume fractions of each species. In bidisperse suspensions, the large particles in the wall layer develop a spatial structure and form chains in the streamwise direction.
A relation between the ensemble average and the nearest particle statistics is derived. The relation accounts for interactions among all the particles, not only the nearest one, and can be used to study long-range particle interactions without the difficulty of a divergent integral. As an example, this relation is applied to calculate the particle sedimentation velocity to the first order of the particle volume fraction. Using the relation, an important particle–fluid–particle stress is introduced for general multiphase flows.
Dual-wavelength (3.2 and 0.32 cm, i.e., X- and W-radar bands) radar ratio (DWR) measurements in ice clouds and precipitation using Canada's National Research Council Institute for Aerospace Research airborne radar are compared to closely collocated particle microphysical in situ sampling data in order to develop relations between DWR and characteristic hydrometeor size. This study uses the radar and in situ data sets collected during the In-Cloud ICing and Large-drop Experiment (ICICLE) campaign in midlatitude frontal clouds. Since atmospheric particle scattering at X band is predominantly in the Rayleigh regime and the W-band frequency is the highest frequency usually used for hydrometeor remote sensing, the X–W-band combination provides a relatively strong dual-wavelength reflectivity difference. This study considers radar and in situ measurements conducted in relatively homogeneous cloud and precipitation conditions. Measurements show that under these conditions, the difference between the X-band radar reflectivities observed with vertical and horizontal pointing of the radar beam are generally small and often negligible. However, W-band reflectivities at vertical beam pointing are, on average, larger than those for horizontal beam pointing by about 4 dB, which is a non-Rayleigh scattering effect from preferentially oriented non-spherical particles. A horizontal radar beam DWR–mean volume particle size relation, D v , provides robust estimates of this characteristic size for populations of particles with different habits. Uncertainties of D v retrievals using DWR are around 0.6 mm when D v is greater than approximately 1 mm. Size estimates using vertical radar beam DWRs have larger uncertainties due to smaller dual-wavelength signals and stronger influences of hydrometeor habits and orientations at this geometry of beam pointing. Mean relations among different characteristic sizes, which describe the entire particle size distribution (PSD), such as D v , and other sizes used in various applications (e.g., the mean, effective, and median sizes) are derived, so the results of this study can be used for estimating different PSD characteristic sizes.
As determined in situ by impactor samplers flown on an ER-2 at 16.5- to 20.7-km pressure altitude and on a DC-8 at 9.5- to 12.6-km pressure altitudes, the 1991 Pinatubo volcanic eruption increased the particle surface area of stratospheric aerosols up to 50-fold and the particle volume up to 2 orders of magnitude. Particle composition was typical of a sulfuric acid-water mixture at ER-2 altitudes. Ash particles coated with sulfuric acid comprised a significant fraction of aerosol at DC-8 altitudes. Mie-computed light extinction increased up to 20-fold at midvisible and greater than 100-fold at near-IR wavelengths. The optical thickness measured through the aerosol layer by an autotracking Sun photometer aboard a DC-8 aircraft at 10.7- to 11.3-km pressure altitudes shows a spectral shape that is similar to the Mie-calculated spectral extinction at ER-2 altitudes. Surface area distributions calculated by inversion of spectral optical depth measurements show characteristics that are similar to the mean surface area distribution resulting from 35 in situ measurements.
During the New England Air Quality Study (NEAQS) in the summer of 2004, airborne measurements were made of the major inorganic ions and the water-soluble organic carbon (WSOC) of the submicron (PM(sub 1.0)) aerosol. These and ancillary data are used to describe the overall aerosol chemical characteristics encountered during the study. Fine particle mass was estimated from particle volume and a calculated density based on measured particle composition. Fine particle organic matter (OM) was estimated from WSOC and a mass balance analysis. The aerosol over the northeastern United States (U.S.) and Canada was predominantly sulfate and associated ammonium, and organic components, although in unique plumes additional ionic components were also periodically above detection limits. In power generation regions, and especially in the Ohio River Valley region, the aerosol tended to be predominantly sulfate (approximately 60% micro gram /micro gram) and apparently acidic, based on an excess of measured anions compared to cations. In all other regions where sulfate concentrations were lower and a smaller fraction of overall mass, the cations and anions were balanced suggesting a more neutral aerosol. In contrast, the WSOC and estimated OM were more spatially uniform and the fraction of OM relative to PM mass was largely influenced by sources of sulfate. The study median OM mass fraction was 40%. Throughout the study region, sulfate and organic aerosol mass were highest near the surface and decreased rapidly with increasing altitude. The relative fraction of organic mass to sulfate was similar throughout all altitudes within the boundary layer (altitude less than 2.5 km), but was significantly higher at altitude layers in the free troposphere (above 2.5 km). A number of distinct biomass burning plumes from fires in Alaska and the Yukon were periodically intercepted, mostly at altitudes between 3 and 4 km. These plumes were associated with highest aerosol concentrations of the study and were largely comprised of organic aerosol components (approximtely 60%).
Data is sub-hourly size distribution, total integrated particle number, and total integrated particle volume from the SMPS. SMPS reports fine (less than 1 um diameter) aerosol. Note, size ranges reported differ slightly between months; read headers. Aerosol is dried to < 50% RH prior to sampling.
The water balance of a hurricane is controlled by boundary layer inflow, near vertical motion in the eyewall causing coalescence precipitation at above and residual ice precipitation at below freezing temperatures, and cirrus outflow at below -40 C aloft. In this paper we address the question of efficiency of water removal by this cirrus outflow which is important for the release of latent heat at high altitudes and its role in the dynamic flow at that level. During NASA's 1998 Convection and Moisture Experiment campaign we acquired microphysical outflow data in order to (1) determine the release and redistribution of latent heat near the top of hurricanes, (2) aid in TRMM algorithm development for remote sensing of precipitation, and (3) determine the optical/radiative characteristics of hurricane outflow. The data were acquired with Particle Measuring Systems two dimensional imaging spectrometers. On 23 August and again during the hurricane's landfall on 26 August, 1998, the NASA DC-8 aircraft penetrated hurricane 'Bonnie' four times each near 200 hPa pressure altitude. The eye crossing times were determined by (1) zero counts of cloud particles, (2) approximately 5 C increases in static and potential temperatures, and (3) minima in speeds and changes of direction of horizontal winds. The vertical winds showed shear between -6 m per second and +4 m per second and tangential winds approached 30 m per second in the eyewall. The particle volumes in the eyewall (determined by the pixels the particles shadowed in the direction of flight [x-direction] and normally to it by the number of diodes that they shadowed [y-direction]) ranged between 0.5 and 5.0 cubic centimeters per cubic meter. With a particle density near 0.2 g per cubic centimeter (determined from in situ melting and evaporation on a surface collector), the 1.0 g per meter corresponding mass of cloud ice ranged between 0.27 and 2.7 g per kilograms yielding horizontal fluxes between 8.1 and 81 g per square meters per second. The outflow ice was concentrated in crystals of a modal size of 190 micrometers. The particle size distributions were heavily skewed toward sizes with 98% of all cirrus particles smaller than the modal size comprising, however, only 20% of the mass. Thus the smaller than modal size particles dominantly affected the optical/radiative characteristics of the cloud, whereas the larger than modal size crystals determined the ice mass, hence dominated latent heating. Questions to be addressed relate to the origin of individual ice particles as the hurricane evolved and the likelihood of pristine and aggregate particle formation under the complicated conditions of rotation and outflow in the eyewall.
A technique for particle size measurement by using Laser Doppler Anemometry is discussed. An additional gate photomultiplier has been introduced at right angles to the optical axis in order to select only those particles passing through the central region of the measurement control volume. Particle sizing measurements have been made in sprays of glass particles using the modified Laser Anemometry system. Measurements in fuel sprays are also reported and compared with the results obtained by a photographic technique. The application of the particle sizing technique to opaque particles is investigated and suitable optical arrangements are suggested. Light scattering characteristics of Laser Anemometry systems for different optical geometries are calculated to select the optimum optical arrangement for the particle sizing measurements.
These data were collected during the Aerosol and Cloud Experiments in the Eastern North Atlantic field campaign as part of ARM Aerial Facility deployment (ACE-ENA, https://www.arm.gov/research/campaigns/aaf2017ace-ena). The ARM Aerial Facility Gulfstream-1 was deployed at Lajes Air Base (IATA: TER, ICAO: LPLA), on Terceira Island in the Azores, Portugal, for the two Intensive Observation Periods from June 20 through July 22, 2017 (IOP#1) and from January 11 through February 22, 2018 (IOP#2). The G-1 aircraft performed 20+19 research flights over the ARM Eastern North Atlantic (ENA) site and Atlantic Ocean to measure atmospheric turbulence, cloud water content and drop size distributions, aerosol precursor gases, aerosol chemical composition and size distributions. The current data set presents re-processed Particle Volume Monitor PVM-100A (aka Gerber probe) data: Liquid Water Content (LWC), Particle Surface Area (PSA), and the effective droplet radius (re) averaged to 50 Hz, 10 Hz, and 1 Hz.
These data were collected during the Cloud, Aerosol, and Complex Terrain Interactions (CACTI; https://www.arm.gov/research/campaigns/amf2018cacti ) field campaign in the Sierras de Córdoba mountain range of north-central Argentina as part of the ARM Aerial Facility (AAF) deployment. The ARM Aerial Facility Gulfstream-1 was operated from Las Higueras Airport (IATA: RCU, ICAO: SAOC), Río Cuarto, Córdoba, Argentina, for the Intensive Observation Period (IOP) from Nov. 1 through Dec. 15, 2018. The G-1 aircraft performed 22 research flights over the first ARM Mobile Facility (AMF1) location in the Sierras de Córdoba mountain range to measure atmospheric state and turbulence, cloud water content and droplet size distributions, aerosol precursor gases, and aerosol chemical composition and size distributions. The current data set presents re-processed Particle Volume Monitor PVM-100A (aka Gerber probe) data: Liquid Water Content (LWC), Particle Surface Area (PSA), and the effective droplet radius (re) averaged to 50Hz, 10Hz, and 1Hz.
The objectives are to use model transparent monotectics to obtain fundamental information applicable to two-phase systems in general, to apply this understanding to materials of interest in the Microgravity Science and Applications program, and to interpret results of flight experimental involving monotectic alloys. A number of model immiscible systems are in use to study various aspects of two-phase behavior within the miscibility gap and during solidification. Particle growth, coalescence and particle motions are under investigation using a holographic microscopy system. The system is capable of working with particle densities up to 10 to the 7th power particles/cubic centimeters through a 100 micron depth and can resolve particles of the order of 2 to 3 micron in diameter throughout the entire cell volume. Particle size, distribution changes with respect to time and temperature are observable from sequential holograms. Initial experiments using diethylene glycol/ethyl salicylate (DEG/ES) have demonstrated the usefulness of the technique. The thermal system controls temperature to at least plus or minus 0.001 K over the course of an experiment. A time-lapse film, made from holograms, of a succinonitrile/water solution shows particle size and number distribution changes with time under isothermal conditions. The observations are consistent with Ostwald ripening theory.
Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.
We report that for a decade, fully ceramic microencapsulated (FCM) fuel, containing tri-structural isotropic (TRISO) fuel particles in a silicon carbide (SiC) matrix, has been investigated as an accident-tolerant fuel for light water reactors (LWRs). Other examples exist of TRISO-based concepts for LWR fuels with different matrix materials. Previous studies assumed TRISO particle volume packing of approximately 0.44 in SiC (or another) matrix, the highest realistic packing fractions possible with conventional manufacturing. Recent advances in advanced manufacturing have yielded the development and demonstration of a fuel form that consists of conventionally manufactured TRISO particles in a 3D-printed SiC matrix with significantly higher possible TRISO packing fractions (0.5–0.7). This increased uranium loading enhances the viability of using TRISO-based particle fuel forms in LWRs. The viability of high-packing-fraction TRISO-based particle fuel forms in LWRs is assessed from the perspective of fuel cycle length, achievable fuel burnup, reactivity coefficients, and fuel cycle performance. Higher-packing-fraction TRISO-based fuel enables either longer cycle lengths (by ~25% at a packing fraction of 0.55 relative to 0.44) at a constant enrichment or decreased enrichments (by ~25% at a packing fraction of 0.55 relative to 0.44) at a constant cycle length. Studies of different fuel kernel types (uranium nitride, uranium oxycarbide, and uranium carbide) yield similar results, although the cycle length of uranium oxycarbide is shorter than for uranium nitride or uranium carbide (due to the lower density of uranium oxide). This work also characterized the production of 14 C resulting from neutron absorption in 14 N during operation for uranium mononitride fuel kernels; the ratio of 14 C/N was 1–2 at. % at discharge. For the fuel cycle evaluation, the activity of spent nuclear fuel and high-level waste at 100 and 100,000 years was lower for high-packing-fraction fuels than for conventional LWR fuel. Environmental impact metrics were similar overall, but higher on the front end of the fuel cycle and lower on the back end of the fuel cycle. Reactivity coefficients of higher-packing-fraction TRISO-based fuel were reasonable compared with those of conventional fuels.