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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Landau ghost pole problem in quantum field theory: From 50th of last century to the present day

In this paper we present our results of the investigation of asymptotical behavior of amplitude at short distances in four-dimensional scalar field theory with ϕ{sup 4} interaction. To formulate of our calculating model – two-particle approximation of the mean-field expansion we have used an Rochev’s iteration scheme of solution of the Schwinger-Dyson equations with the fermion bilocal source. We have considered the nonlinear integral equations in deep-inelastic region of momenta. As result we have a non-trivial behavior of amplitude at large momenta.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Communication: A difference density picture for the self-consistent field ansatz

We formulate self-consistent field (SCF) theory in terms of an interaction picture where the working variable is the difference density matrix between the true system and a corresponding superposition of atomic densities. As the difference density matrix directly represents the electronic deformations inherent in chemical bonding, this “difference self-consistent field (dSCF)” picture provides a number of significant conceptual and computational advantages. We show that this allows for a stable and efficient dSCF iterative procedure with wholly single-precision Coulomb and exchange matrix builds. We also show that the dSCF iterative procedure can be performed with aggressive screening of the pair space. These approximations are tested and found to be accurate for systems with up to 1860 atoms and >10 000 basis functions, providing for immediate overall speedups of up to 70% in the heavily optimized TERACHEM SCF implementation.

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

Communication: Toward an improved control of the fixed-node error in quantum Monte Carlo: The case of the water molecule

All-electron Fixed-node Diffusion Monte Carlo calculations for the nonrelativistic ground-state energy of the water molecule at equilibrium geometry are presented. The determinantal part of the trial wavefunction is obtained from a selected Configuration Interaction calculation [Configuration Interaction using a Perturbative Selection done Iteratively (CIPSI) method] including up to about 1.4 × 10{sup 6} of determinants. Calculations are made using the cc-pCVnZ family of basis sets, with n = 2 to 5. In contrast with most quantum Monte Carlo works no re-optimization of the determinantal part in presence of a Jastrow is performed. For the largest cc-pCV5Z basis set the lowest upper bound for the ground-state energy reported so far of −76.437 44(18) is obtained. The fixed-node energy is found to decrease regularly as a function of the cardinal number n and the Complete Basis Set limit associated with exact nodes is easily extracted. The resulting energy of −76.438 94(12) — in perfect agreement with the best experimentally derived value — is the most accurate theoretical estimate reported so far. We emphasize that employing selected configuration interaction nodes of increasing quality in a given family of basis sets may represent a simple, deterministic, reproducible, and systematic way of controlling the fixed-node error in diffusion Monte Carlo.

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

Communication: An adaptive configuration interaction approach for strongly correlated electrons with tunable accuracy

We introduce a new procedure for iterative selection of determinant spaces capable of describing highly correlated systems. This adaptive configuration interaction (ACI) determines an optimal basis by an iterative procedure in which the determinant space is expanded and coarse grained until self-consistency. Two importance criteria control the selection process and tune the ACI to a user-defined level of accuracy. The ACI is shown to yield potential energy curves of N{sub 2} with nearly constant errors, and it predicts singlet-triplet splittings of acenes up to decacene that are in good agreement with the density matrix renormalization group.

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

Nano-metrology: The art of measuring X-ray mirrors with slope errors <100 nrad

We present a comprehensive investigation of the systematic and random errors of the nano-metrology instruments used to characterize synchrotron X-ray optics at Diamond Light Source. With experimental skill and careful analysis, we show that these instruments used in combination are capable of measuring state-of-the-art X-ray mirrors. Examples are provided of how Diamond metrology data have helped to achieve slope errors of <100 nrad for optical systems installed on synchrotron beamlines, including: iterative correction of substrates using ion beam figuring and optimal clamping of monochromator grating blanks in their holders. Simulations demonstrate how random noise from the Diamond-NOM’s autocollimator adds into the overall measured value of the mirror’s slope error, and thus predict how many averaged scans are required to accurately characterize different grades of mirror.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Full-dimensional vibrational calculations of five-atom molecules using a combination of Radau and Jacobi coordinates: Applications to methane and fluoromethane

Full quantum mechanical calculations of vibrational energies of methane and fluoromethane are carried out using a polyspherical description combining Radau and Jacobi coordinates. The Hamiltonian is built in a potential-optimized discrete variable representation, and vibrational energies are solved using an iterative eigensolver. This new approach can be applied to a large variety of molecules. In particular, we show that it is able to accurately and efficiently compute eigenstates for four different molecules : CH{sub 4}, CHD{sub 3}, CH{sub 2}D{sub 2}, and CH{sub 3}F. Very good agreement is obtained with the results reported previously in the literature with different approaches and with experimental data.

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

Hierarchical models and iterative optimization of hybrid systems

A class of hybrid control systems on the base of two-level discrete-continuous model is considered. The concept of this model was proposed and developed in preceding works as a concretization of the general multi-step system with related optimality conditions. A new iterative optimization procedure for such systems is developed on the base of localization of the global optimality conditions via contraction the control set.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Combined dry plasma etching and online metrology for manufacturing highly focusing x-ray mirrors

A new figuring station was designed and installed at the ESRF beamline BM05. It allows the figuring of mirrors within an iterative process combining the advantage of online metrology with dry etching. The complete process takes place under a vacuum environment to minimize surface contamination while non-contact surfacing tools open up the possibility of performing at-wavelength metrology and eliminating placement errors. The aim is to produce mirrors whose slopes do not deviate from the stigmatic profile by more than 0.1 µrad rms while keeping surface roughness in the acceptable limit of 0.1-0.2 nm rms. The desired elliptical mirror surface shape can be achieved in a few iterations in about a one day time span. This paper describes some of the important aspects of the process regarding both the online metrology and the etching process.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Microtearing modes in tokamak discharges

Microtearing modes (MTMs) have been identified as a source of significant electron thermal transport in tokamak discharges. In order to describe the evolution of these discharges, it is necessary to improve the prediction of electron thermal transport. This can be accomplished by utilizing a model for transport driven by MTMs in whole device predictive modeling codes. The objective of this paper is to develop the dispersion relation that governs the MTM driven transport. A unified fluid/kinetic approach is used in the development of a nonlinear dispersion relation for MTMs. The derivation includes the effects of electrostatic and magnetic fluctuations, arbitrary electron-ion collisionality, electron temperature and density gradients, magnetic curvature, and the effects associated with the parallel propagation vector. An iterative nonlinear approach is used to calculate the distribution function employed in obtaining the nonlinear parallel current and the nonlinear dispersion relation. The third order nonlinear effects in magnetic fluctuations are included, and the influence of third order effects on a multi-wave system is considered. An envelope equation for the nonlinear microtearing modes in the collision dominant limit is introduced in order to obtain the saturation level. In the limit that the mode amplitude does not vary along the field line, slab geometry, and strong collisionality, the fluid dispersion relation for nonlinear microtearing modes is found to agree with the kinetic dispersion relation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Optimization of current waveform tailoring for magnetically driven isentropic compression experiments

The Thor pulsed power generator is being developed at Sandia National Laboratories. The design consists of up to 288 decoupled and transit time isolated capacitor-switch units, called “bricks,” that can be individually triggered to achieve a high degree of pulse tailoring for magnetically driven isentropic compression experiments (ICE) [D. B. Reisman et al., Phys. Rev. Spec. Top.–Accel. Beams 18, 090401 (2015)]. The connecting transmission lines are impedance matched to the bricks, allowing the capacitor energy to be efficiently delivered to an ICE strip-line load with peak pressures of over 100 GPa. Thor will drive experiments to explore equation of state, material strength, and phase transition properties of a wide variety of materials. We present an optimization process for producing tailored current pulses, a requirement for many material studies, on the Thor generator. This technique, which is unique to the novel “current-adder” architecture used by Thor, entirely avoids the iterative use of complex circuit models to converge to the desired electrical pulse. We begin with magnetohydrodynamic simulations for a given material to determine its time dependent pressure and thus the desired strip-line load current and voltage. Because the bricks are connected to a central power flow section through transit-time isolated coaxial cables of constant impedance, the brick forward-going pulses are independent of each other. We observe that the desired equivalent forward-going current driving the pulse must be equal to the sum of the individual brick forward-going currents. We find a set of optimal brick delay times by requiring that the L{sub 2} norm of the difference between the brick-sum current and the desired forward-going current be a minimum. We describe the optimization procedure for the Thor design and show results for various materials of interest.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Assessing the utility of phase-space-localized basis functions: Exploiting direct product structure and a new basis function selection procedure

In this paper we show that it is possible to use an iterative eigensolver in conjunction with Halverson and Poirier’s symmetrized Gaussian (SG) basis [T. Halverson and B. Poirier, J. Chem. Phys. 137, 224101 (2012)] to compute accurate vibrational energy levels of molecules with as many as five atoms. This is done, without storing and manipulating large matrices, by solving a regular eigenvalue problem that makes it possible to exploit direct-product structure. These ideas are combined with a new procedure for selecting which basis functions to use. The SG basis we work with is orders of magnitude smaller than the basis made by using a classical energy criterion. We find significant convergence errors in previous calculations with SG bases. For sum-of-product Hamiltonians, SG bases large enough to compute accurate levels are orders of magnitude larger than even simple pruned bases composed of products of harmonic oscillator functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Approach to the unfolding and folding dynamics of add A-riboswitch upon adenine dissociation using a coarse-grained elastic network model

Riboswitches are noncoding mRNA segments that can regulate the gene expression via altering their structures in response to specific metabolite binding. We proposed a coarse-grained Gaussian network model (GNM) to examine the unfolding and folding dynamics of adenosine deaminase (add) A-riboswitch upon the adenine dissociation, in which the RNA is modeled by a nucleotide chain with interaction networks formed by connecting adjoining atomic contacts. It was shown that the adenine binding is critical to the folding of the add A-riboswitch while the removal of the ligand can result in drastic increase of the thermodynamic fluctuations especially in the junction regions between helix domains. Under the assumption that the native contacts with the highest thermodynamic fluctuations break first, the iterative GNM simulations showed that the unfolding process of the adenine-free add A-riboswitch starts with the denature of the terminal helix stem, followed by the loops and junctions involving ligand binding pocket, and then the central helix domains. Despite the simplified coarse-grained modeling, the unfolding dynamics and pathways are shown in close agreement with the results from atomic-level MD simulations and the NMR and single-molecule force spectroscopy experiments. Overall, the study demonstrates a new avenue to investigate the binding and folding dynamics of add A-riboswitch molecule which can be readily extended for other RNA molecules.

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

A spectrometer for ultrashort gamma-ray pulses with photon energies greater than 10 MeV

Herein we present a design for a pixelated scintillator based gamma-ray spectrometer for non-linear inverse Compton scattering experiments. By colliding a laser wakefield accelerated electron beam with a tightly focused, intense laser pulse, gamma-ray photons up to 100 MeV energies and with few femtosecond duration may be produced. To measure the energy spectrum and angular distribution, a 33 × 47 array of cesium-iodide crystals was oriented such that the 47 crystal length axis was parallel to the gamma-ray beam and the 33 crystal length axis was oriented in the vertical direction. Using an iterative deconvolution method similar to the YOGI code, modeling of the scintillator response using GEANT4 and fitting to a quantum Monte Carlo calculated photon spectrum, we are able to extract the gamma ray spectra generated by the inverse Compton interaction.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Active learning for robust, high-complexity reactive atomistic simulations

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.

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

Recent Advances in Small Angle X-ray Scattering for Superlattice Study

Small-angle x-ray scattering is used for the structure determination of superlattice for its superior resolution, nondestructive nature, and high penetration power of x rays. With the advent of high brilliance x-ray sources and innovative computing algorithms, there have been notable advances in small angle x-ray scattering analysis of superlattices. High brilliance x-ray beams have made data analyses less model-dependent. Additionally, novel data acquisition systems are faster and more competitive than ever before, enabling a more accurate mapping of the superlattices' reciprocal space. Fast and high-throughput computing systems and algorithms also make possible advanced analysis methods, including iterative phasing algorithms, non-parameterized fitting of scattering data with molecular dynamics simulations, and the use of machine learning algorithms. As a result, solving nanoscale structures with high resolutions has become an attainable task. In this review, we highlight new developments in the field and introduce their applications for the analysis of nanoscale ordered structures, including nanoparticle supercrystals, nanoscale lithography patterns, and supramolecular self-assemblies. Particularly, we highlight the reciprocal space mapping techniques and the use of iterative phase retrieval algorithms. We also cover coherent-beam-based small angle x-ray scattering techniques such as ptychography and ptycho-tomography in view of the traditional small angle x-ray scattering perspective.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Drift-ordered fluid vorticity equation with energy consistency

Although drift-ordered fluid models are widely applied in tokamak edge turbulence simulations, the models used are acknowledged not to conserve energy or even electrical charge. The present paper aims to remove many of the existing pitfalls in drift-fluid models, however, with the objective of finding a solution simple enough to be implemented in numerical applications. Our main result is an improved version of the drift-Braginskii equations involving a generalized vorticity function. In the new drift-Braginskii system, the quasi-neutrality condition translates into a transport equation for a generalized vorticity, expressed in conservation form, and related to the total mass-weighted circulation. It is found that kinetic energy conservation can be achieved if the polarization flow is defined recursively. In conclusion, the resulting model conserves the kinetic energy associated with E×B and diamagnetic flows and retains the associated perpendicular kinetic energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Desmearing Bonse–Hart USANS data using Bayesian Gaussian process regression

Ultra-small-angle neutron scattering (USANS) enables access to micrometer-scale structures but is intrinsically affected by strong, anisotropic resolution smearing arising from slit-geometry optics. As a result, recovery of the intrinsic scattering intensity constitutes an ill-posed inverse problem, and commonly used iterative desmearing methods lack rigorous uncertainty quantification. We present a Bayesian desmearing framework for slit-geometry USANS based on Gaussian process regression. In this approach, the scattering intensity is modeled as a smooth random function, and the instrumental point spread function is incorporated explicitly as a forward operator. The resulting formulation yields a closed-form maximum a posteriori solution with well-defined credibility intervals. Computational benchmarks and experimental validation using combined USANS and small-angle neutron scattering (SANS) measurements demonstrate that the framework enables stable desmearing, suppresses experimental noise, and preserves physically meaningful structural features under realistic conditions.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗