Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Neumann Series in MGS-GMRES and Inner-Outer Iterations: Preprint

A low-synchronization MGS-GMRES Krylov solver employing a truncated Neumann series for the inverse compact WY MGS correction matrix T is presented. A corollary to the backward stability result of Paige et al. [1] establishes that T = I - Lk is sufficient for convergence of GMRES when kLkp F = O("p)_p F (B), where the strictly lower triangular matrix L is defined by the inner products of Krylov vectors V T 1:k-2 vk-1. The preconditioner is the classical Ruge-Stuben AMG algorithm with compatible relaxation and inner-outer Gauss-Seidel smoother. This smoother may also be expressed as a truncated Neumann series. Drop tolerances are applied to the lower triangular matrices arising in the smoother in order to reduce the number of non-zeros and accelerate the time to solution. The number of small matrix elements are found to increase from fine to coarse levels and thus the effciency gains are greater for large problems with many levels in the V -cycle. The solver is applied to the pressure continuity equation for the incompressible Navier-Stokes equations. Unlike the inner-outer iteration, the solver convergence rate with the standard Gauss-Seidel smoother deteriorates with dropping. The solver compute time is reduced by up to 50% without a change in the convergence rate.

Gauss-Seidel smoother↗
DOE OSTI2013

Refinement of overlapping local/global iteration method based on Monte Carlo/p-CMFD calculations

In this paper, the overlapping local/global (OLG) iteration method based on Monte Carlo/p-CMFD calculations is refined in two aspects. One is the consistent use of estimators to generate homogenized scattering cross sections. Another is that the incident or exiting angular interval is divided into multi-angular bins to modulate albedo boundary conditions for local problems. Numerical tests show that, compared to the one angle bin case in a previous study, the four angle bin case shows significantly improved results. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Investigations into alternative radiation transport codes for ITER neutronics analysis

As radiation transport models for ITER become increasingly complex, the use of traditional constructive solid geometry (CSG) modelling approaches such as MCWNP presents significant disadvantages during model creation and integration. Studies were performed to assess the capabilities of alternative mesh-based neutronics workflows offered by the MCNP, DAG-MCNP and Serpent-2 codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Iterative sampling of expensive simulations for faster deep surrogate training

Deep neural network (DNN) surrogates of expensive physics simulations are enabling a rapid change in the way that common experimental design and analysis tasks are approached. Surrogate models allow simulations to be performed in parallel and separately from downstream tasks, thereby enabling analyses that would be impossible with the simulation in-the-loop; surrogates based on DNNs can effectively emulate diverse non-scalar data of the types collected in fusion and laboratory-astrophysics experiments. The challenge is in training the surrogate model, for which large ensembles of physics simulations must be run, preferably without wasting computational effort on uninteresting simulations. Here, in this paper, we present an iterative sampling scheme that can preferentially propose simulations in interesting regions of parameter space without neglecting unexplored regions, allowing high-quality and wide-ranging surrogate models to be trained using 2–3 times fewer simulations compare to space-filling designs. Our approach uses an explicit importance function defined on the simulation output space, balanced against a measure of simulation density which serves as a proxy for surrogate accuracy. It is easy to implement and can be tuned to find interesting simulations early in the study, allowing surrogates to be trained quickly and refined as new simulations become available; this represents an important step towards the routine generation of deep surrogate models quickly enough to be truly relevant to experimental work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

UEDGE modeling of plasma detachment of CFETR with ITER ‐like divertor geometry by external impurity seeding

Abstract Efficient handling of high heat flux on the plasma‐facing components, particularly the divertor targets, poses a significant challenge for the Chinese Fusion Engineering Testing Reactor (CFETR) with fusion power of Gigawatt. This work investigates the divertor plasma detachment of CFETR with a standard ITER‐like divertor geometry by neon (Ne) or argon (Ar) impurity seeding using UEDGE code. The cross‐field drifts terms are switched off, and fluid neutral models and a “fixed‐fraction” impurity model are applied to enable efficient simulations for the study of CFETR detachment. In order to reduce the heat load on the divertor targets below the acceptable level (<10 MW/m 2 ), the impurity fraction ( f ), pumping speed ( S ), and upstream density are varied to identify the suitable operations window during Ne seeding. The effects of Ne and Ar impurities on the plasma detachment are compared. It is found that with the power across the core‐edge interface P SOL = 200 MW and separatrix density of 2.8 10 19 , Ne impurity fraction ≥1.7%, and Ar impurity fraction ≥0.24% can achieve the partial detachment. Achieving similar total radiation power (˜148 MW), the Ne fraction is 2.3% and the Ar fraction is 0.24%. Moreover, the simulation results indicate that Ar exhibits better power radiation efficiency and core compatibility compared with Ne.

Zhang, M. Z.↗
DOE OSTI2014

New algorithms for iterative matrix‐free eigensolvers in quantum chemistry

New algorithms for iterative diagonalization procedures that solve for a small set of eigen‐states of a large matrix are described. The performance of the algorithms is illustrated by calculations of low and high‐lying ionized and electronically excited states using equation‐of‐motion coupled‐cluster methods with single and double substitutions (EOM‐IP‐CCSD and EOM‐EE‐CCSD). We present two algorithms suitable for calculating excited states that are close to a specified energy shift (interior eigenvalues). One solver is based on the Davidson algorithm, a diagonalization procedure commonly used in quantum‐chemical calculations. The second is a recently developed solver, called the “Generalized Preconditioned Locally Harmonic Residual (GPLHR) method.” We also present a modification of the Davidson procedure that allows one to solve for a specific transition. The details of the algorithms, their computational scaling, and memory requirements are described. The new algorithms are implemented within the EOM‐CC suite of methods in the Q‐Chem electronic structure program. © 2014 Wiley Periodicals, Inc.

Zuev, Dmitry↗
DOE OSTI2014

Iterative structure of finite loop integrals

In this paper we develop further and refine the method of differential equations for computing Feynman integrals. In particular, we show that an additional iterative structure emerges for finite loop integrals. As a concrete non-trivial example we study planar master integrals of light-by-light scattering to three loops, and derive analytic results for all values of the Mandelstam variables s and t and the mass m. We start with a recent proposal for defining a basis of loop integrals having uniform transcendental weight properties and use this approach to compute all planar two-loop master integrals in dimensional regularization. We then show how this approach can be further simplified when computing finite loop integrals. This allows us to discuss precisely the subset of integrals that are relevant to the problem. We find that this leads to a block triangular structure of the differential equations, where the blocks correspond to integrals of different weight. We explain how this block triangular form is found in an algorithmic way. Another advantage of working in four dimensions is that integrals of different loop orders are interconnected and can be seamlessly discussed within the same formalism. We use this method to compute all finite master integrals needed up to three loops. Finally, we remark that all integrals have simple Mandelstam representations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Iterative Method to Obtain Raman Spectra of Organic Dyes by Their Chaotic Generation in Multiple-Scattering Media

An iterative method is developed to obtain vibrational spectra of laser organic dyes from their secondary radiation spectra in multiple-scattering media, which manifests as a result of coupled process of chaotic generation (CG) and stimulated Raman scattering (SRS). The developed method increases significantly the capability to detect weak Raman lines by using batch processing of a set of experimental SRS-CG spectra, which makes it possible to optimize the deconvolution of these spectra and to average vibrational spectra obtained from the set of experimental SRS-CG spectra. The method is tested on spectra of the SRS-CG laser dye pyrromethene 597 in a vesicular film. The obtained results agree well with Raman spectra calculated by the Hartree–Fock method.

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

Effect of High Magnetic Pulse Iteration on Electrical Property of CNT-Polypyrrole Composite: Alignment Mechanism

In this work, a pulsed magnetic field of 10 T with a duration of 1.5 ms is produced by a resistor-inductor-capacitor (RLC) discharge system. A carbon nanotube-polypyrrole (CNT-PPy) composite (9 wt%) is prepared. The influence of the pulsed magnetic field iteration on the alignment of CNTs embedded in polymer matrix is observed (up to 60 pulses). The microstructure of the magnetically aligned composites is probed by scanning probe microscopy technique. The behavior of a CNT in a magnetic field is theoretically studied, and the equivalent magnetic susceptibility for an anisotropic CNT is calculated. The theory of CNT alignment in a magnetic field is discussed. We found that the perpendicular diamagnetism of a CNT results in its alignment in a magnetic field. We showed that the magnetic torque, which is applied to a CNT, depends on the perpendicular and parallel susceptibilities of the CNT as well as on its initial orientation with respect to the magnetic direction. After exposure to magnetic pulses, the electrical property of the composites is enhanced. The results show that there is a positive correlation between the electrical conductance of CNTs and the number of applied pulses (for 30 pulses). Fourier transform infrared spectroscopy technique is used to control the polymer structure during the experiment. Our theoretical finding reveals the limited capacity of CNT’s alignment in magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Performance of preconditioned iterative solvers in MFiX–Trilinos for fluidized beds

MFiX, a general-purpose Fortran-based suite, simulates the complex flow in fluidized bed applications via BiCGStab and GMRES methods along with plane relaxation preconditioners. Trilinos, an object-oriented framework, contains various first- and second-generation Krylov subspace solvers and preconditioners. We developed a framework to integrate MFiX with Trilinos as MFiX does not possess advanced linear methods. The framework allows MFiX to access advanced linear solvers and preconditioners in Trilinos. The integrated solver is called MFiX–Trilinos, here after. In the present work, we study the performance of variants of GMRES and CGS methods in MFiX–Trilinos and BiCGStab and GMRES solvers in MFiX for a 3D gas–solid fluidized bed problem. Two right preconditioners employed along with various solvers in MFiX–Trilinos are Jacobi and smoothed aggregation. The flow from MFiX–Trilinos is validated against the same from MFiX for BiCGStab and GMRES methods. And, the effect of the preconditioning on the iterative solvers in MFiX–Trilinos is also analyzed. In addition, the effect of left and right smoothed aggregation preconditioning on the solvers is studied. Finally, the performance of the first- and second-generation solver stacks in MFiX–Trilinos is studied as well for two different problem sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Sparse, iterative simulation methods for one-dimensional laminar flames

Sparse, iterative simulation methods for one-dimensional laminar flames are proposed. The resulting steady and unsteady flame solvers exploit approximate Jacobians to greatly reduce the computational cost associated with matrix operations. The constant, non-unity Lewis number assumption is introduced to further reduce the computational cost. The solvers are also parallelized to reduce time-to-solution on distributed memory computer systems. Computed laminar flame speeds and species profiles for a range of chemical mechanisms (from 10 to 2878 species) are compared against a well-validated commercial code and are found to be consistent within solver tolerances. The computation times of both the unsteady and steady solutions increase only linearly with the number of species, which is a significant improvement over the quadratic or cubic scaling of existing steady-state flame solvers. For the largest mechanism tested, the steady-state flame solver is two orders of magnitude faster than commonly-used codes. The use of an approximate Jacobian is shown to reduce the rate of convergence for the steady-state solver, but does not significantly affect the domain of convergence. The steady-state solver with approximate Jacobian is thus well suited for computationally efficient laminar flame speed sweeps with large kinetic mechanisms.

42 ENGINEERING↗
DOE OSTI2021

Irradiation hardening and ductility loss of Eurofer97 steel variants after neutron irradiation to ITER-TBM relevant conditions

Ten Eurofer97 steel variants, produced by non-standard fabrication-processing routes and modified alloying chemistries, were studied by neutron irradiations in the high flux isotope reactor. The irradiations were performed to ITER-TBM relevant conditions of ~255–350 °C, 2.94–3.24 dpa. We quantified the irradiation-induced degradation of the steels using mechanical property tests. All the steels suffered from irradiation hardening, where a significant increase in Vickers microhardness and yield stress (σ YS ) occurred, accompanied with severe loss of tensile elongation. The extent of hardening was material dependent. For T irr = 300±30 °C, most steels showed σ YS increase in the range of ~30% to as high as ~66%, except for a low temperature tempered steel with σYS increase below 15%. Despite large losses in elongation, most failures were ductile. Significant post-necking ductility was retained with reduction in area (RA) between 65–75%, but <50% for low temperature tempered steels. The ultimate tensile stress to yield stress (σ UTS /σ YS ) ratios decreased significantly after irradiation, highlighting irradiation-induced strain hardening capacity reduction. No major effect of irradiation on the plastic instability stress (σ PIS ) and true fracture stress of the steels was observed. By comparing the tensile stresses in true stress units and with literature, the results suggest that RAFM steel designing should target materials with a large separation between σ PIS and σ YS , to ensure the materials can maintain large work hardening and uniform deformation capability after irradiation. The tensile data of the steels additionally revealed a compelling evidence of an inverse trend between the change in RA and increase in σ YS of the neutron irradiated Eurofer97 type steels.

36 MATERIALS SCIENCE↗
DOE OSTI2026

SOLPS-ITER modeling of a dedicated divertor for negative triangularity operation on DIII-D

The design of a new dedicated divertor for negative triangularity (NT) operation on DIII-D with neutral baffles and pumping is informed by SOLPS-ITER transport modeling. This dedicated NT divertor is the latest step in a progression of NT shapes with various divertor characteristics explored on DIII-D, including NT shapes at reduced triangularity and a campaign with stronger shaping that included new armored components on the outboard side. SOLPS simulations played a key role in these divertor designs. Interpretive simulations, using cross-field diffusivities constrained by experimental data in the NT Shelf shape were used to inform the design of the 2023 armor campaign components. A similar procedure used armor campaign data to predict conditions for the dedicated NT divertor. The predictive simulations were used to assess the divertor fluxes, detachment threshold, pumped flux, and neutral leakage. For the dedicated NT divertor, SOLPS simulations and two-point-modeling were used to show the relative impact of magnetic topology (mainly longer connection length) and divertor closure on the divertor conditions relative to the armor campaign. It is predicted that the dedicated NT divertor reaches detachment (measured by target ion flux rollover) at a lower upstream density (≈(1.75−−2.4)×1⁢019m−3) as compared to the armor campaign shape. For the preliminary design geometry, divertor closure reduces the neutral leakage by ≈10%. Parametric optimization indicating further ≈20%–60% improvement in the leakage flux and recycled flux crossing the pump entrance is possible for relatively minor changes to the divertor and baffle layout.

Lore, Jeremy [ORNL] (ORCID:000000029192465X)↗
DOE OSTI2022

An efficient iterative method for dynamical Ginzburg-Landau equations

Here we propose a new finite element approach to solving the time-dependent Ginzburg-Landau equations for superconductivity under the temporal gauge and design an efficient preconditioner for the Newton iteration of the resulting discrete system. It solves the vector magnetic potential in $H$(curl) space by the lowest order of the second kind N´ed´elec element. This approach offers a simple way to deal with the boundary conditions, leading to a stable and reliable performance for superconductor geomtry with reentrant corners. The bench-marking in numerical simulations verifies the efficiency of the proposed preconditioner, which can be employed to significantly speed up in large-scale computations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Tritium distributions on W-coated divertor tiles used in the third JET ITER-like wall campaign

Tritium (T) distributions on tungsten (W)-coated plasma-facing tiles used in the third ITER-like wall campaign (2015–2016) of the Joint European Torus (JET) were examined by means of an imaging plate technique and β-ray induced x-ray spectrometry, and they were compared with the distributions after the second (2013–2014) campaign. Strong enrichment of T in beryllium (Be) deposition layers was observed after the second campaign. In contrast, T distributions after the third campaign was more uniform though Be deposition layers were visually recognized. The one of the possible explanations is enhanced desorption of T from Be deposition layers due to higher tile temperatures caused by higher energy input in the third campaign.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Validation of SOLPS-ITER simulations with kinetic, fluid, and hybrid neutral models for JET-ILW low-confinement mode plasmas

For JET L-mode plasmas in low-recycling conditions (electron temperature at the outer strike point, T e,ot ≳ 30 eV ), SOLPS-ITER simulations agree within the error bars for the experimental profiles at the low-field side (LFS) divertor target. The peak Balmer-α (Dα ) emission in the LFS divertor agrees within the error bars of the KS3 filterscope diagnostic, but is approximately 30% lower than the peak value of the KT1 spectrometer. Simulations have been performed with fluid, kinetic, and hybrid models for the neutrals. The large fluid-kinetic discrepancies of more than a factor 2 are successfully corrected by using a hybrid fluid-kinetic approach, for which kinetic atoms are transferred to the fluid population when the local Knudsen number of the atom becomes smaller than a user-defined transition Knudsen number Kn t . The hybrid-kinetic discrepancies are limited to a few % for Kn t ≤ 100 . When increasing the upstream density to high-recycling conditions, at the onset of detachment ( T e,ot ≈ 5 eV ), the simulations predict more than a factor 2 lower peak ion saturation current to the LFS divertor than the experiments. Also the Dα emission is underpredicted with approximately a factor 2. For these high-recycling conditions, the fluid-kinetic discrepancies are limited to maximum 50%, which are again corrected by using the hybrid approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗