Search NASA⌕ Search

SEARCH · Search NASA

Results for “LOGARITHM”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

The velocity distribution of hypervelocity stars

We consider the process of stellar binaries tidally disrupted by a supermassive black hole (BH). For highly eccentric orbits, as one star is ejected from the three-body system, the companion remains bound to the BH. Hypervelocity stars (HVSs) observed in the Galactic halo and S-stars observed orbiting the central BH may originate from such mechanism. In this paper, we predict the velocity distribution of the ejected stars of a given mass, after they have traveled out of the Galactic potential. We use both analytical methods and Monte Carlo simulations. We find that each part of the velocity distribution encodes different information. At low velocities <800 km s{sup –1}, the Galactic potential universally shapes the observed distribution, which rises toward a peak, related to the Galactic escape velocity. Beyond the peak, the velocity distribution depends on binary mass and separation distributions. Finally, the finite star life introduces a break related to their mass. A qualitative comparison of our models with current observations shows the great potential of HVSs to constrain bulge and Galactic properties. Standard choices for parameter distributions predict velocities below and above ∼800 km s{sup –1} with equal probability, while none are observed beyond ∼700 km s{sup –1} and the current detections are more clustered at low velocities 300-400 km s{sup –1}. These features may indicate that the separation distribution of binaries that reach the tidal sphere is not flat in logarithmic space, as observed in more local massive binaries, but has more power toward larger separations, enhancing smaller velocities. In addition, the binary formation/evolution process or the injection mechanism might also induce a cut-off a {sub min} ∼ 10 R {sub ☉} in the separation distribution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Hubble space telescope emission line galaxies at z ∼ 2: the Lyα escape fraction

We compare the Hβ line strengths of 1.90 < z < 2.35 star-forming galaxies observed with the near-IR grism of the Hubble Space Telescope with ground-based measurements of Lyα from the HETDEX Pilot Survey and narrow-band imaging. By examining the line ratios of 73 galaxies, we show that most star-forming systems at this epoch have a Lyα escape fraction below ∼6%. We confirm this result by using stellar reddening to estimate the effective logarithmic extinction of the Hβ emission line (c {sub Hβ} = 0.5) and measuring both the Hβ and Lyα luminosity functions in a ∼100, 000 Mpc{sup 3} volume of space. We show that in our redshift window, the volumetric Lyα escape fraction is at most 4.4{sub −1.2}{sup +2.1}%, with an additional systematic ∼25% uncertainty associated with our estimate of extinction. Finally, we demonstrate that the bulk of the epoch's star-forming galaxies have Lyα emission line optical depths that are significantly greater than that for the underlying UV continuum. In our predominantly [O III] λ5007-selected sample of galaxies, resonant scattering must be important for the escape of Lyα photons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A Fourier transform method for Vsin i estimations under nonlinear Limb-Darkening laws

Star rotation offers us a large horizon for the study of many important physical issues pertaining to stellar evolution. Currently, four methods are widely used to infer rotation velocities, namely those related to line width calibrations, on the fitting of synthetic spectra, interferometry, and on Fourier transforms (FTs) of line profiles. Almost all of the estimations of stellar projected rotation velocities using the Fourier method in the literature have been addressed with the use of linear limb-darkening (LD) approximations during the evaluation of rotation profiles and their cosine FTs, which in certain cases, lead to discrepant velocity estimates. In this work, we introduce new mathematical expressions of rotation profiles and their Fourier cosine transforms assuming three nonlinear LD laws—quadratic, square-root, and logarithmic—and study their applications with and without gravity-darkening (GD) and geometrical flattening (GF) effects. Through an analysis of He I models in the visible range accounting for both limb and GD, we find out that, for classical models without rotationally driven effects, all the Vsin i values are too close to each other. On the other hand, taking into account GD and GF, the Vsin i obtained with the linear law result in Vsin i values that are systematically smaller than those obtained with the other laws. As a rule of thumb, we apply these expressions to the FT method to evaluate the projected rotation velocity of the emission B-type star Achernar (α Eri).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE HERSCHEL EXPLOITATION OF LOCAL GALAXY ANDROMEDA (HELGA). VI. THE DISTRIBUTION AND PROPERTIES OF MOLECULAR CLOUD ASSOCIATIONS IN M31

In this paper we present a catalog of giant molecular clouds (GMCs) in the Andromeda (M31) galaxy extracted from the Herschel Exploitation of Local Galaxy Andromeda (HELGA) data set. GMCs are identified from the Herschel maps using a hierarchical source extraction algorithm. We present the results of this new catalog and characterize the spatial distribution and spectral energy properties of its clouds based on the radial dust/gas properties found by Smith et al. A total of 326 GMCs in the mass range 10{sup 4}-10{sup 7} M {sub ☉} are identified; their cumulative mass distribution is found to be proportional to M {sup –2.34}, in agreement with earlier studies. The GMCs appear to follow the same correlation of cloud mass to L {sub CO} observed in the Milky Way. However, comparison between this catalog and interferometry studies also shows that the GMCs are substructured below the Herschel resolution limit, suggesting that we are observing associations of GMCs. Following Gordon et al., we study the spatial structure of M31 by splitting the observed structure into a set of spiral arms and offset rings. We fit radii of 10.3 and 15.5 kpc to the two most prominent rings. We then fit a logarithmic spiral with a pitch angle of 8.°9 to the GMCs not associated with either ring. Last, we comment on the effects of deprojection on our results and investigate the effect different models for M31's inclination will have on the projection of an unperturbed spiral arm system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE STRUCTURE OF HE 1104-1805 FROM INFRARED TO X-RAY

The gravitationally lensed quasar HE 1104-1805 has been observed at a variety of wavelengths ranging from the mid-infrared to X-ray for nearly 20 yr. We combine flux ratios from the literature, including recent Chandra data, with new observations from the SMARTS telescope and Hubble Space Telescope, and use them to investigate the spatial structure of the central regions using a Bayesian Monte Carlo analysis of the microlensing variability. The wide wavelength coverage allows us to constrain not only the accretion disk half-light radius r {sub 1/2}, but the power-law slope ξ of the size-wavelength relation r {sub 1/2}∝λ{sup ξ}. With a logarithmic prior on the source size, the (observed-frame) R-band half-light radius log (r {sub 1/2}/cm) is 16.0{sub {sup −}{sup 0}{sup .}{sup 4{sub +}{sub 0}{sub .}{sub 3}}}, and the slope ξ is 1.0{sub {sup −}{sup 0}{sup .}{sup 5}{sup 6{sub +}{sub 0}{sub .}{sub 3}{sub 0}}}. We put upper limits on the source size in soft (0.4–1.2 keV) and hard (1.2–8 keV) X-ray bands, finding 95% upper limits on log (r {sub 1/2}/cm) of 15.33 in both bands. A linear prior yields somewhat larger sizes, particularly in the X-ray bands. For comparison, the gravitational radius, using a black hole mass estimated using the Hβ line, is log (r{sub g} /cm) = 13.94. We find that the accretion disk is probably close to face-on, with cos i = 1.0 being four times more likely than cos i = 0.5. We also find probability distributions for the mean mass of the stars in the foreground lensing galaxy, the direction of the transverse peculiar velocity of the lens, and the position angle of the projected accretion disk's major axis (if not face-on)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A STATISTICAL RECONSTRUCTION OF THE PLANET POPULATION AROUND KEPLER SOLAR-TYPE STARS

Using the cumulative catalog of planets detected by the NASA Kepler mission, we reconstruct the intrinsic occurrence of Earth- to Neptune-size (1-4 R {sub ⊕}) planets and their distributions with radius and orbital period. We analyze 76,711 solar-type (0.8 < R {sub *}/R {sub ☉} < 1.2) stars with 430 planets on 20-200 day orbits, excluding close-in planets that may have been affected by the proximity to the host star. Our analysis considers errors in planet radii and includes an ''iterative simulation'' technique that does not bin the data. We find a radius distribution that peaks at 2-2.8 Earth radii, with lower numbers of smaller and larger planets. These planets are uniformly distributed with logarithmic period, and the mean number of such planets per star is 0.46 ± 0.03. The occurrence is ∼0.66 if planets interior to 20 days are included. We estimate the occurrence of Earth-size planets in the ''habitable zone'' (defined as 1-2 R {sub ⊕}, 0.99-1.7 AU for solar-twin stars) as 6.4{sub −1.1}{sup +3.4}%. Our results largely agree with those of Petigura et al., although we find a higher occurrence of 2.8-4 Earth-radii planets. The reasons for this excess are the inclusion of errors in planet radius, updated Huber et al. stellar parameters, and also the exclusion of planets that may have been affected by proximity to the host star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Does the presence of planets affect the frequency and properties of extrasolar Kuiper belts? Results from the HErschel DEBRIS and DUNES surveys

The study of the planet–debris disk connection can shed light on the formation and evolution of planetary systems and may help “predict” the presence of planets around stars with certain disk characteristics. In preliminary analyses of subsamples of the Herschel DEBRIS and DUNES surveys, Wyatt et al. and Marshall et al. identified a tentative correlation between debris and the presence of low-mass planets. Here we use the cleanest possible sample out of these Herschel surveys to assess the presence of such a correlation, discarding stars without known ages, with ages <1 Gyr, and with binary companions <100 AU to rule out possible correlations due to effects other than planet presence. In our resulting subsample of 204 FGK stars, we do not find evidence that debris disks are more common or more dusty around stars harboring high-mass or low-mass planets compared to a control sample without identified planets. There is no evidence either that the characteristic dust temperature of the debris disks around planet-bearing stars is any different from that in debris disks without identified planets, nor that debris disks are more or less common (or more or less dusty) around stars harboring multiple planets compared to single-planet systems. Diverse dynamical histories may account for the lack of correlations. The data show a correlation between the presence of high-mass planets and stellar metallicity, but no correlation between the presence of low-mass planets or debris and stellar metallicity. Comparing the observed cumulative distribution of fractional luminosity to those expected from a Gaussian distribution in logarithmic scale, we find that a distribution centered on the solar system’s value fits the data well, while one centered at 10 times this value can be rejected. This is of interest in the context of future terrestrial planet detection and characterization because it indicates that there are good prospects for finding a large number of debris disk systems (i.e., with evidence of harboring planetesimals, the building blocks of planets) with exozodiacal emission low enough to be appropriate targets for an ATLAST-type mission to search for biosignatures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

RESIK solar X-ray flare element abundances on a non-isothermal assumption

Solar X-ray spectra from the REntgenovsky Spektrometr s Izognutymi Kristalami (RESIK) crystal spectrometer on the CORONAS-F spacecraft (spectral range 3.3–6.1 Å) are analyzed for 33 flares using a method to derive abundances of Si, S, Ar, and K, emission lines of which feature prominently in the spectra. For each spectrum, the method first optimizes element abundances and then derives the differential emission measure as a function of temperature based on a procedure given by Sylwester et al. and Withbroe. This contrasts with our previous analyses of RESIK spectra in which an isothermal assumption was used. The revised abundances (on a logarithmic scale with A(H)=12) averaged for all the flares in the analysis are A(Si)=7.53±0.08 (previously 7.89 ± 0.13), A(S)=6.91±0.07 (7.16 ± 0.17), A(Ar)=6.47±0.08 (6.45 ± 0.07), and A(K)=5.73±0.19 (5.86 ± 0.20), with little evidence for time variations of abundances within the evolution of each flare. Our previous estimates of the Ar and K flare abundances are thus confirmed by this analysis, but those for Si and S are reduced. This suggests that the flare abundances of Si and Ar are very close to the photospheric abundance or solar proxies, while S is significantly less than photospheric and the K abundance is much higher than photospheric. These estimates differ to some extent from those in which a single enhancement factor applies to elements with first ionization potential less than 10 eV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A new model for the radio emission from SN 1994I and an associated search for radio transients in M51

We revisit the exquisite archival radio data for the Type Ic supernova SN 1994I and present a revised model for the supernova (SN) radio emission and a pilot study that aims to constrain the rate of C-band radio transients within the face-on host galaxy, M51 (NGC 5194). We find that the temporal and spectral evolution of the SN 1994I radio emission are well fit by a synchrotron self-absorption model and use this to estimate physical parameters. We compute a pre-explosion mass loss rate of M-dot =3.0×10{sup −5} M{sub ⊙} yr{sup −1} for the progenitor, consistent with those observed from galactic Wolf–Rayet stars. Our model makes different assumptions for the dynamical model for the shockwave interaction than the model previously published by Weiler et al., but our M-dot is consistent with theirs to within errors and assumptions. Drawing from a subset of the archival radio observations from the Very Large Array collected for the monitoring of SN 1994I, we conduct a pilot study to search for previously unidentified transients. Data were primarily taken at a frequency of 4.9 GHz and are logarithmic in cadence, enabling sensitivity to transients with variability timescales ranging from days to months. We find no new transient detections in 31 epochs of data, allowing us to place a 2σ upper limit of 17 deg{sup −2} for the source density of radio transients above 0.5 mJy (L ≳ 4×10{sup 25} erg s{sup −1} Hz{sup −1} at the distance of M51). This study highlights the feasibility of utilizing archival high-cadence radio studies of SN host galaxies to place constraints on the radio transient rate as a function of luminosity in the local universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The energetics of cusp destruction

We present a new analytic estimate for the energy required to create a constant density core within a dark matter halo that, based on more realistic assumptions, leads to demands that are orders of magnitude lower than claimed in earlier works. We define a core size based on the logarithmic slope of the dark matter density profile as it is insensitive to the functional form used to fit observed data. The energy required to form a core sensitively depends on the radial scale over which dark matter within the cusp is redistributed within the halo. Simulations indicate that within a region in size comparable to the active star forming regions of the central galaxy inhabiting a halo, dark matter particles have their orbits radially increased by a factor of 2–3 during core formation. Thus, the inner properties of the dark matter halo set the energy requirements. The energy cost increases slowly with halo mass as M{sub h} {sup 0.3−0.7} for core sizes ≲1 kpc. We use the expected star formation history for a given halo mass to predict dwarf galaxy core sizes. We find that supernovae alone would create well over 4 kpc cores in 10{sup 10} M{sub ⊙} galaxies if 100% of the energy were transferred to dark matter particle orbits. We can directly constrain the efficiency factor by studying galaxies with known stellar content and core size. We find that the efficiency of coupling between stellar feedback and dark matter orbital energy need only be ≲1% to explain Fornax’s 1 kpc core.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Physical properties of tidal features of interacting disk galaxies: Three-dimensional self-consistent models

Using self-consistent three-dimensional (3D) N-body simulations, we investigate the physical properties of nonaxisymmetric features in a disk galaxy created by a tidal interaction with its companion. The primary galaxy consists of a stellar disk, a bulge, and a live halo, corresponding to Milky-Way-type galaxies, while the companion is represented by a halo alone. We vary the companion mass and the pericenter distance to explore situations with differing tidal strength parameterized by either the relative tidal force P or the relative imparted momentum S. We find that the formation of a tidal tail in the outer parts requires P≳0.05 or S≳0.07. A stronger interaction results in a stronger, less wound tail that forms earlier. Similarly, a stronger tidal forcing produces stronger, more loosely wound spiral arms in the inner parts. The arms are approximately logarithmic in shape, with both amplitude and pitch angle decaying with time. The derived pattern speed decreases with radius and is close to the Ω−κ/2 curve at late time, with Ω and κ denoting the angular and epicycle frequencies, respectively. This suggests that the tidally induced spiral arms are most likely kinematic density waves weakly modified by self-gravity. Compared to the razor-thin counterparts, arms in the 3D models are weaker, have a smaller pitch angle, and wind and decay more rapidly. The 3D density structure of the arms is well described by the concentrated and sinusoidal models when the arms are in the nonlinear and linear regimes, respectively. We demonstrate that dynamical friction between interacting galaxies transfers the orbital angular momentum of one galaxy to the spin angular momentum of the companion halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE MASS DISTRIBUTION FUNCTION OF PLANETS

The distribution of orbital period ratios of adjacent planets in extrasolar planetary systems discovered by the Kepler space telescope exhibits a peak near ∼1.5–2, a long tail of larger period ratios, and a steep drop-off in the number of systems with period ratios below ∼1.5. We find from these data that the dimensionless orbital separations have an approximately log-normal distribution. Using Hill’s criterion for the dynamical stability of two planets, we find an upper bound on planet masses such that the most common planet mass does not exceed 10{sup −3.2}m{sub ∗}, or about two-thirds of Jupiter’s mass for solar-mass stars. Assuming that the mass ratio and the dynamical separation (orbital spacings in units of mutual Hill radius) of adjacent planets are independent random variates, and adopting empirical distributions for these, we use Hill’s criterion in a statistical way to estimate the planet mass distribution function from the observed distribution of orbital separations. We find that the planet mass function is peaked in logarithm of mass, with a peak value and standard deviation of logm/M{sub ⊕} of ∼(0.6−1.0) and ∼(1.1−1.2), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Non-universality of dark-matter halos: Cusps, cores, and the central potential

Dark-matter halos grown in cosmological simulations appear to have central Navarro–Frenk–White-like density cusps with mean values of dlogρ/dlogr≈−1, and some dispersion, which is generally parametrized by the varying index α in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp—which has only one shape parameter, ϕ{sub 0}—fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp ϕ{sub 0}, Einasto α, or the central logarithmic slope in the Dehnen–Tremaine analytical γ-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

ADIABATIC MASS LOSS IN BINARY STARS. II. FROM ZERO-AGE MAIN SEQUENCE TO THE BASE OF THE GIANT BRANCH

In the limit of extremely rapid mass transfer, the response of a donor star in an interacting binary becomes asymptotically one of adiabatic expansion. We survey here adiabatic mass loss from Population I stars (Z = 0.02) of mass 0.10 M{sub ⊙}–100 M{sub ⊙} from the zero-age main sequence to the base of the giant branch, or to central hydrogen exhaustion for lower main sequence stars. The logarithmic derivatives of radius with respect to mass along adiabatic mass-loss sequences translate into critical mass ratios for runaway (dynamical timescale) mass transfer, evaluated here under the assumption of conservative mass transfer. For intermediate- and high-mass stars, dynamical mass transfer is preceded by an extended phase of thermal timescale mass transfer as the star is stripped of most of its envelope mass. The critical mass ratio q{sub ad} (throughout this paper, we follow the convention of defining the binary mass ratio as q ≡ M{sub donor}/M{sub accretor}) above which this delayed dynamical instability occurs increases with advancing evolutionary age of the donor star, by ever-increasing factors for more massive donors. Most intermediate- or high-mass binaries with nondegenerate accretors probably evolve into contact before manifesting this instability. As they approach the base of the giant branch, however, and begin developing a convective envelope, q{sub ad} plummets dramatically among intermediate-mass stars, to values of order unity, and a prompt dynamical instability occurs. Among low-mass stars, the prompt instability prevails throughout main sequence evolution, with q{sub ad} declining with decreasing mass, and asymptotically approaching q{sub ad} = 2/3, appropriate to a classical isentropic n = 3/2 polytrope. Our calculated q{sub ad} values agree well with the behavior of time-dependent models by Chen and Han of intermediate-mass stars initiating mass transfer in the Hertzsprung gap. Application of our results to cataclysmic variables, as systems that must be stable against rapid mass transfer, nicely circumscribes the range in q{sub ad} as a function of the orbital period in which they are found. These results are intended to advance the verisimilitude of population synthesis models of close binary evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

DECONSTRUCTING THERMAL SUNYAEV–ZEL’DOVICH—GRAVITATIONAL LENSING CROSS-CORRELATIONS: IMPLICATIONS FOR THE INTRACLUSTER MEDIUM

Recent first detections of the cross-correlation of the thermal Sunyaev–Zel’dovich (tSZ) signal in Planck cosmic microwave background (CMB) temperature maps with gravitational lensing maps inferred from the Planck CMB data and the CFHTLenS galaxy survey provide new probes of the relationship between baryons and dark matter. Using cosmological hydrodynamics simulations, we show that these cross-correlation signals are dominated by contributions from hot gas in the intracluster medium (ICM), rather than diffuse, unbound gas located beyond the virial radius (the “missing baryons”). Thus, these cross-correlations offer a tool with which to study the ICM over a wide range of halo masses and redshifts. In particular, we show that the tSZ—CMB lensing cross-correlation is more sensitive to gas in lower-mass, higher-redshift halos and gas at larger cluster-centric radii than the tSZ—galaxy lensing cross-correlation. Combining these measurements with primary CMB data will constrain feedback models through their signatures in the ICM pressure profile. We forecast the ability of ongoing and future experiments to constrain the parameters of a phenomenological ICM model, including the mean amplitude of the pressure–mass relation, the redshift evolution of this amplitude, and the mean outer logarithmic slope of the pressure profile. The results are promising, with ≈5%–20% precision constraints achievable with upcoming experiments, even after marginalizing over cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SHARDS: A global view of the star formation activity at z ∼ 0.84 and z ∼ 1.23

In this paper, we present a comprehensive analysis of star-forming galaxies (SFGs) at intermediate redshifts (z ∼ 1). We combine the ultra-deep optical spectro-photometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) with deep UV-to-FIR observations in the GOODS-N field. Exploiting two of the 25 SHARDS medium-band filters, F687W17 and F823W17, we select [O ii] emission line galaxies at z ∼ 0.84 and z ∼ 1.23 and characterize their physical properties. Their rest-frame equivalent widths (EW{sub rf}([O ii])), line fluxes, luminosities, star formation rates (SFRs), and dust attenuation properties are investigated. The evolution of EW{sub rf}([O ii]) closely follows the SFR density evolution of the universe, with a trend of EW{sub rf}([O ii]) ∝ (1 + z){sup 3} up to redshift z ≃ 1, followed by a possible flattening. The SF properties of the galaxies selected on the basis of their [O ii] emission are compared with complementary samples of SFGs selected by their MIR and FIR emission, and also with a general mass-selected sample of galaxies at the same redshifts. We demonstrate observationally that the UVJ diagram (or, similarly, a cut in the specific SFR) is only partially able to distinguish the quiescent galaxies from the SFGs. The SFR–M{sub *} relation is investigated for the different samples, yielding a logarithmic slope ∼1, in good agreement with previous results. The dust attenuations derived from different SFR indicators (UV(1600), UV(2800), [O ii], IR) are compared and show clear trends with respect to both the stellar mass and total SFR, with more massive and highly star-forming galaxies being affected by stronger dust attenuation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Monte Carlo uncertainty estimation for an oscillating-vessel viscosity measurement

This paper discusses the initial design and evaluation of a high temperature viscosity measurement system with the focus on the uncertainty assessment. Numerical simulation of the viscometer is used to estimate viscosity uncertainties through the Monte Carlo method. The simulation computes the system response for a particular set of inputs (viscosity, moment of inertia, spring constant and hysteretic damping), and the viscosity is calculated using two methods: the Roscoe approximate solution and a numerical-fit method. For numerical fitting, a residual function of the logarithmic decay of oscillation amplitude and oscillation period is developed to replace the residual function of angular oscillation, which is mathematically stiff. The results of this study indicate that the method using computational solution of the equations and fitting for the parameters should be used, since it almost always out-performs the Roscoe approximation in uncertainty. The hysteretic damping and spring stiffness uncertainties translate into viscosity uncertainties almost directly, whereas the moment of inertial and vessel-height uncertainties are magnified approximately two-fold. As the hysteretic damping increases, so does the magnification of its uncertainty, therefore it should be minimized in the system design. The result of this study provides a general guide for the design and application of all oscillation-vessel viscosity measurement systems.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2015

Impact of a narrow limiter SOL heat flux channel on the ITER first wall panel shaping

We report the inboard limiters for ITER were initially designed on the assumption that the parallel heat flux density in the scrape-off layer (SOL) could be approximated by a single exponential with decay length λ q . This assumption was found not to be adequate in 2012, when infra-red (IR) thermography measurements on the inner column during JET limiter discharges clearly revealed the presence of a narrow heat flux channel adjacent to the last closed flux surface. This near-SOL decay occurs with λ q ~ few mm, much shorter than the main SOL λ q , and can raise the heat flux at the limiter apex a factor up to ~4 above the value expected from a single, broader exponential. The original logarithmically shaped ITER inner wall first wall panels (FWPs) would be unsuited to handling the power loads produced by such a narrow feature. A multi-machine study involving the C-Mod, COMPASS, DIII-D and TCV tokamaks, employing inner wall IR measurements and/or inner wall reciprocating probes, was initiated to investigate the narrow limiter SOL heat flux channel. This paper describes the new results which have provided an experimental database for the narrow feature and presents an ITER inner wall FWP toroidal shape optimized for a double-exponential profile with λ q = 4 (narrow feature) and 50 mm (main-SOL), the latter also derived from a separate multi-machine database constituted recently within the International Tokamak Physics Activity. It is shown that the new shape allows the power handling capability of the original shape design to be completely recovered for a wide variety of limiter start-up equilibria in the presence of a narrow feature, even taking assembly tolerances into account. Additionally, it is further shown that the new shape has the interesting property of both mitigating the impact of the narrow feature and resulting in only a very modest increase in heat load, compared to the current design, if the narrow feature is not eventually found on ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗