Search NASA⌕ Search

SEARCH · Search NASA

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.

15 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A second-order distributed memory parallel fast sweeping method for the Eikonal equation

The Eikonal equation is used to calculate wave propagation and distance fields, and due to its complexity requires numerical treatment for its solution. In this work, we present a second-order distributed memory parallel fast sweeping method. The second-order solution switches on a two-point stencil when two upwind points are available, and reverts to first-order otherwise. In all examples, the second-order method improves the solution over the first-order, allowing for significant savings in memory while achieving the same accuracy. Parallelization over distributed memory saw good weak scaling with optimal convergence. The computational time for second-order was approximately 2.5 times slower than first-order, where the largest amount of mesh points ran on 144 cores (512 GB) was ≈20 billion. The savings in memory from the second-order method combined with the distributed memory algorithm result in the ability to solve problems much larger than are possible with the serial first-order method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Exact signed distance fields using parallel Fast Sweeping Method

Signed distance fields are often used in multiphysics simulations to track material interfaces. We present a simple methodology based on the fast sweeping method to generate the exact signed distance from triangular meshes and linear paths on Cartesian grids. The methodology propagates the closest primitive to the boundary to the rest of the domain following the characteristics. A local upwind criterion is used to decide between the new and existing closest primitive at each grid point while capturing the correct sign of the global function. The methodology has optimal computational complexity and runs efficiently in distributed-memory architectures. We include 2D and 3D test cases along with a resolution study up to 0.512 trillion zones and 1,000 computer cores. The solution strategy can also be applied to other types of meshes or collections of primitives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

High Order Anchoring and Reinitialization of Level Set Function for Simulating Interface Motion

A second order interface anchoring method has been developed and used with fast sweeping algorithm for reinitialization of a level set function. The algebraic anchoring formulation ensures that the location of the actual interface is preserved, leading to better mass conservation property. It also provides high order accurate algebraic constraint for solving the Eikonal equation on a finite difference grid. Geometric properties of the interface such as surface normal and curvature are subsequently computed from the reinitialized distance function. Various analytical functions for modeling distortion in level set field are considered and accuracy of reinitialization is evaluated using first and second order anchoring schemes. It is also shown that accurate computation of interface curvature requires a high order anchor in addition to a high order fast sweeping method. Mass conservation property of reinitialization is also analyzed by considering test problems from literature including the classic Rider–Kothe single vortex problem. The formulation is suitable for efficient parallelization for both distributed memory and on-node shared memory parallel systems. In conclusion, scalability and performance of the reinitialization scheme on multiple architectures are demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Low-Field EMR Studies of Permalloy Films and Gratings

Flat and profile-modulated permalloy films have been studied by the electron magnetic resonance (EMR) method. In addition to ferromagnetic and spin-wave resonances, the structures demonstrate low-field EMR signals of an unusual shape, which form a hysteresis loop in sweeping fields. The low-field signals are attributed to a fast reorientation of magnetic domains. The low-field EMR behavior is comparable to the behavior in magneto-dependent photovoltage previously observed in the optical experiments. The shapes of the loops and typical values of the switching fields depend on the profile modulation parameters confirming the possibility of controlling magnetic properties and the coupling of magnetic and optical effects with nanoscale geometry.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Fast sweeping probe system for characterization of spokes in E × B discharges

In this work, we developed a rapidly swept, back-to-back 100 kHz Langmuir probe system using a tunable compensating network to study the temporal evolution of low frequency oscillations in Penning discharges, Hall Thrusters, and other E × B discharges. Experimental validation of the probe system is performed at low and high sweeping frequencies in a stable Penning discharge. Then application of the probe system to measurements of plasma parameter fluctuations in a low frequency (4 kHz) rotating spoke and an analysis method using the Hilbert transform are shown. We proposs that the rotating spoke oscillation conducts approximately a third of the cross field current in our Penning device.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Fast yaw optimization for wind plant wake steering using Boolean yaw angles

Abstract. In wind plants, turbines can be yawed into the wind to steer their wakes away from downstream turbines and achieve an overall increase in plant power. Mathematical optimization is typically used to determine the best yaw angles at which to operate the turbines in a plant. In this paper, we present a new heuristic to rapidly determine the yaw angles in a wind plant. In this method, we define the turbine yaw angles as Boolean – either yawed at a predefined angle or nonyawed – as opposed to the typical methods of defining yaw angles as continuous or with fine discretizations. We then optimize which turbines should be yawed with an algorithm that sweeps through the turbines from the most upstream to the most downstream. We demonstrate that our new Boolean optimization method can find turbine yaw angles that perform well compared to a traditionally used gradient-based optimizer for which the yaw angles are defined as continuous. There is less than 0.6 % difference in the optimized power between the two optimization methods for randomly placed turbine layouts and less than a 0.6 % difference in the optimal annual energy production between the two optimization methods for a real wind farm. Additionally, we show that our new method is much more computationally efficient than the traditional method. For plants with nonzero optimal yaw angles, our new method is generally able to solve for the turbine yaw angles 50–150 times faster, and in some extreme cases up to 500 times faster, than the traditional method.

17 WIND ENERGY↗
DOE OSTI2020

A fast implicit solver for semiconductor models in one space dimension

Several different approaches are proposed for solving fully implicit discretizations of a simplified Boltzmann-Poisson system with a linear relaxation-type collision kernel. This system models the evolution of free electrons in semiconductor devices under a low-density assumption. At each implicit time step, the discretized system is formulated as a fixed-point problem, which can then be solved with a variety of methods. A key algorithmic component in all the approaches considered here is a recently developed sweeping algorithm for Vlasov-Poisson systems. A synthetic acceleration scheme has been implemented to accelerate the convergence of iterative solvers by using the solution to a drift-diffusion equation as a preconditioner. The performance of four iterative solvers and their accelerated variants has been compared on problems modeling semiconductor devices with various electron mean-free-path.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Optimizing Process Conditions during Catalytic Fast Pyrolysis of Pine with Pt/TiO 2 —Improving the Viability of a Multiple-Fixed-Bed Configuration

Catalytic fast pyrolysis (CFP) has been identified as a promising pathway for the production of renewable fuels and co-products. However, continued technology development is needed to increase process efficiency and reduce process costs. This report builds upon previous research in which a bifunctional metal–acid Pt/TiO 2 catalyst was utilized in a fixed-bed reactor operated with co-fed H 2 to improve product yield and reduce coke generation compared to conventional CFP methods. Here, we report further process optimization, in which we achieved similar CFP oil carbon efficiency (>35%) and CFP oil oxygen content (<20 wt %) to our previous report while reducing catalyst and equipment costs by increasing time-on-stream between regenerations by 40–95% and decreasing required regeneration time by more than a factor of 2. These process improvements were achieved by conducting parameter sweeps to determine optimum conditions for CFP and regeneration with key variables including pyrolysis temperature, catalytic upgrading temperature, hydrogen partial pressure, and regeneration oxygen concentration. Coupled with comprehensive oil analyses, these data provide foundational insight into the deoxygenation and coking chemistries for CFP under realistic process conditions while also advancing the technology through applied engineering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Cell dark current–voltage from non-calibrated module electroluminescence image analysis

Here, we present a fast, accurate, and reliable method of obtaining cell dark current–voltage (I–V) curves from module electroluminescence (EL) images without requiring calibration or correction. For a pristine module, EL-derived dark I–V are compared to directly probed data for a variety of changing imaging parameters: camera sensor, lens, filter, aperture width, exposure time (level of sensor saturation), number of images used, and various combinations of these. Pristine modules and those experiencing different modes and degrees of degradation are examined. A recent study of modules using five different cell technologies demonstrates the practicality of our “EL sweep” technique for performance and degradation studies.

14 SOLAR ENERGY↗
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↗
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↗