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,495 records · Page 83

Alternative optical concept for electron cyclotron emission imaging

The implementation of advanced electron cyclotron emission imaging (ECEI) systems on tokamak experiments has revolutionized the diagnosis of magnetohydrodynamic (MHD) activities and improved our understanding of instabilities, which lead to disruptions. It is therefore desirable to have an ECEI system on the ITER tokamak. However, the large size of optical components in presently used ECEI systems have, up to now, precluded the implementation of an ECEI system on ITER. This paper describes a new optical ECEI concept that employs a single spherical mirror as the only optical component and exploits the astigmatism of such a mirror to produce an image with one-dimensional spatial resolution on the detector. Since this alternative approach would only require a thin slit as the viewing port to the plasma, it would make the implementation of an ECEI system on ITER feasible. The results obtained from proof-of-principle experiments with a 125 GHz microwave system are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Electron kinetic effects on interferometry, polarimetry and Thomson scattering measurements in burning plasmas (invited)

At anticipated high electron temperatures in ITER, the effects of electron thermal motion on Thomson scattering (TS), toroidal interferometer/polarimeter (TIP), and poloidal polarimeter (PoPola) diagnostics will be significant and must be accurately treated. The precision of the previous lowest order linear in τ = T{sub e}/m{sub e}c{sup 2} model may be insufficient; we present a more precise model with τ{sup 2}-order corrections to satisfy the high accuracy required for ITER TIP and PoPola diagnostics. The linear model is extended from Maxwellian to a more general class of anisotropic electron distributions that allows us to take into account distortions caused by equilibrium current, ECRH, and RF current drive effects. The classical problem of the degree of polarization of incoherent Thomson scattered radiation is solved analytically exactly without any approximations for the full range of incident polarizations, scattering angles, and electron thermal motion from non-relativistic to ultra-relativistic. The results are discussed in the context of the possible use of the polarization properties of Thomson scattered light as a method of T{sub e} measurement relevant to ITER operational scenarios.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Diagnostics and required R and D for control of DEMO grade plasmas

Even if the diagnostics of ITER performs as expected, installation and operation of the diagnostic systems in Demo device will be much harsher than those of the present ITER device. In order to operate the Demo grade plasmas, which may have a higher beta limit, safely with very limited number of simple diagnostic system, it requires a well defined predictable plasma modelling in conjunction with the reliable control system for burn control and potential harmful instabilities. Development of such modelling in ITER is too risky and the logical choice would be utilization of the present day steady state capable devices such as KSTAR and EAST. In order to fulfill this mission, sophisticated diagnostic systems such as 2D/3D imaging systems can validate the physics in the theoretical modeling and challenge the predictable capability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron kinetic effects on optical diagnostics in fusion plasmas

At anticipated high electron temperatures in ITER, the effects of electron thermal motion on Thomson scattering (TS), toroidal interferometer/polarimeter (TIP) and poloidal polarimeter (PoPola) diagnostics will be significant and must be accurately treated. We calculate electron thermal corrections to the interferometric phase and polarization state of an EM wave propagating along tangential and poloidal chords (Faraday and Cotton-Mouton polarimetry) and perform analysis of the degree of polarization for incoherent TS. The precision of the previous lowest order linear in τ = T{sub e}/m{sub e}c{sup 2} model may be insufficient; we present a more precise model with τ{sup 2}-order corrections to satisfy the high accuracy required for ITER TIP and PoPola diagnostics. The linear model is extended from Maxwellian to a more general class of anisotropic electron distributions that allows us to take into account distortions caused by equilibrium current, ECRH and RF current drive effects. The classical problem of degree of polarization of incoherent Thomson scattered radiation is solved analytically exactly without any approximations for the full range of incident polarizations, scattering angles, and electron thermal motion from non-relativistic to ultra-relativistic. The results are discussed in the context of the possible use of the polarization properties of Thomson scattered light as a method of T{sup e} measurement relevant to ITER operational scenarios.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Can we estimate plasma density in ICP driver through electrical parameters in RF circuit?

To avoid regular maintenance, invasive plasma diagnostics with probes are not included in the inductively coupled plasma (ICP) based ITER Neutral Beam (NB) source design. Even non-invasive probes like optical emission spectroscopic diagnostics are also not included in the present ITER NB design due to overall system design and interface issues. As a result, negative ion beam current through the extraction system in the ITER NB negative ion source is the only measurement which indicates plasma condition inside the ion source. However, beam current not only depends on the plasma condition near the extraction region but also on the perveance condition of the ion extractor system and negative ion stripping. Nevertheless, inductively coupled plasma production region (RF driver region) is placed at distance (∼ 30cm) from the extraction region. Due to that, some uncertainties are expected to be involved if one tries to link beam current with plasma properties inside the RF driver. Plasma characterization in source RF driver region is utmost necessary to maintain the optimum condition for source operation. In this paper, a method of plasma density estimation is described, based on density dependent plasma load calculation.

43 PARTICLE ACCELERATORS↗

Status of the ELISE test facility

The test facility ELISE, equipped with a large radio frequency (RF) driven ion source (1×0.9 m2) of half the size of the ion source for the ITER neutral beam injection (NBI) system, is operational since beginning of 2013. The first experimental campaign was dedicated to a thorough qualification of the test facility and its diagnostic tools at low RF power (80 kW in total, i.e. 20 kW per driver) in volume operation, i.e. operation without cesium, where the negative hydrogen ion production is done in the plasma volume only. This paper reports on the main results of the second and third experimental campaigns, where Cs was inserted in the ion source for an enhancement of the negative ion production by the surface process. The second experimental campaign was done still with low RF power, both for hydrogen and deuterium, with pulse lengths of up to 500 s. The results of this campaign are rather encouraging, especially in hydrogen, where large current densities with respect to the low RF power could be achieved at a ratio of co-extracted electrons to extracted ions of 0.5-0.6 at the relevant source pressure of 0.3 Pa. Similar large extracted ion currents could be achieved also in deuterium, but with larger amounts of co-extracted electrons. The required ratio of co-extracted electrons to extracted ions of one could be achieved only in short pulses. The third experimental campaign aimed then for approaching the required ITER NBI parameters with respect to the ion and electron extracted currents, both for hydrogen and deuterium, by increasing the RF power with short pulses, i.e. beam-on times of up to 10 s and RF-on time up to 20 s. Current densities near the ITER NBI requirements could be achieved in hydrogen at a ratio of co-extracted electrons to extracted ions of 0.5-0.6 at the relevant source pressure of 0.3 Pa. As it was the case for the low RF operation, the required filter field was significantly lower than expected from the experience with the small prototype RF source. Similar large extracted ion currents could be achieved also in deuterium, but with larger amounts of co-extracted electrons; the main problem in deuterium operation are the non-stable currents of the co-extracted electrons, most probably caused by the high dynamic of the Cs redistribution in the system. The reasons for this larger instability are still under investigation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fusion power production in International Thermonuclear Experimental Reactor baseline H-mode scenarios

Self-consistent simulations of 15 MA ITER H-mode DT scenarios, from ramp-up through flat-top, are carried out. Electron and ion temperatures, toroidal angular frequency, and currents are evolved, in simulations carried out using the predictive TRANSPort and integrated modeling code starting with initial profiles and equilibria obtained from tokamak simulation code studies. Studies are carried out examining the dependence and sensitivity of fusion power production on electron density, argon impurity concentration, choice of radio frequency heating, pedestal temperature without and with E × B flow shear effects included, and the degree of plasma rotation. The goal of these whole-device ITER simulations is to identify dependencies that might impact ITER fusion performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of long pulse RF heating and current drive for H-mode scenarios with metallic walls in WEST

The longstanding expertise of the Tore Supra team in long pulse heating and current drive with radiofrequency (RF) systems will now be exploited in the WEST device (tungsten-W Environment in Steady-state Tokamak) [1]. WEST will allow an integrated long pulse tokamak programme for testing W-divertor components at ITER-relevant heat flux (10-20 MW/m{sup 2}), while treating crucial aspects for ITER-operation, such as avoidance of W-accumulation in long discharges, monitoring and control of heat fluxes on the metallic plasma facing components (PFCs) and coupling of RF waves in H-mode plasmas. Scenario modelling using the METIS-code shows that ITER-relevant heat fluxes are compatible with the sustainment of long pulse H-mode discharges, at high power (up to 15 MW / 30 s at I{sub P} = 0.8 MA) or high fluence (up to 10 MW / 1000 s at I{sub P} = 0.6 MA) [2], all based on RF heating and current drive using Ion Cyclotron Resonance Heating (ICRH) and Lower Hybrid Current Drive (LHCD). This paper gives a description of the ICRH and LHCD systems in WEST, together with the modelling of the power deposition of the RF waves in the WEST-scenarios.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dependence of intrinsic torque and momentum confinement on normalized gyroradius and collisionality in the DIII-D tokamak

We measured the dependence of intrinsic torque and momentum confinement time on normalized gyroradius (ρ * ) and collisionality (v * ) in the DIII-D tokamak. The intrinsic torque normalized to temperature is found to have ρ * and v * dependencies of ρ * -1.5 ± 0.8 and v * -0.26 ± 0.04 . This dependence on ρ * is unexpectedly favorable (increasing as ρ * decreases). The choice of normalization is important, and the implications are discussed. The unexpected dependence on ρ * is found to be robust, despite some uncertainty in the choice of normalization. Furthermore, the dependence of momentum confinement on ρ * does not clearly demonstrate Bohm or gyro-Bohm like scaling, and a weaker dependence on v * is found. The calculations required to use these dependencies to determine the intrinsic torque in future tokamaks such as ITER are presented, and the importance of the normalization is explained. Based on the currently available information, the intrinsic torque predicted for ITER is 33 N m, comparable to the expected torque available from neutral beam injection. The expected average intrinsic rotation associated with this intrinsic torque is small compared to current tokamaks, but it may still aid stability and performance in ITER. Published by AIP Publishing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Benchmarking the completely renormalised equation-of-motion coupled-cluster approaches for vertical excitation energies

When vertical energies are excited for a comprehensive test set of about 150 singlet excited states of 28 medium-sized organic molecules computed using two variants of the completely renormalised (CR) equation-of-motion (EOM) coupled-cluster (CC) method with singles, doubles, and non-iterative triples, abbreviated as δ-CR-EOMCCSD(T), and the analogous two variants of the newer, left-eigenstate δ-CR-EOMCC(2,3) approach are benchmarked against the previously published CASPT2, CC3, and EOMCCSDT-3 results, as well as the suggested theoretical best estimate (TBE) values. The δ-CR-EOMCC approaches are also used to identify and characterise about 50 additional excited states, including several states having substantial two-electron excitation components, which have not been found in the previous work and which can be used in future benchmark studies. We demonstrated that the non-iterative triples corrections to the EOMCCSD excitation energies defining the relatively inexpensive, single-reference, black-box δ-CR-EOMCC approaches provide significant improvements in the EOMCCSD data, while closely matching the results of the iterative and considerably more expensive CC3 and EOMCCSDT-3 calculations and their CASPT2 and TBE counterparts. It is also shown that the δ-CR-EOMCC methods, especially δ-CR-EOMCC(2,3), are capable of bringing the results of the CC3 and EOMCCSDT-3 calculations to a closer agreement with the CASPT2 and TBE data, demonstrating the utility of the cost-effective δ-CR-EOMCC methods in applications involving molecular electronic spectra. Finally, we show that there may exist a relationship between the magnitude of the triples corrections defining δ-CR-EOMCC approaches and the reduced excitation level diagnostic resulting from EOMCCSD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ON THE GRAVITATIONAL FIELDS OF MACLAURIN SPHEROID MODELS OF ROTATING FLUID PLANETS

Hubbard recently derived an important iterative equation for calculating the gravitational coefficients of a Maclaurin spheroid that does not require an expansion in a small distortion parameter. We show that this iterative equation, which is based on an incomplete solution of the Poisson equation, diverges when the distortion parameter is not sufficiently small. We derive a new iterative equation that is based on a complete solution of the Poisson equation and, hence, always converges when calculating the gravitational coefficients of a Maclaurin spheroid.

79 ASTRONOMY AND ASTROPHYSICS↗

Full waveform inversion for time-distance helioseismology

Inferring interior properties of the Sun from photospheric measurements of the seismic wavefield constitutes the helioseismic inverse problem. Deviations in seismic measurements (such as wave travel times) from their fiducial values estimated for a given model of the solar interior imply that the model is inaccurate. Contemporary inversions in local helioseismology assume that properties of the solar interior are linearly related to measured travel-time deviations. It is widely known, however, that this assumption is invalid for sunspots and active regions and is likely for supergranular flows. Here, we introduce nonlinear optimization, executed iteratively, as a means of inverting for the subsurface structure of large-amplitude perturbations. Defining the penalty functional as the L {sub 2} norm of wave travel-time deviations, we compute the total misfit gradient of this functional with respect to the relevant model parameters at each iteration around the corresponding model. The model is successively improved using either steepest descent, conjugate gradient, or the quasi-Newton limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm. Performing nonlinear iterations requires privileging pixels (such as those in the near field of the scatterer), a practice that is not compliant with the standard assumption of translational invariance. Measurements for these inversions, although similar in principle to those used in time-distance helioseismology, require some retooling. For the sake of simplicity in illustrating the method, we consider a two-dimensional inverse problem with only a sound-speed perturbation.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code

The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Upgrades of Diagnostic Techniques and Technologies for JET next D-T Campaigns

JET next D-T campaign is presently scheduled for the year 2017. The main scientific objectives include the assessment of the isotopic effects on various plasma aspects: mainly on confinement, on the threshold to access the H mode and on ELM behaviour. From a technical point of view, the total yield of the entire D-T phase is expected to be 1.7 1021 neutrons, about a factor of six higher than the previous main D-T campaign on JET, DTE1. Therefore the radiation field will be quite relevant for next step devices, since the neutron flux at the first wall (∼10{sup 16} n/cm{sup 2}s), for example, will be comparable to the one in ITER behind the blanket. From a point of view of diagnostics developments, for many years JET diagnostics have been upgraded in order to provide adequate support for the scientific exploitation of a D-T campaign. The main efforts have concentrated on improving three main aspects of JET measuring capability: 1) the quality of the measurements of the electron and ion fluids to support the plasma physics programme 2) the diagnostic for the fusion products 3) diagnostic technologies for ITER. In terms of general diagnostic capability, compared to the previous DTE1, JET diagnostics have a much better spatial and temporal resolution of both the ion and electron fluid (about one order of magnitude improvement for each parameter). The consistency of the various independent measurements of the same parameters has also increased significantly; the three independent measurements of the electron temperature, for example, agree now within 5%. Moreover, solutions are being addressed to operate some cameras, both visible and IR, even during the full D-T phase to provide imaging of the plasma and the first wall. Various upgrades of neutral particle analysis are being considered, mainly to measure the isotopic composition. A new set of reflectometers is expected to provide valuable information about the changes in the turbulence with the different fuel mixtures. With regard to the fusion products, JET now can deploy a consistent set of techniques to measure the neutron yield and neutron spectra and to diagnose the fast particles. A full calibration of the neutron diagnostics with a 14 MeV source is being considered, after the recent very successful calibration for the 2.45 MeV neutrons. Vertical and horizontal lines of sight are foreseen for neutron and gamma spectrometry, in order to better determine the thermal neutron yield and to separate the trapped and passing components of the alphas. Various gamma ray spectrometers are being developed to cover all the various operational scenarios, from trace tritium to 50-50 D-T operation. The redistribution of the alphas will be measured with the gamma ray cameras, recently upgraded with full digital electronics; new detectors are being considered to bring the time resolution of the system in the ten of ms range. The lost alphas will also be diagnosed with improved spatial and temporal resolution, using Faraday cups and a scintillator probe. From a technological perspective, the D-T campaign will provide a unique opportunity to test ITER relevant technologies. From radiation hard detectors, for example Hall probes, to neutron absorbers and to shielding concepts, the potential of various solutions in a realistic 14 MeV radiation field will be assessed. The effects of neutrons and gamma on ancillary technologies and systems, such as fibre optics and electronics circuits, are also expected to be sufficiently high to derive useful information about the competitive advantage of various alternatives.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A methodology for accident analysis of fusion breeder blankets and its application to helium-cooled lead–lithium blanket

'Fusion for Energy' (F4E) is designing, developing, and implementing the European Helium-Cooled Lead-Lithium (HCLL) and Helium-Cooled Pebble-Bed (HCPB) Test Blanket Systems (TBSs) for ITER (Nuclear Facility INB-174). Safety demonstration is an essential element for the integration of these TBSs into ITER and accident analysis is one of its critical components. A systematic approach to accident analysis has been developed under the F4E contract on TBS safety analyses. F4E technical requirements, together with Amec Foster Wheeler and INL efforts, have resulted in a comprehensive methodology for fusion breeding blanket accident analysis that addresses the specificity of the breeding blanket designs, materials, and phenomena while remaining consistent with the approach already applied to ITER accident analyses. Furthermore, the methodology phases are illustrated in the paper by its application to the EU HCLL TBS using both MELCOR and RELAP5 codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diagnostics for Burning Plasmas

Here, recent progress in the design and manufacturing development of in-vessel diagnostic components by the U.S. ITER team is reported. Design and manufacturing choices are driven by engineering challenges stemming from ITER’s burning plasma mission, which entails levels of plasma neutron emission, pulse lengths, and size scales well beyond the existing database. Challenges include port access limitations, neutron irradiation and activation risks, nuclear heating of components, complex interfaces, remote handling compatibility, and constraints on component and processing materials. Responses to these challenges are presented using examples from the low field side reflectometer, toroidal interferometer polarimeter, and electron cyclotron emission diagnostics. Relevance of the experience from ITER diagnostics engineering to the design of future integrated fusion facilities based on burning plasmas is addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Network-Level Optimization for Unbalanced Power Distribution System: Approximation and Relaxation

The nonlinear programming (NLP) problem to solve distribution-level optimal power flow (D-OPF) poses convergence issues and does not scale well for unbalanced distribution systems. The existing scalable D-OPF algorithms either use approximations that are not valid for an unbalanced power distribution system, or apply relaxation techniques to the nonlinear power flow equations that do not guarantee a feasible power flow solution. In this paper, we propose scalable D-OPF algorithms that simultaneously achieve optimal and feasible solutions by solving multiple iterations of approximate, or relaxed, D-OPF subproblems of low complexity. The first algorithm is based on a successive linear approximation of the nonlinear power flow equations around the current operating point, where the D-OPF solution is obtained by solving multiple iterations of a linear programming (LP) problem. The second algorithm is based on the relaxation of the nonlinear power flow equations as conic constraints together with directional constraints, which achieves optimal and feasible solutions over multiple iterations of a second-order cone programming (SOCP) problem. Finally, it is demonstrated that the proposed algorithms are able to reach an optimal and feasible solution while significantly reducing the computation time as compared to an equivalent NLPD-OPF model for the same distribution system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Temporal resolution and motion artifacts in single-source and dual-source cardiac CT

Purpose: The temporal resolution of a given image in cardiac computed tomography (CT) has so far mostly been determined from the amount of CT data employed for the reconstruction of that image. The purpose of this paper is to examine the applicability of such measures to the newly introduced modality of dual-source CT as well as to methods aiming to provide improved temporal resolution by means of an advanced image reconstruction algorithm. Methods: To provide a solid base for the examinations described in this paper, an extensive review of temporal resolution in conventional single-source CT is given first. Two different measures for assessing temporal resolution with respect to the amount of data involved are introduced, namely, either taking the full width at half maximum of the respective data weighting function (FWHM-TR) or the total width of the weighting function (total TR) as a base of the assessment. Image reconstruction using both a direct fan-beam filtered backprojection with Parker weighting as well as using a parallel-beam rebinning step are considered. The theory of assessing temporal resolution by means of the data involved is then extended to dual-source CT. Finally, three different advanced iterative reconstruction methods that all use the same input data are compared with respect to the resulting motion artifact level. For brevity and simplicity, the examinations are limited to two-dimensional data acquisition and reconstruction. However, all results and conclusions presented in this paper are also directly applicable to both circular and helical cone-beam CT. Results: While the concept of total TR can directly be applied to dual-source CT, the definition of the FWHM of a weighting function needs to be slightly extended to be applicable to this modality. The three different advanced iterative reconstruction methods examined in this paper result in significantly different images with respect to their motion artifact level, despite exactly the same amount of data being used in the reconstruction process. Conclusions: The concept of assessing temporal resolution by means of the data employed for reconstruction can nicely be extended from single-source to dual-source CT. However, for advanced (possibly nonlinear iterative) reconstruction algorithms the examined approach fails to deliver accurate results. New methods and measures to assess the temporal resolution of CT images need to be developed to be able to accurately compare the performance of such algorithms.

62 RADIOLOGY AND NUCLEAR MEDICINE↗