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

Secondary electron emission and self-consistent charge transport in semi-insulating samples

Electron beam induced self-consistent charge transport and secondary electron emission (SEE) in insulators are described by means of an electron-hole flight-drift model (FDM) now extended by a certain intrinsic conductivity (c) and are implemented by an iterative computer simulation. Ballistic secondary electrons (SE) and holes, their attenuation to drifting charge carriers, and their recombination, trapping, and field- and temperature-dependent detrapping are included. As a main result the time dependent ''true'' secondary electron emission rate {delta}(t) released from the target material and based on ballistic electrons and the spatial distributions of currents j(x,t), charges {rho}(x,t), field F(x,t), and potential V(x,t) are obtained where V{sub 0} = V(0,t) presents the surface potential. The intrinsic electronic conductivity limits the charging process and leads to a conduction sample current to the support. In that case the steady-state total SE yield will be fixed below the unit: i.e., {sigma} {eta} + {delta} < 1.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Communication: Multistate quantum dynamics of photodissociation of carbon dioxide between 120 nm and 160 nm

UV absorption cross section of CO{sub 2} is studied using high level ab initio quantum chemistry for electrons and iterative quantum dynamics for nuclear motion on interacting global full dimensional potential energy surfaces. Six electronic states-1, 2, 3{sup 1}A{sup Prime} and 1, 2, 3{sup 1}A{sup Double-Prime }-are considered. At linearity, they correspond to the ground electronic state X(tilde sign){sup 1}{Sigma}{sub g}{sup +} and the optically forbidden but vibronically allowed valence states 1{sup 1}{Delta}{sub u}, 1{sup 1}{Sigma}{sub u}{sup -}, and 1{sup 1}{Pi}{sub g}. In the Franck-Condon region, these states interact via Renner-Teller and conical intersections and are simultaneously involved in an intricate network of non-adiabatic couplings. The absorption spectrum, calculated for many rotational states, reproduces the distinct two-band shape of the experimental spectrum measured at 190 K and the characteristic patterns of the diffuse structures in each band. Quantum dynamics unravel the relative importance of different vibronic mechanisms, while metastable resonance states, underlying the diffuse structures, provide dynamically based vibronic assignments of individual lines.

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

Pulsed laser deposition to synthesize the bridge structure of artificial nacre: Comparison of nano- and femtosecond lasers

Nature offers inspiration to new adaptive technologies that allow us to build amazing shapes and structures such as nacre using synthetic materials. Consequently, we have designed a pulsed laser ablation manufacturing process involving thin film deposition and micro-machining to create hard/soft layered 'brick-bridge-mortar' nacre of AlMgB{sub 14} (hard phase) with Ti (soft phase). In this paper, we report pulsed laser deposition (PLD) to mimic brick and bridge structures of natural nacre in AlMgB{sub 14}. Particulate formation inherent in PLD is exploited to develop the bridge structure. Mechanical behavior analysis of the AlMgB{sub 14}/Ti system revealed that the brick is to be 250 nm thick, 9 {mu}m lateral dimensions while the bridge (particle) is to have a diameter of 500 nm for a performance equivalent to natural nacre. Both nanosecond (ns) and femtosecond (fs) pulsed lasers were employed for PLD in an iterative approach that involves varying pulse energy, pulse repetition rate, and target-to-substrate distance to achieve the desired brick and bridge characteristics. Scanning electron microscopy, x-ray photoelectron spectroscopy, and optical profilometer were used to evaluate the film thickness, particle size and density, stoichiometry, and surface roughness of thin films. Results indicated that both ns-pulsed and fs-pulsed lasers produce the desired nacre features. However, each laser may be chosen for different reasons: fs-pulsed laser is preferred for much shorter deposition time, better stoichiometry, uniform-sized particles, and uniform film thickness, while ns-pulsed laser is favored for industrial acceptance, reliability, ease of handling, and low cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Analytic and algorithmic aspects of generalized harmonic sums and polylogarithms

In recent three-loop calculations of massive Feynman integrals within Quantum Chromodynamics (QCD) and, e.g., in recent combinatorial problems the so-called generalized harmonic sums (in short S-sums) arise. They are characterized by rational (or real) numerator weights also different from ±1. In this article we explore the algorithmic and analytic properties of these sums systematically. We work out the Mellin and inverse Mellin transform which connects the sums under consideration with the associated Poincaré iterated integrals, also called generalized harmonic polylogarithms. In this regard, we obtain explicit analytic continuations by means of asymptotic expansions of the S-sums which started to occur frequently in current QCD calculations. In addition, we derive algebraic and structural relations, like differentiation with respect to the external summation index and different multi-argument relations, for the compactification of S-sum expressions. Finally, we calculate algebraic relations for infinite S-sums, or equivalently for generalized harmonic polylogarithms evaluated at special values. The corresponding algorithms and relations are encoded in the computer algebra package HarmonicSums.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Evaluating the performance of two neutron spectrum unfolding codes based on iterative procedures and artificial neural networks

In this work the performance of two neutron spectrum unfolding codes based on iterative procedures and artificial neural networks is evaluated. The first one code based on traditional iterative procedures and called Neutron spectrometry and dosimetry from the Universidad Autonoma de Zacatecas (NSDUAZ) use the SPUNIT iterative algorithm and was designed to unfold neutron spectrum and calculate 15 dosimetric quantities and 7 IAEA survey meters. The main feature of this code is the automated selection of the initial guess spectrum trough a compendium of neutron spectrum compiled by the IAEA. The second one code known as Neutron spectrometry and dosimetry with artificial neural networks (NDSann) is a code designed using neural nets technology. The artificial intelligence approach of neural net does not solve mathematical equations. By using the knowledge stored at synaptic weights on a neural net properly trained, the code is capable to unfold neutron spectrum and to simultaneously calculate 15 dosimetric quantities, needing as entrance data, only the rate counts measured with a Bonner spheres system. Similarities of both NSDUAZ and NSDann codes are: they follow the same easy and intuitive user's philosophy and were designed in a graphical interface under the LabVIEW programming environment. Both codes unfold the neutron spectrum expressed in 60 energy bins, calculate 15 dosimetric quantities and generate a full report in HTML format. Differences of these codes are: NSDUAZ code was designed using classical iterative approaches and needs an initial guess spectrum in order to initiate the iterative procedure. In NSDUAZ, a programming routine was designed to calculate 7 IAEA instrument survey meters using the fluence-dose conversion coefficients. NSDann code use artificial neural networks for solving the ill-conditioned equation system of neutron spectrometry problem through synaptic weights of a properly trained neural network. Contrary to iterative procedures, in neural net approach it is possible to reduce the rate counts used to unfold the neutron spectrum. To evaluate these codes a computer tool called Neutron Spectrometry and dosimetry computer tool was designed. The results obtained with this package are showed. The codes here mentioned are freely available upon request to the authors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Photoionization cross section by Stieltjes imaging applied to coupled cluster Lanczos pseudo-spectra

A recently implemented asymmetric Lanczos algorithm for computing (complex) linear response functions within the coupled cluster singles (CCS), coupled cluster singles and iterative approximate doubles (CC2), and coupled cluster singles and doubles (CCSD) is coupled to a Stieltjes imaging technique in order to describe the photoionization cross section of atoms and molecules, in the spirit of a similar procedure recently proposed by Averbukh and co-workers within the Algebraic Diagrammatic Construction approach. Pilot results are reported for the atoms He, Ne, and Ar and for the molecules H{sub 2}, H{sub 2}O, NH{sub 3}, HF, CO, and CO{sub 2}.

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

Response of a partial wall to an external perturbation of rotating plasma

In this paper, we present the response of a 3D thin multiply connected wall to an external kink mode perturbation in axisymmetric tokamak configurations. To calculate the contribution of the plasma perturbed magnetic field in the vacuum region, we have made use of the concept of surface currents [following C. V. Atanasiu, A. H. Boozer, L. E. Zakharov, and A. A. Subbotin, Phys. Plasmas 6, 2781 (1999)]. The wall response is expressed in terms of a stream function of the wall surface currents, which are obtained by solving a diffusion type equation, taking into account the contribution of the wall currents themselves iteratively. The use of stream function makes the approach applicable for both well-studied earlier Resistive Wall Modes and for Wall Touching Kink Modes, which were discovered recently as a key phenomenon in disruptions [L. E. Zakharov, S. A. Galkin, and S. N. Gerasimov, Phys. Plasmas 19, 055703 (2012)]. New analytical expressions, suitable for numerical calculations of toroidal harmonics of the vacuum magnetic fields from the surface currents on axisymmetric shells, are derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Mechanical-plowing-based high-speed patterning on hard material via advanced-control and ultrasonic probe vibration

In this paper, we present a high-speed direct pattern fabrication on hard materials (e.g., a tungsten-coated quartz substrate) via mechanical plowing. Compared to other probe-based nanolithography techniques based on chemical- and/or physical-reactions (e.g., the Dip-pen technique), mechanical plowing is meritorious for its low cost, ease of process control, and capability of working with a wide variety of materials beyond conductive and/or soft materials. However, direct patterning on hard material faces two daunting challenges. First, the patterning throughput is ultimately hindered by the “writing” (plowing) speed, which, in turn, is limited by the adverse effects that can be excited/induced during high-speed, and/or large-range plowing, including the vibrational dynamics of the actuation system (the piezoelectric actuator, the cantilever, and the mechanical fixture connecting the cantilever to the actuator), the dynamic cross-axis coupling between different axes of motion, and the hysteresis and the drift effects related to the piezoelectric actuators. Secondly, it is very challenging to directly pattern on ultra-hard materials via plowing. Even with a diamond probe, the line depth of the pattern via continuous plowing on ultra-hard materials such as tungsten, is still rather small (<0.5 nm), particularly when the “writing” speed becomes high. To overcome these two challenges, we propose to utilize a novel iterative learning control technique to achieve precision tracking of the desired pattern during high-speed, large-range plowing, and introduce ultrasonic vibration of the probe in the normal (vertical) direction during the plowing process to enable direct patterning on ultra hard materials. The proposed approach was implemented to directly fabricate patterns on a mask with tungsten coating and quartz substrate. The experimental results demonstrated that a large-size pattern of four grooves (20 μm in length with 300 nm spacing between lines) can be fabricated at a high speed of ∼5 mm/s, with the line width and the line depth at ∼95 nm and 2 nm, respectively. A fine pattern of the word “NANO” is also fabricated at the speed of ∼5 mm/s.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Momentum space iterative solution of the time-dependent Schrödinger equation

We present a novel approach, the iterative solution of the time-dependent Schrödinger equation (iTDSE model), for the investigation of atomic systems interacting with external laser fields. This model is the extension of the momentum-space strong-field approximation (MSSFA) [1], in which the Coulomb potential was considered only as a first order perturbation. In the iTDSE approach higher order terms were gradually introduced until convergence was achieved. Benchmark calculations were done on the hydrogen atom, and the obtained results were compared to the direct numerical solution [2].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Effects of pressure anisotropy on magnetospheric magnetohydrodynamics equilibrium of an internal ring current system

Effects of pressure anisotropy on magnetospheric magnetohydrodynamics equilibrium are studied analytically and numerically, where the plasma is confined by only poloidal magnetic field generated by an internal ring current. The plasma current due to finite pressure can be divided into two components; one remains at isotropic pressure and the other arises from pressure anisotropy. When p{sub ⊥}, the pressure perpendicular to the magnetic field, is larger than p{sub ∥}, the pressure parallel to the magnetic field, those two components of plasma current tend to cancel each other to reduce the total amount of plasma current. Equilibrium beta limit is also examined, where the beta is a ratio of the plasma pressure to the magnetic pressure. The equilibrium beta limit decreases as the pressure anisotropy becomes strong. The beta value is strictly limited by ellipticity of the equilibrium equation when p{sub ∥}>p{sub ⊥}. On the other hand, when p{sub ⊥}>p{sub ∥}, although the tendency of the beta limit agrees with the ellipticity condition of the equilibrium equation, equilibria with a hyperbolic region can be obtained by iterative procedure with practically reasonable convergence criteria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Neutron source reconstruction from pinhole imaging at National Ignition Facility

The neutron imaging system at the National Ignition Facility (NIF) is an important diagnostic tool for measuring the two-dimensional size and shape of the neutrons produced in the burning deuterium-tritium plasma during the ignition stage of inertial confinement fusion (ICF) implosions at NIF. Since the neutron source is small (∼100 μm) and neutrons are deeply penetrating (>3 cm) in all materials, the apertures used to achieve the desired 10-μm resolution are 20-cm long, single-sided tapers in gold. These apertures, which have triangular cross sections, produce distortions in the image, and the extended nature of the pinhole results in a non-stationary or spatially varying point spread function across the pinhole field of view. In this work, we have used iterative Maximum Likelihood techniques to remove the non-stationary distortions introduced by the aperture to reconstruct the underlying neutron source distributions. We present the detailed algorithms used for these reconstructions, the stopping criteria used and reconstructed sources from data collected at NIF with a discussion of the neutron imaging performance in light of other diagnostics.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

Approximate equations of state in two-temperature plasma mixtures

Approximate thermodynamic state relations for multicomponent atomic and molecular gas mixtures are often constructed by artificially partitioning the mixture into its constituent materials and requiring the separated materials to be in temperature and pressure equilibrium. Iterative numerical algorithms have been employed to enforce this equilibration and compute the resulting approximate state relations in single-temperature mixtures. In partially ionized gas mixtures, there is both theoretical and empirical evidence that equilibrating the chemical potentials, number densities, or partial pressures of the free electrons is likely to produce more accurate results than equilibrating the total pressures. Moreover, in many situations of practical interest the free electrons and heavy particles have different temperatures. In this paper, we present a generalized algorithm for equilibrating the heavy-particle and electron temperatures and a third user-specified independent thermodynamic variable in a two-temperature plasma mixture. Test calculations based on the equilibration of total pressure vs. electron pressure are presented for three different mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

The chemical abundances of the Ap star HD94660

In this work I present the determination of chemical abundances of the Ap star HD94660, a possible rapid oscillating star. As all the magnetic chemically peculiar objects, it presents CNO underabundance and overabundance of iron peak elements of ∼100 times and of rare earths up to 4 dex with respect to the Sun. The determination was based on the conversion of the observed equivalent widths into abundances simultaneously to the determination of effective temperature and gravity. Since the Balmer lines of early type stars are very sensitive to the surface gravity while the flux distribution is sensitive to the effective temperature, I have adopted an iterative procedure to match the H{sub α} line profile and the observed UV-Vis-NIR magnitudes of HD94660 looking for a consistency between the metallicity of the atmosphere model and the derived abundances. From my spectroscopic analysis, this star belongs to the no-rapid oscillating class.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Continuous-scanning laser Doppler vibrometry: Extensions to arbitrary areas, multi-frequency and 3D capture

To date, differing implementations of continuous scan laser Doppler vibrometry have been demonstrated by various academic institutions, but since the scan paths were defined using step or sine functions from function generators, the paths were typically limited to 1D line scans or 2D areas such as raster paths or Lissajous trajectories. The excitation was previously often limited to a single frequency due to the specific signal processing performed to convert the scan data into an ODS. In this paper, a configuration of continuous-scan laser Doppler vibrometry is demonstrated which permits scanning of arbitrary areas, with the benefit of allowing multi-frequency/broadband excitation. Various means of generating scan paths to inspect arbitrary areas are discussed and demonstrated. Further, full 3D vibration capture is demonstrated by the addition of a range-finding facility to the described configuration, and iteratively relocating a single scanning laser head. Here, the range-finding facility was provided by a Microsoft Kinect, an inexpensive piece of consumer electronics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Iterative reconstruction using a Monte Carlo based system transfer matrix for dedicated breast positron emission tomography

To maximize sensitivity, it is desirable that ring Positron Emission Tomography (PET) systems dedicated for imaging the breast have a small bore. Unfortunately, due to parallax error this causes substantial degradation in spatial resolution for objects near the periphery of the breast. In this work, a framework for computing and incorporating an accurate system matrix into iterative reconstruction is presented in an effort to reduce spatial resolution degradation towards the periphery of the breast. The GATE Monte Carlo Simulation software was utilized to accurately model the system matrix for a breast PET system. A strategy for increasing the count statistics in the system matrix computation and for reducing the system element storage space was used by calculating only a subset of matrix elements and then estimating the rest of the elements by using the geometric symmetry of the cylindrical scanner. To implement this strategy, polar voxel basis functions were used to represent the object, resulting in a block-circulant system matrix. Simulation studies using a breast PET scanner model with ring geometry demonstrated improved contrast at 45% reduced noise level and 1.5 to 3 times resolution performance improvement when compared to MLEM reconstruction using a simple line-integral model. The GATE based system matrix reconstruction technique promises to improve resolution and noise performance and reduce image distortion at FOV periphery compared to line-integral based system matrix reconstruction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Phonon transport in single-layer transition metal dichalcogenides: A first-principles study

Two-dimensional transition metal dichalcogenides (TMDCs) are finding promising electronic and optical applications due to their unique properties. In this letter, we systematically study the phonon transport and thermal conductivity of eight semiconducting single-layer TMDCs, MX{sub 2} (M = Mo, W, Zr, and Hf, X = S and Se), by using the first-principles-driven phonon Boltzmann transport equation approach. The validity of the single-mode relaxation time approximation to predict the thermal conductivity of TMDCs is assessed by comparing the results with the iterative solution of the phonon Boltzmann transport equation. We find that the phononic thermal conductivities of 2H-type TMDCs are above 50 W/mK at room temperature while the thermal conductivity values of the 1T-type TMDCs are much lower, when the size of the sample is 1 μm. A very high thermal conductivity value of 142 W/mK was found in single-layer WS{sub 2}. The large atomic weight difference between W and S leads to a very large phonon bandgap which in turn forbids the scattering between acoustic and optical phonon modes and thus resulting in very long phonon relaxation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Blip decomposition of the path integral: Exponential acceleration of real-time calculations on quantum dissipative systems

The real-time path integral representation of the reduced density matrix for a discrete system in contact with a dissipative medium is rewritten in terms of the number of blips, i.e., elementary time intervals over which the forward and backward paths are not identical. For a given set of blips, it is shown that the path sum with respect to the coordinates of all remaining time points is isomorphic to that for the wavefunction of a system subject to an external driving term and thus can be summed by an inexpensive iterative procedure. This exact decomposition reduces the number of terms by a factor that increases exponentially with propagation time. Further, under conditions (moderately high temperature and/or dissipation strength) that lead primarily to incoherent dynamics, the “fully incoherent limit” zero-blip term of the series provides a reasonable approximation to the dynamics, and the blip series converges rapidly to the exact result. Retention of only the blips required for satisfactory convergence leads to speedup of full-memory path integral calculations by many orders of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A coupled-cluster study of photodetachment cross sections of closed-shell anions

We investigate the performance of Stieltjes Imaging applied to Lanczos pseudo-spectra generated at the coupled cluster singles and doubles, coupled cluster singles and approximate iterative doubles and coupled cluster singles levels of theory in modeling the photodetachment cross sections of the closed shell anions H{sup −}, Li{sup −}, Na{sup −}, F{sup −}, Cl{sup −}, and OH{sup −}. The accurate description of double excitations is found to play a much more important role than in the case of photoionization of neutral species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗