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

Coherent diffractive imaging using randomly coded masks

We experimentally demonstrate an extension to coherent diffractive imaging that encodes additional information through the use of a series of randomly coded masks, removing the need for typical object-domain constraints while guaranteeing a unique solution to the phase retrieval problem. Phase retrieval is performed using a numerical convex relaxation routine known as “PhaseCut,” an iterative algorithm known for its stability and for its ability to find the global solution, which can be found efficiently and which is robust to noise. The experiment is performed using a laser diode at 532.2 nm, enabling rapid prototyping for future X-ray synchrotron and even free electron laser experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Polarization effects on the electric properties of urea and thiourea molecules in solid phase

We present theoretical results for the dipole moment, linear polarizability, and first hyperpolarizability of the urea and thiourea molecules in solid phase. The in-crystal electric properties were determined by applying a supermolecule approach in combination with an iterative electrostatic scheme, in which the surrounding molecules are represented by point charges. It is found for both urea and thiourea molecules that the influence of the polarization effects is mild for the linear polarizability, but it is marked for the dipole moment and first hyperpolarizability. The replacement of oxygen atoms by sulfur atoms increases, in general, the electric responses. Our second-order Møller–Plesset perturbation theory based iterative scheme predicts for the in-crystal dipole moment of urea and thiourea the values of 7.54 and 9.19 D which are, respectively, increased by 61% and 58%, in comparison with the corresponding isolated values. The result for urea is in agreement with the available experimental result of 6.56 D. In addition, we present an estimate of macroscopic quantities considering explicit unit cells of urea and thiourea crystals including environment polarization effects. These supermolecule calculations take into account partially the exchange and dispersion effects. The results illustrate the role played by the electrostatic interactions on the static second-order nonlinear susceptibility of the urea crystal.

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

The glass transition in cured epoxy thermosets: A comparative molecular dynamics study in coarse-grained and atomistic resolution

We investigate the volumetric glass transition temperature T{sub g} in epoxy thermosets by means of molecular dynamics simulations. The epoxy thermosets consist of the resin bisphenol A diglycidyl ether and the hardener diethylenetriamine. A structure based coarse-grained (CG) force field has been derived using iterative Boltzmann inversion in order to facilitate simulations of larger length scales. We observe that T{sub g} increases clearly with the degree of cross-linking for all-atomistic (AA) and CG simulations. The transition T{sub g} in CG simulations of uncured mixtures is much lower than in AA-simulations due to the soft nature of the CG potentials, but increases all the more with the formation of rigid cross-links. Additional simulations of the CG mixtures in contact with a surface show the existence of an interphase region of about 3 nm thickness in which the network properties deviate significantly from the bulk. In accordance to experimental studies, we observe that T{sub g} is reduced in this interphase region and gradually increases to its bulk value with distance from the surface. The present study shows that the glass transition is a local phenomenon that depends on the network structure in the immediate environment.

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

Iterative quantum-classical path integral with dynamically consistent state hopping

We investigate the convergence of iterative quantum-classical path integral calculations in sluggish environments strongly coupled to a quantum system. The number of classical trajectories, thus the computational cost, grows rapidly (exponentially, unless filtering techniques are employed) with the memory length included in the calculation. We argue that the choice of the (single) trajectory branch during the time preceding the memory interval can significantly affect the memory length required for convergence. At short times, the trajectory branch associated with the reactant state improves convergence by eliminating spurious memory. We also introduce an instantaneous population-based probabilistic scheme which introduces state-to-state hops in the retained pre-memory trajectory branch, and which is designed to choose primarily the trajectory branch associated with the reactant at early times, but to favor the product state more as the reaction progresses to completion. Test calculations show that the dynamically consistent state hopping scheme leads to accelerated convergence and a dramatic reduction of computational effort.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Local conditions for the Pauli potential in order to yield self-consistent electron densities exhibiting proper atomic shell structure

The local conditions for the Pauli potential that are necessary in order to yield self-consistent electron densities from orbital-free calculations are investigated for approximations that are expressed with the help of a local position variable. It is shown that those local conditions also apply when the Pauli potential is given in terms of the electron density. An explicit formula for the Ne atom is given, preserving the local conditions during the iterative procedure. The resulting orbital-free electron density exhibits proper shell structure behavior and is in close agreement with the Kohn-Sham electron density. This study demonstrates that it is possible to obtain self-consistent orbital-free electron densities with proper atomic shell structure from simple one-point approximations for the Pauli potential at local density level.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Absorption spectroscopy of heavy alkaline earth metals Ba and Sr in rare gas matrices—CCSD(T) calculations and atomic site occupancies

Isolation of the heavier alkaline earth metals Ba and Sr in the solid rare gases (RGs) Ar, Kr, and Xe is analysed with absorption spectroscopy and interpreted partly with the assistance of ab initio calculations of the diatomic M ⋅ RG ground state interaction potentials. The y{sup 1}P←a{sup 1}S resonance transitions in the visible spectral region are used to compare the isolation conditions of these two metal atom systems and calcium. Complex absorption bands were recorded in all three metal atom systems even after extensive sample annealing. Coupled cluster calculations conducted on the ground states of the nine M ⋅ RG diatomics (M = Ca, Sr, and Ba; RG = Ar, Kr, and Xe) at the coupled cluster single, double, and non-iterative triple level of theory revealed long bond lengths (>5 Å) and shallow bound regions (<130 cm{sup −1}). All of the M ⋅ RG diatomics have bond lengths considerably longer than those of the rare gas dimers, with the consequence that isolation of these metal atoms in a single substitutional site of the solid rare gas is unlikely, with the possible exception of Ca/Xe. The luminescence of metal dimer bands has been recorded for Ba and Sr revealing very different behaviours. Resonance fluorescence with a lifetime of 15 ns is observed for the lowest energy transition of Sr{sub 2} while this transition is quenched in Ba{sub 2}. This behaviour is consistent with the absence of vibrational structure on the dimer absorption band in Ba{sub 2} indicating lifetime broadening arising from efficient relaxation to low-lying molecular states. More extensive 2D excitation-emission data recorded for the complex site structures present on the absorption bands of the atomic Ba and Sr systems will be presented in future publications.

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

Design of a multi beam klystron cavity from its single beam parameters

The klystron is a well-known microwave amplifier which uses kinetic energy of an electron beam for amplification of the RF signal. There are some limitations of conventional single beam klystron such as high operating voltage, low efficiency and bulky size at higher power levels, which are very effectively handled in Multi Beam Klystron (MBK) that uses multiple low purveyance electron beams for RF interaction. Each beam propagates along its individual transit path through a resonant cavity structure. Multi-Beam klystron cavity design is a critical task due to asymmetric cavity structure and can be simulated by 3D code only. The present paper shall discuss the design of multi beam RF cavities for klystrons operating at 2856 MHz (S-band) and 5 GHz (C-band) respectively. The design approach uses some scaling laws for finding the electron beam parameters of the multi beam device from their single beam counter parts. The scaled beam parameters are then used for finding the design parameters of the multi beam cavities. Design of the desired multi beam cavity can be optimized through iterative simulations in CST Microwave Studio.

42 ENGINEERING↗
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↗