Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Athena”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

The THESEUS space mission concept: science case, design and expected performances

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5–1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift ~10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late ’20s/early ’30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. In conclusion, THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A fourth-order accurate finite volume method for ideal MHD via upwind constrained transport

Here we present a fourth-order accurate finite volume method for the solution of ideal magnetohydrodynamics (MHD). The numerical method combines high-order quadrature rules in the solution of semi-discrete formulations of hyperbolic conservation laws with the upwind constrained transport (UCT) framework to ensure that the divergence-free constraint of the magnetic field is satisfied. A novel implementation of UCT that uses the piecewise parabolic method (PPM) for the reconstruction of magnetic fields at cell corners in 2D is introduced. The resulting scheme can be expressed as the extension of the second-order accurate constrained transport (CT) Godunov-type scheme that is currently used in the Athena astrophysics code. After validating the base algorithm on a series of hydrodynamics test problems, we present the results of multidimensional MHD test problems which demonstrate formal fourth-order convergence for smooth problems, robustness for discontinuous problems, and improved accuracy relative to the second-order scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Simulation studies related to the particle identification by the forward and backward RICH detectors at Electron Ion Collider

The Electron-Ion collider (EIC) will be the ultimate facility to study the dynamics played by colored quarks and gluons in the phenomenology of nucleons and nuclei, described by Quantum Chromodynamics. The physics programs will greatly rely on efficient particle identification (PID) in both the forward and the backward regions. The forward and the backward RICHes of the EIC have to be able to cover wide acceptance and momentum ranges; in the forward region a dual radiator RICH (dRICH) is foreseen and in the backward region a proximity-focusing RICH can be foreseen to be employed. The geometry and the performance studies of the dRICH have been performed as prescribed in the EIC Yellow Report using the ATHENA software framework. Furthermore, this part of our work reports the effort following the call for EIC detector proposal and the studies related to the forward and the backward RICH performance. In the forward region, the dRICH performance showed a pion- kaon separation from around 1 GeV/c to 50 GeV/c at a three sigma level; the proximity focusing RICH (pfRICH) foreseen for the backward region can reach three sigma separation up to 3 GeV/c for e/$π$ and up to 10 GeV/c for $π$/K mass hypotheses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Stable and unstable supersonic stagnation of an axisymmetric rotating magnetized plasma

The Naval Research Laboratory ‘Mag Noh problem’, described in this paper, is a self-similar magnetized implosion flow, which contains a fast magnetohydrodynamic (MHD) outward propagating shock of constant velocity. We generalize the classic Noh ( OSTI Tech. Rep. 577058, 1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five parametric, each solution having its own self-similarity index, gas gamma, magnetization, the ratio of axial to the azimuthal field and rotation. While the classic Noh problem must have a supersonic implosion velocity to create a shock, our solutions have an interesting three-parametric special case with zero initial velocity in which magnetic tension, instead of implosion flow, creates the shock at $t=0+$ . Our self-similar solutions are indeed realized when we solve the initial value MHD problem with the finite volume MHD code Athena. We numerically investigated the stability of these solutions and found both stable and unstable regions in parameter space. Stable solutions can be used to test the accuracy of numerical codes. Unstable solutions have also been widely used to test how codes reproduce linear growth, transition to turbulence and the practically important effects of mixing. Now we offer a family of unstable solutions featuring all three elements relevant to magnetically driven implosions: convergent flow, magnetic field and a shock wave.

Mechanics↗
DOE OSTI2020

High-resolution X-ray spectroscopy of the stellar wind in Vela X-1 during a flare

In this work, we present a ~130 ks observation of the prototypical wind-accreting, high-mass X-ray binary Vela X-1 collected with XMM-Newton at orbital phases between 0.12 and 0.28. A strong flare took place during the observation that allows us to investigate the reaction of the clumpy stellar wind to the increased X-ray irradiation. To examine the wind’s reaction to the flare, we performed both time-averaged and time-resolved analyses of the RGS spectrum and examined potential spectral changes. We focused on the high-resolution XMM-Newton RGS spectra and divided the observation into pre-flare, flare, and post-flare phases. We modeled the time-averaged and time-resolved spectra with phenomenological components and with the self-consistent photoionization models calculated via CLOUDY and XSTAR in the pre-flare phase, where strong emission lines due to resonant transitions of highly ionized ions are seen. In the spectra, we find emission lines corresponding to K-shell transitions in highly charged ions of oxygen, neon, magnesium, and silicon as well as radiative recombination continua (RRC) of oxygen. Additionally, we observe potential absorption lines of magnesium at a lower ionization stage and features identified as iron L lines. The, CLOUDY and XSTAR photoionization models provide contradictory results, either pointing towards uncertainties in theory or possibly a more complex multi-phase plasma, or both. We are able to demonstrate the existence of a plethora of variable narrow features, including the firm detection of oxygen lines and RRC that RGS enables to observe in this source for the first time. We show that Vela X-1 is an ideal source for future high-resolution missions, such as XRISM and Athena.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Looking through the photoionisation wake: Vela X-1 at φ orb ≈ 0.75 with Chandra /HETG

In this work, we detect and identify five narrow radiative recombination continua (Mg XI-XII, Ne IX-X, O VIII) and several emission lines from Fe, S, Si, Mg, Ne, Al, and Na, including four He-like triplets (S XV, Si XIII, Mg XI, and Ne IX). Photoionisation models reproduce the overall spectrum well, except for the near-neutral fluorescence lines of Fe, S, and Si. We conclude that the plasma is mainly photoionised, but more than one component is most likely present, which is consistent with a multi-phase plasma scenario, where denser and colder clumps of matter are embedded in the hot, photoionised wind of the companion star. Simulations with the future X-ray satellites Athena and XRISM show that a few hundred seconds of exposure is sufficient to disentangle the lines of the Fe Kα doublet and the He-like Fe XXV, improving, in general, the determination of the plasma parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Framework for custom event sample augmentations for ATLAS analysis data

For HEP event processing, data is typically stored in column-wise synchronized containers, such as most prominently ROOT’s TTree, which have been used for several decades to store by now over 1 exabyte. These containers can combine row-wise association capabilities needed by most HEP event processing frameworks (e.g. Athena for ATLAS) with column-wise storage, which typically results in better compression and more efficient support for many analysis use-cases. One disadvantage is that these containers, TTree in the HEP use-case, require to contain the same attributes for each entry/row (representing events), which can make extending the list of attributes very costly in storage, even if those are only required for a small subsample of events. Since the initial design, the ATLAS software framework features powerful navigational infrastructure to allow storing custom data extensions for subsamples of events in separate, but synchronized containers. This allows adding event augmentations to ATLAS standard data products (such as DAOD-PHYS or PHYSLITE) avoiding duplication of those core data products, while limiting their size increase. For this functionality, the framework does not rely on any associations made by the I/O technology (i.e. ROOT), however it supports TTree friends and builds the associated index to allow for analysis outside of the ATLAS framework. A prototype based on the Long-Lived Particle search is implemented and preliminary results with this prototype will be presented. At this point, augmented data are stored within the same file as the core data. Storing them in separate files will be investigated in future, as this could provide more flexibility, e.g. certain sites may only want a subset of several augmentations or augmentations can be archived to tape once their analysis is complete.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Multithreading ATLAS offline software: A retrospective

For Run 3, ATLAS redesigned its offline software, Athena, so that the main workflows run completely multithreaded. The resulting substantial reduction in the overall memory requirements allows for better use of machines with many cores. This note will discuss the performance achieved by the multithreaded reconstruction, the process of migrating the large ATLAS code base, and tools and techniques that were useful in debugging threading-related problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Accelerating science: The usage of commercial clouds in ATLAS Distributed Computing

The ATLAS experiment at CERN is one of the largest scientific machines built to date and will have ever growing computing needs as the Large Hadron Collider collects an increasingly larger volume of data over the next 20 years. ATLAS is conducting R&D projects on Amazon Web Services and Google Cloud as complementary resources for distributed computing, focusing on some of the key features of commercial clouds: lightweight operation, elasticity and availability of multiple chip architectures. The proof of concept phases have concluded with the cloud-native, vendoragnostic integration with the experiment’s data and workload management frameworks. Google Cloud has been used to evaluate elastic batch computing, ramping up ephemeral clusters of up to O(100k) cores to process tasks requiring quick turnaround. Amazon Web Services has been exploited for the successful physics validation of the Athena simulation software on ARM processors. We have also set up an interactive facility for physics analysis allowing endusers to spin up private, on-demand clusters for parallel computing with up to 4 000 cores, or run GPU enabled notebooks and jobs for machine learning applications. The success of the proof of concept phases has led to the extension of the Google Cloud project, where ATLAS will study the total cost of ownership of a production cloud site during 15 months with 10k cores on average, fully integrated with distributed grid computing resources and continue the R&D projects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Low temperature high frequency coaxial pulse tube for space application

The 4K stage is a critical step for space missions. The Hershel mission is using a helium bath, which is consumed day by day (after depletion, the space mission is over) while the Plank mission is equipped with one He4 Joule-Thomson cooler. Cryogenic chain without helium bath is a challenge for space missions and 4.2K Pulse-Tube working at high frequency (around 30Hz) is one option to take it up. A low temperature Pulse-Tube would be suitable for the ESA space mission EChO (Exoplanet Characterisation Observatory, expected launch in 2022), which requires around 30mW cooling power at 6K; and for the ESA space mission ATHENA (Advanced Telescope for High ENergy Astrophysics), to pre-cool the sub-kelvin cooler (few hundreds of mW at 15K). The test bench described in this paper combines a Gifford-McMahon with a coaxial Pulse-Tube. A thermal link is joining the intercept of the Pulse-Tube and the second stage of the Gifford-McMahon. This intercept is a separator between the hot and the cold regenerators of the Pulse-Tube. The work has been focused on the cold part of this cold finger. Coupled with an active phase shifter, this Pulse-Tube has been tested and optimized and temperatures as low as 6K have been obtained at 30Hz with an intercept temperature at 20K.

42 ENGINEERING↗
DOE OSTI2015

Study of ultrasonic characterization and propagation in austenitic welds: The MOSAICS project

Regulatory requirements enforce a volumetric inspection of welded components of nuclear equipments. However, the multi-pass austenitic welds are characterized by anisotropic and heterogeneous structures which lead to numerous disturbances of the ultrasonic beam. The MOSAICS project supported by the ANR (French National Research Agency) aims at matching various approaches to improve the prediction of the ultrasonic testing in those welds. The first stage consists in characterizing the weld structure (determination of the columnar grain orientation and measurements of elastic constants and attenuation coefficients). The techniques of characterization provide input data for the modeling codes developed in another task of the project. For example, a 3D version of the finite elements code ATHENA is developed by EDF R and D to take into account anisotropic texture in any direction. Semi-analytical models included in CIVA software are also improved to better predict the ultrasonic propagation in highly anisotropic and heterogeneous structures. The last stage deals with modeling codes validation based on experimental inspections on representative mock-ups containing calibrated defects. The objective of this paper is to give an overview of the MOSAICS project and to present specific results illustrating the various tasks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Application of x-ray absorption fine structure (XAFS) to local-order analysis in Fe-Cr maghemite-like materials

The maghemite-like oxide system γ-Fe{sub 2-x}Cr{sub x}O{sub 3} (x=0.75, 1 and 1.25) was studied by X-ray absorption fine structure (XAFS) and by synchrotron radiation X-ray diffraction (XRD). Measurements were performed at the Stanford Synchrotron Radiation Lightsource at room temperature, at beamlines 2-1, 2-3 and 4-3. High-resolution XRD patterns were processed by means of the Rietveld method. In cases of atoms being neighbors in the Periodic Table, the order/disorder degree of the considered solutions is indiscernible by “normal” (absence of “anomalous scattering”) diffraction experiments. Thus, maghemite-like materials were investigated by XAFS in both Fe and Cr K-edges to clarify, via short-range structure characterization, the local ordering of the investigated system. Athena and Artemis graphic user interfaces for IFEFFIT and FEFF8.4 codes were employed for XAFS spectra interpretation. Pre-edge decomposition and theoretical modeling of X-ray absorption near edge structure (XANES) transitions were performed. By analysis of the Cr K-edge XANES, it has been confirmed that Cr is located in an octahedral environment. Fitting of the extended X-ray absorption fine structure (EXAFS) spectra was performed under the consideration that the central atom of Fe is allowed to occupy octa- and tetrahedral positions, while Cr occupies only octahedral ones. Coordination number of neighboring atoms, interatomic distances and their quadratic deviation average were determined for x=1, by fitting simultaneously the EXAFS spectra of both Fe and Cr K-edges. The results of fitting the experimental spectra with theoretical standards showed that the cation vacancies tend to follow a regular pattern within the structure of the iron-chromium maghemite (FeCrO{sub 3})

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

NONLINEAR EVOLUTION OF RAYLEIGH-TAYLOR INSTABILITY IN A RADIATION-SUPPORTED ATMOSPHERE

The nonlinear regime of Rayleigh-Taylor instability (RTI) in a radiation supported atmosphere, consisting of two uniform fluids with different densities, is studied numerically. We perform simulations using our recently developed numerical algorithm for multi-dimensional radiation hydrodynamics based on a variable Eddington tensor (VET) as implemented in Athena, focusing on the regime where scattering opacity greatly exceeds absorption opacity. We find that the radiation field can reduce the growth and mixing rate of RTI, but this reduction is only significant when radiation pressure significantly exceeds gas pressure. Small-scale structures are also suppressed in this case. In the nonlinear regime, dense fingers sink faster than rarefied bubbles can rise, leading to asymmetric structures about the interface. By comparing the calculations that use a VET versus the Eddington approximation, we demonstrate that anisotropy in the radiation field can affect the nonlinear development of RTI significantly. We also examine the disruption of a shell of cold gas being accelerated by strong radiation pressure, motivated by models of radiation driven outflows in ultraluminous infrared galaxies. We find that when the growth timescale of RTI is smaller than acceleration timescale, the amount of gas that would be pushed away by the radiation field is reduced due to RTI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

TESTING CONVERGENCE FOR GLOBAL ACCRETION DISKS

Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by magneto-rotational instability (MRI). Using them to accurately predict quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: (Q{sub z} ) and (Q{sub {phi}}), the ratio of the grid zone size to the characteristic MRI wavelength, {alpha}{sub mag}, the ratio of the Maxwell stress to the magnetic pressure, and / , the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by (Q{sub z} ) {>=} 15, (Q{sub {phi}}) {>=} 20, {alpha}{sub mag} Almost-Equal-To 0.45, and / {approx}0.2. These values are associated with {>=}35 zones per scaleheight H, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the Harten-Lax-van Leer-Einfeldt flux solver or second-order interpolation can significantly degrade the effective resolution compared to the Harten-Lax-van Leer discontinuities flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation. Without it, however, quantitative measures important to predictions of observables are subject to large systematic errors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

TESTING THE NO-HAIR THEOREM WITH OBSERVATIONS IN THE ELECTROMAGNETIC SPECTRUM. IV. RELATIVISTICALLY BROADENED IRON LINES

According to the no-hair theorem, astrophysical black holes are fully characterized by their masses and spins and are described by the Kerr metric. This theorem can be tested observationally by measuring (at least) three different multipole moments of the spacetimes of black holes. In this paper, we calculate the profiles of fluorescent iron lines emitted from the accretion flows around black hole candidates within a framework that allows us to perform the calculation as a function of its mass and spin as well as of a free parameter that measures potential deviations from the Kerr metric. We show that such deviations lead to line profiles that are significantly altered and may exhibit a modified flux ratio of the two peaks in their characteristic double-peaked shape. We also show that the disk inclination can be measured independently of the spin and the deviation parameter at low to intermediate inclination angles, as in the case of Kerr black holes. We estimate the precision that near-future X-ray missions such as Astro-H and ATHENA+ are required to achieve in order to resolve deviations from the Kerr metric in iron line profiles and show that constraints on such deviations will be strongest for rapidly spinning black holes. More generally, we show that measuring the line profile with a precision of {approx}5% at disk inclinations of 30 Degree-Sign or 60 Degree-Sign constrains the deviation parameter to order unity for values of the spin a {approx}> 0.5M.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

ROLE OF MAGNETIC FIELD STRENGTH AND NUMERICAL RESOLUTION IN SIMULATIONS OF THE HEAT-FLUX-DRIVEN BUOYANCY INSTABILITY

The role played by magnetic fields in the intracluster medium (ICM) of galaxy clusters is complex. The weakly collisional nature of the ICM leads to thermal conduction that is channeled along field lines. This anisotropic heat conduction profoundly changes the instabilities of the ICM atmosphere, with convective stabilities being driven by temperature gradients of either sign. Here, we employ the Athena magnetohydrodynamic code to investigate the local non-linear behavior of the heat-flux-driven buoyancy instability (HBI) relevant in the cores of cooling-core clusters where the temperature increases with radius. We study a grid of two-dimensional simulations that span a large range of initial magnetic field strengths and numerical resolutions. For very weak initial fields, we recover the previously known result that the HBI wraps the field in the horizontal direction, thereby shutting off the heat flux. However, we find that simulations that begin with intermediate initial field strengths have a qualitatively different behavior, forming HBI-stable filaments that resist field-line wrapping and enable sustained vertical conductive heat flux at a level of 10%-25% of the Spitzer value. While astrophysical conclusions regarding the role of conduction in cooling cores require detailed global models, our local study proves that systems dominated by the HBI do not necessarily quench the conductive heat flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

On the thermal stability of radiation-dominated accretion disks

We study the long-term thermal stability of radiation-dominated disks in which the vertical structure is determined self-consistently by the balance of heating due to the dissipation of MHD turbulence driven by magneto-rotational instability (MRI) and cooling due to radiation emitted at the photosphere. The calculations adopt the local shearing box approximation and utilize the recently developed radiation transfer module in the Athena MHD code based on a variable Eddington tensor rather than an assumed local closure. After saturation of the MRI, in many cases the disk maintains a steady vertical structure for many thermal times. However, in every case in which the box size in the horizontal directions are at least one pressure scale height, fluctuations associated with MRI turbulence and dynamo action in the disk eventually trigger a thermal runaway that causes the disk to either expand or contract until the calculation must be terminated. During runaway, the dependence of the heating and cooling rates on total pressure satisfy the simplest criterion for classical thermal instability. We identify several physical reasons why the thermal runaway observed in our simulations differ from the standard α disk model; for example, the advection of radiation contributes a non-negligible fraction to the vertical energy flux at the largest radiation pressure, most of the dissipation does not happen in the disk mid-plane, and the change of dissipation scale height with mid-plane pressure is slower than the change of density scale height. We discuss how and why our results differ from those published previously. Such thermal runaway behavior might have important implications for interpreting temporal variability in observed systems, but fully global simulations are required to study the saturated state before detailed predictions can be made.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

An analytic model for buoyancy resonances in protoplanetary disks

Zhu et al. found in three-dimensional shearing box simulations a new form of planet-disk interaction that they attributed to a vertical buoyancy resonance in the disk. We describe an analytic linear model for this interaction. We adopt a simplified model involving azimuthal forcing that produces the resonance and permits an analytic description of its structure. We derive an analytic expression for the buoyancy torque and show that the vertical torque distribution agrees well with the results of the Athena simulations and a Fourier method for linear numerical calculations carried out with the same forcing. The buoyancy resonance differs from the classic Lindblad and corotation resonances in that the resonance lies along tilted planes. Its width depends on damping effects and is independent of the gas sound speed. The resonance does not excite propagating waves. At a given large azimuthal wavenumber k{sub y} > h {sup –1} (for disk thickness h), the buoyancy resonance exerts a torque over a region that lies radially closer to the corotation radius than the Lindblad resonance. Because the torque is localized to the region of excitation, it is potentially subject to the effects of nonlinear saturation. In addition, the torque can be reduced by the effects of radiative heat transfer between the resonant region and its surroundings. For each azimuthal wavenumber, the resonance establishes a large scale density wave pattern in a plane within the disk.

79 ASTRONOMY AND ASTROPHYSICS↗