Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

WE-G-207-04: Non-Local Total-Variation (NLTV) Combined with Reweighted L1-Norm for Compressed Sensing Based CT Reconstruction

Purpose: Compressed sensing (CS) has been used for CT (4DCT/CBCT) reconstruction with few projections to reduce dose of radiation. Total-variation (TV) in L1-minimization (min.) with local information is the prevalent technique in CS, while it can be prone to noise. To address the problem, this work proposes to apply a new image processing technique, called non-local TV (NLTV), to CS based CT reconstruction, and incorporate reweighted L1-norm into it for more precise reconstruction. Methods: TV minimizes intensity variations by considering two local neighboring voxels, which can be prone to noise, possibly damaging the reconstructed CT image. NLTV, contrarily, utilizes more global information by computing a weight function of current voxel relative to surrounding search area. In fact, it might be challenging to obtain an optimal solution due to difficulty in defining the weight function with appropriate parameters. Introducing reweighted L1-min., designed for approximation to ideal L0-min., can reduce the dependence on defining the weight function, therefore improving accuracy of the solution. This work implemented the NLTV combined with reweighted L1-min. by Split Bregman Iterative method. For evaluation, a noisy digital phantom and a pelvic CT images are employed to compare the quality of images reconstructed by TV, NLTV and reweighted NLTV. Results: In both cases, conventional and reweighted NLTV outperform TV min. in signal-to-noise ratio (SNR) and root-mean squared errors of the reconstructed images. Relative to conventional NLTV, NLTV with reweighted L1-norm was able to slightly improve SNR, while greatly increasing the contrast between tissues due to additional iterative reweighting process. Conclusion: NLTV min. can provide more precise compressed sensing based CT image reconstruction by incorporating the reweighted L1-norm, while maintaining greater robustness to the noise effect than TV min.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

A special generalized HSS method for discrete ill-posed problems

In this paper, we use a new special generalized Hermitian and skew-Hermitian splitting (SGHSS) method for solving ill-posed inverse problems. Based on an augmented system formulation, we apply a new splitting for the Hermitian part of the coefficient matrix and discuss on the convergence of this iterative method. Moreover, the optimal parameter minimizing the spectral radius of iteration matrix is derived. Finally, we present some numerical examples to demonstrate the performance of the presented method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Seven open problems in applied combinatorics

We present and discuss seven different open problems in applied combinatorics. Additionally, the application areas relevant to this compilation include quantum computing, algorithmic differentiation, topological data analysis, iterative methods, hypergraph cut algorithms, and power systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Performance evaluation of OpenFOAM on many-core architectures

In this article application of Open Source Field Operation and Manipulation (OpenFOAM) C++ libraries for solving engineering problems on many-core architectures is presented. Objective of this article is to present scalability of OpenFOAM on parallel platforms solving real engineering problems of fluid dynamics. Scalability test of OpenFOAM is performed using various hardware and different implementation of standard PCG and PBiCG Krylov iterative methods. Speed up of various implementations of linear solvers using GPU and MIC accelerators are presented in this paper. Numerical experiments of 3D lid-driven cavity flow for several cases with various number of cells are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

A nearly-linear computational-cost scheme for the forward dynamics of an N-body pendulum

The dynamic equations of motion of an n-body pendulum with spherical joints are derived to be a mixed system of differential and algebraic equations (DAE's). The DAE's are kept in implicit form to save arithmetic and preserve the sparsity of the system and are solved by the robust implicit integration method. At each solution point, the predicted solution is corrected to its exact solution within given tolerance using Newton's iterative method. For each iteration, a linear system of the form J delta X = E has to be solved. The computational cost for solving this linear system directly by LU factorization is O(n exp 3), and it can be reduced significantly by exploring the structure of J. It is shown that by recognizing the recursive patterns and exploiting the sparsity of the system the multiplicative and additive computational costs for solving J delta X = E are O(n) and O(n exp 2), respectively. The formulation and solution method for an n-body pendulum is presented. The computational cost is shown to be nearly linearly proportional to the number of bodies.

Chou, Jack C. K.↗
DOE OSTI2015

Final Report, DE-FG01-06ER25718 Domain Decomposition and Parallel Computing

The goal of this project is to develop and improve domain decomposition algorithms for a variety of partial differential equations such as those of linear elasticity and electro-magnetics.These iterative methods are designed for massively parallel computing systems and allow the fast solution of the very large systems of algebraic equations that arise in large scale and complicated simulations. A special emphasis is placed on problems arising from Maxwell's equation. The approximate solvers, the preconditioners, are combined with the conjugate gradient method and must always include a solver of a coarse model in order to have a performance which is independent of the number of processors used in the computer simulation. A recent development allows for an adaptive construction of this coarse component of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Absorption effects in electron-sulfur-dioxide collisions

A joint experimental-theoretical study on electron-SO{sub 2} collisions in the low and intermediate energy range is reported. More specifically, experimental elastic differential, integral, and momentum transfer cross sections in absolute scale are measured in the 100-1000 eV energy range using the relative-flow technique. Calculated elastic differential, integral, and momentum transfer cross sections as well as grand-total and total absorption cross sections are also presented in the 1-1000 eV energy range. A complex optical potential is used to represent the electron-molecule interaction dynamics, whereas the Schwinger variational iterative method combined with the distorted-wave approximation is used to solve the scattering equations. Comparison of the present results is made with the theoretical and experimental results available in the literature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2011

Thermophysical properties of nitrogen plasmas under thermal equilibrium and non-equilibrium conditions

Calculated thermophysical properties of nitrogen plasmas in and out of thermal equilibrium are presented. The cut-off of the partition functions due to the lowering of the ionization potential has been taken into account, together with the contributions from different core excited electronic states. The species composition and thermodynamic properties are determined numerically using the Newton-Raphson iterative method, taking into account the corrections due to Coulomb interactions. The transport properties including diffusion coefficient, viscosity, thermal conductivity, and electrical conductivity are calculated using the most recent collision interaction potentials by adopting Devoto's electron and heavy particle decoupling approach, expanded to the third-order approximation (second-order for viscosity) in the framework of Chapman-Enskog method. Results are presented in the pressure range of 0.1 atm-10 atm and in electron temperature range from 300 to 40 000 K, with the ratio of electron temperature to heavy-particle temperature varied from 1 to 20. Results are compared with those from previous works, and the influences of different definitions of the Debye length are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

The multigrid preconditioned conjugate gradient method

A multigrid preconditioned conjugate gradient method (MGCG method), which uses the multigrid method as a preconditioner of the PCG method, is proposed. The multigrid method has inherent high parallelism and improves convergence of long wavelength components, which is important in iterative methods. By using this method as a preconditioner of the PCG method, an efficient method with high parallelism and fast convergence is obtained. First, it is considered a necessary condition of the multigrid preconditioner in order to satisfy requirements of a preconditioner of the PCG method. Next numerical experiments show a behavior of the MGCG method and that the MGCG method is superior to both the ICCG method and the multigrid method in point of fast convergence and high parallelism. This fast convergence is understood in terms of the eigenvalue analysis of the preconditioned matrix. From this observation of the multigrid preconditioner, it is realized that the MGCG method converges in very few iterations and the multigrid preconditioner is a desirable preconditioner of the conjugate gradient method.

Tatebe, Osamu↗
NASA NTRS2016

Assessment of Preconditioner for a USM3D Hierarchical Adaptive Nonlinear Method (HANIM) (Invited)

Enhancements to the previously reported mixed-element USM3D Hierarchical Adaptive Nonlinear Iteration Method (HANIM) framework have been made to further improve robustness, efficiency, and accuracy of computational fluid dynamic simulations. The key enhancements include a multi-color line-implicit preconditioner, a discretely consistent symmetry boundary condition, and a line-mapping method for the turbulence source term discretization. The USM3D iterative convergence for the turbulent flows is assessed on four configurations. The configurations include a two-dimensional (2D) bump-in-channel, the 2D NACA 0012 airfoil, a three-dimensional (3D) bump-in-channel, and a 3D hemisphere cylinder. The Reynolds Averaged Navier Stokes (RANS) solutions have been obtained using a Spalart-Allmaras turbulence model and families of uniformly refined nested grids. Two types of HANIM solutions using line- and point-implicit preconditioners have been computed. Additional solutions using the point-implicit preconditioner alone (PA) method that broadly represents the baseline solver technology have also been computed. The line-implicit HANIM shows superior iterative convergence in most cases with progressively increasing benefits on finer grids.

Pandya, Mohagna J.↗
DOE OSTI2013

Polymer/metal hybrid multilayers modified Schottky devices

Insulating, polymethylmethacrylate (PMMA), and semiconducting, poly(3-hexylthiophene) (P3HT), nanometer thick polymers/Au nanoparticles based hybrid multilayers (HyMLs) were fabricated on p-Si single-crystal substrate. An iterative method, which involves, respectively, spin-coating (PMMA and P3HT deposition) and sputtering (Au nanoparticles deposition) techniques to prepare Au/HyMLs/p-Si Schottky device, was used. The barrier height and the ideality factor of the Au/HyMLs/p-Si Schottky devices were investigated by current-voltage measurements in the thickness range of 1–5 bilayers. It was observed that the barrier height of such hybrid layered systems can be tuned as a function of bilayers number and its evolution was quantified and analyzed.

36 MATERIALS SCIENCE↗
DOE OSTI2015

TU-EF-204-04: Feasibility Study of a Novel Stationary Spectral CT Scanner

Purpose: To evaluate a cadmium zinc telluride (CZT) detector module and to develop a proof-of-concept bench-top system to examine the concept and feasibility of a stationary spectral scanner that would consist of a CZT detector ring and a distributed source for rotationless imaging. Methods: The bench-top system employs a double-rotation set-up with two concentric but independent rotating stages, one for turning the object at sparse angles to mimic the imaging from a sequentially activated source distributed on a ring, while the other for translating the detector module to simulate a projection arc along a virtual detector ring. A source of 30-degree fan is used in the set-up, and the ring is designed to be 85cm in diameter in order to provide a 22cm field-of-view that is large enough for an average human head. The CZT detector module has 16×16 0.5mm pixels and eight global energy thresholds, and is attached to an experimental data acquisition board from the same vendor (Nova R&D). Its pixel gain is first calibrated using radioisotope Tc-99m to align the photo-peaks. Count rate performance is assessed using the X-ray source. Air and offset corrections are done prior to each CT scan. A 10cm acrylic phantom with 5 different inserts is imaged at 100kVp, 2mA, with an effective exposure window of 150ms per projection. Reconstruction is done using a total variation-type iterative method. Results: The detector module exhibited large gain variations among pixels. Average count rate was over 1Mcps per pixel. Parallax effect was observed along the arcs but was corrected. Reconstructed images from two energy bins demonstrated good contrast, spatial resolution and image uniformity, with little artifacts. Conclusion: The CZT detector tested is suitable for imaging relatively small objects. A stationary spectral CT appears feasible should we have enough buttable modules and a distributed source. This work was supported by Ministry of Science and Technology of China through grants 11YQ040082 and 2012BAI13B04.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

The Transfer of Resonance Line Polarization with Partial Frequency Redistribution in the General Hanle–Zeeman Regime

The spectral line polarization encodes a wealth of information about the thermal and magnetic properties of the solar atmosphere. Modeling the Stokes profiles of strong resonance lines is, however, a complex problem both from a theoretical and computational point of view, especially when partial frequency redistribution (PRD) effects need to be taken into account. In this work, we consider a two-level atom in the presence of magnetic fields of arbitrary intensity (Hanle–Zeeman regime) and orientation, both deterministic and micro-structured. Working within the framework of a rigorous PRD theoretical approach, we have developed a numerical code that solves the full non-LTE radiative transfer problem for polarized radiation, in one-dimensional models of the solar atmosphere, accounting for the combined action of the Hanle and Zeeman effects, as well as for PRD phenomena. After briefly discussing the relevant equations, we describe the iterative method of solution of the problem and the numerical tools that we have developed and implemented. We finally present some illustrative applications to two resonance lines that form at different heights in the solar atmosphere, and provide a detailed physical interpretation of the calculated Stokes profiles. We find that magneto-optical effects have a strong impact on the linear polarization signals that PRD effects produce in the wings of strong resonance lines. We also show that the weak-field approximation has to be used with caution when PRD effects are considered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Thermoelectric power of supported graphene - An iterative approach

Thermoelectric power, S, of graphene supported on SiO{sub 2} substrate is studied for 10K < T < 400K within the framework of Boltzmann transport formalism by an iterative method. Numerical calculations of diffusion thermopower, S{sub d}, as a function of temperature, are presented assuming the electrons to be scattered by impurities, vacancies, surface roughness, acoustic phonons, inelastic optical phonons and surface polar optical phonons. For the range of temperatures considered, S{sub d} is found to be dominated by impurities for T < 40K and by acoustic phonon and vacancy scatterings for T > 40K. The optical phonons are found to influence S{sub d} for T > 300K. Our calculations, assuming the drag component to be negligible, obtain good agreement with the recent experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

2022 AI Testbed Expeditions Report

By exploiting the coherent properties of a light source, coherent diffraction imaging (CDI) is able to obtain the sample image at a nanoscale resolution using the measured diffraction pattern. Bragg Coherent Diffraction Imaging (BCDI) has become valuable for recovering the displacement and strain field of crystals, providing a valuable tool in material science and solid-state physics. X-ray ptychography is another emerging CDI technique that can produce a high-resolution image of the extended sample and has become popular in many research areas (e.g., materials science, biology, electronics, and optics characterization). CDI including BCDI and ptychography has become an established technique in Synchrotron Facilities including the Advanced Photon Source (APS) and will greatly benefit from the 100x coherent flux increase of the upcoming APS Upgrade (APSU). The current image formation process in CDI employs iterative phase retrieval algorithms, which is a time-consuming and computationally expensive process. Especially after APSU, the traditional iterative methods will not be able to match the experimental data acquisition speed. We employ deep learning (DL) approach to replace the iterative approaches, therefore allowing hundreds of times faster recovery of the object. We developed AutoPhaseNN, a DL-based approach which learns to solve the inverse problem without labeled data. Taking 3D BCDI as a representative technique, AutoPhaseNN has been demonstrated to be one hundred times faster than traditional iterative phase retrieval methods while providing comparable image quality. The current network is trained with 64 x 64 x 64 data size, to achieve higher resolution imaging, we will need to scale the network to input and train/infer 3D arrays of size 256 x 256 x 256 (today) and of size 2560x2560x2560 (APSU). However, the scalability of the network is restricted due to the memory-intensive training process. To perform the training for a 256 x 256 x 256 data size, the required memory exceeds the capacity of the current machine. In this project, we explore using Sambanova system to train the network for the direct data inversion for CDI.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Attenuation correction in emission tomography using the emission data—A review

The problem of attenuation correction (AC) for quantitative positron emission tomography (PET) had been considered solved to a large extent after the commercial availability of devices combining PET with computed tomography (CT) in 2001; single photon emission computed tomography (SPECT) has seen a similar development. However, stimulated in particular by technical advances toward clinical systems combining PET and magnetic resonance imaging (MRI), research interest in alternative approaches for PET AC has grown substantially in the last years. In this comprehensive literature review, the authors first present theoretical results with relevance to simultaneous reconstruction of attenuation and activity. The authors then look back at the early history of this research area especially in PET; since this history is closely interwoven with that of similar approaches in SPECT, these will also be covered. We then review algorithmic advances in PET, including analytic and iterative algorithms. The analytic approaches are either based on the Helgason–Ludwig data consistency conditions of the Radon transform, or generalizations of John’s partial differential equation; with respect to iterative methods, we discuss maximum likelihood reconstruction of attenuation and activity (MLAA), the maximum likelihood attenuation correction factors (MLACF) algorithm, and their offspring. The description of methods is followed by a structured account of applications for simultaneous reconstruction techniques: this discussion covers organ-specific applications, applications specific to PET/MRI, applications using supplemental transmission information, and motion-aware applications. After briefly summarizing SPECT applications, we consider recent developments using emission data other than unscattered photons. In summary, developments using time-of-flight (TOF) PET emission data for AC have shown promising advances and open a wide range of applications. These techniques may both remedy deficiencies of purely MRI-based AC approaches in PET/MRI and improve standalone PET imaging.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

RADIATIVE TRANSFER MODELING OF THE ENIGMATIC SCATTERING POLARIZATION IN THE SOLAR Na i D{sub 1} LINE

The modeling of the peculiar scattering polarization signals observed in some diagnostically important solar resonance lines requires the consideration of the detailed spectral structure of the incident radiation field as well as the possibility of ground level polarization, along with the atom's hyperfine structure and quantum interference between hyperfine F-levels pertaining either to the same fine structure J-level, or to different J-levels of the same term. Here we present a theoretical and numerical approach suitable for solving this complex non-LTE radiative transfer problem. This approach is based on the density-matrix metalevel theory (where each level is viewed as a continuous distribution of sublevels) and on accurate formal solvers of the transfer equations and efficient iterative methods. We show an application to the D-lines of Na i, with emphasis on the enigmatic D{sub 1} line, pointing out the observable signatures of the various physical mechanisms considered. We demonstrate that the linear polarization observed in the core of the D{sub 1} line may be explained by the effect that one gets when the detailed spectral structure of the anisotropic radiation responsible for the optical pumping is taken into account. This physical ingredient is capable of introducing significant scattering polarization in the core of the Na i D{sub 1} line without the need for ground-level polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

A numerical method for the solution of three dimensional, incompressible, viscous flows over slender bodies

A marching iterative method for the solution of the three dimensional, incompressibhle, steady and parabolized Navier-Stokes equations is described. The equations are written in primitive variables and discretized in general axisymmetric orthogonal coordinate systems. The coupled set of finite-difference equations are solved without any splitting or factorization errors. Moreover, the continuity equation and the two crossflow momentum equations are exactly satisfied at every step of the iterative process. The solution scheme is equivalent to the solution of one Poisson equation by the Successive Plane Over Relaxation method and has good convergence properties. Other existing solution methods resemble a Jacobi-type iterative scheme and therefore are less efficient. Numerical experiments include the laminar, incompressible flow over prolate spheroids at incidence.

Rosenfeld, Moshe↗