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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE 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 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↗

NASA NTRS1966

The potential-sweep method - A theoretical analysis.

Potential sweep method for organic and adsorption analyses, showing relation between galvanostatic and fast potential sweep transients, noting mathematical analysis of current-potential transients

ELECTROCHEMISTRY↗

NASA NTRS1985

Noise constraints effecting optimal propeller designs

A preliminary design tool for advanced propellers was developed combining a fast vortex lattice aerodynamic analysis, a fast subsonic point source noise analysis, and an optimization scheme using a conjugate directions method. Twist, chord and sweep distributions are optimized to simultaneously improve both the aerodynamic performance and the noise observed at a fixed relative position. The optimal noise/performance tradeoffs for straight and advanced concept blades are presented. The techniques used include increasing the blade number, blade sweep, reducing the rotational speed, shifting the spanwise loading and diameter changes.

Miller, C. J.↗

NASA NTRS1985

Noise constraints effecting optimal propeller designs

Miller, C. J.↗

NASA NTRS1984

Special Test Methods for Batteries

Various methods are described for measuring heat generation in primary and secondary batteries as well as the specific heat of batteries and cell thermal conductance. Problems associated with determining heat generation in large batteries are examined. Special attention is given to monitoring temperature gradients in nickel cadmium cells, the use of auxiliary electrodes for conducting tests on battery charge control, evaluating the linear sweep of current from charge to discharge, and determining zero current voltage. The fast transient behavior of batteries in the microsecond range, and the electrical conductance of nickel sinters in the thickness direction are also considered. Mechanical problems experienced in the vibration of Ni-Cd batteries and tests to simulate cyclic fatigue of the steel table connecting the plates to the comb are considered. Methods of defining the distribution of forces when cells are compressed during battery packaging are also explored.

Gross, S.↗

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↗

NASA NTRS1979

Multi-level adaptive finite element methods. 1: Variation problems

A general numerical strategy for solving partial differential equations and other functional problems by cycling between coarser and finer levels of discretization is described. Optimal discretization schemes are provided together with very fast general solvers. It is described in terms of finite element discretizations of general nonlinear minimization problems. The basic processes (relaxation sweeps, fine-grid-to-coarse-grid transfers of residuals, coarse-to-fine interpolations of corrections) are directly and naturally determined by the objective functional and the sequence of approximation spaces. The natural processes, however, are not always optimal. Concrete examples are given and some new techniques are reviewed. Including the local truncation extrapolation and a multilevel procedure for inexpensively solving chains of many boundary value problems, such as those arising in the solution of time-dependent problems.

Brandt, A.↗

NASA NTRS1986

Drag-reduction characteristics of aft-swept wing tips

Results are presented of a study conducted to analyze the potential aerodynamic improvements which can be obtained with aft-swept wing tips at subsonic speeds. A review and summary of previous swept-tip applications and studies are given. Additionally, the through evolution optimized tail-fin shapes of fast-cruising aquatic animals are examined as well as the wing planform shapes of sea birds. These superior long-distance travellers both display lifting-surface planform shapes with increased leading-edge sweep towards the tip. These observations and analytical results obtained with a nonlinear surface-panel method demonstrate that substantial induced-drag improvements can be obtained through the adoption of planform shapes with highly-tapered aft-swept tips.

Van Dam, C. P.↗

NASA NTRS1990

Novel laser gain and time-resolved FTIR studies of photochemistry

Several techniques are discussed which can be used to explore laboratory photochemical processes and kinetics relevant to planetary atmospheres; these include time-resolved laser gain-versus-absorption spectroscopy and time-resolved Fourier transform infrared (FTIR) emission studies. The laser gain-versus-absorption method employed tunable diode and F-center lasers to determine the yields of excited photofragments and their kinetics. The time-resolved FTIR technique synchronizes the sweep of a commercial FTIR with a pulsed source of light to obtain emission spectra of novel transient species in the infrared. These methods are presently being employed to investigate molecular photodissociation, the yields of excited states of fragments, their subsequent reaction kinetics, Doppler velocity distributions, and velocity-changing collisions of translationally fast atoms. Such techniques may be employed in future investigations of planetary atmospheres, for example to study polycyclic aromatic hydrocarbons related to cometary emissions, to analyze acetylene decomposition products and reactions, and to determine spectral features in the near infrared and infrared wavelength regions for planetary molecules and clusters.

Leone, Stephen R.↗

NASA NTRS1990

Implicit methods for the Navier-Stokes equations

Numerical solutions of the Navier-Stokes equations using explicit schemes can be obtained at the expense of efficiency. Conventional implicit methods which often achieve fast convergence rates suffer high cost per iteration. A new implicit scheme based on lower-upper factorization and symmetric Gauss-Seidel relaxation offers very low cost per iteration as well as fast convergence. High efficiency is achieved by accomplishing the complete vectorizability of the algorithm on oblique planes of sweep in three dimensions.

Yoon, S.↗

NASA NTRS1980

Data acquisition techniques for exploiting the uniqueness of the time-of-flight mass spectrometer: Application to sampling pulsed gas systems

Mass spectra are produced in most mass spectrometers by sweeping some parameter within the instrument as the sampled gases flow into the ion source. It is evident that any fluctuation in the gas during the sweep (mass scan) of the instrument causes the output spectrum to be skewed in its mass peak intensities. The time of flight mass spectrometer (TOFMS) with its fast, repetitive mode of operation produces spectra without skewing or varying instrument parameters and because all ion species are ejected from the ion source simultaneously, the spectra are inherently not skewed despite rapidly changing gas pressure or composition in the source. Methods of exploiting this feature by utilizing fast digital data acquisition systems, such as transient recorders and signal averagers which are commercially available are described. Applications of this technique are presented including TOFMS sampling of vapors produced by both pulsed and continuous laser heating of materials.

Lincoln, K. A.↗

NASA NTRS1990

Plane waves and structures in turbulent channel flow

A direct simulation of turbulent flow in a channel is analyzed by the method of empirical eigenfunctions (Karhunen-Loeve procedure, proper orthogonal decomposition). This analysis reveals the presence of propagating plane waves in the turbulent flow. The velocity of propagation is determined by the flow velocity at the location of maximal Reynolds stress. The analysis further suggests that the interaction of these waves appears to be essential to the local production of turbulence via bursting or sweeping events in the turbulent boundary layer, with the additional suggestion that the fast acting plane waves act as triggers.

Sirovich, L.↗

NASA NTRS1996

Towards and FVE-FAC Method for Determining Thermocapillary Effects on Weld Pool Shape

Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.

Canright, David↗

NASA NTRS1989

Operator induced multigrid algorithms using semirefinement

A variant of multigrid, based on zebra relaxation, and a new family of restriction/prolongation operators is described. Using zebra relaxation in combination with an operator-induced prolongation leads to fast convergence, since the coarse grid can correct all error components. The resulting algorithms are not only fast, but are also robust, in the sense that the convergence rate is insensitive to the mesh aspect ratio. This is true even though line relaxation is performed in only one direction. Multigrid becomes a direct method if an operator-induced prolongation is used, together with the induced coarse grid operators. Unfortunately, this approach leads to stencils which double in size on each coarser grid. The use of an implicit three point restriction can be used to factor these large stencils, in order to retain the usual five or nine point stencils, while still achieving fast convergence. This algorithm achieves a V-cycle convergence rate of 0.03 on Poisson's equation, using 1.5 zebra sweeps per level, while the convergence rate improves to 0.003 if optimal nine point stencils are used. Numerical results for two and three dimensional model problems are presented, together with a two level analysis explaining these results.

Decker, Naomi↗

NASA NTRS1989

Operator induced multigrid algorithms using semirefinement

A variant of multigrid, based on zebra relaxation, and a new family of restriction/prolongation operators is described. Using zebra relaxation in combination with an operator-induced prolongation leads to fast convergence, since the coarse grid can correct all error components. The resulting algorithms are not only fast, but are also robust, in the sense that the convergence rate is insensitive to the mesh aspect ratio. This is true even though line relaxation is performed in only one direction. Multigrid becomes a direct method if an operator-induced prolongation is used, together with the induced coarse grid operators. Unfortunately, this approach leads to stencils which double in size on each coarser grid. The use of an implicit three point restriction can be used to factor these large stencils, in order to retain the usual five or nine point stencils, while still achieving fast convergence. This algorithm achieves a V-cycle convergence rate of 0.03 on Poisson's equation, using 1.5 zebra sweeps per level, while the convergence rate improves to 0.003 if optimal nine point stencils are used. Numerical results for two- and three-dimensional model problems are presented, together with a two level analysis explaining these results.

Decker, Naomi Henderson↗

NASA NTRS1981

Determination of wind tunnel constraint effects by a unified pressure signature method. Part 1: Applications to winged configurations

A new, fast, non-iterative version of the "Wall Pressure Signature Method" is described and used to determine blockage and angle-of-attack wind tunnel corrections for highly-powered jet-flap models. The correction method is complemented by the application of tangential blowing at the tunnel floor to suppress flow breakdown there, using feedback from measured floor pressures. This tangential blowing technique was substantiated by subsequent flow investigations using an LV. The basic tests on an unswept, knee-blown, jet flapped wing were supplemented to include the effects of slat-removal, sweep and the addition of unflapped tips. C sub mu values were varied from 0 to 10 free-air C sub l's in excess of 18 were measured in some cases. Application of the new methods yielded corrected data which agreed with corresponding large tunnel "free air" resuls to within the limits of experimental accuracy in almost all cases. A program listing is provided, with sample cases.

Hackett, J. E.↗