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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Construction of fractal nanostructures based on Kepler-Shubnikov nets

A system of information codes for deterministic fractal lattices and sets of multifractal curves is proposed. An iterative modular design was used to obtain a series of deterministic fractal lattices with generators in the form of fragments of 2D structures and a series of multifractal curves (based on some Kepler-Shubnikov nets) having Cantor set properties. The main characteristics of fractal structures and their lacunar spectra are determined. A hierarchical principle is formulated for modules of regular fractal structures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2013

Effective regularized algebraic reconstruction technique for computed tomography

A new fast version of the reconstruction algorithm for computed tomography based on the simultaneous algebraic reconstruction technique (SART) is proposed. The algorithm iteration is asymptotically accelerated using the fast Hough transform from O(n{sup 3}) to O(n{sup 2}logn). Similarly to the algebraic reconstruction technique (RegART), which was proposed by us previously, the regularization operator is applied after each iteration. A bilateral filter plays the role of this operator. The algorithm behavior is investigated using the model experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Dynamics of an electron in a relativistically intense laser field including radiaion reaction

The dynamics of an electron in a relativistically intense laser pulse field is described with the radiation reaction being taken into account. The study is based on solving the Newton equation with the Lorentz and the radiation reaction forces. Validation is provided for an iteration technique which makes it possible to remove the discrepancies found in the theoretical models of radiation reaction. It is demonstrated that an electron having a high initial velocity and colliding head-on with a laser pulse sheds a considerable part of its kinetic energy due to the radiation reaction. A broadening of the electromagnetic pulse emitted by the electron occurs as a result of the same effect. The findings obtained can be used to experimentally verify the effect of radiation reaction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Helium trimer calculations with a public quantum three-body code

We present an illustration of using a quantumthree-body code being prepared for public release. The code is based on iterative solving of the three-dimensional Faddeev equations. The code is easy to use and allows users to perform highly-accurate calculations of quantum three-body systems. The previously known results for He{sub 3} ground state are well reproduced by the code.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Elastic scattering of hadrons outside the diffraction cone

The model-independent solution of the s-channel unitarity condition for the imaginary part of the hadronic elastic scattering amplitude outside the diffraction peak allows to make conclusions about the shape of the differential cross sections in this region. The eigensolution of the homogeneous linear integral equation and the iterative solutions of the inhomogeneous equation with different assumptions about the behavior of the overlap function and the real part of the amplitude are presented. The properties of the overlap function and of the ratio of the real-to-imaginary part of the amplitude are discussed. In particular, it is explicitly shown that the ratio changes its sign at a definite value of the transferred momentum. Some comments concerning the present day experimental results about the behavior of the differential cross section of elastic scattering outside the diffraction cone are given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Investigation of heat transfer in liquid-metal flows under fusion-reactor conditions

The effect discovered in studying a downward liquid-metal flow in vertical pipe and in a channel of rectangular cross section in, respectively, a transverse and a coplanar magnetic field is analyzed. In test blanket modules (TBM), which are prototypes of a blanket for a demonstration fusion reactor (DEMO) and which are intended for experimental investigations at the International Thermonuclear Experimental Reactor (ITER), liquid metals are assumed to fulfil simultaneously the functions of (i) a tritium breeder, (ii) a coolant, and (iii) neutron moderator and multiplier. This approach to testing experimentally design solutions is motivated by plans to employ, in the majority of the currently developed DEMO blanket projects, liquid metals pumped through pipes and/or rectangular channels in a transvers magnetic field. At the present time, experiments that would directly simulate liquid-metal flows under conditions of ITER TBM and/or DEMO blanket operation (irradiation with thermonuclear neutrons, a cyclic temperature regime, and a magnetic-field strength of about 4 to 10 T) are not implementable for want of equipment that could reproduce simultaneously the aforementioned effects exerted by thermonuclear plasmas. This is the reason why use is made of an iterative approach to experimentally estimating the performance of design solutions for liquid-metal channels via simulating one or simultaneously two of the aforementioned factors. Therefore, the investigations reported in the present article are of considerable topical interest. The respective experiments were performed on the basis of the mercury magneto hydrodynamic (MHD) loop that is included in the structure of the MPEI—JIHT MHD experimental facility. Temperature fields were measured under conditions of two- and one-sided heating, and data on averaged-temperature fields, distributions of the wall temperature, and statistical fluctuation features were obtained. A substantial effect of counter thermo gravitational convection (TGC) on averaged and fluctuating quantities were found. The development of TGC in the presence of a magnetic field leads to the appearance of low-frequency fluctuations whose anomalously high intensity exceeds severalfold the level of turbulence fluctuations. This effect manifest itself over a broad region of regime parameters. It was confirmed that low-energy fluctuations penetrate readily through the wall; therefore, it is necessary to study this effect further—in particular, from the point of view of the fatigue strength of the walls of liquid-metal channels.

42 ENGINEERING↗
DOE OSTI2016

Algorithm for solving the linear Cauchy problem for large systems of ordinary differential equations with the use of parallel computations

An algorithm for solving the linear Cauchy problem for large systems of ordinary differential equations is presented. The algorithm for systems of first-order differential equations is implemented in the EDELWEISS code with the possibility of parallel computations on supercomputers employing the MPI (Message Passing Interface) standard for the data exchange between parallel processes. The solution is represented by a series of orthogonal polynomials on the interval [0, 1]. The algorithm is characterized by simplicity and the possibility to solve nonlinear problems with a correction of the operator in accordance with the solution obtained in the previous iterative process.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Conditioning Analysis of Incomplete Cholesky Factorizations with Orthogonal Dropping

The analysis of preconditioners based on incomplete Cholesky factorization in which the neglected (dropped) components are orthogonal to the approximations being kept is presented. General estimate for the condition number of the preconditioned system is given which only depends on the accuracy of individual approximations. The estimate is further improved if, for instance, only the newly computed rows of the factor are modified during each approximation step. In this latter case it is further shown to be sharp. The analysis is illustrated with some existing factorizations in the context of discretized elliptic partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Generic Discretization Library

The GenDiL library is a collection of C++ software abstractions designed to discretize and solve partial differential equations (PDEs) for high-performance computing (HPC) applications. Its primary focus is on modern C++ generic programming, which helps ensure portability across various hardware architectures. The central idea behind the library is to provide building blocks for numerical algorithms-such as discretization methods and iteration patterns-so that domain experts can focus on the math, rather than the low-level details of hardware or implementation. By defining abstractions for data types, iteration over computational grids, and scheduling of operations, the library isolates the high-level PDE algorithms from the platform-specific optimizations needed to achieve efficient performance.

Dudouit, Yohann [Lawrence Livermore National Labor↗
DOE OSTI2024

Dynamic Graph Sequence Data from Simulated Neutron Reflectometry Measurements

This dataset comprises dynamic graph sequences derived from simulated in-situ neutron reflectometry measurements, capturing the gradual evolution of a layer structure over time. Each graph sequence represents a synthetic sample, with node features detailing the scattering vector and corresponding reflectivity measurements, while adjacency matrices have corresponding reference material parameters attached as metadata. The dataset spans multiple sets, each with a different number of sequences, offering a comprehensive basis for training models that handle dynamic input sequences with embedded physics. This dataset is particularly suited for tackling inverse problems with hidden physical states that evolve over time, challenges that are typically difficult to address using conventional iterative fitting methods.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Zeolite membrane cascade for tritium extraction and recovery systems

Membrane separation by zeolite membranes has been proposed as a pre-concentration stage for the tritium extraction from the purge helium of the breeding blanket combined with a final recovery by the catalytic membrane reactor PERMCAT. This fully continuous operation improves the tritium management in fusion machines, minimizing the tritium inventory. For the first time, the permeation measurements for H{sub 2} - He mixtures through a MFI-alumina hollow fibre membrane has been measured for different compositions (0.1 - 20% H{sub 2}) and temperatures. Such a highly permeable membrane, although it shows a limited selectivity, appears attractive for tritium recovery in the blanket. This will imply its operation in a membrane cascade, for which simulation work is ongoing. Mathematically the process is modeled using mass balance equations that can be transformed into the matrix form and solved iteratively assuming a permeate concentration on the first step of iteration, until the separation requirements are fulfilled.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Dynamic tomography with a priori information

We present ''dynamic tomography'' algorithms that allow for the high-resolution, time-resolved imaging of dynamic (i.e., continuously time evolving) complex systems at existing x-ray micro-CT facilities. The behavior of complex systems is constrained by the underlying physics. By exploiting a priori knowledge of the geometry of the physical process being studied to allow the use of sophisticated iterative reconstruction techniques that incorporate constraints, we improve on current frame rates by at least an order of magnitude. This allows time-resolved imaging of previously intractable processes, such as two-phase fluid flow. We present reconstructions from experimental data collected at the Australian National University x-ray micro-CT facility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Information content of and the ability to reconstruct dichroic X-ray tomography and laminography

Dichroic tomography is a 3D imaging technique in which the polarization of the incident beam is used to induce contrast due to the magnetization or orientation of a sample. The aim is to reconstruct not only the optical density but the dichroism of the sample. The theory of dichroic tomographic and laminographic imaging in the parallel-beam case is discussed as well as the problem of reconstruction of the sample’s optical properties. The set of projections resulting from a single tomographic/laminographic measurement is not sufficient to reconstruct the magnetic moment for magnetic circular dichroism unless additional constraints are applied or data are taken at two or more tilt angles. For linear dichroism, three polarizations at a common tilt angle are insufficient for unconstrained reconstruction. However, if one of the measurements is done at a different tilt angle than the other, or the measurements are done at a common polarization but at three distinct tilt angles, then there is enough information to reconstruct without constraints. Possible means of applying constraints are discussed. Furthermore, it is shown that for linear dichroism, the basic assumption that the absorption through a ray path is the integral of the absorption coefficient, defined on the volume of the sample, along the ray path, is not correct when dichroism or birefringence is strong. This assumption is fundamental to tomographic methods. An iterative algorithm for reconstruction of linear dichroism is demonstrated on simulated data.

Marcus, Matthew A. (ORCID:0000000325277586)↗
DOE OSTI2015

Argonne Bubble Experiment Thermal Model Development

This report will describe the Computational Fluid Dynamics (CFD) model that was developed to calculate the temperatures and gas volume fractions in the solution vessel during the irradiation. It is based on the model used to calculate temperatures and volume fractions in an annular vessel containing an aqueous solution of uranium . The experiment was repeated at several electron beam power levels, but the CFD analysis was performed only for the 12 kW irradiation, because this experiment came the closest to reaching a steady-state condition. The aim of the study is to compare results of the calculation with experimental measurements to determine the validity of the CFD model.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2015

Verification and Validation of Digitally Upgraded Control Rooms

As nuclear power plants undertake main control room modernization, a challenge is the lack of a clearly defined human factors process to follow. Verification and validation (V&V) as applied in the nuclear power community has tended to involve efforts such as integrated system validation, which comes at the tail end of the design stage. To fill in guidance gaps and create a step-by-step process for control room modernization, we have developed the Guideline for Operational Nuclear Usability and Knowledge Elicitation (GONUKE). This approach builds on best practices in the software industry, which prescribe an iterative user-centered approach featuring multiple cycles of design and evaluation. Nuclear regulatory guidance for control room design emphasizes summative evaluation—which occurs after the design is complete. In the GONUKE approach, evaluation is also performed at the formative stage of design—early in the design cycle using mockups and prototypes for evaluation. The evaluation may involve expert review (e.g., software heuristic evaluation at the formative stage and design verification against human factors standards like NUREG-0700 at the summative stage). The evaluation may also involve user testing (e.g., usability testing at the formative stage and integrated system validation at the summative stage). An additional, often overlooked component of evaluation is knowledge elicitation, which captures operator insights into the system. In this report we outline these evaluation types across design phases that support the overall modernization process. The objective is to provide industry-suitable guidance for steps to be taken in support of the design and evaluation of a new human-machine interface (HMI) in the control room. We suggest the value of early-stage V&V and highlight how this early-stage V&V can help improve the design process for control room modernization. We argue that there is a need to overcome two shortcomings of V&V in current practice—the propensity for late-stage V&V and the use of increasingly complex psychological assessment measures for V&V.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

MOC Efficiency Improvements Using a Jacobi Inscatter Approximation

In recent weeks, attention has been given to resolving the convergence issues encountered with TCP 0 by trying a Jacobi (J) inscatter approach when group sweeping, where the inscatter source is constructed using the previous iteration flux. This is in contrast to a Gauss-Seidel (GS) approach, which has been the default to-date, where the scattering source uses the most up-to-date flux values. The former is consistent with CASMO, which has no issues with TCP 0 convergence. Testing this out on a variety of problems has demonstrated that the Jacobi approach does indeed provide substantially more stability, though can take more outer iterations to converge. While this is not surprising, there are improvements that can be made to the MOC sweeper to capitalize on the Jacobi approximation and provide substantial speedup. For example, the loop over groups, which has traditionally been the outermost loop in MPACT, can be moved to the interior, avoiding duplicate modular ray trace and coarse ray trace setup (mapping coarse mesh surface indexes), which needs to be performed repeatedly when group is outermost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Progress Update on Iterative Reconstruction of Neutron Tomographic Images

This report satisfies the fiscal year 2016 technical deliverable to report on progress in development of fast iterative reconstruction algorithms for project OR16-3DTomography-PD2Jb, "3D Tomography and Image Processing Using Fast Neutrons." This project has two overall goals. The first of these goals is to extend associated-particle fast neutron transmission and, particularly, induced-reaction tomographic imaging algorithms to three dimensions. The second of these goals is to automatically segment the resultant tomographic images into constituent parts, and then extract information about the parts, such as the class of shape and potentially shape parameters. This report addresses of the component of the project concerned with three-dimensional (3D) image reconstruction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗