Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fast sweeping method”

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.

50 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A Fast Solver for Implicit Integration of the Vlasov--Poisson System in the Eulerian Framework

In this paper, we present a domain decomposition algorithm to accelerate the solution of Eulerian-type discretizations of the linear, steady-state Vlasov equation. The steady-state solver then forms a key component in the implementation of fully implicit or nearly fully implicit temporal integrators for the nonlinear Vlasov--Poisson system. The solver relies on a particular decomposition of phase space that enables the use of sweeping techniques commonly used in radiation transport applications. The original linear system for the phase space unknowns is then replaced by a smaller linear system involving only unknowns on the boundary between subdomains, which can then be solved efficiently with Krylov methods such as GMRES. Steady-state solves are combined to form an implicit Runge--Kutta time integrator, and the Vlasov equation is coupled self-consistently to the Poisson equation via a linearized procedure or a nonlinear fixed-point method for the electric field. Finally, numerical results for standard test problems demonstrate the efficiency of the domain decomposition approach when compared to the direct application of an iterative solver to the original linear system.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

Power System Optimization Codes Modified

A major modification of and addition to existing Closed Brayton Cycle (CBC) space power system optimization codes was completed. These modifications relate to the global minimum mass search driver programs containing three nested iteration loops comprising iterations on cycle temperature ratio, and three separate pressure ratio iteration loops--one loop for maximizing thermodynamic efficiency, one for minimizing radiator area, and a final loop for minimizing overall power system mass. Using the method of steepest ascent, the code sweeps through the pressure ratio space repeatedly, each time with smaller iteration step sizes, so that the three optimum pressure ratios can be obtained to any desired accuracy for each of the objective functions referred to above (i.e., maximum thermodynamic efficiency, minimum radiator area, and minimum system mass). Two separate options for the power system heat source are available: 1. A nuclear fission reactor can be used. It is provided with a radiation shield 1. (composed of a lithium hydride (LiH) neutron shield and tungsten (W) gamma shield). Suboptions can be used to select the type of reactor (i.e., fast spectrum liquid metal cooled or epithermal high-temperature gas reactor (HTGR)). 2. A solar heat source can be used. This option includes a parabolic concentrator and heat receiver for raising the temperature of the recirculating working fluid. A useful feature of the code modifications is that key cycle parameters are displayed, including the overall system specific mass in kilograms per kilowatt and the system specific power in watts per kilogram, as the results for each temperature ratio are computed. As the minimum mass temperature ratio is encountered, a message is printed out. Several levels of detailed information on cycle state points, subsystem mass results, and radiator temperature profiles are stored for this temperature ratio condition and can be displayed or printed by users.

Juhasz, Albert J.↗
DOE OSTI2019

ZIF‐8 Membrane Separation Performance Tuning by Vapor Phase Ligand Treatment

Abstract Vapor phase ligand treatment (VPLT) of 2‐aminobenzimidazole (2abIm) for 2‐methylimidazole (2mIm) in ZIF‐8 membranes prepared by two different methods (LIPS: ligand induced permselectivation and RTD: rapid thermal deposition) results in a notable shift of the molecular level cut‐off to smaller molecules establishing selectivity improvements from ca. 1.8 to 5 for O 2 /N 2 ; 2.2 to 32 for CO 2 /CH 4 ; 2.4 to 24 for CO 2 /N 2 ; 4.8 to 140 for H 2 /CH 4 and 5.2 to 126 for H 2 /N 2 . Stable (based on a one‐week test) oxygen‐selective air separation performance at ambient temperature, 7 bar(a) feed, and 1 bar(a) sweep‐free permeate with a mixture separation factor of 4.5 and oxygen flux of 2.6×10 −3 mol m −2 s −1 is established. LIPS and RTD membranes exhibit fast and gradual evolution upon a 2abIm‐VPLT, respectively, reflecting differences in their thickness and microstructure. Functional reversibility is demonstrated by showing that the original permeation properties of the VPLT‐LIPS membranes can be recovered upon 2mIm‐VPLT.

Eum, Kiwon↗
DOE OSTI2019

ZIF-8 Membrane Separation Performance Tuning by Vapor Phase Ligand Treatment

Vapor phase ligand treatment (VPLT) of 2-aminobenzimidazole (2abIm) for 2-methylimidazole (2mIm) in ZIF-8 membranes prepared by two different methods (LIPS: ligand induced permselectivation and RTD: rapid thermal deposition) results in a notable shift of the molecular level cut-off to smaller molecules establishing selectivity improvements from ca. 1.8 to 5 for O 2 /N 2 ; 2.2 to 32 for CO 2 /CH 4 ; 2.4 to 24 for CO 2 /N 2 ; 4.8 to 140 for H 2 /CH 4 and 5.2 to 126 for H 2 /N 2 . Furthermore, stable (based on a one-week test) oxygen-selective air separation performance at ambient temperature, 7 bar(a) feed, and 1 bar(a) sweep-free permeate with a mixture separation factor of 4.5 and oxygen flux of 2.6×10 -3 mol m -2 s -1 is established. LIPS and RTD membranes exhibit fast and gradual evolution upon a 2abIm-VPLT, respectively, reflecting differences in their thickness and microstructure. Functional reversibility is demonstrated by showing that the original permeation properties of the VPLT-LIPS membranes can be recovered upon 2mIm-VPLT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1971

The stopping rate of negative cosmic-ray muons near sea level

A production rate of 0.065 + or - 0.003 Ar-37 atom/kg min of K-39 at 2-mwe depth below sea level was measured by sweeping argon from potassium solutions. This rate is unaffected by surrounding the solution by paraffin and is attributed to negative muon captures and the electromagnetic interaction of fast muons, and not to nucleonic cosmic ray component. The Ar-37 yield from K-39 by the stopping of negative muons in a muon beam of a synchrocyclotron was measured to be 8.5 + or - 1.7%. The stopping rate of negative cosmic ray muons at 2-mwe depth below sea level from these measurements and an estimated 17% electromagnetic production is 0.63 + or - 0.13 muon(-)/kg min. Previous measurements on the muon stopping rate vary by a factor of 5. Our value is slightly higher but is consistent with two previous high values. The sensitivity of the Ar-37 radiochemical method for the detection of muons is considerably higher than that of the previous radiochemical methods and could be used to measure the negative muon capture rates at greater depths.

Spannagel, G.↗
DOE OSTI2017

Model Reduction for Fast Assessment of Grid Impact of High Penetration PV

A key challenge in properly analyzing the impact of solar photovoltaics (PV) during planning, interconnection and operation in the distribution grid is computational complexity. Such studies require the use of Quasi-Steady-Time-Series (QSTS) power flow to capture the complex interaction of PV with voltage regulating equipment and protection devices and a large number of scenario studies are needed to capture location specific constraints and uncertainties. In addition, a typical distribution feeder can have thousands of nodes and there can be thousands of such feeders in a given utility system. To solve this problem, a new and systematic method for distribution network reduction is proposed. The objective is to reduce the size of distribution systems by more than 80% while preserving the accuracy of power flow results. To this end, an extended Ward-equivalent method has been proposed to create an equivalent circuit. Moreover, load and PV aggregation based on the backward sweep method has been developed to further simplify the feeder. The proposed method is tested with a modified IEEE 123-node distribution test system. The simulation results demonstrate the high efficiency and accuracy of the proposed approach for distribution network reduction.

reduced order systems, load flow, load modeling, v↗
DOE OSTI2016

Extreme high temperature redox kinetics in ceria: exploration of the transition from gas-phase to material-kinetic limitations

The redox kinetics of undoped ceria (CeO 2-δ ) are investigated by the electrical conductivity relaxation method in the oxygen partial pressure range of -4.3 ≤ log(pO 2 /atm) ≤ -2.0 at 1400 °C. It is demonstrated that extremely large gas flow rates, relative to the mass of the oxide, are required in order to overcome gas phase limitations and access the material kinetic properties. Using these high flow rate conditions, the surface reaction rate constant k chem is found to obey the correlation log(k chem /cm s -1 ) = (0.84 ± 0.02) × log(pO 2 /atm) - (0.99 ± 0.05) and increases with oxygen partial pressure. This increase contrasts the known behavior of the dominant defect species, oxygen vacancies and free electrons, which decrease in concentration with increasing oxygen partial pressure. For the sample geometries employed, diffusion was too fast to be detected. At low gas flow rates, the relaxation process becomes limited by the capacity of the sweep gas to supply/remove oxygen to/from the oxide. An analytical expression is derived for the relaxation in the gas-phase limited regime, and the result reveals an exponential decay profile, identical in form to that known for a surface reaction limited process. Thus, measurements under varied gas flow rates are required to differentiate between surface reaction limited and gas flow limited behavior.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Dual Ion Spectrometers and Their Calibration for the Fast Plasma Investigation on NASA's Magnetospheric Multiscale Mission

The scientific target of NASA's Magnetospheric Multiscale (MMS) mission is to study the fundamentally important phenomenon of magnetic reconnection. Theoretical models of this process predict a small size, on the order of hundred kilometers, for the ion diffusion region where ions are demagnetized at the dayside magnetopause. This region may typically sweep over the spacecraft at relatively high speeds of 50 km/s, requiring the fast plasma investigation (FPI) instrument suite to have an extremely high time resolution for measurements of the 3D particle distribution functions. As part of the FPI on MMS, the 16 dual ion spectrometers (DIS) will provide fast (150 ms) 3D ion velocity distributions, from 10 to 30,000 eV/q, by combining the measurements from four dual spectrometers on each of four MMS spacecraft. For any multispacecraft mission, the response uniformity among the spectrometer set assumes an enhanced importance. Due to these demanding instrument requirements and the effort of calibrating more than 32 sensors (16 × 2) within a tight schedule, a highly systematic and precise calibration was required for measurement repeatability. To illustrate how this challenge was met, a brief overview of the FPI DIS was presented with a detailed discussion of the calibration method of approach and implementation. Finally, a discussion of DIS performance results, their unit-to-unit variation, and the lessons learned from this calibration effort are presented.

Coffey, V. N.↗
DOE OSTI2016

Analysis of the transient response of nuclear spins in GaAs with/without nuclear magnetic resonance

As an alternative to studying the steady-state responses of nuclear spins in solid state systems, working within a transient-state framework can reveal interesting phenomena. The response of nuclear spins in GaAs to a changing magnetic field was analyzed based on the time evolution of nuclear spin temperature. Simulation results well reproduced our experimental results for the transient oblique Hanle signals observed in an all-electrical spin injection device. The analysis showed that the so called dynamic nuclear polarization can be treated as a cooling tool for the nuclear spins: It works as a provider to exchange spin angular momentum between polarized electron spins and nuclear spins through the hyperfine interaction, leading to an increase in the nuclear polarization. In addition, a time-delay of the nuclear spin temperature with a fast sweep of the external magnetic field produces a possible transient state for the nuclear spin polarization. On the other hand, the nuclear magnetic resonance acts as a heating tool for a nuclear spin system. This causes the nuclear spin temperature to jump to infinity: i.e., the average nuclear spins along with the nuclear field vanish at resonant fields of {sup 75}As, {sup 69}Ga and {sup 71}Ga, showing an interesting step-dip structure in the oblique Hanle signals. These analyses provide a quantitative understanding of nuclear spin dynamics in semiconductors for application in future computation processing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2015

In Flight Calibration of the Magnetospheric Multiscale Mission Fast Plasma Investigation

The Fast Plasma Investigation (FPI) on the Magnetospheric Multiscale mission (MMS) combines data from eight spectrometers, each with four deflection states, into a single map of the sky. Any systematic discontinuity, artifact, noise source, etc. present in this map may be incorrectly interpreted as legitimate data and incorrect conclusions reached. For this reason it is desirable to have all spectrometers return the same output for a given input, and for this output to be low in noise sources or other errors. While many missions use statistical analyses of data to calibrate instruments in flight, this process is difficult with FPI for two reasons: 1. Only a small fraction of high resolution data is downloaded to the ground due to bandwidth limitations and 2: The data that is downloaded is, by definition, scientifically interesting and therefore not ideal for calibration. FPI uses a suite of new tools to calibrate in flight. A new method for detection system ground calibration has been developed involving sweeping the detection threshold to fully define the pulse height distribution. This method has now been extended for use in flight as a means to calibrate MCP voltage and threshold (together forming the operating point) of the Dual Electron Spectrometers (DES) and Dual Ion Spectrometers (DIS). A method of comparing higher energy data (which has low fractional voltage error) to lower energy data (which has a higher fractional voltage error) will be used to calibrate the high voltage outputs. Finally, a comparison of pitch angle distributions will be used to find remaining discrepancies among sensors.

Detection-system↗
NASA NTRS2015

In Flight Calibration of the Magnetospheric Multiscale Mission Fast Plasma Investigation

The Fast Plasma Investigation (FPI) on the Magnetospheric Multiscale mission (MMS) combines data from eight spectrometers, each with four deflection states, into a single map of the sky. Any systematic discontinuity, artifact, noise source, etc. present in this map may be incorrectly interpreted as legitimate data and incorrect conclusions reached. For this reason it is desirable to have all spectrometers return the same output for a given input, and for this output to be low in noise sources or other errors. While many missions use statistical analyses of data to calibrate instruments in flight, this process is insufficient with FPI for two reasons: 1. Only a small fraction of high resolution data is downloaded to the ground due to bandwidth limitations and 2: The data that is downloaded is, by definition, scientifically interesting and therefore not ideal for calibration. FPI uses a suite of new tools to calibrate in flight. A new method for detection system ground calibration has been developed involving sweeping the detection threshold to fully define the pulse height distribution. This method has now been extended for use in flight as a means to calibrate MCP voltage and threshold (together forming the operating point) of the Dual Electron Spectrometers (DES) and Dual Ion Spectrometers (DIS). A method of comparing higher energy data (which has low fractional voltage error) to lower energy data (which has a higher fractional voltage error) will be used to calibrate the high voltage outputs. Finally, a comparison of pitch angle distributions will be used to find remaining discrepancies among sensors.

Heliophysics↗
DOE OSTI2015

Influence of Microstructure and Surface Activation of Dual-Phase Membrane Ce 0.8 Gd 0.2 O 2-δ -FeCo 2 O 4 on Oxygen Permeation

In dual-phase oxygen transport membranes we noticed that there is fast-growing interest in research for oxyfuel combustion process application. One such potential candidate is CGO-FCO (60wt% Ce 0.8 Gd 0.2 O 2-δ -40wt% FeCo 2 O4) identified to provide good oxygen permeation flux with substantial stability in harsh atmosphere. Dense CGO-FCO membranes of 1mm thickness were fabricated by sintering dry pellets pressed from powders synthesized by one-pot method (modified Pechini process) at 1200 degrees C for 10h. Microstructure analysis indicates presence of a third orthorhombic perovskite phase in the sintered composite. We also identified that the spinel phase tends to form an oxygen deficient phase at the grain boundary of spinel and CGO phases. Surface exchange limitation of the membranes was overcome by La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ (LSCF) porous layer coating over the composite. Moreover, the oxygen permeation flux of the CGO-FCO screen printed with a porous layer of 10mthick LSCF is 0.11mL/cm 2 per minute at 850 degrees C with argon as sweep and air as feed gas at the rates of 50 and 250mL/min.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Diffusion Synthetic Acceleration for Heterogeneous Domains, Compatible with Voids

A standard approach to solving the S N transport equations is to use source iteration with diffusion synthetic acceleration (DSA). Although this approach is widely used and effective on many problems, there remain some practical issues with DSA preconditioning, particularly on highly heterogeneous domains. For large-scale parallel simulation, it is critical that both (a) preconditioned source iteration converges rapidly and (b) the action of the DSA preconditioner can be applied using fast, scalable solvers, such as algebraic multigrid (AMG). For heterogeneous domains, these two interests can be at odds. In particular, there exist DSA diffusion discretizations that can be solved rapidly using AMG, but they do not always yield robust/fast convergence of the larger source iteration. Conversely, there exist robust DSA discretizations where source iteration converges rapidly on difficult heterogeneous problems, but fast parallel solvers like AMG tend to struggle applying the action of such operators. Moreover, very few current methods for the solution of deterministic transport are compatible with voids. This paper develops a new heterogeneous DSA preconditioner based on only preconditioning the optically thick subdomains. The resulting method proves robust on a variety of heterogeneous transport problems, including a linearized hohlraum mesh related to inertial confinement fusion. Moreover, the action of the preconditioner is easily computed using O(1) AMG iterations, convergence of the transport iteration typically requires 2 to 5× fewer iterations than current state-of-the-art “full” DSA, and the proposed method is trivially compatible with voids. On the hohlraum problem, rapid convergence is obtained by preconditioning less than 3% of the mesh elements with five to ten AMG iterations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska’s North Slope – Polymer Injection Performance

Alaska North Slope (ANS) holds an estimated 20-30+ billion barrel heavy oil resources, yet the development pace has been very slow due to high development costs and low oil recovery using conventional waterflood and EOR methods. The objective of this pilot is to perform a field experiment to validate the use of an advanced polymer flooding technology to unlock the vast heavy oil resources on ANS. The advanced polymer flooding technology combines polymer flooding, low salinity water flooding, horizontal wells, and if necessary, injection conformance control treatments into one integrated process to significantly improve oil recovery from heavy oil reservoirs. Two pairs of horizontal injection and production wells have been deployed in an isolated fault block of the Schrader Bluff heavy oil reservoir at the Milne Point Field to conduct a polymer flood pilot. The pilot will acquire scientific knowledge and field performance data to optimize polymer flood design in the Schrader Bluff heavy oil reservoirs on ANS. Polymer injection started on August 28, 2018 using a custom made polymer blending and pumping unit. This paper focuses on the facility setup and polymer injection performance into the horizontal injectors drilled and completed in the Schrader Bluff heavy oil reservoir. Partially hydrolyzed polyacrylamide (HPAM) polymer was selected and the initial target viscosity was set at 45 centipoise. Polymer injection rate was set at 2200 bbl/day for one injector (J-23A) and 1200 bbl/day for the other (J-24A) based on production voidage. Injection pressure was controlled at below fracture pressure to prevent fracturing the reservoir and causing fast breakthroughs. Step rate and pressure falloff tests indicate that short term polymer injectivity is similar to water injectivity, which means that injectivity is mostly controlled by fluid mobility deep in the reservoir rather than that in the vicinity of the injection wellbore. Long term injection data indicate that polymer injectivity has been decreasing in both injectors as the reservoir is filled by polymer. No polymer has been observed in the production stream 7 months after the start of polymer injection compared with a 3-month breakthrough time with waterflood. This indicates that polymer significantly delays breakthrough time which will lead to increased sweep efficiency.

Ning, Samson↗