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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Modeling Impacts of Dust Mineralogy on Fast Climate Response

Mineralogical composition drives dust impacts on Earth's climate systems. However, most climate models still use homogeneous dust, without accounting for the temporal and spatial variation in mineralogy. To quantify the radiative impact of resolving dust mineralogy on Earth's climate, we implement and simulate the distribution of dust minerals (i.e., illite, kaolinite, smectite, hematite, calcite, feldspar, quartz, and gypsum) from Claquin et al. (1999) (C1999) and activate their interaction with radiation in the GFDL AM4.0 model. Resolving mineralogy reduces dust absorption compared to the homogeneous dust used in the standard GFDL AM4.0 model that assumes a globally uniform hematite volume content of 2.7 % (HD27). The reduction in dust absorption results in improved agreement with observation-based single-scattering albedo (SSA), radiative fluxes from CERES (the Clouds and the Earth's Radiant Energy System), and land surface temperature from the CRU (Climatic Research Unit) compared to the baseline HD27 model version. It also results in distinct radiative impacts on Earth's climate over North Africa. Over the 19-year (from 2001 to 2019) modeled period during JJA (June–July–August), the reduction in dust absorption in AM4.0 leads to a reduction of over 50 % in net downward radiation across the Sahara and approximately 20 % over the Sahel at the top of the atmosphere (TOA) compared to the baseline HD27 model version. The reduced dust absorption weakens the atmospheric warming effect of dust aerosols and leads to an alteration in land surface temperature, resulting in a decrease of 0.66 K over the Sahara and an increase of 0.7 K over the Sahel. The less warming in the atmosphere suppresses ascent and weakens the monsoon inflow from the Gulf of Guinea. This brings less moisture to the Sahel, which combined with decreased ascent induces a reduction of precipitation. To isolate the effect of reduced absorption compared to resolving spatial and temporal mineralogy, we carry out a simulation where the hematite volume content of homogeneous dust is reduced from 2.7 % to 0.9 % (HD09). The dust absorption (e.g., single-scattering albedo) of HD09 is comparable to that of the mineralogically speciated model on a global mean scale, albeit with a lower spatial variation that arises solely from particle size. Comparison of the two models indicates that the spatial inhomogeneity in dust absorption resulting from resolving mineralogy does not have significant impacts on Earth's radiation and climate, provided there is a similar level of dust absorption on a global mean scale before and after resolving dust mineralogy. However, uncertainties related to emission and distribution of minerals may blur the advantages of resolving minerals to study their impact on radiation, cloud properties, ocean biogeochemistry, air quality, and photochemistry. On the other hand, lumping together clay minerals (i.e., illite, kaolinite, and smectite), but excluding externally mixed hematite and gypsum, appears to provide both computational efficiency and relative accuracy. Nevertheless, for specific research, it may be necessary to fully resolve mineralogy to achieve accuracy.

Soil dust↗
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↗
NASA NTRSDate not supplied

Evolving the SpaceFOM: Lessons Learned and Future Development

For over fifty years, simulation has been a cornerstone technology for space missions. In many cases, one single simulator, developed by one team, is not enough to meet the simulation objectives. Astronaut crew training, docking, vehicle testing and complex mission planning all require several simulations to be connected and to interoperate, in order to achieve the goal of the simulation. To facilitate the development of interoperable simulations, an open standard, the SISO SpaceFOM was developed in collaboration between government, industry and academia. The standard was released in 2020. The initial release of SpaceFOM focuses on the core interoperability topics for space simulation: (i) specifying spatial data using several, well defined reference frames, (ii) managing time (including real-time, hard real-time, and as-fast-as-possible execution), (iii) composition of systems with subsystems such as space vehicles, and (iv) execution control, such as initializing, running, pausing, and terminating simulations. A number of requirements for the next version of SpaceFOM have been identified and proposed solutions are already under development. One key area is the object classes used for data exchange, where a richer set of classes have been proposed, based on practical use cases. For time management, the use of variable time resolutions and variable time steps is also proposed. For execution control, save and restore of simulations is considered as well as improved robustness, to handle failing components. A number of clarifications are also under consideration. This paper provides both a forward-looking view of the technical developments underway and a retrospective on the standard's evolution. The authors emphasize that prioritizing core challenges and delivering a practical first version were crucial steps. While not every idea was implemented in the initial release, the timely delivery of a functional standard enabled real-world applications, the growth of a user community, and the collection of valuable insights to shape future versions. The authors look forward to working with the next version of SpaceFOM together with the growing user community.

FOMs↗
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↗
NASA NTRS2001

Models of Small-Scale Patchiness

Patchiness is perhaps the most salient characteristic of plankton populations in the ocean. The scale of this heterogeneity spans many orders of magnitude in its spatial extent, ranging from planetary down to microscale. It has been argued that patchiness plays a fundamental role in the functioning of marine ecosystems, insofar as the mean conditions may not reflect the environment to which organisms are adapted. Understanding the nature of this patchiness is thus one of the major challenges of oceanographic ecology. The patchiness problem is fundamentally one of physical-biological-chemical interactions. This interconnection arises from three basic sources: (1) ocean currents continually redistribute dissolved and suspended constituents by advection; (2) space-time fluctuations in the flows themselves impact biological and chemical processes, and (3) organisms are capable of directed motion through the water. This tripartite linkage poses a difficult challenge to understanding oceanic ecosystems: differentiation between the three sources of variability requires accurate assessment of property distributions in space and time, in addition to detailed knowledge of organismal repertoires and the processes by which ambient conditions control the rates of biological and chemical reactions. Various methods of observing the ocean tend to lie parallel to the axes of the space/time domain in which these physical-biological-chemical interactions take place. Given that a purely observational approach to the patchiness problem is not tractable with finite resources, the coupling of models with observations offers an alternative which provides a context for synthesis of sparse data with articulations of fundamental principles assumed to govern functionality of the system. In a sense, models can be used to fill the gaps in the space/time domain, yielding a framework for exploring the controls on spatially and temporally intermittent processes. The following discussion highlights only a few of the multitude of models which have yielded insight into the dynamics of plankton patchiness. In addition, this particular collection of examples is intended to furnish some exposure to the diversity of modeling approaches which can be brought to bear on the problem. These approaches range from abstract theoretical models intended to elucidate specific processes, to complex numerical formulations which can be used to actually simulate observed distributions in detail.

McGillicuddy, D. J.↗
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↗
DOE OSTI2014

Orbital perturbations of the Galilean satellites during planetary encounters

The Nice model of the dynamical instability and migration of the giant planets can explain many properties of the present solar system, and can be used to constrain its early architecture. In the jumping-Jupiter version of the Nice model, required from the terrestrial planet constraint and dynamical structure of the asteroid belt, Jupiter has encounters with an ice giant. Here, we study the survival of the Galilean satellites in the jumping-Jupiter model. This is an important concern because the ice-giant encounters, if deep enough, could dynamically perturb the orbits of the Galilean satellites and lead to implausible results. We performed numerical integrations where we tracked the effect of planetary encounters on the Galilean moons. We considered three instability cases from Nesvorný and Morbidelli that differed in the number and distribution of encounters. We found that in one case, where the number of close encounters was relatively small, the Galilean satellite orbits were not significantly affected. In the other two, the orbital eccentricities of all moons were excited by encounters, Callisto's semimajor axis changed, and, in a large fraction of trials, the Laplace resonance of the inner three moons was disrupted. The subsequent evolution by tides damps eccentricities and can recapture the moons in the Laplace resonance. A more important constraint is represented by the orbital inclinations of the moons, which can be excited during the encounters and not appreciably damped by tides. We find that one instability case taken from Nesvorný and Morbidelli clearly does not meet this constraint. This shows how the regular satellites of Jupiter can be used to set limits on the properties of encounters in the jumping-Jupiter model, and help us to better understand how the early solar system evolved.

79 ASTRONOMY AND ASTROPHYSICS↗