Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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,225 records · Page 68

Direct force wall shear measurements in pressure-driven three-dimensional turbulent boundary layers

Unique, simultaneous direct measurements of the magnitude and direction of the local wall shear stress in a pressure-driven three-dimensional turbulent boundary layer are presented. The flow is also described with an oil streak wall flow pattern, a map of the wall shear stress-wall pressure gradient orientations, a comparison of the wall shear stress directions relative to the directions of the nearest wall velocity as measured with a typical, small boundary layer directionally sensitive claw probe, as well as limiting wall streamline directions from the oil streak patterns, and a comparison of the freestream streamlines and the wall flow streamlines. A review of corrections for direct force sensing shear meters for two-dimensional flows is presented with a brief discussion of their applicability to three-dimensional devices.

Mcallister, J. E.↗

A viscous shock-layer flowfield analysis by an explicit-implicit method

An implicit analogue of a widely used explicit method to external axisymmetric laminar flows with strong entropy gradients is extended. The details of the "numerics" of the implicit part are provided in a body oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analog may appeal to researchers for using it in the analysis of 3-D flow problems.

Gupta, R. M.↗

Advanced Containerless Processing Technology

Breadboards for high temperature containerless processing systems are to be developed, the principles of operation are to be studied, the performance is to be characterized, the limitations identified, and the influence of the acoustic field on the samples established. The subjects to be addressed are experimental and theoretical studies of: (1) acoustic positioning and manipulation capabilities in a high temperature gradient environment (from 26 C to 990 C); (2) acoustic waveforms, harmonic contents, power transfer, sample transport and stability associated with high temperature gradient system; (3) high temperature ground based levitation systems which will allow melting, processing, and solidifying samples without crucibles in the laboratory, (4) KC-135 and laboratory tests of various acoustic geometries which may have spherical applications in the Materials Processing in Space Program; and (5) provide technical information to Acoustics Containerless Experimental System (ACES) engineering team and establish the operation conditions for ACES.

Wang, T. G.↗

Effect of vertical magnetic field on convection and segregation in vertical Bridgman crystal growth

A previous finite-element analysis of vertical Bridgman growth for dilute and nondilute alloys is extended to include the effect of a vertically-aligned magnetic field in the limit of zero magnetic Reynolds number. Calculations are presented for growth of a dilute gallium-germanium alloy in a vertically stabilized Bridgman-Stockbarger system and in a furnace with a uniform temperature gradient imposed along the ampoule. Steady cellular convection driven by radial temperature gradients causes good axial and radial mixing in both systems without a magnetic field. A weak magnetic field decreases the intensity of convection and the effectiveness of solute mixing. The radial nonuniformity is greatest for an intermediate field strength. Stronger fields suppress flow recirculation completely, and lead to uniform solute segregation across the crystal and to diffusion-controlled axial segregation.

Kim, Do Hyun↗

A theory of nonlocal mixing-length convection. I - The moment formalism

A flexible and potentially powerful theory of convection, based on the mixing length picture, is developed to make unbiased self-consistent predictions about overshooting and other complicated phenomena in convection. The basic formalism is set up, and the method's power is demonstrated by showing that a simplified version of the theory reproduces all the standard results of local convection. The second-order equations of the theory are considered in the limit of a steady state and vanishing third moments, and it is shown that they reproduce all the standard results of local mixing-length convection. There is a particular value of the superadiabatic gradient, below which the only possible steady state of a fluid is nonconvecting. Above this critical value, a fluid is convectively unstable. Two distinct regimes of convection, which are identified as efficient and inefficient convection, are determined.

Grossman, Scott A.↗

Surface Penetrating Radar Simulations for Jupiter's Icy Moons

The icy moons of Jupiter (Europa, Callisto, and Ganymede) are of similar overall composition but show different surface features as a result of different sub-surface processes. Furthermore, each of these moons could have a liquid ocean of water buried underneath the icy crust, but their depth can only be speculated. For Europa, estimates put the thickness of the ice shell anywhere between 2-30 km, with'a few models predicting up to 100 km. Much of the uncertainties are due to the largely unknown temperature gradients and levels of water impurities across different surface layers. One of the most important geological processes is the possible transportation of heat by ice convection. If the ice is convecting, then an upper limit of about 20 km is set for the depth of the ocean underneath. Convection leads to a sharp increase in temperature followed by a thick region of nearly constant temperature. If ice is not convecting, then an exponentially increasing temperature profile is expected. The crust is thought to be a mixture of ice and rock, and although the exact percentage of rock is not known, it is expected to be low. Additionally, the ice crust could contain salt, similar to sea ice on Earth. The exact amount of salt and how that amount changes with depth is also unknown. In preparation for the Jupiter Icy Moons Orbiter (JIMO) mission, we performed simulations for a surface-penetrating radar investigating signatures for different possible surface and sub-surface structures of these moons in order to estimate the applicability of using radar with a frequency range between 1 and 50 MHz. This includes simulations of power requirements, attenuation losses, layer resolutions for scenarios with and without the presence of a liquid ocean underneath the ice, cases of convecting and non-convecting ice, different impurities within the ice, and different surface roughnesses.

Markus, Thorsten↗

Gravitational Influences on Flame Propagation Through Non-Uniform, Premixed Gas Systems

Flame propagation through non-uniformly premixed (or layered) gases has importance both in useful combustion systems and in unintentional fires. As summarized recently and in previous Microgravity Workshop papers, non-uniform premixed gas combustion receives scant attention compared to the more usual limiting cases of diffusion or uniformly premixed flames, especially regarding the role gravity plays. This paper summarizes our recent findings on gravitational effects on layered combustion along a floor, in which the fuel concentration gradient exists normal to the direction of flame spread. In an effort to understand the mechanism by which the flames spread faster in microgravity (and much faster, in laboratory coordinates, than the laminar burning velocity for uniform mixtures), we have begun making pressure measurements across the spreading flame front that are described here. Earlier researchers, testing in 1g, claimed that hydrostatic pressure differences could account for the rapid spread rates. Additionally, we present the development of a new apparatus to study flame spread in free (i.e., far from walls), non-homogeneous fuel layers formed in a flow tunnel behind an airfoil that has been tested in normal gravity.

Miller, Fletcher J.↗

The Use of Subsurface Barriers to Support Treatment of Metals and Reduce the Flux of Tritium to Fourmile Branch at the Savannah River Site in South Carolina - 13358

The Savannah River Site (SRS) produced tritium, plutonium, and special nuclear materials for national defense, medicine, and the space programs. Acidic groundwater plumes containing metals, metallic radionuclides, non-metallic radionuclides and tritium sourced from the F and H Area Seepage Basins have impacted the surface water of Fourmile Branch on SRS. Tritium releases from Fourmile Branch have impacted the water quality within areas of the Savannah River adjacent to the SRS, and this circumstance has been an ongoing regulatory concern. The F and H Area Seepage Basins operated until 1988 for the disposition of deionized acidic waste water from the F and H Separations Facilities. The waste water contained dilute nitric acid and low concentrations of non-radioactive metals, and radionuclides, with the major isotopes being Cs-137, Sr-90, U-235, U-238, Pu-239, Tc-99, I-129, and tritium. The tritium concentration in the waste water was relatively elevated because there is not a practicable removal method in water. The acid content of the waste water during the operational period of the basins was equal to 12 billion liters of nitric acid. The seepage basins were closed in 1988 and backfilled and capped by 1991. The plumes associated with the F and H basins cover an area of nearly 2.4 square kilometers (600 acres) and discharge along ∼2,600 meters of Fourmile Branch. The acidic nature of the plumes and their overall discharge extent along the branch represent a large challenge with respect to reducing contaminant flux to Fourmile Branch. The introduction of nitric acid into the groundwater over a long time effectively reduced the retardation of metal migration from the basins to the groundwater and in the groundwater to Fourmile Branch, because most negatively charged surfaces on the aquifer materials were filled with hydrogen ion. Two large pump and treat systems were constructed in 1997 and operated until 2003 in an attempt to capture and control the releases to Fourmile Branch. The operating cost, including waste disposal, for the two systems was ∼$1.3 M/month. Both systems employed reinjection of tritiated water up gradient of the extraction, and produced large quantities of waste from non-tritium isotopes and metals removal prior to reinjection. Both systems were determined to be ineffective and potentially detrimental with respect to limiting the flux of contaminants to Fourmile Branch. After it became apparent that there was very little benefit to continued operation of the systems, and the staggering cost of operations was recognized by the SRS and regulators, a new remedy was developed. The new system uses vertical subsurface barriers to redirect groundwater flow to limit the transport of contaminants to the stream. The barriers were constructed of acid resistant grout using deep soil mixing techniques. The grout mixture used low swelling clay, fly ash, and sodium hydroxide to form a pozzolana material with low permeability and low strength. The SRS and regulators agreed to a series of remedial goals, with the first goal to reduce tritium flux to the stream by 70% and bring constituents other than tritium to groundwater protection standards. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hydride Microstructure at the Metal-Oxide Interface of a Zircaloy-4 Fuel Clad from the H.B. Robinson Nuclear Reactor

During nuclear reactor operation, some of the hydrogen that is produced by the water side corrosion is uptaken by cladding. When the uptaken hydrogen content exceeds its solubility limit in the cladding at a specific temperature, it precipitates as hydride particles in the cladding. This process occurs immediately at the metal-oxide interface forming a hydride rim because of the temperature gradient in the cladding. The presence of hydride rim and its microstructure severely reduces the cladding strength and can induce crack initiation and propagation during design basis accidents, transportation, and wet/dry storage. Because hydrides platelets are one of the potent crack initiation sites, microstructural characterization of hydride rim becomes essential to the safe operation of nuclear reactor and spent nuclear fuel retrievability. Therefore, we investigate the hydride rim microstructure at the metal oxide interface using scanning transmission electron microscopy (S/TEM) with Energy Dispersive X-ray Spectroscopy (EDS) that allows for the microstructural and local chemistry characterization of hydride precipitates to be determined, and enables a relationship to be derived on how the precipitate sites are affected by the oxide and their complex stacking behavior of the precipitates.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Seawater pH reconstruction using boron isotopes in multiple planktonic foraminifera species with different depth habitats and their potential to constrain pH and <i>p</i>CO<sub>2</sub> gradients

Boron isotope systematics of planktonic foraminifera from core-top sediments and culture experiments have been studied to investigate the sensitivity ofδ11B of calcite tests to seawater pH. However, our knowledge of the relationship between δ11B and pH remains incomplete for many taxa. Thus, to expand the potential scope of application of this proxy, we report δ11B data for seven different species of planktonic foraminifera from sediment core tops. We utilize a method for the measurement of small samples of foraminifera and calculate the δ11B-calcite sensitivity to pH for Globigerinoides ruber, Trilobus sacculifer (sacc or without sacc), Orbulina universa, Pulleniatina obliquiloculata, Neogloboquadrina dutertrei, Globorotalia menardii, and Globorotalia tumida, including for unstudied core tops and species. These taxa have diverse ecological preferences and are from sites that span a range of oceanographic regimes, including some that are in regions of air–sea equilibrium and others that are out of equilibrium with the atmosphere. The sensitivity of δ11Bcarbonate to δ11Bborate (e.g., Δδ11Bcarbonate/Δδ11Bborate) in core tops is consistent with previous studies for T. sacculifer and G. ruber and close to unity for N. dutertrei, O. universa, and combined deep-dwelling species. Deep-dwelling species closely follow the core-top calibration fo rO. universa, which is attributed to respiration-driven microenvironments likely caused by light limitation and/or symbiont–host interactions. Our data support the premise that utilizing boron isotope measurements of multiple species within a sediment core can be utilized to constrain vertical profiles of pH andpCO2 at sites spanning different oceanic regimes, thereby constraining changes in vertical pH gradients and yielding insights into the past behavior of the oceanic carbon pumps.

54 ENVIRONMENTAL SCIENCES↗

Final Report for contract DE-SC0010151, "A Cost-Effective Oxygen Separation System Based on Open Gradient Magnetic Field by Polymer Beads"

At present much of the oxygen enrichment market is based primarily on cryogenic systems, ion transport membranes made of ceramics, or polymer membranes technologies. All these methods have some positive and some negative attributes, but remain too expensive for the cost effective production of pure oxygen needed for clean renewable energy conversion projects. In a very successful Phase I period, the ITN-TAMU team accomplished the completion of most of the sub-tasks on oxygen separation based on our proposed magnetic field gradients. Theoretical and computer simulation data under different conditions of flow rates, magnetic flux, magnetic beads packing density, pressure and temperature have resulted very valuable observations. Chitosan based polymer magnetic beads showed the best results for uniformity and column packing. Due to limited magnetic force of 1 Tesla, we were restricted to use only single small column of magnetic beads and limited air pressure with controlled flow rate. We also fabricated smaller sized beads column with metal and clear plastic frames to be used for more concentrated packed capacity. Beads made with rare earth Neodymium were not uniform and were quite brittle. Different surfactants have minimized the effect of the de-aggregation of polymer magnetic beads. This Phase I development work also included overcoming key packaging hurdles such as high packing capacity of magnetic polymer beads and formation of consistent batch of magnetic beads. Market and commercialization studies were also initiated to show that there is a growing market for the proposed cost-effective oxygen separation system for energy conversion projects.

20 FOSSIL-FUELED POWER PLANTS↗

Degradation Reduces Microbial Richness and Alters Microbial Functions in an Australian Peatland

Abstract Peatland ecosystems cover only 3% of the world’s land area; however, they store one-third of the global soil carbon (C). Microbial communities are the main drivers of C decomposition in peatlands, yet we have limited knowledge of their structure and function. While the microbial communities in the Northern Hemisphere peatlands are well documented, we have limited understanding of microbial community composition and function in the Southern Hemisphere peatlands, especially in Australia. We investigated the vertical stratification of prokaryote and fungal communities from Wellington Plains peatland in the Australian Alps. Within the peatland complex, bog peat was sampled from the intact peatland and dried peat from the degraded peatland along a vertical soil depth gradient (i.e., acrotelm, mesotelm, and catotelm). We analyzed the prokaryote and fungal community structure, predicted functional profiles of prokaryotes using PICRUSt, and assigned soil fungal guilds using FUNGuild. We found that the structure and function of prokaryotes were vertically stratified in the intact bog. Soil carbon, manganese, nitrogen, lead, and sodium content best explained the prokaryote composition. Prokaryote richness was significantly higher in the intact bog acrotelm compared to degraded bog acrotelm. Fungal composition remained similar across the soil depth gradient; however, there was a considerable increase in saprotroph abundance and decrease in endophyte abundance along the vertical soil depth gradient. The abundance of saprotrophs and plant pathogens was two-fold higher in the degraded bog acrotelm. Soil manganese and nitrogen content, electrical conductivity, and water table level (cm) best explained the fungal composition. Our results demonstrate that both fungal and prokaryote communities are shaped by soil abiotic factors and that peatland degradation reduces microbial richness and alters microbial functions. Thus, current and future changes to the environmental conditions in these peatlands may lead to altered microbial community structures and associated functions which may have implications for broader ecosystem function changes in peatlands.

Environmental Sciences & Ecology↗

Aspects of linear Landau damping in discretized systems

Basic linear eigenmode spectra for electrostatic Langmuir waves and drift-kinetic slab ion temperature gradient modes are examined in a series of scenarios. Collisions are modeled via a Lenard-Bernstein collision operator which fundamentally alters the linear spectrum even for infinitesimal collisionality [Ng et al., Phys. Rev. Lett. 83, 1974 (1999)]. A comparison between different discretization schemes reveals that a Hermite representation is superior for accurately resolving the spectra compared to a finite differences scheme using an equidistant velocity grid. Additionally, it is shown analytically that any even power of velocity space hyperdiffusion also produces a Case-Van Kampen spectrum which, in the limit of zero hyperdiffusivity, matches the collisionless Landau solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Avalanches driven by pressure gradients in a magnetized plasma

The results are presented for a basic heat transport experiment involving an off-axis heat source in which avalanche events occur. The configuration consists of a long, hollow, cylindrical region of elevated electron temperature embedded in a colder plasma, and far from the device walls [Van Compernolle et al. Phys. Rev. E 91, 031102(R) (2015)]. The avalanche events are identified as sudden rearrangements of the pressure profile following the growth of fluctuations from ambient noise. The intermittent collapses of the plasma pressure profile are associated with unstable drift-Alfvén waves and exhibit both radial and poloidal dynamics. After each collapse, the plasma enters a quiescent phase in which the pressure profile slowly recovers and steepens until a threshold is exceeded, and the process repeats. The use of reference probes as time markers allows for the visualization of the 2D spatio-temporal evolution of the avalanche events. Avalanches are observed only for a limited combination of heating powers and magnetic fields. At higher heating powers, the system transits from the avalanche regime into a regime dominated by sustained drift-Alfvén wave activity. Here, this manuscript focuses on new results that illustrate the individual contributions to the avalanche process from density and temperature gradients in the presence of zero-order, sheared flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cryogenic testing of the 2.1 GHz five-cell superconducting RF cavity with a photonic band gap coupler cell

We present results from cryogenic tests of the multi-cell superconducting radio frequency (SRF) cavity with a photonic band gap (PBG) coupler cell. Achieving high average beam currents is particularly desirable for future light sources and particle colliders based on SRF energy-recovery-linacs (ERLs). Beam current in ERLs is limited by the beam break-up instability, caused by parasitic higher order modes (HOMs) interacting with the beam in accelerating cavities. A PBG cell incorporated in an accelerating cavity can reduce the negative effect of HOMs by providing a frequency selective damping mechanism, thus allowing significantly higher beam currents. The multi-cell cavity was designed and fabricated of niobium. Two cryogenic (vertical) tests were conducted. The high unloaded Q-factor was demonstrated at a temperature of 4.2 K at accelerating gradients up to 3 MV/m. The measured value of the unloaded Q-factor was 1.55 × 10 8 , in agreement with prediction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A predictive discrete-continuum multiscale model of plasticity with quantified uncertainty

Multiscale models of materials, consisting of upscaling discrete simulations to continuum models, are unique in their capability to simulate complex materials behavior. The fundamental limitation in multiscale models is the presence of uncertainty in the computational predictions delivered by them. In this work, a sequential multiscale model has been developed, incorporating discrete dislocation dynamics (DDD) simulations and a strain gradient plasticity (SGP) model to predict the size effect in plastic deformations of metallic micro-pillars. The DDD simulations include uniaxial compression of micro-pillars with different sizes and over a wide range of initial dislocation densities and spatial distributions of dislocations. An SGP model is employed at the continuum level that accounts for the size-dependency of flow stress and hardening rate. Sequences of uncertainty analyses have been performed to assess the predictive capability of the multiscale model. The variance-based global sensitivity analysis determines the effect of parameter uncertainty on the SGP model prediction. The multiscale model is then constructed by calibrating the continuum model using the data furnished by the DDD simulations. A Bayesian calibration method is implemented to quantify the uncertainty due to microstructural randomness in discrete dislocation simulations (density and spatial distribution of dislocations) on the macroscopic continuum model prediction (size effect in plastic deformation). Here, the outcomes of this study indicate that the discrete-continuum multiscale model can accurately simulate the plastic deformation of micro-pillars, despite the significant uncertainty in the DDD results. Additionally, depending on the macroscopic features represented by the DDD simulations, the SGP model can reliably predict the size effect in plasticity responses of the micropillars with below 10% of error.

36 MATERIALS SCIENCE↗

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗

Hydrocarbon Recovery from Williston Basin Shale and Mudrock Cores with Supercritical CO 2 : 2. Mechanisms That Control Oil Recovery Rates and CO 2 Permeation

In part 1 (10.1021/acs.energyfuels.9b01177), CO2 was used to recover oil from 51 source shale and reservoir mudrock cores collected from the Bakken Petroleum System in the Williston Basin. Oil hydrocarbon recoveries after 24 hexposures to CO2 at reservoir pressure and temperature were >94% for the reservoir mudrock cores and ranged from a few percent to as much as 80% for the source shales depending on the well location. In part 2, the experimental parameters that control oil recoveries were investigated, and the results show that exposed rock surface areas and CO2 contact times are primary factors controlling oil recovery at reservoir temperature and pressure. Compared to the 11.2 mm diameter rods used in part 1,increasing the surface area by splitting the rods into a stack of 2−3 mm thick “coins” doubled the oil recovery rates from both reservoir and source rocks, while further increasing the surface area using ground and sieved 1−3.4 mm samples at least tripledrecovery rates. In addition, extending the exposure time for the 11.2 mm diameter source shale rods to 96 h yielded nearly complete oil recoveries (as did the 24 h exposures with the 1−3.4 mm samples), indicating that pore spaces in the source shales as well as the more permeable reservoir rocks can be accessed by the CO2 and the associated oil hydrocarbons recovered. Higher CO 2 pressures yielded higher oil recoveries from both the reservoir rocks and source shales regardless of whether the exposure pressure was at or a little below, somewhat above, or substantially above the minimum miscibility pressure (MMP).Laboratory experiments also demonstrated that crude oil recoveries are based primarily on the ability of the CO 2 to penetrate the rock matrix and dissolve the oil hydrocarbons via hydrocarbon vaporization into the CO 2 phase rather than bulk physical processes (e.g., swelling, lowered viscosity, the physical “sweeping” effect) that are important in conventional CO 2 floods. Lighter hydrocarbons (e.g., C7 to C14) were recovered at much faster rates than heavier hydrocarbons (e.g., >C20) from all rock samples and geometries as well as for all CO2 pressures tested, as might be expected because lighter hydrocarbons have both higher diffusion coefficients and higher solubilities in CO 2 than heavier hydrocarbons. Limiting the amount of CO2 had little or no effect on the recovery rates of the lighter hydrocarbons but greatly reduced those of the heavier hydrocarbons. These two observations are consistent with a concentration-gradient-driven diffusion recovery mechanism. Laboratory results suggest that oil recovery in the Bakken play with CO2 will be enhanced by longer soak times, larger exposed rock surface areas, and higher pressures.

Energy & Fuels↗