Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solids flow”

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 1,783 records · Page 99

Radiant Extinction Of Gaseous Diffusion Flames

The absence of buoyancy-induced flows in microgravity (mu:g) and the resulting increase in the reactant residence time significantly alters the fundamentals of many combustion processes. Substantial differences between normal gravity (ng) and :g flames have been reported in experiments on candle flames [1, 2], flame spread over solids [3, 4], droplet combustion [5,6], and others. These differences are more basic than just in the visible flame shape. Longer residence times and higher concentration of combustion products in the flame zone create a thermochemical environment that changes the flame chemistry and the heat and mass transfer processes. Processes such as flame radiation, that are often ignored in ng, become very important and sometimes even controlling. Furthermore, microgravity conditions considerably enhance flame radiation by: (i) the build-up of combustion products in the high-temperature reaction zone which increases the gas radiation, and (ii) longer residence times make conditions appropriate for substantial amounts of soot to form which is also responsible for radiative heat loss. Thus, it is anticipated that radiative heat loss may eventually extinguish the Aweak@ (low burning rate per unit flame area) :g diffusion flame. Yet, space shuttle experiments on candle flames show that in an infinite ambient atmosphere, the hemispherical candle flame in :g will burn indefinitely [1]. This may be because of the coupling between the fuel production rate and the flame via the heat-feedback mechanism for candle flames, flames over solids and fuel droplet flames. Thus, to focus only on the gas-phase phenomena leading to radiative extinction, aerodynamically stabilized gaseous diffusion flames are examined. This enables independent control of the fuel flow rate to help identify conditions under which radiative extinction occurs. Also, spherical geometry is chosen for the :g experiments and modeling because: (i) It reduces the complexity by making the problem one-dimensional. (ii) The spherical diffusion flame completely encloses the soot which is formed on the fuel rich side of the reaction zone. This increases the importance of flame radiation because now both soot and gaseous combustion products co-exist inside the high temperature spherical diffusion flame. (iii) For small fuel injection velocities, as is usually the case for a pyrolyzing solid, the diffusion flame in :g around the solid naturally develops spherical symmetry. Thus, spherical diffusion flames are of interest to fires in :g and identifying conditions that lead to radiation-induced extinction is important for spacecraft fire safety.

Berhan, S.↗

Multiphase Model of Semisolid Slurry Generation and Isothermal Holding During Cooling Slope Rheoprocessing of A356 Al Alloy

In the present paper, we present an experimentally validated 3D multiphase and multiscale solidification model to understand the transport processes involved during slurry generation with a cooling slope. In this process, superheated liquid alloy is poured at the top of the cooling slope and allowed to flow along the slope under the influence of gravity. As the melt flows down the slope, it progressively loses its superheat, starts solidifying at the melt/slope interface with formation of solid crystals, and eventually exits the slope as semisolid slurry. In the present simulation, the three phases considered are the parent melt as the primary phase, and the solid grains and air as secondary phases. The air phase forms a definable air/liquid melt interface as the free surface. After exiting the slope, the slurry fills an isothermal holding bath maintained at the slope exit temperature, which promotes further globularization of microstructure. The outcomes of the present model include prediction of volume fractions of the three different phases considered, grain evolution, grain growth, size, sphericity and distribution of solid grains, temperature field, velocity field, macrosegregation and microsegregation. In addition, the model is found to be capable of making predictions of morphological evolution of primary grains at the onset of isothermal coarsening. The results obtained from the present simulations are validated by performing quantitative image analysis of micrographs of the rapidly oil-quenched semisolid slurry samples, collected from strategic locations along the slope and from the isothermal slurry holding bath.

36 MATERIALS SCIENCE↗

Topology optimization of additively manufactured fluidic components free of internal support structures

This report integrates projection-based approaches for implementing overhang constraints with fluid topology optimization to design additively manufactured fluidic components that do not require internal support structures. Internal support structures are challenging, potentially impossible to remove in complex fluid networks, yet would degrade fluidic performance if left within the part. This presents a major challenge in coupling topology optimization with additive manufacturing in the design of fluid components. The proposed approach leverages past work in overhang constraints for solid mechanics, including projection formulations and adjoint sensitivity analysis, and considers incompressible Navier–Stokes equations for laminar fluid flow. The approach is demonstrated on 2D and 3D fluid topology optimization problems to minimize pressure drop and constrain mass flow rates in pipes and manifolds. Resulting designs are crisp, logical, satisfy performance constraints, and satisfy overhang constraints, eliminating the need for internal support structures.

42 ENGINEERING↗

A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations

We introduce an extension of the divide-and-conquer (DC) algorithmic paradigm called divide-conquer-recombine (DCR) to perform large quantum molecular dynamics (QMD) simulations on massively parallel supercomputers, in which interatomic forces are computed quantum mechanically in the framework of density functional theory (DFT). In DCR, the DC phase constructs globally informed, overlapping local-domain solutions, which in the recombine phase are synthesized into a global solution encompassing large spatiotemporal scales. For the DC phase, we design a lean divide-and-conquer (LDC) DFT algorithm, which significantly reduces the prefactor of the O(N) computational cost for N electrons by applying a density-adaptive boundary condition at the peripheries of the DC domains. Our globally scalable and locally efficient solver is based on a hybrid real-reciprocal space approach that combines: (1) a highly scalable real-space multigrid to represent the global charge density; and (2) a numerically efficient plane-wave basis for local electronic wave functions and charge density within each domain. Hybrid space-band decomposition is used to implement the LDC-DFT algorithm on parallel computers. A benchmark test on an IBM Blue Gene/Q computer exhibits an isogranular parallel efficiency of 0.984 on 786 432 cores for a 50.3 × 10{sup 6}-atom SiC system. As a test of production runs, LDC-DFT-based QMD simulation involving 16 661 atoms is performed on the Blue Gene/Q to study on-demand production of hydrogen gas from water using LiAl alloy particles. As an example of the recombine phase, LDC-DFT electronic structures are used as a basis set to describe global photoexcitation dynamics with nonadiabatic QMD (NAQMD) and kinetic Monte Carlo (KMC) methods. The NAQMD simulations are based on the linear response time-dependent density functional theory to describe electronic excited states and a surface-hopping approach to describe transitions between the excited states. A series of techniques are employed for efficiently calculating the long-range exact exchange correction and excited-state forces. The NAQMD trajectories are analyzed to extract the rates of various excitonic processes, which are then used in KMC simulation to study the dynamics of the global exciton flow network. This has allowed the study of large-scale photoexcitation dynamics in 6400-atom amorphous molecular solid, reaching the experimental time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neutron tomography of axially symmetric objects using 14 MeV neutrons from a portable neutron generator

In nuclear boiling water reactor cores, the distribution of water and steam (void) is essential for both safety and efficiency reasons. In order to enhance predictive capabilities, void distribution assessment is performed in two-phase test-loops under reactor-relevant conditions. This article proposes the novel technique of fast-neutron tomography using a portable deuterium-tritium neutron generator to determine the time-averaged void distribution in these loops. Fast neutrons have the advantage of high transmission through the metallic structures and pipes typically concealing a thermal-hydraulic test loop, while still being fairly sensitive to the water/void content. However, commercially available fast-neutron generators also have the disadvantage of a relatively low yield and fast-neutron detection also suffers from relatively low detection efficiency. Fortunately, some loops are axially symmetric, a property which can be exploited to reduce the amount of data needed for tomographic measurement, thus limiting the interrogation time needed. In this article, three axially symmetric test objects depicting a thermal-hydraulic test loop have been examined; steel pipes with outer diameter 24 mm, thickness 1.5 mm, and with three different distributions of the plastic material POM inside the pipes. Data recorded with the FANTOM fast-neutron tomography instrument have been used to perform tomographic reconstructions to assess their radial material distribution. Here, a dedicated tomographic algorithm that exploits the symmetry of these objects has been applied, which is described in the paper. Results are demonstrated in 20 rixel (radial pixel) reconstructions of the interior constitution and 2D visualization of the pipe interior is demonstrated. The local POM attenuation coefficients in the rixels were measured with errors (RMS) of 0.025, 0.020, and 0.022 cm{sup −1}, solid POM attenuation coefficient. The accuracy and precision is high enough to provide a useful indication on the flow mode, and a visualization of the radial material distribution can be obtained. A benefit of this system is its potential to be mounted at any axial height of a two-phase test section without requirements for pre-fabricated entrances or windows. This could mean a significant increase in flexibility of the void distribution assessment capability at many existing two-phase test loops.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Droplet and Particle Methods to Investigate Turbulent Particle Laden Jets

The SARS-CoV-2 pandemic has heightened the interest in particle-laden turbulent jets generated by breathing, talking, coughing and sneezing, and how these can contribute to disease transmission. We present quantitative measurement methods for such flows, while exploring and offering improvements for common shortcomings. We generate jets consisting of either liquid droplets or solid particles in an isothermal, quiescent and electrically isopotential experimental chamber that was constructed to control the effects of ambient forcing on jet behavior. For liquid droplets, we find promise in surface deposition analysis based on fluorescent tracer use. For particles, we explore the performance of commercially available adhesive sampling strips and develop conductive grounded carbon tape based sampling strips. We explore ways in which the smallest of thermal gradients or electrostatic charge issues can affect particle dispersion, and suggest practical methods to address these issues. The developed methods are applied to study the simultaneous deposition of (Formula presented.) 25, 50 and 200 μm solid particles from a particle laden turbulent jet with a mean velocity of 33.2 m/s. The deposition location as a function of particle size was compared to results from a simple numerical RANS model, and illustrates ways in which imprecise initial or boundary conditions can lead to a notable deviation from experimental results. The differences in deposition pattern seen in experimental and numerical results despite a carefully controlled environment and characterized particle ejection indicate the need for a more stringent numerical model validation, especially when studying fate and transport of mid-range (neither purely aerosol or ballistic) sized particles.

42 ENGINEERING↗

Universal behavior in fragmenting brittle, isotropic solids across material properties

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. Importantly, this evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

36 MATERIALS SCIENCE↗

Technetium Reduction and Permanent Sequestration by Abiotic and Biotic Formation of Low-Solubility Sulfide Mineral Phases

One way to minimize the mobility of the Tc VII oxyanion pertechnetate (TcO 4 - ) is to effect reduction under sulfidogenic conditions (generated abiotically by Fe 0 or biotically) to form TcS x , which is significantly slower to oxidize than Tc IV O 2 . In sediment systems, TcS x and other precipitates may oxidize more slowly due to oxygen diffusion limitations to these low permeability precipitate zones. In addition, the TcO 4 - reduction rate may be more rapid in the presence of sediment because of additional reductive surface phases. This project aims to provide a fundamental understanding of the feasibility of immobilization of TcO 4 - as TcS x in the vadose zone or groundwater by application nano zero-valent iron (nZVI), and sulfide or sulfate. Biotic batch experiments have used the sulfate-reducing bacterium (SRB) Desulfotomaculum reducens . The iron sulfide mineral mackinawite was generated under these conditions, while vivianite was formed in nZVI only controls. The sulfide/bacteria-containing system consistently reduced aqueous pertechnetate rapidly (> 95% in the first hour), a rate similar to that for the sulfide-free, nZVI only system. Reduced Tc (aged for 3 months) generated in both SRB/nZVI systems was highly resistant to reoxidation. In reduced samples, Tc was found associated with solid phases containing Fe and S (D. reducens/nZVI) or Fe (nZVI only). Experiments using D. reducens without nZVI provided some additional insights. Firstly, stationary phase cultures were able to slowly reduce pertechnetate. Secondly, addition of pertechnetate at the beginning of cell growth (lag phase) resulted in a faster rate of Tc reduction, possibly indicating a direct (e.g. enzymatic) role for D. reducens in Tc reduction. Abiotic batch experiments were conducted with Na 2 S as the sulfide source. Pertechnetate reduction was rapid in the presence of sulfide and nZVI, although the rate was suppressed at the higher S/Fe ratios tested. This suppression appeared to be due to the formation of Tc-containing colloids. As with the biotic experiments, pertechnetate reduced under sulfidic conditions was highly resistant to reoxidation. The microscopic morphology of abiotically-transformed nZVI particles varied significantly with those in the biotic experiment, although mackinawite was formed in both systems (as indicated by μXRD and Mössbauer spectroscopy). Preliminary XAS analysis pointed to a mixture of Tc-O and Tc-S binding in the abiotic sulfide/nZVI system, while the major reduced solids under non-sulfidic conditions were TcO 2 •nH 2 O. The presence of sediment and advective flow to the TcO 4 - /nZVI/sulfide system results in additional processes occurring. Although the natural Hanford sediment used has sufficient available ferrous iron to slowly reduce TcO 4 - , under anaerobic conditions, that rate is orders of magnitude slower than reduction by nZVI/sulfide. Batch and 1-D column experiments showed that the TcO 4 - reduction rate increased with the sediment surface area (with the same nZVI mass). As in batch systems, column studies showed that the presence of sulfide with TcO 4 - at low (2-5 mM) concentrations increased the TcO 4 - reduction rate and high (10-30 mM) sulfide decreased the rate. This change is attributed to the formation of sulfide precipitates on the nZVI and sediment surfaces. Injection of low and high sulfide (i.e. pretreatment) prior to TcO 4 - /sulfide injection also greatly decreased the TcO 4 - reduction rate, likely decreasing the generation of ferrous iron from the nZVI. Although the high sulfide systems have slower Tc reduction rates, 190 times more Tc mass precipitated than in the low sulfide systems and the highest fraction of Tc mass remained immobilized.

54 ENVIRONMENTAL SCIENCES↗

Novel Zoned Waste-forms for High-Priority Radionuclide Waste Streams

This report describes the potential of a novel class of materials—α-ZrW 2 O 8 , Zr 2 WP 2 O 12 , and related compounds that contract upon amorphization as possible radionuclide waste-forms. The proposed ceramic waste-forms would consist of zoned grains, or sintered ceramics with center- loaded radionuclides and barren shells. Radiation-induced amorphization would result in core shrinkage but would not fracture the shells or overgrowths, maintaining isolation of the radionuclide. In this report, we have described synthesis techniques to produce phase-pure forms of the materials, and how to fully densify those materials. Structural models for the materials were developed and validated using DFPT approaches, and radionuclide substitution was evaluated; U(IV), Pu(IV), Tc(IV) and Tc(VII) all readily substitute into the material structures. MD modeling indicated that strain associated with radiation-induced amorphization would not affect the integrity of surrounding crystalline materials, and these results were validated via ion beam experimental studies. Finally, we have evaluated the leach rates of the barren materials, as determined by batch and flow-through reactor experiments. ZrW 2 O 8 leaches rapidly, releasing tungstate while Zr is retained as a solid oxide or hydroxide. Tungsten release rates remain elevated over time and are highly sensitive to contact times, suggesting that this material will not be an effective waste-form. Conversely, tungsten releases rates from Zr 2 WP 2 O 12 rapidly drop, show little dependence on short-term changes in fluid contact time, and in over time, become tied to P release rates. The results presented here suggest that this material may be a viable waste-form for some hard-to-handle radionuclides such as Pu and Tc.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Model Validation and Uncertainty Quantification on the KRUSTY Microreactor Design Using GRIFFIN Neutron Transport Code [Poster]

Argonne National Laboratory (ANL) and INL have developed a GRIFFIN steady state neutronics model for the multiphysics simulations of the Kilopower Reactor Using Sterling TechnologY (KRUSTY) microreactor in the Multiphysics Object Oriented Simulation Environment (MOOSE). The reliability of such deterministic neutronics models can be validated by comparing with computations from Monte Carlo codes (e.g. MCNP, SERPENT, OpenMC, Shift, etc). Furthermore, potential modeling/design improvements can be identified by incorporating uncertainty quantification (UQ), which can be performed by MOOSE’s Stochastic Tools Module (STM). KRUSTY is a prototype for a 5-kW thermal nuclear-powered space reactor. Its primary components consist of nuclear fuel, heat pipes, a control rod, a reflector, and the shielding. The fuel consists of 3 stacked U-7.65Mo cylinders with a hole in the center for the control rod. 8 liquid sodium heat pipes transfer fission energy from the solid fuel block to the Sterling power conversion system where the energy is extracted, and the cooled sodium flows back to the core via capillary action . The movable Boron Carbide control rod regulates the neutron population during startup or when a reactor temperature boost is needed . The beryllium oxide reflector is in 3 places in the reactor; it surrounds the core axially, it lies beneath the core on a platen, and it is present in the shim. The axial and lower reflectors rest on an adjustable stainless-steel platen that moves upward to cover the fuel and help the reactor reach criticality. Lastly, radial stainless steel surrounds the core offering protection from radiation exposure .

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Phase Field Simulations of Microstructure Evolution in IN718 using a Surrogate Ni–Fe–Nb Alloy during Laser Powder Bed Fusion

The solidification microstructure in IN718 during additive manufacturing was modeled using phase field simulations. The novelty of the research includes the use of a surrogate Ni–Fe–Nb alloy that has the same equilibrium solidification range as IN718 as the model system for phase field simulations, the integration of the model alloy thermodynamics with the phase field simulations, and the use of high-performance computing tools to perform the simulations with a high enough spatial resolution for realistically capturing the dendrite morphology and the level of microsegregation seen under additive manufacturing conditions. Heat transfer and fluid flow models were used to compute the steady state temperature gradient and an average value of the solid-liquid (s-l) interface velocity that were used as input for the phase field simulations. The simulations show that the solidification morphology is sensitive to the spacing between the columnar structures. Spacing narrower than a critical value results in continued growth of a columnar microstructure, while above a critical value the columnar structure evolves into a columnar dendritic structure through the formation of secondary arms. These results are discussed in terms of the existing columnar to dendritic transition (CDT) theories. The measured interdendritic Nb concentration, the primary and secondary arm spacing is in reasonable agreement with experimental measurements performed on the nickel-base superalloy IN718.

36 MATERIALS SCIENCE↗

Determination of Plasma Screening Effects for Thermonuclear Reactions in Laser-generated Plasmas

Due to screening effects, nuclear reactions in astrophysical plasmas may behave differently than in the laboratory. The possibility to determine the magnitude of these screening effects in colliding laser-generated plasmas is investigated theoretically, having as a starting point a proposed experimental setup with two laser beams at the Extreme Light Infrastructure facility. A laser pulse interacting with a solid target produces a plasma through the Target Normal Sheath Acceleration scheme, and this rapidly streaming plasma (ion flow) impacts a secondary plasma created by the interaction of a second laser pulse on a gas jet target. We model this scenario here and calculate the reaction events for the astrophysically relevant reaction {sup 13}C({sup 4}He, n ){sup 16}O. We find that it should be experimentally possible to determine the plasma screening enhancement factor for fusion reactions by detecting the difference in reaction events between two scenarios of ion flow interacting with the plasma target and a simple gas target. This provides a way to evaluate nuclear reaction cross-sections in stellar environments and can significantly advance the field of nuclear astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

High-Resolution X-Ray and Neutron Computed Tomography of an Engine Combustion Network Spray G Gasoline Injector

Given the importance of the fuel-injection process on the combustion and emissions performance of gasoline direct injected engines, there has been significant recent interest in understanding the fluid dynamics within the injector, particularly around the needle and through the nozzles. Furthermore, the pressure losses and transients that occur in the flow passages above the needle are also of interest. Simulations of these injectors typically use the nominal design geometry, which does not always match the production geometry. Computed tomography (CT) using x-ray and neutron sources can be used to obtain the real geometry from production injectors, but there are trade-offs in using these techniques. X-ray CT provides high resolution, but cannot penetrate through the thicker parts of the injector. Neutron CT has excellent penetrating power but lower resolution. Here, we present results from a joint effort to characterize a gasoline direct injector representative of the Spray G injector as defined by the Engine Combustion Network. High-resolution (1.2 to 3 µm) x-ray CT measurements from the Advanced Photon Source at Argonne National Laboratory were combined with moderate-resolution (40 µm) neutron CT measurements from the High Flux Isotope Reactor at Oak Ridge National Laboratory to generate a complete internal geometry for the injector. This effort combined the strengths of both facilities’ capabilities, with extremely fine spatially resolved features in the nozzles and injector tips and fine resolution of internal features of the needle along the length of injector. Analysis of the resulting surface model of the internal fluid flow volumes of the injector reveals how the internal cross-sectional area and nozzle hole geometry differs slightly from the design dimensions. A simplified numerical simulation of the internal flow shows how deviations from the design geometry can alter the flow inside the sac and holes. Our results of this study will provide computational modelers with very accurate solid and surface models for use in computational fluid dynamics studies and experimentalists with increased insight into the operating characteristics of their injectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Aerothermodynamic technology for Space Shuttle - and beyond.

Wind-tunnel testing ground rules are discussed together with questions of orbiter heating, orbiter aerodynamics, flow fields, and ascent configurations. Aerodynamic heating is not expected to be of concern during ascent prior to solid rocket booster staging because of the relatively low staging velocity. It is pointed out that aerothermodynamics technology for the Space Shuttle orbiter has a three-year data base and is in hand. As shuttle support requirements diminish, and more effort can be brought to bear on advanced studies, advanced aerospace vehicle types, and their potential mission and operational requirements, will be explored in increasing depth.

Henderson, A., Jr.↗