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 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Steady state thermal-hydraulic analyses of the MITICA cooling circuits

Megavolt ITER Injector Concept Advancement is the full scale prototype of the heating and current drive neutral beam injectors for ITER, to be built at Consorzio RFX (Padova). The engineering design of its components is challenging: the total heat loads they will be subjected to (expected between 2 and 19 MW), the high heat fluxes (up to 20 MW/m{sup 2}), and the beam pulse duration up to 1 h, set demanding requirements for reliable active cooling circuits. In support of the design, the thermo-hydraulic behavior of each cooling circuit under steady state condition has been investigated by using one-dimensional models. The final results, obtained considering a number of optimizations for the cooling circuits, show that all the requirements in terms of flow rate, temperature, and pressure drop are properly fulfilled.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Quantum learning of classical stochastic processes: The completely positive realization problem

Among several tasks in Machine Learning, a specially important one is the problem of inferring the latent variables of a system and their causal relations with the observed behavior. A paradigmatic instance of this is the task of inferring the hidden Markov model underlying a given stochastic process. This is known as the positive realization problem (PRP), [L. Benvenuti and L. Farina, IEEE Trans. Autom. Control 49(5), 651–664 (2004)] and constitutes a central problem in machine learning. The PRP and its solutions have far-reaching consequences in many areas of systems and control theory, and is nowadays an important piece in the broad field of positive systems theory. We consider the scenario where the latent variables are quantum (i.e., quantum states of a finite-dimensional system) and the system dynamics is constrained only by physical transformations on the quantum system. The observable dynamics is then described by a quantum instrument, and the task is to determine which quantum instrument — if any — yields the process at hand by iterative application. We take as a starting point the theory of quasi-realizations, whence a description of the dynamics of the process is given in terms of linear maps on state vectors and probabilities are given by linear functionals on the state vectors. This description, despite its remarkable resemblance with the hidden Markov model, or the iterated quantum instrument, is however devoid of any stochastic or quantum mechanical interpretation, as said maps fail to satisfy any positivity conditions. The completely positive realization problem then consists in determining whether an equivalent quantum mechanical description of the same process exists. We generalize some key results of stochastic realization theory, and show that the problem has deep connections with operator systems theory, giving possible insight to the lifting problem in quotient operator systems. Our results have potential applications in quantum machine learning, device-independent characterization and reverse-engineering of stochastic processes and quantum processors, and more generally, of dynamical processes with quantum memory [M. Guţă, Phys. Rev. A 83(6), 062324 (2011); M. Guţă and N. Yamamoto, e-print http://arxiv.org/abs/1303.3771 (2013)].

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

FJK - A cap-free fragment approach with embedding Fock potentials

We present a cap-free fragment-based approach for Hartree-Fock calculations where full Fock potentials are used for monomer embedding. In this FJK (Fragmentation with Coulomb (J) and exchange (K) embedding) method the system is divided into fragments by heterolytic fission of covalent bonds. The one-electron operators for each fragment are built including the mean-field Coulomb and exchange potential of all the remaining monomers. The description is improved by iterative updates of the monomer potentials within which the method converges to the full system Hartree-Fock result. It is shown that a small number of iterations is sufficient to yield errors in the absolute energy comparable to those usually observed in density fitting procedures. Conformational energies of organic molecules are used as benchmark, demonstrating the robustness of the procedure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes.

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

Vapor shielding models and the energy absorbed by divertor targets during transient events

The erosion of divertor targets caused by high heat fluxes during transients is a serious threat to ITER operation, as it is going to be the main factor determining the divertor lifetime. Under the influence of extreme heat fluxes, the surface temperature of plasma facing components can reach some certain threshold, leading to an onset of intense material evaporation. The latter results in formation of cold dense vapor and secondary plasma cloud. This layer effectively absorbs the energy of the incident plasma flow, turning it into its own kinetic and internal energy and radiating it. This so called vapor shielding is a phenomenon that may help mitigating the erosion during transient events. In particular, the vapor shielding results in saturation of energy (per unit surface area) accumulated by the target during single pulse of heat load at some level E{sub max}. Matching this value is one of the possible tests to verify complicated numerical codes, developed to calculate the erosion rate during abnormal events in tokamaks. The paper presents three very different models of vapor shielding, demonstrating that E{sub max} depends strongly on the heat pulse duration, thermodynamic properties, and evaporation energy of the irradiated target material. While its dependence on the other shielding details such as radiation capabilities of material and dynamics of the vapor cloud is logarithmically weak. The reason for this is a strong (exponential) dependence of the target material evaporation rate, and therefore the “strength” of vapor shield on the target surface temperature. As a result, the influence of the vapor shielding phenomena details, such as radiation transport in the vapor cloud and evaporated material dynamics, on the E{sub max} is virtually completely masked by the strong dependence of the evaporation rate on the target surface temperature. However, the very same details define the amount of evaporated particles, needed to provide an effective shielding to the target, and, therefore, strongly influence resulting erosion rate. Thus, E{sub max} cannot be used for validation of shielding models and codes, aimed at the target material erosion calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Impurity-induced divertor plasma oscillations

Two different oscillatory plasma regimes induced by seeding the plasma with high- and low-Z impurities are found for ITER-like divertor plasmas, using computer modeling with the DUSTT/UEDGE and SOLPS4.3 plasma-impurity transport codes. The oscillations are characterized by significant variations of the impurity-radiated power and of the peak heat load on the divertor targets. Qualitative analysis of the divertor plasma oscillations reveals different mechanisms driving the oscillations in the cases of high- and low-Z impurity seeding. The oscillations caused by the high-Z impurities are excited near the X-point by an impurity-related instability of the radiation-condensation type, accompanied by parallel impurity ion transport affected by the thermal and plasma friction forces. The driving mechanism of the oscillations induced by the low-Z impurities is related to the cross-field transport of the impurity atoms, causing alteration between the high and low plasma temperature regimes in the plasma recycling region near the divertor targets. The implications of the impurity-induced plasma oscillations for divertor operation in the next generation tokamaks are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Probing off-Hugoniot states in Ta, Cu, and Al to 1000 GPa compression with magnetically driven liner implosions

We report on a new technique for obtaining off-Hugoniot pressure vs. density data for solid metals compressed to extreme pressure by a magnetically driven liner implosion on the Z-machine (Z) at Sandia National Laboratories. In our experiments, the liner comprises inner and outer metal tubes. The inner tube is composed of a sample material (e.g., Ta and Cu) whose compressed state is to be inferred. The outer tube is composed of Al and serves as the current carrying cathode. Another aluminum liner at much larger radius serves as the anode. A shaped current pulse quasi-isentropically compresses the sample as it implodes. The iterative method used to infer pressure vs. density requires two velocity measurements. Photonic Doppler velocimetry probes measure the implosion velocity of the free (inner) surface of the sample material and the explosion velocity of the anode free (outer) surface. These two velocities are used in conjunction with magnetohydrodynamic simulation and mathematical optimization to obtain the current driving the liner implosion, and to infer pressure and density in the sample through maximum compression. This new equation of state calibration technique is illustrated using a simulated experiment with a Cu sample. Monte Carlo uncertainty quantification of synthetic data establishes convergence criteria for experiments. Results are presented from experiments with Al/Ta, Al/Cu, and Al liners. Symmetric liner implosion with quasi-isentropic compression to peak pressure ∼1000 GPa is achieved in all cases. These experiments exhibit unexpectedly softer behavior above 200 GPa, which we conjecture is related to differences in the actual and modeled properties of aluminum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Self-consistent second-order Green’s function perturbation theory for periodic systems

Despite recent advances, systematic quantitative treatment of the electron correlation problem in extended systems remains a formidable task. Systematically improvable Green’s function methods capable of quantitatively describing weak and at least qualitatively strong correlations appear as promising candidates for computational treatment of periodic systems. We present a periodic implementation of temperature-dependent self-consistent 2nd-order Green’s function (GF2) method, where the self-energy is evaluated in the basis of atomic orbitals. Evaluating the real-space self-energy in atomic orbitals and solving the Dyson equation in k-space are the key components of a computationally feasible algorithm. We apply this technique to the one-dimensional hydrogen lattice — a prototypical crystalline system with a realistic Hamiltonian. By analyzing the behavior of the spectral functions, natural occupations, and self-energies, we claim that GF2 is able to recover metallic, band insulating, and at least qualitatively Mott regimes. We observe that the iterative nature of GF2 is essential to the emergence of the metallic and Mott phases.

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

Exploration of the Super H-mode regime on DIII-D and potential advantages for burning plasma devices

A new high pedestal regime (“Super H-mode”) has been predicted and accessed on DIII-D. Super H-mode was first achieved on DIII-D using a quiescent H-mode edge, enabling a smooth trajectory through pedestal parameter space. By exploiting Super H-mode, it has been possible to access high pedestal pressures at high normalized densities. While elimination of Edge localized modes (ELMs) is beneficial for Super H-mode, it may not be a requirement, as recent experiments have maintained high pedestals with ELMs triggered by lithium granule injection. Simulations using TGLF for core transport and the EPED model for the pedestal find that ITER can benefit from the improved performance associated with Super H-mode, with increased values of fusion power and gain possible. Similar studies demonstrate that the Super H-mode pedestal can be advantageous for a steady-state power plant, by providing a path to increasing the bootstrap current while simultaneously reducing the demands on the core physics performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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

Rotation profile flattening and toroidal flow shear reversal due to the coupling of magnetic islands in tokamaks

The electromagnetic coupling of helical modes, even those having different toroidal mode numbers, modifies the distribution of toroidal angular momentum in tokamak discharges. This can have deleterious effects on other transport channels as well as on magnetohydrodynamic (MHD) stability and disruptivity. At low levels of externally injected momentum, the coupling of core-localized modes initiates a chain of events, whereby flattening of the core rotation profile inside successive rational surfaces leads to the onset of a large m/n = 2/1 tearing mode and locked-mode disruption. With increased torque from neutral beam injection, neoclassical tearing modes in the core may phase-lock to each other without locking to external fields or structures that are stationary in the laboratory frame. The dynamic processes observed in these cases are in general agreement with theory, and detailed diagnosis allows for momentum transport analysis to be performed, revealing a significant torque density that peaks near the 2/1 rational surface. However, as the coupled rational surfaces are brought closer together by reducing q{sub 95}, additional momentum transport in excess of that required to attain a phase-locked state is sometimes observed. Rather than maintaining zero differential rotation (as is predicted to be dynamically stable by single-fluid, resistive MHD theory), these discharges develop hollow toroidal plasma fluid rotation profiles with reversed plasma flow shear in the region between the m/n = 3/2 and 2/1 islands. The additional forces expressed in this state are not readily accounted for, and therefore, analysis of these data highlights the impact of mode coupling on torque balance and the challenges associated with predicting the rotation dynamics of a fusion reactor—a key issue for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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

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

Transport and deceleration of fusion products in microturbulence

The velocity-space distribution of alpha particles born in fusion devices is subject to modification at moderate energies due to turbulent transport. Therefore, one must calculate the evolution of an equilibrium distribution whose functional form is not known a priori. Using a novel technique, applicable to any trace impurity, we have made this calculation for fully nonlinear gyrokinetic simulations not only possible but also particularly efficient. We demonstrate a microturbulence-induced departure from the local slowing-down distribution, an inversion of the energy distribution, and associated modifications to the alpha heating and pressure profiles in an ITER-like scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Thermal quench mitigation and current quench control by injection of mixed species shattered pellets in DIII-D

Injection of large shattered pellets composed of variable quantities of the main ion species (deuterium) and high-Z impurities (neon) in the DIII-D tokamak demonstrates control of thermal quench (TQ) and current quench (CQ) properties in mitigated disruptions. As the pellet composition is varied, TQ radiation fractions increase continuously with the quantity of radiating impurity in the pellet, with a corresponding decrease in divertor heating. Post-TQ plasma resistivities increase as a result of the higher radiation fraction, allowing control of current decay timescales based on the pellet composition. Magnetic reconstructions during the CQ show that control of the current decay rate allows continuous variation of the minimum safety factor during the vertically unstable disruption, reducing the halo current fraction and resulting vessel displacement. Both TQ and CQ characteristics are observed to saturate at relatively low quantities of neon, indicating that effective mitigation of disruption loads by shattered pellet injection (SPI) can be achieved with modest impurity quantities, within injection quantities anticipated for ITER. This mixed species SPI technique provides a possible approach for tuning disruption properties to remain within the limited ranges allowed in the ITER design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Exploration of magnetic perturbation effects on advanced divertor configurations in NSTX-U

The control of divertor heat loads—both steady state and transient—remains a key challenge for the successful operation of ITER and FNSF. Magnetic perturbations provide a promising technique to control ELMs (Edge Localized Modes) (transients), but understanding their detailed impact is difficult due to their symmetry breaking nature. One approach for reducing steady state heat loads is so called “advanced divertors” which aim at optimizing the magnetic field configuration: the snowflake and the (super-)X-divertor. It is likely that both concepts—magnetic perturbations and advanced divertors—will have to work together, and we explore their interaction based on the NSTX-U setup. An overview of different divertor configurations under the impact of magnetic perturbations is presented, and the resulting impact on plasma edge transport is investigated with the EMC3-EIRENE code. Variations in size of the magnetic footprint of the perturbed separatrix are found, which are related to the level of flux expansion on the divertor target. Non-axisymmetric peaking of the heat flux related to the perturbed separatrix is found at the outer strike point, but only in locations where flux expansion is not too large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Analysis of non-equilibrium phenomena in inductively coupled plasma generators

This article addresses the modeling of non-equilibrium phenomena in inductively coupled plasma discharges. In the proposed computational model, the electromagnetic induction equation is solved together with the set of Navier-Stokes equations in order to compute the electromagnetic and flow fields, accounting for their mutual interaction. Semi-classical statistical thermodynamics is used to determine the plasma thermodynamic properties, while transport properties are obtained from kinetic principles, with the method of Chapman and Enskog. Particle ambipolar diffusive fluxes are found by solving the Stefan-Maxwell equations with a simple iterative method. Two physico-mathematical formulations are used to model the chemical reaction processes: (1) A Local Thermodynamics Equilibrium (LTE) formulation and (2) a thermo-chemical non-equilibrium (TCNEQ) formulation. In the TCNEQ model, thermal non-equilibrium between the translational energy mode of the gas and the vibrational energy mode of individual molecules is accounted for. The electronic states of the chemical species are assumed in equilibrium with the vibrational temperature, whereas the rotational energy mode is assumed to be equilibrated with translation. Three different physical models are used to account for the coupling of chemistry and energy transfer processes. Numerical simulations obtained with the LTE and TCNEQ formulations are used to characterize the extent of non-equilibrium of the flow inside the Plasmatron facility at the von Karman Institute. Each model was tested using different kinetic mechanisms to assess the sensitivity of the results to variations in the reaction parameters. A comparison of temperatures and composition profiles at the outlet of the torch demonstrates that the flow is in non-equilibrium for operating conditions characterized by pressures below 30 000 Pa, frequency 0.37 MHz, input power 80 kW, and mass flow 8 g/s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Refine your search to explore more results.