Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STAR MODELS”

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 199 records · Page 11

COMPARING NUMERICAL METHODS FOR ISOTHERMAL MAGNETIZED SUPERSONIC TURBULENCE

Many astrophysical applications involve magnetized turbulent flows with shock waves. Ab initio star formation simulations require a robust representation of supersonic turbulence in molecular clouds on a wide range of scales imposing stringent demands on the quality of numerical algorithms. We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine popular astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. These applications employ a variety of numerical approaches, including both split and unsplit, finite difference and finite volume, divergence preserving and divergence cleaning, a variety of Riemann solvers, and a range of spatial reconstruction and time integration techniques. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss the convergence of various characteristics for the turbulence decay test and the impact of various components of numerical schemes on the accuracy of solutions. The nine codes gave qualitatively the same results, implying that they are all performing reasonably well and are useful for scientific applications. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method and using little to no explicit artificial viscosity. Codes that fall short in one or more of these areas are still useful, but they must compensate for higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices about choosing optimal applications for their specific astrophysical problems.

79 ASTRONOMY AND ASTROPHYSICS↗

A MAGNETIC RIBBON MODEL FOR STAR-FORMING FILAMENTS

We develop a magnetic ribbon model for molecular cloud filaments. These result from turbulent compression in a molecular cloud in which the background magnetic field sets a preferred direction. We argue that this is a natural model for filaments and is based on the interplay between turbulence, strong magnetic fields, and gravitationally driven ambipolar diffusion, rather than pure gravity and thermal pressure. An analytic model for the formation of magnetic ribbons that is based on numerical simulations is used to derive a lateral width of a magnetic ribbon. This differs from the thickness along the magnetic field direction, which is essentially the Jeans scale. We use our model to calculate a synthetic observed relation between apparent width in projection versus observed column density. The relationship is relatively flat, similar to observations, and unlike the simple expectation based on a Jeans length argument.

79 ASTRONOMY AND ASTROPHYSICS↗

UH cosmic rays: Possible origin in massive stars

The origin of the Z greater than 28, ultraheavy, cosmic rays in supernova explosions of massive stars is considered. For Z greater than 70, the UH data is dominated by an r-process source distribution, but for the elements just beyond iron, 29 or = Z less than 36, the data cannot be explained by any single process of nucleosynthesis. This problem is solved naturally in a massive star model by secondary neutron capture reactions occuring during core helium burning and during explosive carbon burning. Interstellar propagation calculations were performed with these episodes of synthesis as source distributions, and the results offer an explanation for the current UH cosmic-ray data. The heavy element synthesis during explosive carbon burning is reexamined using more realistic initial conditions given by the post-helium-burning configuration of the star. Effects of preferential acceleration are considered, and experimental tests are discussed.

Wefel, J. P.↗

X-RAY OBSERVATIONS OF RADIO TRANSIENTS WITHOUT OPTICAL HOSTS

We present a 50 ks Chandra/ACIS-I X-ray observation of the Bower et al. Very Large Array archival field. The observations reach a limiting sensitivity of {approx}10{sup -4} counts s{sup -1}, corresponding to a flux of a few times 10{sup -15} erg s{sup -1} cm{sup -2} for the models we explore. The Chandra observations were undertaken to search for X-ray counterparts to the eight transient sources without optical counterparts, and the two transient sources with optical counterparts seen by Bower et al. Neither of the sources with optical counterparts was detected in X-rays. One of the eight optical non-detections is associated with a marginal (2.4{sigma}) X-ray detection in our Chandra image. A second optically undetected Bower et al. transient may be associated with a z = 1.29 X-ray-detected quasar or its host galaxy, or alternatively is undetected in X-rays and is a chance association with the nearby X-ray source. The X-ray flux upper limits, and the one marginal detection, are consistent with the interpretation of Ofek et al. that the optically undetected radio transients are flares from isolated old Galactic neutron stars. The marginal X-ray detection has a hardness ratio that implies a temperature too high for a simple one-temperature neutron star model, but plausible multi-component fits are not excluded, and in any case the marginal X-ray detection may be due to cosmic rays or particle background. The X-ray flux upper limits are also consistent with flare star progenitors at {approx}> 1 kpc (which would require the radio luminosity of the transient to be unusually high for such an object) or less extreme flares from brown dwarfs at distances of around 100 pc.

79 ASTRONOMY AND ASTROPHYSICS↗

Star forming regions in gas-rich SO galaxies

The first results of an H alpha imaging survey of HI rich SO galaxies, which were searched for HII regions and other sources of emission, are presented. The charge coupled device H alpha interference filter images were made of 16 galaxies. Eight of these galaxies show evidence for on-going star formation, one has nuclear emission but no HII regions, and the remaining seven have no emissions detected within well defined upper limits. With the exception of one notably peculiar galaxy in which the emission from HII regions appears pervasive, the HII regions are either organized into inner-disk rings or randomly distributed throughout the disk. A few of these galaxies are found to be clearly not SO's; or peculiar objects atypical of the SO class. Using simple models star formation rates (SFRs) and gas depletion times from the observed H alpha fluxes were estimated. In general, the derived SFRs are much lower than those found in isolated field spiral galaxies and the corresponding gas depletion time scales are also longer.

Pogge, Richard W.↗

Coordinated X-ray optical and radio observations of YZ Canis Minoris

Coordinated X ray, optical, and radio observations of the flare star YZ CMi are reported. Twenty-two minor optical flares and twelve radio events were recorded. No major optical flares, greater than 3 magnitudes, were observed. Although no flare related X ray emission was observed, the measured upper limits in this band enable meaningful comparisons with published flare star models. Three of the five models predicting the relative X ray to optical or radio flare luminosities are in serious disagreement with the observations. For the largest optical flare with coincident X ray coverage, the 3 sigma upper limit on X ray emission in the 0.15 to 0.8 keV band is 8.7 x 10 to the 28th power erg/s, corresponding to a ratio of X ray to B-band luminosity of less than 0.3. Based on the present results, the contribution of the flares UV Ceti flare stars to the galactic component of the diffuse soft X ray background is less than 0.2 percent.

Karpen, J. T.↗

Interpreting short gamma-ray burst progenitor kicks and time delays using the host galaxy-dark matter halo connection

Nearly 20% of short gamma-ray bursts (sGRBs) have no observed host galaxies. Combining this finding with constraints on galaxies' dark matter halo potential wells gives strong limits on the natal kick velocity distribution for sGRB progenitors. For the best-fitting velocity distribution, one in five sGRB progenitors receives a natal kick above 150 km s{sup –1}, consistent with merging neutron star models but not with merging white dwarf binary models. This progenitor model constraint is robust to a wide variety of systematic uncertainties, including the sGRB progenitor time-delay model, the Swift redshift sensitivity, and the shape of the natal kick velocity distribution. We also use constraints on the galaxy-halo connection to determine the host halo and host galaxy demographics for sGRBs, which match extremely well with available data. Most sGRBs are expected to occur in halos near 10{sup 12} M {sub ☉} and in galaxies near 5 × 10{sup 10} M {sub ☉} (L {sub *}); unobserved faint and high-redshift host galaxies contribute a small minority of the observed hostless sGRB fraction. We find that sGRB redshift distributions and host galaxy stellar masses weakly constrain the progenitor time-delay model; the active versus passive fraction of sGRB host galaxies may offer a stronger constraint. Finally, we discuss how searches for gravitational wave optical counterparts in the local universe can reduce follow-up times using these findings.

79 ASTRONOMY AND ASTROPHYSICS↗

Analytic closures for M1 neutrino transport

Carefully accounting for neutrino transport is an essential component of many astrophysical studies. Solving the full transport equation is too expensive for most realistic applications, especially those involving multiple spatial dimensions. For such cases, resorting to approximations is often the only viable option for obtaining solutions. One such approximation, which recently became popular, is the M1 method. It utilizes the system of the lowest two moments of the transport equation and closes the system with an ad hoc closure relation. The accuracy of the M1 solution depends on the quality of the closure. Several closures have been proposed in the literature and have been used in various studies. We carry out an extensive study of these closures by comparing the results of M1 calculations with precise Monte Carlo calculations of the radiation field around spherically symmetric protoneutron star models. We find that no closure performs consistently better or worse than others in all cases. The level of accuracy that a given closure yields depends on the matter configuration, neutrino type and neutrino energy. As a result, given this limitation, the maximum entropy closure by Minerbo on average yields relatively accurate results in the broadest set of cases considered in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

A computer program for modeling non-spherical eclipsing binary star systems

The accurate analysis of eclipsing binary light curves is fundamental to obtaining information on the physical properties of stars. The model described accounts for the important geometric and photometric distortions such as rotational and tidal distortion, gravity brightening, and reflection effect. This permits a more accurate analysis of interacting eclipsing star systems. The model is designed to be useful to anyone with moderate computing resources. The programs, written in FORTRAN 4 for the IBM 360, consume about 80k bytes of core. The FORTRAN program listings are provided, and the computational aspects are described in some detail.

Wood, D. B.↗

Synthetic Stromgren photometry for F dwarf stars

Recent synthetic spectrum and color calculations for cool dwarf star models are tested by comparison with observation. The accuracy of the computed dependence of the thermal colors B-V and b-y on effective temperature is examined, and H-beta indices are presented and compared with observed values. The accuracy of the predictions of the Stromgren uvby system metal-abundance indicator m1 and luminosity indicator c1 are tested. A new calibration of the c1, b-y diagram in terms of absolute magnitudes is given, making use of recent calculations of stellar isochrones. Observations of very metal-poor subdwarfs are used to study the accuracy of the isochrones. The c1, b-y diagram of the subdwarfs is compared with that of the turnoff-region stars in the very metal-poor globular cluster NGC 6397.

Bell, R. A.↗

Universal Scaling Laws for Solar and Stellar Atmospheric Heating: Catalog of Power-law Index between Solar Activity Proxies and Various Spectral Irradiances

The formation of extremely hot outer atmospheres is one of the most prominent manifestations of magnetic activity common to late-type dwarf stars, including the Sun. It is widely believed that these atmospheric layers, the corona, transition region, and chromosphere, are heated by the dissipation of energy transported upwards from the stellar surface by the magnetic field. This is signified by the spectral line fluxes at various wavelengths, scaled with power-law relationships against the surface magnetic flux over a wide range of formation temperatures, which are universal to the Sun and Sunlike stars of different ages and activity levels. This study describes a catalog of power-law indices between solar activity proxies and various spectral line fluxes. Compared to previous studies, we expanded the number of proxies, which now includes the total magnetic flux, total sunspot number, total sunspot area, and the F10.7 cm radio flux, and further enhanced the number of spectral lines by a factor of 2. This provides the data to study in detail the flux–flux scaling laws from the regions specified by the temperatures of the corona (log(T/K) = 6–7) to those of the chromosphere (log(T/K) ∼ 4), as well as the reconstruction of various spectral line fluxes of the Sun in the past, F-, G-, and K-type dwarfs, and the modeled stars.

Stellar coronae↗

Evaluation of DNB Predictive Capability in Multiphase Boiling Model of STAR-CCM+

The current milestone report is to evaluate the predictive capability of the existing boiling model in the commercial Computational Fluid Dynamics (CFD) package (so called, STARCCM+ 11.04.012), and validate the numerical calculation of a boiling curve with identification of departure from nucleate boiling (DNB) point by comparing with experimental measurements in the literatures. The current work directly demonstrates the DNB predictive capability in multiphase boiling model for nuclear reactor applications, which is tightly coupled with the CASL mission goal and evolving need from the nuclear industry. The boiling model tested here is based on the heat flux partitioning RPI approach (so called Generation-I boiling model in STAR-CCM+). A traditional experimental DNB test with a vertical pipe flow performed by Argonne National Laboratory (ANL), and a realistic 5x5-bundle fuel bundle with non-mixing vane DNB test conducted by Westinghouse Electric Company (WEC) are selected as reference experimental measurement. A baseline DNB simulation methodology is reported and its predictive capability is evaluated qualitatively as well as quantitatively with existing experiments. A parametric study of the microscopic boiling closure relations is investigated. For the pipe flow DNB validation study, the calculated DNB from the current study demonstrates good agreement with measured DNB data. The maximum deviation is observed less than 20%. The trend behaviors of DNB in varying operating conditions are similar with the observation reported in the literature. In addition, a meshing strategy for 5x5 bundle case is reported, some challenges on 5x5 bundle DNB model are also discussed. The current study provide a fundamental insight on the usefulness and limitations of the generation-I boiling model approach and may contribute to an advanced boiling model development for a practical nuclear thermal hydraulic application.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Constraints for models of Be stars derived from UV and IRAS observations

The effect of different geometrical models (disk and circumstellar shell) on the mass loss determinations of Be stars is investigated. The characteristics of the disk and circumstellar shell models are described. The models are applied to the analysis of the column density derived from the UV lines and emission measurements derived from IRAS observations of 10 Be stars, and the mass loss rate for these stars are calculated. It is observed that the mass loss rates obtained from the circumstellar shell model are a factor of 0.01 smaller than those of the disk model; the mass of the IR emitting region is three times larger in the shell model than that in the disk model; and the radial outflow velocity of the IR emitting region for the shell model is of the order of 0.01 km/sec and for the disk model of the order of 10 km/sec. It is noted that the UV and IRAS observations can be analyzed using the two models if the proper geometry is considered.

Lamers, H. J. G. L. M.↗

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. Furthermore, to improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 to study the dependence of the results on these parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Rotation in solar-type stars. I - Evolutionary models for the spin-down of the sun

Models of rotating, 1 solar mass stars through the pre-main sequence and main sequence stages have been calculated. The calculations begin on the Hayashi track with rapid, rigid-body rotation, with angular momentum removed from the convective envelope to simulate the effect of the solar wind as the models evolve. Interior angular momentum redistribution by circulation flows and rotational instabilities is computed by a diffusion technique. It is found that at the present age of solar evolution, the models employed have internal rotation rates significantly larger than those allowed by measured solar oblateness, suggesting that the effects of circulation currents near the outer boundary of the radiative interior have been underestimated. An important result of this study is the implication that circulation currents produce deep mixing of such minor constituents as He-3 and C-13, perhaps also explaining the observed steady depletion of lithium during the main-sequence stage.

Endal, A. S.↗

PROBING PRE-GALACTIC METAL ENRICHMENT WITH HIGH-REDSHIFT GAMMA-RAY BURSTS

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Population III (Pop III) GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature {approx}> 10{sup 4} K. We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1 + z) with values of about a few cm{sup -3}. In more massive halos, corresponding to the first galaxies, the density may be larger, n {approx}> 100 cm{sup -3}. The resulting afterglow fluxes are weakly dependent on redshift at a fixed observed time, and may be detectable with the James Webb Space Telescope and Very Large Array in the near-IR and radio wavebands, respectively, out to redshift z {approx}> 20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. The metal absorption line signature is expected to be detectable in the near future. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages (C II, O I, Si II and Fe II). We show that, if the afterglow can be observed sufficiently early, analysis of the metal lines may distinguish whether the first heavy elements were produced in a pair-instability supernova or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing Pre-Galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Population III (Pop III) GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature approximately greater than10(exp 4) K.We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1 + z) with values of about a few cm(exp −3). In more massive halos, corresponding to the first galaxies, the density may be larger, n approximately greater than100 cm(exp −3). The resulting afterglow fluxes are weakly dependent on redshift at a fixed observed time, and may be detectable with the James Webb Space Telescope and Very Large Array in the near-IR and radio wavebands, respectively, out to redshift z approximately greater than 20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. The metal absorption line signature is expected to be detectable in the near future. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages (C II, O I, Si II and Fe II). We show that, if the afterglow can be observed sufficiently early, analysis of the metal lines may distinguish whether the first heavy elements were produced in a pair-instability supernova or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

Metal↗