Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed simulation”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Automated evaluation of setup errors in carbon ion therapy using PET: Feasibility study

Purpose: To investigate the possibility of detecting patient mispositioning in carbon-ion therapy with particle therapy positron emission tomography (PET) in an automated image registration based manner. Methods: Tumors in the head and neck (H and N), pelvic, lung, and brain region were investigated. Biologically optimized carbon ion treatment plans were created with TRiP98. From these treatment plans, the reference β{sup +}-activity distributions were calculated using a Monte Carlo simulation. Setup errors were simulated by shifting or rotating the computed tomography (CT). The expected β{sup +} activity was calculated for each plan with shifts. Finally, the reference particle therapy PET images were compared to the “shifted” β{sup +}-activity distribution simulations using the Pearson's correlation coefficient (PCC). To account for different PET monitoring options the inbeam PET was compared to three different inroom scenarios. Additionally, the dosimetric effects of the CT misalignments were investigated. Results: The automated PCC detection of patient mispositioning was possible in the investigated indications for cranio-caudal shifts of 4 mm and more, except for prostate tumors. In the rather homogeneous pelvic region, the generated β{sup +}-activity distribution of the reference and compared PET image were too much alike. Thus, setup errors in this region could not be detected. Regarding lung lesions the detection strongly depended on the exact tumor location: in the center of the lung tumor misalignments could be detected down to 2 mm shifts while resolving shifts of tumors close to the thoracic wall was more challenging. Rotational shifts in the H and N and lung region of +6° and more could be detected using inroom PET and partly using inbeam PET. Comparing inroom PET to inbeam PET no obvious trend was found. However, among the inroom scenarios a longer measurement time was found to be advantageous. Conclusions: This study scopes the use of various particle therapy PET verification techniques in four indications. The automated detection of patients' setup errors was investigated in a broad accumulation of data sets. The evaluation of introduced setup errors is performed automatically, which is of utmost importance to introduce highly required particle therapy monitoring devices into the clinical routine.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Improving the efficiency of configurational-bias Monte Carlo: A density-guided method for generating bending angle trials for linear and branched molecules

A new method has been developed to generate bending angle trials to improve the acceptance rate and the speed of configurational-bias Monte Carlo. Whereas traditionally the trial geometries are generated from a uniform distribution, in this method we attempt to use the exact probability density function so that each geometry generated is likely to be accepted. In actual practice, due to the complexity of this probability density function, a numerical representation of this distribution function would be required. This numerical table can be generated a priori from the distribution function. This method has been tested on a united-atom model of alkanes including propane, 2-methylpropane, and 2,2-dimethylpropane, that are good representatives of both linear and branched molecules. It has been shown from these test cases that reasonable approximations can be made especially for the highly branched molecules to reduce drastically the dimensionality and correspondingly the amount of the tabulated data that is needed to be stored. Despite these approximations, the dependencies between the various geometrical variables can be still well considered, as evident from a nearly perfect acceptance rate achieved. For all cases, the bending angles were shown to be sampled correctly by this method with an acceptance rate of at least 96% for 2,2-dimethylpropane to more than 99% for propane. Since only one trial is required to be generated for each bending angle (instead of thousands of trials required by the conventional algorithm), this method can dramatically reduce the simulation time. The profiling results of our Monte Carlo simulation code show that trial generation, which used to be the most time consuming process, is no longer the time dominating component of the simulation.

2-2-DIMETHYLPROPANE↗
DOE OSTI2014

A spectroscopic census in young stellar regions: the σ Orionis cluster

We present a spectroscopic survey of the stellar population of the σ Orionis cluster. We have obtained spectral types for 340 stars. Spectroscopic data for spectral typing come from several spectrographs with similar spectroscopic coverage and resolution. More than half of the stars in our sample are members confirmed by the presence of lithium in absorption, strong Hα in emission or weak gravity-sensitive features. In addition, we have obtained high-resolution (R ∼ 34,000) spectra in the Hα region for 169 stars in the region. Radial velocities were calculated from this data set. The radial velocity distribution for members of the cluster is in agreement with previous work. Analysis of the profile of the Hα line and infrared observations reveals two binary systems or fast rotators that mimic the Hα width expected in stars with accretion disks. On the other hand, there are stars with optically thick disks and narrow Hα profiles not expected in stars with accretion disks. This contribution constitutes the largest homogeneous spectroscopic data set of the σ Orionis cluster to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Ensuring Workforce Readiness for the Energy Industry Through 3D Simulations

In February 2009, during his joint address to Congress, President Obama challenged every American to commit to at least one year of college education and set a goal that the United States should once again have the highest proportion of college graduates in the world by the year 2020. This push by the Obama administration reflected the understanding that while a college degree was becoming a prerequisite to obtain many of the jobs in the 21. century economy; the United States still lagged behind other countries with respect to higher education attainment. Because of this, the U.S. is facing a skills gap between the credentials required for the jobs of today and tomorrow, and the number of skilled workers eligible to fill those jobs. Within the energy industry, the skills gap is compounded by two factors: an expanding domestic industry in energy generation and transmission, and an aging workforce set to retire in large numbers. In fact, a 2013 survey administered by Center for Energy Workforce Development (CEWD) estimates that 36% of the workforce may need to be replaced in the next five years, and nearly 50% will need to be replaced by the year 2023. As a 2014 report by CEWD notes, 'The message for workforce development is clear: qualified and skilled utility technicians and engineers will continue to be in high demand in all types of energy generation, transmission, and distribution jobs'. Simulations provide students the ability to practice proper procedure in a safe environment, test out multiple scenarios to determine an optimal outcome, and repeat practice to solidify learning. Industry currently applies simulation (i.e., 3D Electrical Substation Simulation) as it affords an opportunity for technicians to train in a safe, virtual environment, thus eliminating risks associated with arc flash hazards. The intent of 3D simulation is to enhance the students' understanding of the concepts learned in the course and better prepare students for field work in a related technical discipline. The 3D simulation will reinforce the concepts learned during a specific course and are intended to add real world context to the curriculum. The 3D simulation will be designed such that a student will perform specific job activities incorporating concepts learned during the course. This method of instruction ensures students have an understanding of the concepts and are able to apply the concepts in a real world context. Different scenarios will be developed to ensure a thorough understanding of the concepts with real time feedback when there is a misunderstanding of the concepts. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Virtual Time III, Part 2: Combining Conservative and Optimistic Synchronization

This is Part 2 of a trio of works intended to provide a unifying framework in which conservative and optimistic synchronization for parallel discrete event simulations can be freely and transparently combined in the same logical process on an event-by-event basis. Here, in this article, we continue the outline of an approach called Unified Virtual Time (UVT) that was introduced in Part 1, showing in detail via two extended examples how conservative synchronization can be refactored and combined with optimistic synchronization in the UVT framework. We describe UVT versions of both a basic time windowing algorithm called Unified Simple Time Windows and a refactored version of the Chandy-Misra-Bryant Null Message algorithm called Unified CMB.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

On lithium-rich red giants. I. Engulfment of substellar companions

A small fraction of red giants are known to be lithium (Li) rich, in contradiction with expectations from stellar evolutionary theory. A possible explanation for these atypical giants is the engulfment of an Li-rich planet or brown dwarf by the star. In this work, we model the evolution of Li abundance in canonical red giants including the accretion of a substellar mass companion. We consider a wide range of stellar and companion masses, Li abundances, stellar metallicities, and planetary orbital periods. Based on our calculations, companions with masses lower than 15 M{sub J} dissolve in the convective envelope and can induce Li enrichment in regimes where extra mixing does not operate. Our models indicate that the accretion of a substellar companion can explain abundances up to A(Li) ≈ 2.2, setting an upper limit for Li-rich giants formed by this mechanism. Giants with higher abundances need another mechanism to be explained. For reasonable planetary distributions, we predict the Li abundance distribution of low-mass giants undergoing planet engulfment, finding that between 1% and 3% of them should have A(Li)⩾1.5. We show that depending on the stellar mass range, this traditional definition of Li-rich giants is misleading, as isolated massive stars would be considered anomalous while giants engulfing a companion would be set aside, flagged as normal. We explore the detectability of companion engulfment, finding that planets with masses higher than ∼7 M{sub J} produce a distinct signature, and that descendants of stars originating in the Li dip and low-luminosity red giants are ideal tests of this channel.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The O{sub 2} A-Band in the Fluxes and Polarization of Starlight Reflected by Earth-Like Exoplanets

Earth-like, potentially habitable exoplanets are prime targets in the search for extraterrestrial life. Information about their atmospheres and surfaces can be derived by analyzing the light of the parent star reflected by the planet. We investigate the influence of the surface albedo A {sub s}, the optical thickness b {sub cloud}, the altitude of water clouds, and the mixing ratio of biosignature O{sub 2} on the strength of the O{sub 2} A-band (around 760 nm) in the flux and polarization spectra of starlight reflected by Earth-like exoplanets. Our computations for horizontally homogeneous planets show that small mixing ratios ( η < 0.4) will yield moderately deep bands in flux and moderate-to-small band strengths in polarization, and that clouds will usually decrease the band depth in flux and the band strength in polarization. However, cloud influence will be strongly dependent on properties such as optical thickness, top altitude, particle phase, coverage fraction, and horizontal distribution. Depending on the surface albedo and cloud properties, different O{sub 2} mixing ratios η can give similar absorption-band depths in flux and band strengths in polarization, especially if the clouds have moderate-to-high optical thicknesses. Measuring both the flux and the polarization is essential to reduce the degeneracies, although it will not solve them, especially not for horizontally inhomogeneous planets. Observations at a wide range of phase angles and with a high temporal resolution could help to derive cloud properties and, once those are known, the mixing ratio of O{sub 2} or any other absorbing gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Predictive coarse-graining

We propose a data-driven, coarse-graining formulation in the context of equilibrium statistical mechanics. In contrast to existing techniques which are based on a fine-to-coarse map, we adopt the opposite strategy by prescribing a probabilistic coarse-to-fine map. This corresponds to a directed probabilistic model where the coarse variables play the role of latent generators of the fine scale (all-atom) data. From an information-theoretic perspective, the framework proposed provides an improvement upon the relative entropy method and is capable of quantifying the uncertainty due to the information loss that unavoidably takes place during the coarse-graining process. Furthermore, it can be readily extended to a fully Bayesian model where various sources of uncertainties are reflected in the posterior of the model parameters. The latter can be used to produce not only point estimates of fine-scale reconstructions or macroscopic observables, but more importantly, predictive posterior distributions on these quantities. Predictive posterior distributions reflect the confidence of the model as a function of the amount of data and the level of coarse-graining. The issues of model complexity and model selection are seamlessly addressed by employing a hierarchical prior that favors the discovery of sparse solutions, revealing the most prominent features in the coarse-grained model. A flexible and parallelizable Monte Carlo – Expectation–Maximization (MC-EM) scheme is proposed for carrying out inference and learning tasks. A comparative assessment of the proposed methodology is presented for a lattice spin system and the SPC/E water model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Hydrological analysis in watersheds with a variable-resolution global climate model (VR-CESM)

Traditionally, watershed-scale hydrology is simulated by distributed hydrological models with offline meteorological forcing data, or by regional regional climate models that link atmospheric and land hydrology interactions. Global climate models (GCMs) are rarely used to study watershed-scale hydrology due to the relatively coarse grid resolution, computationally expensive downscaling, and simplified physical processes. Recently, however, watershed-scale hydrology analysis has become possible in GCMs due to the development of variable-resolution GCMs that dynamically couple the hydrological processes between atmosphere and land systems at fine resolutions in selected regions and coarse resolution elsewhere. In this study, we used the variable-resolution Community Earth System Model (VR-CESM) with refined-resolution (14 km) in the western U.S. and eastern China to evaluate watershed-scale hydrology. Furthermore, we compared the historical VR-CESM outputs with gauge measurements and other hydrological models (e.g., National Water Model in the U.S.) and calibrated the subsurface runoff capacities in four mountainous watersheds. An RCP8.5 projection from 2007 to 2050 is used to estimate the impact of changing precipitation and snow climatology on watershed hydrology. We also analyzed the long-term runoff variability and the possibility of extreme runoff events as simulated by the VR-CESM. Although calibration is not possible in larger-scale watersheds, VR-CESM simulates the long-term annual variability of watersheds and provides insights on climate change impacts on hydrology. We conclude that refined-resolution VR-CESM can be used for watershed-scale hydrology analysis to understand water resources and water balance, in addition to traditional watershed-scale hydrological models. It enables hydrological analysis at multiple watersheds in one simulation and can help understand the two-way dynamics between land surface hydrology and atmospheric processes, and is especially practical for projecting climate change impacts. However, it is challenging to apply VR-CESM for hydrologic analysis in regulated watersheds as human factors (e.g., pumping, irrigation, water diversion) have not been fully addressed in VR-CESM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Early-type eclipsing binaries with intermediate orbital periods

We analyze 221 eclipsing binaries (EBs) in the Large Magellanic Cloud with B-type main-sequence (MS) primaries (M{sub 1} ≈ 4–14 M{sub ⊙}) and orbital periods P = 20–50 days that were photometrically monitored by the Optical Gravitational Lensing Experiment. We utilize our three-stage automated pipeline to (1) classify all 221 EBs, (2) fit physical models to the light curves of 130 detached well-defined EBs from which unique parameters can be determined, and (3) recover the intrinsic binary statistics by correcting for selection effects. We uncover two statistically significant trends with age. First, younger EBs tend to reside in dustier environments with larger photometric extinctions, an empirical relation that can be implemented when modeling stellar populations. Second, younger EBs generally have large eccentricities. This demonstrates that massive binaries at moderate orbital periods are born with a Maxwellian “thermal” orbital velocity distribution, which indicates they formed via dynamical interactions. In addition, the age-eccentricity anticorrelation provides a direct constraint for tidal evolution in highly eccentric binaries containing hot MS stars with radiative envelopes. The intrinsic fraction of B-type MS stars with stellar companions q = M{sub 2}/M{sub 1} > 0.2 and orbital periods P = 20–50 days is (7 ± 2)%. We find early-type binaries at P = 20–50 days are weighted significantly toward small mass ratios q ≈ 0.2–0.3, which is different than the results from previous observations of closer binaries with P < 20 days. This indicates that early-type binaries at slightly wider orbital separations have experienced substantially less competitive accretion and coevolution during their formation in the circumbinary disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

HATS-7b: A HOT SUPER NEPTUNE TRANSITING A QUIET K DWARF STAR

We report the discovery by the HATSouth network of HATS-7b, a transiting Super-Neptune with a mass of 0.120 ± 0.012 M{sub J}, a radius of 0.563{sub −0.034}{sup +0.046} R{sub J}, and an orbital period of 3.1853 days. The host star is a moderately bright (V=13.340 ± 0.010 mag, K{sub S}=10.976 ± 0.026 mag) K dwarf star with a mass of 0.849 ± 0.027 M{sub ⊙}, a radius of 0.815{sub −0.035}{sup +0.049} R{sub ⊙}, and a metallicity of [Fe/H] =+0.250 ± 0.080. The star is photometrically quiet to within the precision of the HATSouth measurements, has low RV jitter, and shows no evidence for chromospheric activity in its spectrum. HATS-7b is the second smallest radius planet discovered by a wide-field ground-based transit survey, and one of only a handful of Neptune-size planets with mass and radius determined to 10% precision. Theoretical modeling of HATS-7b yields a hydrogen–helium fraction of 18 ± 4% (rock-iron core and H{sub 2}–He envelope), or 9 ± 4% (ice core and H{sub 2}–He envelope), i.e., it has a composition broadly similar to that of Uranus and Neptune, and very different from that of Saturn, which has 75% of its mass in H{sub 2}–He. Based on a sample of transiting exoplanets with accurately (<20%) determined parameters, we establish approximate power-law relations for the envelopes of the mass–density distribution of exoplanets. HATS-7b, which, together with the recently discovered HATS-8b, is one of the first two transiting super-Neptunes discovered in the Southern sky, is a prime target for additional follow-up observations with Southern hemisphere facilities to characterize the atmospheres of Super-Neptunes (which we define as objects with mass greater than that of Neptune, and smaller than halfway between that of Neptune and Saturn, i.e., 0.054 M{sub J}<M{sub p}<0.18 M{sub J}).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

The statistical multifragmentation model: Origins and recent advances

We review the Statistical Multifragmentation Model (SMM) which considers a generalization of the liquid-drop model for hot nuclei and allows one to calculate thermodynamic quantities characterizing the nuclear ensemble at the disassembly stage. We show how to determine probabilities of definite partitions of finite nuclei and how to determine, through Monte Carlo calculations, observables such as the caloric curve, multiplicity distributions, heat capacity, among others. Some experimental measurements of the caloric curve confirmed the SMM predictions of over 10 years before, leading to a surge in the interest in the model. However, the experimental determination of the fragmentation temperatures relies on the yields of different isotopic species, which were not correctly calculated in the schematic, liquid-drop picture, employed in the SMM. This led to a series of improvements in the SMM, in particular to the more careful choice of nuclear masses and energy densities, specially for the lighter nuclei. With these improvements the SMM is able to make quantitative determinations of isotope production. We show the application of SMM to the production of exotic nuclei through multifragmentation. These preliminary calculations demonstrate the need for a careful choice of the system size and excitation energy to attain maximum yields.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Fast Quasi-Static Time-Series (QSTS) for yearlong PV impact studies using vector quantization

The rapidly growing penetration levels of distributed photovoltaic (PV) systems requires more comprehensive studies to understand their impact on distribution feeders. IEEE P.1547 highlights the need for Quasi-Static Time Series (QSTS) simulation in conducting distribution impact studies for distributed resource interconnection. Unlike conventional scenario-based simulation, the time series simulation can realistically assess time-dependent impacts such as the operation of various controllable elements (e.g. voltage regulating tap changers) or impacts of power fluctuations. However, QSTS simulations are still not widely used in the industry because of the computational burden associated with running yearlong simulations at a 1-s granularity, which is needed to capture device controller effects responding to PV variability. This study presents a novel algorithm that reduces the number of times that the non-linear 3-phase unbalanced AC power flow must be solved by storing and reassigning power flow solutions as it progresses through the simulation. Each unique power flow solution is defined by a set of factors affecting the solution that can easily be queried. We demonstrate a computational time reduction of 98.9% for a yearlong simulation at 1-s resolution with minimal errors for metrics including: number of tap changes, capacitor actions, highest and lowest voltage on the feeder, line losses, and ANSI voltage violations. Finally, the key contribution of this work is the formulation of an algorithm capable of: (i) drastically reducing the computational time of QSTS simulations, (ii) accurately modeling distribution system voltage-control elements with hysteresis, and (iii) efficiently compressing result time series data for post-simulation analysis.

14 SOLAR ENERGY↗
DOE OSTI2015

Real-time sensing of epithelial cell-cell and cell-substrate interactions by impedance spectroscopy on porous substrates

Oxidized porous silicon (PS) is a common topographical biocompatible substrate that potentially provides a distinct in vitro environment for better understanding of in vivo behavior. But in the reported studies on oxidized PS, cell-cell and cell-substrate interactions have been detected only by fluorescent labeling. This paper is the first attempt to investigate real-time sensing of these interactions on HaCaT cells by label-free impedance spectroscopy on oxidized PS of two pore diameters (50 and 500 nm). One of the major requirements for successful impedance spectroscopy measurement is to restrict the channeling of electric field lines through the pores. To satisfy this criterion, we have designed the pore depths after analyzing the penetration of the medium by using computational fluid dynamics simulation. A distributed electrical model was also developed for estimating the various cellular attributes by considering a pseudorandom distribution of pores. It is observed from the impedance measurements and from the model that the proliferation rate increases for 50 nm pores but decreases for 500 nm pores compared to that for planar substrates. The rate of decrease in cell substrate separation (h) in the initial stage is more than the rate of increase in cell-cell junction resistance (R{sub b}) corresponding to the initial adhesion phase of cells. It is observed that R{sub b} and h are higher for 50 nm pores than those for planar substrates, corresponding to the fact that substrates more conducive toward cell adhesion encourage cell-cell interactions than direct cell-substrate interactions. Thus, the impedance spectroscopy coupled with the proposed theoretical framework for PS substrates can sense and quantify the cellular interactions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

THERMAL CHARACTERISTICS AND THE DIFFERENTIAL EMISSION MEASURE DISTRIBUTION DURING A B8.3 FLARE ON 2009 JULY 4

We investigate the evolution of the differential emission measure distribution (DEM[ T ]) in various phases of a B8.3 flare which occurred on 2009 July 04. We analyze the soft X-ray (SXR) emission in the 1.6–8.0 keV range, recorded collectively by the Solar Photometer in X-rays (SphinX; Polish) and the Solar X-ray Spectrometer (Indian) instruments. We conduct a comparative investigation of the best-fit DEM[ T ] distributions derived by employing various inversion schemes, namely, single Gaussian, power-law functions and a Withbroe–Sylwester (W–S) maximum likelihood algorithm. In addition, the SXR spectrum in three different energy bands, that is, 1.6–5.0 keV (low), 5.0–8.0 keV (high), and 1.6–8.0 keV (combined), is analyzed to determine the dependence of the best-fit DEM[ T ] distribution on the selection of the energy interval. The evolution of the DEM[ T ] distribution, derived using a W–S algorithm, reveals multi-thermal plasma during the rise to the maximum phase of the flare, and isothermal plasma in the post-maximum phase of the flare. The thermal energy content is estimated by considering the flare plasma to be (1) isothermal and (2) multi-thermal in nature. We find that the energy content during the flare, estimated using the multi-thermal approach, is in good agreement with that derived using the isothermal assumption, except during the flare maximum. Furthermore, the (multi-) thermal energy estimated while employing the low-energy band of the SXR spectrum results in higher values than that derived from the combined energy band. On the contrary, the analysis of the high-energy band of the SXR spectrum leads to lower thermal energy than that estimated from the combined energy band.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Stochastic acceleration of electrons by fast magnetosonic waves in solar flares: the effects of anisotropy in velocity and wavenumber space

We develop a model for stochastic acceleration of electrons in solar flares. As in several previous models, the electrons are accelerated by turbulent fast magnetosonic waves ({sup f}ast waves{sup )} via transit-time-damping (TTD) interactions. (In TTD interactions, fast waves act like moving magnetic mirrors that push the electrons parallel or anti-parallel to the magnetic field). We also include the effects of Coulomb collisions and the waves' parallel electric fields. Unlike previous models, our model is two-dimensional in both momentum space and wavenumber space and takes into account the anisotropy of the wave power spectrum F{sub k} and electron distribution function f {sub e}. We use weak turbulence theory and quasilinear theory to obtain a set of equations that describes the coupled evolution of F{sub k} and f {sub e}. We solve these equations numerically and find that the electron distribution function develops a power-law-like non-thermal tail within a restricted range of energies E in (E {sub nt}, E {sub max}). We obtain approximate analytic expressions for E {sub nt} and E {sub max}, which describe how these minimum and maximum energies depend upon parameters such as the electron number density and the rate at which fast-wave energy is injected into the acceleration region at large scales. We contrast our results with previous studies that assume that F{sub k} and f {sub e} are isotropic, and we compare one of our numerical calculations with the time-dependent hard-X-ray spectrum observed during the 1980 June 27 flare. In our numerical calculations, the electron energy spectra are softer (steeper) than in models with isotropic F{sub k} and f {sub e} and closer to the values inferred from observations of solar flares.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Sizing up the milky way: a bayesian mixture model meta-analysis of photometric scale length measurements

The exponential scale length (L {sub d}) of the Milky Way’s (MW’s) disk is a critical parameter for describing the global physical size of our Galaxy, important both for interpreting other Galactic measurements and helping us to understand how our Galaxy fits into extragalactic contexts. Unfortunately, current estimates span a wide range of values and are often statistically incompatible with one another. Here, we perform a Bayesian meta-analysis to determine an improved, aggregate estimate for L {sub d}, utilizing a mixture-model approach to account for the possibility that any one measurement has not properly accounted for all statistical or systematic errors. Within this machinery, we explore a variety of ways of modeling the nature of problematic measurements, and then employ a Bayesian model averaging technique to derive net posterior distributions that incorporate any model-selection uncertainty. Our meta-analysis combines 29 different (15 visible and 14 infrared) photometric measurements of L {sub d} available in the literature; these involve a broad assortment of observational data sets, MW models and assumptions, and methodologies, all tabulated herein. Analyzing the visible and infrared measurements separately yields estimates for L {sub d} of 2.71{sub −0.20}{sup +0.22} kpc and 2.51{sub −0.13}{sup +0.15} kpc, respectively, whereas considering them all combined yields 2.64 ± 0.13 kpc. The ratio between the visible and infrared scale lengths determined here is very similar to that measured in external spiral galaxies. We use these results to update the model of the Galactic disk from our previous work, constraining its stellar mass to be 4.8{sub −1.1}{sup +1.5}×10{sup 10} M {sub ⊙}, and the MW’s total stellar mass to be 5.7{sub −1.1}{sup +1.5}×10{sup 10} M {sub ⊙}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

On extreme-ultraviolet helium line intensity enhancement factors on the sun

Helium lines in the solar EUV spectrum provide useful diagnostics of the solar atmosphere plasma. Helium is one of the few elements that exhibits strong emission lines formed in the lower transition region, and it is the second most abundant element in the Sun. However, the analysis of helium lines is complicated by their optical thickness and the unusual behavior of their intensities, with enhancements by a factor of up to 15 reported in the literature. Detailed study requires spatially and spectrally resolved observations in the EUV range, as well as sophisticated atomic modeling. The present work focuses on the application of the differential emission measure distribution to reproduce the observed fluxes of the He i and He ii lines observed by the Solar and Heliospheric Observatory (SOHO)/Coronal Diagnostic Spectrometer and Hinode/EIS spectrometers, using the latest atomic data. It is found that the comparison between observed and reconstructed intensities for He i resonance lines λλ537.03, 522.21, and 515.62 and the intercombination line λ591.41 does not show a real enhancement. By contrast, He i λ584.33, the first line of the 1s{sup 2} {sup 1}S−1s np {sup 1}P resonance series, shows a depletion of a factor ∼2, due to the opacity effect, as supported by non-LTE radiative transfer calculations. For single ionized helium lines 303.78 and 256.32 Å, the enhancement factors obtained are higher and agree with those of previous work. The different behavior of He i and He ii lines suggests a mechanism that affects ionized helium only.

79 ASTRONOMY AND ASTROPHYSICS↗