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 343 records · Page 19

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗

Stellar occultation spikes as probes of atmospheric structure and composition

The characteristics of spikes observed in occultation light curves of Beta Scorpii by Jupiter are discussed in terms of the gravity-gradient model. The occultation of Beta Sco by Jupiter on May 13, 1971, is reviewed, and the gravity-gradient model is defined as an isothermal atmosphere of constant composition in which the refractivity is a function only of the radial coordinate from the center of refraction, which is assumed to lie parallel to the local gravity gradient. The derivation of the occultation light curve in terms of the atmosphere, the angular diameter of the occulted star, and the occultation geometry is outlined. It is shown that analysis of the light-curve spikes can yield the He/H2 concentration ratio in a well-mixed atmosphere, information on fine-scale atmospheric structure, high-resolution images of the occulted star, and information on ray crossing. Observational limits are placed on the magnitude of horizontal refractivity gradients, and it is concluded that the spikes are the result of local atmospheric density variations: atmospheric layers, density waves, or turbulence.

Elliot, J. L.↗

High-sensitivity silicon ultraviolet p+-i-n avalanche photodiode using ultra-shallow boron gradient doping

Photo detection of ultraviolet (UV) light remains a challenge since the penetration depth of UV light is limited to the nanometer scale. Therefore, the doping profile and electric field in the top nanometer range of the photo detection devices become critical. Traditional UV photodetectors usually use a constant doping profile near the semiconductor surface, resulting in a negligible electric field, which limits the photo-generated carrier collection efficiency of the photodetector. Here, we demonstrate, via the use of an optimized gradient boron doping technique, that the carrier collection efficiency and photo responsivity under the UV wavelength region have been enhanced. Moreover, the ultrathin p+-i-n junction shows an avalanche gain of 2800 and an ultra-low junction capacitance (sub pico-farad), indicating potential applications in the low timing jitter single photon detection area.

36 MATERIALS SCIENCE↗

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Oxygen Electrocatalysis on Mixed-Metal Oxides/Oxyhydroxides: From Fundamentals to Membrane Electrolyzer Technology

Catalyzing the oxygen evolution reaction (OER) is important for key energy-storage technologies, particularly water electrolysis and photoelectrolysis for hydrogen fuel production. Under neutral-to-alkaline conditions, first-row transition-metal oxides/(oxy)hydroxides are the fastest-known OER catalysts and have been the subject of intense study for the past decade. Critical to their high performance is the intentional or accidental addition of Fe to Ni/Co oxides that convert to layered (oxy)hydroxide structures during the OER. Unraveling the role that Fe plays in the catalysis and the molecular identity of the true “active site” has proved challenging, however, due to the dynamics of the host structure and absorbed Fe sites as well as the diversity of local structures in these disordered active phases. In this Account, we highlight our work to understand the role of Fe in Ni/Co (oxy)hydroxide OER catalysts. We first discuss how we characterize the intrinsic activity of the first-row transition-metal (oxy)hydroxide catalysts as thin films by accounting for the contributions of the catalyst-layer thickness (mass loading) and electrical conductivity as well as the underlying substrate’s chemical interactions with the catalyst and the presence of Fe species in the electrolyte. We show how Fe-doped Ni/Co (oxy)hydroxides restructure during catalysis, absorb/desorb Fe, and in some cases degrade or regenerate their activity during electrochemical testing. We highlight the relevant techniques and procedures that allowed us to better understand the role of Fe in activating other first-row transition metals for OER. We find several modes of Fe incorporation in Ni/Co (oxy)hydroxides and show how those modes correlate with activity and durability. We also discuss how this understanding informs the incorporation of earth-abundant transition-metal OER catalysts in anion-exchange-membrane water electrolyzers (AEMWE) that provide a locally basic anode environment but run on pure water and have advantages over the more-developed proton-exchange-membrane water electrolyzers (PEMWE) that use platinum-group-metal (PGM) catalysts. We outline the key issues of introducing Fe-doped Ni/Co (oxy)hydroxide catalysts at the anode of the AEMWE, such as the oxidative processes triggered by Fe species traveling through the polymer membrane, pH-gradient effects on the catalyst stability, and possibly limited catalyst utilization in the compressed stack configuration. We also suggest possible mitigation strategies for these issues. Lastly, we summarize remaining challenges including the long-term stability of Fe-doped Ni/Co (oxy)hydroxides under OER conditions and the lack of accurate models of the dynamic active surface that hinder our understanding of, and thus ability to design, these catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Soil gas probes for monitoring trace gas messengers of microbial activity

Abstract Soil microbes vigorously produce and consume gases that reflect active soil biogeochemical processes. Soil gas measurements are therefore a powerful tool to monitor microbial activity. Yet, the majority of soil gases lack non-disruptive subsurface measurement methods at spatiotemporal scales relevant to microbial processes and soil structure. To address this need, we developed a soil gas sampling system that uses novel diffusive soil probes and sample transfer approaches for high-resolution sampling from discrete subsurface regions. Probe sampling requires transferring soil gas samples to above-ground gas analyzers where concentrations and isotopologues are measured. Obtaining representative soil gas samples has historically required balancing disruption to soil gas composition with measurement frequency and analyzer volume demand. These considerations have limited attempts to quantify trace gas spatial concentration gradients and heterogeneity at scales relevant to the soil microbiome. Here, we describe our new flexible diffusive probe sampling system integrated with a modified, reduced volume trace gas analyzer and demonstrate its application for subsurface monitoring of biogeochemical cycling of nitrous oxide (N 2 O) and its site-specific isotopologues, methane, carbon dioxide, and nitric oxide in controlled soil columns. The sampling system observed reproducible responses of soil gas concentrations to manipulations of soil nutrients and redox state, providing a new window into the microbial response to these key environmental forcings. Using site-specific N 2 O isotopologues as indicators of microbial processes, we constrain the dynamics of in situ microbial activity. Unlocking trace gas messengers of microbial activity will complement -omics approaches, challenge subsurface models, and improve understanding of soil heterogeneity to disentangle interactive processes in the subsurface biome.

04 OIL SHALES AND TAR SANDS↗

Communication: A new class of non-empirical explicit density functionals on the third rung of Jacob’s ladder

We construct an orbital-free non-empirical meta-generalized gradient approximation (GGA) functional, which depends explicitly on density through the density overlap regions indicator [P. de Silva and C. Corminboeuf, J. Chem. Theory Comput. 10, 3745 (2014)]. The functional does not depend on either the kinetic energy density or the density Laplacian; therefore, it opens a new class of meta-GGA functionals. By construction, our meta-GGA yields exact exchange and correlation energy for the hydrogen atom and recovers the second order gradient expansion for exchange in the slowly varying limit. We show that for molecular systems, overall performance is better than non-empirical GGAs. For atomization energies, performance is on par with revTPSS, without any dependence on Kohn-Sham orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Destabilization of counter-propagating Alfvénic instabilities by tangential, co-current neutral beam injection

Injection of high-energy neutrals is a common tool to heat the plasma and drive current non-inductively in fusion devices. Once neutrals ionize, the resulting energetic particles can drive instabilities that are detrimental for the performance and the predictability of plasma discharges. A broad deposition profile of neutrals from neutral beam injection, e.g. by aiming the beam tangentially on the outboard midplane (i.e. off-axis), is often assumed to limit those undesired effects by reducing the radial gradient of the EP density, thus reducing the drive for instabilities. However, this paper presents new evidence that tangential neutral beam injection, including off-axis injection near the plasma mid-radius, can also lead to undesired effects such as the destabilization of Alfvénic instabilities. Time-dependent analysis with the TRANSP code indicates that instabilities are driven by a combination of radial and energy gradients in the distribution function of the energetic particles. Finally, the mechanisms for wave-particle interaction revealed by the energetic particle phase space resolved analysis are the basis to identify strategies to mitigate or suppress the observed instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-F-T-449: Dosimetric Comparison of Acuros XB, Adaptive Convolve in Intensity Modulated Radiotherapy for Head and Neck Cancer

Purpose: There have been several publications focusing on dose calculation in lung for a new dose calculation algorithm of Acuros XB (AXB). AXB could contribute to dose calculation for high-density media for bone and dental prosthesis rather than in lung. We compared the dosimetric performance of AXB, Adaptive Convolve (AC) in head and neck IMRT plans. Methods: In a phantom study, the difference in depth profile between AXB and AC was evaluated using Kodak EDR2 film sandwiched with tough water phantoms. 6 MV x-ray using the TrueBeam was irradiated. In a patient study, 20 head and neck IMRT plans had been clinically approved in Pinnacle3 and were transferred to Eclipse. Dose distribution was recalculated using AXB in Eclipse while maintaining AC-calculated monitor units and MLC sequence planned in Pinnacle. Subsequently, both the dose-volumetric data obtained using the two different calculation algorithms were compared. Results: The results in the phantom evaluation for the shallow area ahead of the build-up region shows over-dose for AXB and under-dose for AC, respectively. In the patient plans, AXB shows more hot spots especially around the high-density media than AC in terms of PTV (Max difference: 4.0%) and OAR (Max. difference: 1.9%). Compared to AC, there were larger dose deviations in steep dose gradient region and higher skin-dose. Conclusion: In head and neck IMRT plans, AXB and AC show different dosimetric performance for the regions inside the target volume around high-density media, steep dose gradient regions and skin-surface. There are limitations in skin-dose and complex anatomic condition using even inhomogeneous anthropomorphic phantom Thus, there is the potential for an increase of hot-spot in AXB, and an underestimation of dose in substance boundaries and skin regions in AC.

60 APPLIED LIFE SCIENCES↗

Thermal / structural analysis of the HB 650 thermal shield

Fermilab’s PIP-II project’s superconducting linear accelerator will drive the next generation of particle accelerators through a revolution in beam intensity. Key to beam intensification are the high efficiencies of niobium-tin superconducting radiofrequency (SRF) cavities operating at cryogenic temperatures near 5 Kelvin. A multifaceted approach is employed to achieve and maintain the extreme temperature. Vacuum provides the first barrier to thermal convection. Physical thermal intercept zones further isolate the exterior shell vacuum vessel at 300 Kelvin from the supercooled beamline. The first thermal intercept lies just inside the exterior vacuum vessel forming the 40 Kelvin zone. This 40K thermal shield is the focus of this investigation. In preparation for operation, the thermal shield is cooled from 300K to 40K in a slow process over 2-1/2 days. The factors limiting the cooling rate are a high thermal gradient which produces thermal strain and high mechanical stress. D ecreasing the time required to reach operational temperatures is the desired research goal. Physical testing data and finite element method (FEM) computer models from the Single Spoke Resonator – 1 (SSR1) and Linac Coherent Light Source (LCLS-II) cryomodules will be evaluated to determine needed design changes. The thermal shield design for the High Beta 650 MHz (HB650) will be modified and analyzed using FEM computer models to determine maximum cooling rate with a target of 50 thermal cycles of life. Maximum allowable stress will be determined through low cycle fatigue calculations.

43 PARTICLE ACCELERATORS↗

Structure development in silicon sheet by shaped crystallization

Models are presented for the development of a parallel twinned structure of the 110 plane type and the 112 line type in silicon ribbons. The models are believed to be mutually compatible and operable. The first model relates the requirements for super-cooling during crystallization. The existence of reentrant angles associated with the twin structure is proposed to provide a rough interface to reduce super-cooling. The spacing of the twins is proposed to be limited by the geometrical relationship between the thermal gradient in the liquid and the dimensions of the twinned crystallization front. The second model relates the thermal stress configuration to detail dislocation reactions which would be expected to develop twins. While a specific dislocation mechanism cannot yet be defined, a number of alternatives are presented. All of these various dislocation mechanisms would result in the observed crystalline configuration and the choice among them is not critical.

Leipold, M. H.↗

Theoretical analysis of heat flow in horizontal ribbon growth from a melt

A theoretical heat flow analysis for horizontalribbon growth is presented. Equations are derived relating pull speed, ribbon thickness, thermal gradient in the melt, and melt temperature for limiting cases of heat removal by radiation only and isothermal heat removal from the solid surface over the melt. Geometrical cross sections of the growth zone are shown to be triangular and nearly parabolic for the two respective cases. Theoretical pull speed for silicon ribbon 0.01 cm thick, where the loss of latent heat of fusion is by radiation to ambient temperature (300 K) only, is shown to be 1 cm/sec for horizontal growth extending 2 cm over the melt and with no heat conduction either to or from the melt. Further enhancement of ribbon growth rate by placing cooling blocks adjacent to the top surface is shown to be theoretically possible.

Zoutendyk, J. A.↗

Interaction effects in thermocapillary bubble migration

Two bubbles migrating along their line of centers under the influence of an imposed thermal gradient are considered in the quasi-static limit. Results are reported for representative values of the governing parameters.

Meyyappan, M.↗

Parabolic Dish Concentrator (PDC-1)

The design, construction, and installation of the Parabolic Dish Concentrator, Type 1 (PDC-1) has been one of the most significant JPL concentrator projects because of the knowledge gained about this type of concentrator and the development of design, testing, and analysis procedures which are applicable to all solar concentrator projects. The need for these procedures was more clearly understood during the testing period which started with the prototype panel evaluation and ended with the performance characterization of the completed concentrator. For each phase of the test program, practical test procedures were required and these procedures defined the mathematical analysis which was essential for successful concentrator development. The concentrator performance appears to be limited only by the distortions resulting from thermal gradients through the reflecting panels. Simple optical testing can be extremely effective, but comprehensive mechanical and optical analysis is essential for cost effective solar concentrator development.

Dennison, E. W.↗

Thermostructural Analysis of Carbon Cloth Phenolic Material Tested at the Laser Hardened Material Evaluation Laboratory

RSRM nozzle liner components have been analyzed and tested to explore the occurrence of anomalous material performance known as pocketing erosion. Primary physical factors that contribute to pocketing seem to include the geometric permeability, which governs pore pressure magnitudes and hence load, and carbon fiber high temperature tensile strength, which defines a material limiting capability. The study reports on the results of a coupled thermostructural finite element analysis of Carbon Cloth Phenolic (CCP) material tested at the Laser Hardened Material Evaluation Laboratory (the LHMEL facility). Modeled test configurations will be limited to the special case of where temperature gradients are oriented perpendicular to the composite material ply angle. Analyses were conducted using a transient, one-dimensional flow/thermal finite element code that models pore pressure and temperature distributions and in an explicitly coupled formulation, passes this information to a 2-dimensional finite element structural model for determination of the stress/deformation behavior of the orthotropic fiber/matrix CCP. Pore pressures are generated by thermal decomposition of the phenolic resin which evolve as a multi-component gas phase which is partially trapped in the porous microstructure of the composite. The nature of resultant pressures are described by using the Darcy relationships which have been modified to permit a multi-specie mass and momentum balance including water vapor condensation. Solution to the conjugate flow/thermal equations were performed using the SINDA code. Of particular importance to this problem was the implementation of a char and deformation state dependent (geometric) permeability as describing a first order interaction between the flow/thermal and structural models. Material property models are used to characterize the solid phase mechanical stiffness and failure. Structural calculations were performed using the ABAQUS code. Iterations were made between the two codes involving the dependent variables temperature, pressure and across-ply strain level. Model results comparisons are made for three different surface heat rates and dependent variable sensitivities discussed for the various cases.

Clayton, J. Louie↗

Addressing Technical and Policy Challenges for Remediation of Metals and Radionuclides within Complex Deep Vadose Zone Environments - 17028

Deep vadose zone contamination is a significant issue in many regions of the world. Although remediation of deep vadose zone contamination is important, direct contact with the contamination does not threaten human health or the environment due to its depth in the vadose zone. Yet, movement of contamination from the deep vadose zone to the groundwater creates the potential for exposure and risk to receptors. Therefore, the challenge of deep vadose zone remediation is to maintain a minimal rate of contaminant transport such that flux from the contaminated vadose zone is limited to protect groundwater resources and down-gradient receptors. We herein review major processes for viable remediation of deep vadose zone metal and radionuclides that form the practical constraints on remedial actions and provide an overview of how these principles are being practiced to address key contaminants at the U.S. Department of Energy Hanford Site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of Hydrogen Production Feasibility for a Light Water Reactor in the Midwest

Increased electricity production from renewable energy resources coupled with low natural gas prices has caused existing light-water reactor (LWRs) to experience ever diminishing returns from the electricity market. Via partnership among Idaho National Laboratory (INL), The National Renewable Energy Laboratory (NREL), Argonne National Laboratory (ANL), Exelon, and Fuel Cell Energy a techno-economic analysis on the viability of retrofitting existing pressurized water reactors (PWRs) to produce hydrogen (H2) via high temperature steam electrolysis (HTSE) has been conducted. Such integration would allow nuclear facilities to expand into additional markets that may be more profitable in the long-term. To accommodate such an integration, a detailed analysis of HTSE process operation, requirements, and flexibility was conducted. Technical analysis includes proposed nuclear system control scheme modifications to allow dynamic operation of the HTSE via both thermal and electrical connection to the nuclear plant. High fidelity Modelica simulations showcase the viability of such control schemes. However, due to limited knowledge of solid oxide fuel cell (SOFC) stack degradation due to thermal gradients; thermal cycling of the HTSE was not included. Therefore, the control schemes proposed are only utilized to re-distribute steam at startup and only the portion of electricity utilized in the electrolyzers is cycled. From the detailed analysis of the nuclear integration and the HTSE process design, a comprehensive cost estimation was conducted in the APEA and H2A models to elucidate capital and operational costs associated with the production, compression, and distribution of hydrogen from a nuclear facility. Alongside this costing analysis, market analyses were conducted by NREL and ANL on the electric and hydrogen markets respectively in the PJM interconnect. Utilizing the electricity data market projections in the PJM interconnect from NREL and hydrogen demand/pricing projections from ANL a five variable sweep over component capacities, discount rates, and hydrogen pricing was completed using the stochastic framework RAVEN (Risk Analysis Virtual ENVironment) through its resource dispatch plugin HERON (Heuristic Energy Resource Optimization Network). Each combination of variables was evaluated stochastically four times over a seventeen year timespan from 2026-2042 (inclusive) to determine the most economically advantageous solution. Following the five variable sweep an optimization was conducted about the best sweep point to determine optimal component sizing and setpoints. Results suggest positive gain is achievable at all projected hydrogen market pricing levels and at all discount rates. However, exact component sizing and net returns vary based on these values and if incorrect sizing is selected; major net losses can occur. The optimal result occurred with set points as follows: high hydrogen prices, the largest possible HTSE unit in the sweep set at 7.47 kg/sec [645.4 tpd], a contractual hydrogen market agreement 7.29 kg/sec [629.8 tpd], and a hydrogen storage size 115,188 kg. The analysis suggested that with a discount rate of 8%, a ΔNPV = 1.2 billion over the seventeen year span can be achieved. The results illuminate that by operating in multiple markets the nuclear facility is capable of avoiding sale of electricity during times of low electricity market pricing, while maintaining the ability to capitalize on the high electricity market pricing. It should be noted that the analysis conducted in this report is a differential cash flow analysis and as such does not present profit levels. Instead it highlights the net benefit between building and competing in the hydrogen market utilizing nuclear facilities and conducting business as usual in the electricity markets. Additionally, results presented in this report exhibit conservatism due to four key assumptions: Given the limited knowledge on SOFC stack degradation due to thermal gradients the high temperature steam electrolysis plant is not allowed to thermally cycle. This limitation decreases electricity generation capacity in the nuclear plant, reducing capacity payments and the maximum output to the electrical grid. This analysis considers building a separate hydrogen pipeline for use in the nuclear facility. Research is currently being conducted to determine if existing natural gas pipelines can accommodate direct hydrogen injection. If hydrogen can be integrated into existing natural gas pipelines then a capital savings of ~$19,000,000 per kg/sec [1 kg/sec = 86.4 tpd] of installed HTSE capacity could be realized. Fast transients in the electric grid are typically served by flywheels and electric batteries. With electrolysis only operation the nuclear plant can flex its electricity just as quickly. However, for this analysis, the ancillary services market has been neglected. Should the nuclear plant have the ability to operate in the ancillary services market during periods when it is producing hydrogen, an additional substantial income stream that would increase overall system profitability could be considered. Subsidies for non-emitting generation technologies are not introduced in the cash flow analysis, but instead are completed assuming zero subsidies. Profitability is based solely on electricity pricing. Electric pricing is assumed to never go negative. Instead, it is assumed overall system curtailment by the Independent System Operator (ISO) will occur in such a manner that pricing remains at or above zero. This assumption inflates the “business as usual” scenario, thus decreasing the ΔNPV calculated. If negative pricing scenarios continue to occur as they do today, one would expect a further increase in ΔNPV. Overall, results advocate that, through market diversification, nuclear power plants have the potential to substantially increase current profit margins, increase market penetration, and ultimately solidify their place as a mainstay in energy production.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Pioneer-10 measurements of the differential and integral cosmic-ray gradient between 1 and 3 AU

The cosmic-ray experiment on Pioneer has provided measurements of the differential radial gradient of cosmic-ray helium nuclei between 1 and 3 AU. Upper limits are quoted between 10 and 500 MeV/nuc which are everywhere 25%/AU and in some cases substantially smaller. The integral proton gradient ( 56 MeV) was also measured and found to have the following values: 2.4 + or - 0.3%/AU (1-2 AU), 7.6 + or - 0.7%/AU (2-3 AU).

Teegarden, B. J.↗