Search NASA⌕ Search

SEARCH · Search NASA

Results for “MCNP Model”

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 OSTI2017

Calculation of Activated Water Radiation Maps in ITER

ITER is a nuclear machine designed for 500 MW DT fusion power emitting an immense power in the form of 14 MeV neutrons (∼400 MW). These high energy neutrons activate the water used to cool the first wall, blanket modules, vacuum vessel etc. The principal isotopes generated in water coolant are {sup 16}N and {sup 17}'N from oxygen isotopes. Both these radioisotopes are generated by neutrons having energy more than 10 MeV. Even though these isotopes are short-lived (T{sub 1/2} {sup 16}N: 7.13 s and T{sub 1/2} {sup 17}N: 4.17 s) they generate significant radiation fields outside the bio-shield due to high neutron fluence at irradiation end, large volume of the coolant, high speed of the coolant, emission of high energy nuclear radiations from {sup 16}N (gammas) and {sup 17}N (neutrons). Therefore, radiation fields generated by the activated water need to be calculated in order to design the plant systems to withstand radiation dosage and for the safety of the workers as per regulatory limits. This paper presents the calculation methodology of the {sup 16}N and {sup 17}N sources in the coolant lines, manifolds, heat exchangers etc., which are distributed over the wide area outside the bio-shield. The calculated source is modelled in the MCNP building model which is approximately 120 m x 80 m x 50 m. The building model has significant amount of shielding walls (∼3 m thickness). This presents a challenge to get the sufficient statistics everywhere in the global model. The radiation transport is performed with the code MCNP to calculate the radiation fields (biological dose rate, absorbed dose in silicon, 1 MeV silicon equivalent fluence, neutron flux) throughout the tokamak complex.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

NCERC-CL TREAT Free Field Characterization Foil Reaction Rate Results

This report summarizes the analysis and preliminary MCNP reaction rate modeling of free field characterization measurements performed in the Transient Reaction Test Facility (TREAT) at Idaho National Laboratory (INL) between 12/18/2023 and 2/22/2024. A series of seven irradiations were performed using the Big-BUSTER (Broad Use Specimen Transient Experiment Rig) core configuration. Foil sets including 19.5% LEU-Zr alloy wire, S, Au, Fe, Ni, Co, Ti, and Zr were deployed in each irradiation. Identical foil sets were supplied to Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL), and INL. All foils were supplied by INL. Additional details on the experiment can be found in [2]. The foil sets were shipped to the National Criticality Experiments Research Center Counting Laboratory (NCERC-CL) branch at LANL (NCERC-CL NISC) and were received on 3/25/2024. Additional details on receipt and counting are reported in [3]. All reactor dosimetry measurements were performed in adherence to the ASTM standards applicable to reactor dosimetry. Measured reaction rates from Au, Co, Fe, Ni, and Ti are reported. Select reaction rate ratios and simulated reaction rates and ratios are discussed. Reaction rates for all measured reaction products were modeled with MCNP 6.3.1 and a TREAT input deck supplied by Edward Lum. Neutron emission estimates calculated using the measured and simulated reaction rates show excellent agreement between the three fast-threshold n-p reactions with the percent differences being < 9% between the three reactions. The capture reactions had poorer, but still reasonable, agreement with < 20% percent differences between the neutron emission estimates of the three capture reactions. The percent difference between the fast-threshold and capture reactions is ∼ 100% for each reaction. This supports the TREAT model is not accurately modeling the neutron spectrum and additional measurements are required to characterize TREAT and match measurements to simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2021

Reactor cell neutron dose for the molten salt breeder reactor conceptual design

The private sector’s interest in the active development of molten salt reactors has led to the need to develop and test advanced modeling and simulation tools to analyze various advanced reactor types under numerous conditions. This paper discusses the effort undertaken to model the Oak Ridge National Laboratory (ORNL) Molten Salt Breeder Reactor (MSBR) design using ORNL’s Shift Monte Carlo code. The MSBR model integrates a Monte Carlo N-Particle (MCNP) MSBR core model with an MCNP model that was generated from a CAD model of the external components and the reactor building, which was subsequently run in Shift. This paper focuses on development of the fully integrated model and its use in performing neutron transport calculations in the reactor cell area. This model is intended to aid in understanding radiological dose conditions during operation, as well as the iron dpa rates in the reactor vessel. The neutron biological dose rates and flux calculated in the reactor cell are much higher in the MSBR than in typical light-water reactors. The implications of these results and future work are also discussed in this paper.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Validation of the Fast and Thermal Neutron Flux Profiles in the Advanced Test Reactor for the AGR-3/4 Experiment

The AGR-3/4 experiment was designed to irradiate TRISO particle fuel in the Advanced Test Reactor (ATR) at the Idaho National Laboratory. Temperature control of the test capsules containing the TRISO particles was important since the goal was to measure diffusion rates of fission product elements through graphite as a function of temperature. During the experiment, differences were observed between temperatures measured by the thermocouples and temperatures calculated with the thermal analysis model. To rule out the possibility that the used ATR driver fuel elements were the root cause of these temperature differences, an experiment was devised to measure the axial fast and thermal neutron flux profiles in the northeast flux trap containing the AGR-3/4 experiment. An important input parameter to the thermal analysis model was the calculated compact heat rates from the MCNP physics model. To verify these calculated heat rates were consistent with the ATR driver fuel element loading and operating lobe power, axial flux wires were inserted into ATR northeast flux trap. The measured fast and thermal neutron axial flux profiles derived from the flux wire activation showed excellent agreement with the calculated flux profiles. This agreement dispelled any concerns that the used ATR fuel elements in the northeast lobe might have introduced unusual axial flux variation. In addition, the agreement provided a validation basis for the physics calculation, specifically, (1) MCNP physics model, (2) Monte Carlo depletion methodology, (3) method used to load the ATR driver fuel, and (4) ATR measured lobe power data used in the calculation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

NEET Enhanced Micro-Pocket Fission Detector for High Temperature Reactors - FY16 Status Report

A collaboration between the Idaho National Laboratory (INL), the Kansas State University (KSU), and the French Atomic Energy Agency, Commissariat à l'Énergie Atomique et aux Energies Alternatives, (CEA), has been initiated by the Nuclear Energy Enabling Technologies (NEET) Advanced Sensors and Instrumentation (ASI) program for developing and testing High Temperature Micro-Pocket Fission Detectors (HT MPFD), which are compact fission chambers capable of simultaneously measuring thermal neutron flux, fast neutron flux and temperature within a single package for temperatures up to 800 °C. The MPFD technology utilizes a small, multi-purpose, robust, in-core fission chambers and thermocouple. As discussed within this report, the small size, variable sensitivity, and increased accuracy of the MPFD technology represent a revolutionary improvement over current methods used to support irradiations in US Material Test Reactors (MTRs). Previous research conducted through NEET ASI1-3 has shown that the MPFD technology could be made robust and was successfully tested in a reactor core. This new project will further the MPFD technology for higher temperature regimes and other reactor applications by developing a HT MPFD suitable for temperatures up to 800 °C. This report summarizes the research progress for year two of this three year project. Highlights from research accomplishments include: • Continuation of a joint collaboration between INL, KSU, and CEA. Note that CEA is participating at their own expense because of interest in this unique new sensor. • An updated parallel wire HT MPFD design was developed. • Program support for HT MPFD deployments was given to Accident Tolerant Fuels (ATF) and Advanced Gas-cooled Reactor (AGR) irradiation test programs. • Quality approved materials for HT MPFD construction were procured by irradiation test programs for upcoming deployments. • KSU improved and performed electrical contact and fissile material plating. • KSU delivered fissile HT MPFD parts to INL for final construction of HT MPFD prototype. • A prototype HT MPFD was constructed and analyzed at INL. • The HT MPFD has been modeled in MCNP to optimize the amount of fissile material deposition. • The HT MPFD has been modeled in MCNP to optimize the sensor location in the irradiation test. • The fissile material deposition is undergoing independent verifications. • Detector amplifier electronics have been revised and tested by KSU. • Several project meetings were held at INL and KSU to discuss the roles and responsibilities between INL, KSU, and CEA for development and deployment of the HT MPFDs. As documented in this report, FY16 funding has allowed the project to meet year two planned accomplishments to develop a HT MPFD. In addition, the accomplishments of this project have attracted independent funding from other Department of Energy Office of Nuclear Energy (DOE-NE) programs for MTR irradiations of the MPFD technology. These are significant opportunities for this NEET Enhanced Micro-Pocket Fission Detector for High Temperature Reactors project because the irradiation expense of these experiments could not be included in the original project scope.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

SPND Sensitivity Calculations Using MCNP and Experimental Data from ACRR

The use of the Monte Carlo N-Particle Transport Code (MCNP) to calculate detector sensitivity for Self-Powered Neutron Detectors (SPNDs) in the Annular Core Research Reactor (ACRR) could be a vital tool in the effort to optimize the design of next-generation SPNDs. Next-generation SPND designs, which consider specific materials and geometry, may provide experimenters with capabilities for advanced mixed field dosimetry. These detectors will need to be optimized for configuration, materials, and geometries and the ability to model and iterate must be available in order to decide on the ideal. SPND design. SPNDs were modeled in MCNP which closely resembled the dimensions and location of actual detectors used in the ACRR. Tallies were used to calculate detector sensitivity. Using metrics from a previous report, oscilloscope data from pulses were manipulated in a Matrix Laboratory computing environment (MATLAB) script to calculate experimental detector sensitivity. This report outlines the process in which experimental data from ACRR pulses verified results from tallies in an MCNP ACRR model. The sensitivity values from experiments and MCNP calculations agreed within one standard deviation. Parametric studies were also performed with MCNP to investigate the effects of materials and dimensions of different SPNDs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Characterization of Filters Loaded With Reactor Strontium Carbonate - 13203

A collection of three highly radioactive filters containing reactor strontium carbonate were being prepared for disposal. All three filters were approximately characterized at the time of manufacture by gravimetric methods. The first filter had been partially emptied, and the quantity of residual activity was uncertain. Dose rate to activity modeling using the Monte-Carlo N Particle (MCNP) code was selected to confirm the gravimetric characterization of the full filters, and to fully characterize the partially emptied filter. Although dose rate to activity modeling using MCNP is a common technique, it is not often used for Bremsstrahlung-dominant materials such as reactor strontium. As a result, different MCNP modeling options were compared to determine the optimum approach. This comparison indicated that the accuracy of the results were heavily dependent on the MCNP modeling details and the location of the dose rate measurement point. The optimum model utilized a photon spectrum generated by the Oak Ridge Isotope Generation and Depletion (ORIGEN) code and dose rates measured at 30 cm. Results from the optimum model agreed with the gravimetric estimates within 15%. It was demonstrated that dose rate to activity modeling can be successful for Bremsstrahlung-dominant radioactive materials. However, the degree of success is heavily dependent on the choice of modeling techniques. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Validation of nuclear data based on the ASPIS experiments from the SINBAD database - 25428

SINBAD, the Shielding Integral Benchmark Archive and Database, includes the experimental data relative to integral benchmark experiments relevant for radiation shielding applications. The ASPIS experiments, included in the Reactor Shielding section of SINBAD, provide neutron activation foils and spectra measurements performed before/between shielding slabs of Graphite, Sodium, Iron, Steels, Water and combinations of them. Benchmark models for the measurement simulations with the MCNP code are usually provided in SINBAD. Nonetheless, a user-friendly tool is presented, which would be useful for an independent production of MCNP benchmark models. It is based on GPL C++ code called CombLayer. The reports issued for SINBAD on the quality assessment of the ASPIS experiments are also described. The discussion on the quality of the ASPIS experiments for their use in nuclear data evaluations had a positive outcome, except for two experiments out of ten. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Processing MCNP Elemental Edit Outputs

The Los Alamos National Laboratory's (LANL) Monte Carlo N-Particle® (MCNP®) transport code version 6 has the capability for tracking particles on unstructured mesh (UM) geometry models. The MCNP UM feature was originally designed for the UM models created by the Abaqus/CAE software product. The MCNP version 6.2.0 and later recognizes the UM models read from Abaqus INP files or MCNPUM files converted from Abaqus INP files. Other finite element analysis software packages may generate UM models and then convert these models into Abaqus INP formats.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Generating MCNP Input Files for Unstructured Mesh Geometries

Los Alamos National Laboratory's (LANL) Monte Carlo N-Particle (MCNP) transport code version 6 has the capability for tracking particles on unstructured mesh (UM) geometry models. The MCNP UM feature has been developed for performing calculations of complex geometry models. This capability tracks particles on hybrid geometries where finite element meshes are embedded into constructive solid geometry (CSG) cells. The MCNP UM feature was originally designed to read UM models created by Abaqus/CAE software suite. MCNP versions 6.2.0 and later can process UM models read from Abaqus input les or MCNPUM les converted from Abaqus input les. Sandia National Laboratory's Cubit Toolkit and other finite element analysis software packages may generate UM models and then convert these models into Abaqus input formats.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Validating Intra-Plate Power Distributions in the Advanced Test Reactor

We report that the neutronics software, HELIOS, was validated in 2015 for performing core reload design and safety analysis of the Advanced Test Reactor. However, when HELIOS was benchmarked against historic fission-wire measurements (i.e., zero-power full-core measurements) a statistically resolved calculation-to-measurement bias was discovered. The azimuthal power along each fuel plate computed by HELIOS has consistently shown to under-predict measurements made by fission-wires in historic zero-power tests near the fuel element side-plates. It was hypothesized during the HELIOS software validation work, that this bias is attributable to local moderation in coolant vents in the side-plates axially just above and below the fission-wires on the fuel-plate edges. This work used detailed MCNP and MC21 models of the side-plate vents to test this hypothesis. By comparing the average azimuthal biases between HELIOS, and 2D and 3D MCNP models, and a 3D MC21 model, it was found that the HELIOS azimuthal bias is not due to the measurement.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Integration of the ITER Tokamak C-model with the Building Model

The pellet injector system is designed to accommodate the requirement of continuous plasma refueling during ITER Tokamak operations. Through pellet flight tubes/injection lines, the pellets can directly reach the interior of the tokamak vacuum vessel. The pellet injector system is located at a port cell outside of the Tokamak bio-shield. The plasma radiation source in the tokamak vacuum vessel, together with the highly activated surrounding water in port cells, present major radiation challenges for the pellet injector system during ITER operation and after shutdown. These sources produce a radiation field encompassing the entire Tokamak complex. The ITER neutronics reference model C-model is an ITER radiation transport model. It represents a 40 deg. regular torus sector out to the bio-shield, which consists of 38 individual systems. The Monte Carlo N-Particle (MCNP) input of the C-model consists of more than a hundred thousand cell definitions and surface definitions, and over three hundred material compositions. The latest tokamak complex MCNP model is significantly more detailed than any previous releases. It consists of 36,863 cell definitions, 57,105 surface definitions, and a total of 21 different types of materials. The C-model and tokamak building model share the bio-shield as the boundary between them. In general, the radiation fields in the building due to the plasma source during tokamak operation are determined by generating the radiation source in the inner faces of the bio-shield. An inherent limitation of this approach is that it can only approximate the directions in which streaming paths are produced through the bio-shield plugs. To better estimate the shutdown dose rate for the pellet injector system due to the plasma source, we have decided to integrate the tokamak model (C-model) with the complex building model. Both MCNP models are provided by the ITER International Organization (IO) to the Oak Ridge National Laboratory (ORNL) Radiation Transport (RT) group.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Cubit for MCNP Unstructured Mesh Analysis of Oktavian Benchmarks

The Monte Carlo N-Particle (MCNP) transport code developed by Los Alamos National Laboratory (LANL) can be used to transport various particles across user defined three dimensional (3D) geometries. Traditionally, these geometries are created as constructive solid geometry (CSG), involving the use of Boolean operators on defined surfaces to create 3D regions known as cells. A newer method of geometry definition in the MCNP code is unstructured mesh (UM) embedded within a CSG cell using the universe and fill repeated-structure features. An MCNP UM calculation requires an MCNP input file and mesh geometry file. The MCNP code cannot be used to generate UM geometry models. A computer-aided design (CAD) software is typically used to construct a solid geometry model, and a CAD file is then imported into a mesh generation software to prepare and generate a mesh model. Some mesh generation software packages have the ability to create a solid geometry and thus a separate CAD software for creating a CAD model is not needed. In this work, Cubit, a geometry creation and meshing software developed by Sandia National Laboratories, is used to construct UM models for MCNP simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Optimization of the moderators in the STS preliminary design

This report details the results for an optimization of the dimensions of the moderators in the preliminary design of the Spallation Neutron Source Second Target Station (STS). This study uses the optimization algorithms of Dakota and an unstructured mesh model for the moderators in MCNP. More details on the unstructured mesh model and the automated mesh generation can be found in [3]. Parallel to this effort, the same moderator geometries have been optimized using a constructive solid geometry (CSG) MCNP model. More details on this model and its results can be found in [4]. Three optimal designs are selected for each moderator: one that is optimized for maximum peak brightness, one for maximum time-integrated brightness, and one for a combination of peak and time-integrated brightness. The backbone of the optimization work flow is provided by Dakota. For each set of design parameters requested by Dakota, a new solid geometry is automatically built in Creo and SpaceClaim, and subsequently exported to Attila4MC to generate an unstructured mesh geometry for MCNP. After the MCNP calculation is finished, the objective function (e.g., brightness metric) is returned to Dakota. After the new design has been evaluated, a result-file is written, and Dakota proposes the next set of design parameters to be evaluated. The loop continues until a specified convergence criterion has been met. The design parameters of the cylindrical (upper) moderator include the hydrogen radius, the premoderator thickness (top, bottom, radial), the beryllium radius and the horizontal position of the moderator. The crucial design choice is the hydrogen radius. A radius of 62 mm is shown to provide the maximum time-integrated brightness. The maximum peak brightness occurs with a radius of 40 mm. A combined (middle) design, which balances peak and time-integrated brightnesses, is obtained with a hydrogen radius of 50 mm. The premoderator thicknesses and the beryllium radius are slightly larger in the design optimized for time-integrated brightness than in the design optimized for peak brightness. The sensitivity to these two parameters is relatively small close to the optimal configurations. The hydrogen vessel and vacuum vessel wall thicknesses are dependent on the radius of the liquid hydrogen due to structural integrity requirements. The increased wall thicknesses for larger vessels significantly penalize the time-integrated brightness, with the maximum obtainable value reduced by more than 10% relative to earlier studies which used fixed vessel wall thicknesses. The impact of the variable wall thicknesses is much less for the peak brightness and combined brightness designs. The design parameters of the tube (lower) moderator selected for the optimization are the tube length, the annular premoderator thickness, the beryllium radius and the horizontal position of the moderator. The tube length is the crucial parameter and is chosen large (210 mm) and small (125 mm) in the designs optimized for time-integrated and peak brightness respectively. A combined optimal design has a tube length of 170 mm. The premoderator thickness and the beryllium radius are chosen larger in the design optimized for time-integrated brightness.

42 ENGINEERING↗
DOE OSTI2021

Analysis of Approximations in Modeling of BWR Bundle Void Distributions

In boiling water reactors, complex heterogeneous bundle designs, control blades adjacent to the corner of bundles, and the presence of boiling can lead to complex internal void distributions. A few approximations exist to model these void distributions. They could be modeled using a 1D axial solver in which each axial node is assumed to be at an average void, or each pin cell could be modeled with its own void concentration. In the latter case, the void could be discretized in pin-centered or coolant-centered channels. The goal of this project was to quantify the effect of using the different approximations for modeling internal void distributions on neutronics calculations. Using 3D void distributions calculated with CTF, Monte Carlo Neutral Particle (MCNP) transport code models were created for GE-9 and GE-14 lattices. For each model, the internal void distribution from CTF at a given axial node was selected, and a lattice calculation was carried out with MCNP. Comparisons between models using a lattice-averaged void, or using a void distribution in coolant-centered channels, showed large differences in reactivity which in some cases were well above 1,000 pcm, and it also showed differences in normalized fission rates greater than 20%. It was also found that using a lattice average void can lead to a significant difference in the estimation of the worth of a control blade. The differences found when comparing results from models using pin-centered and coolant-centered channels were up to 200 pcm in reactivity and up to 1.4% in the normalized fission rates. In addition to these two sets of comparisons, MCNP models were set up so that each subchannel had a saturated liquid component around the fuel pins and a saturated vapor component in the center to approximate annular flow. In comparison to the models using coolant-centered subchannels, up to 1–3% differences in normalized fission rates could be found.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Methods for Estimating Dose from Galactic Cosmic Rays and Jovian Radiation Belts on Reduced-Order Monte Carlo Geometries

This document outlines the methods and tools developed to estimate dose rates from galactic cosmic rays and the trapped radiation belts surrounding Jupiter using Monte Carlo particle transport software. Furthermore, the techniques described here are specific for translating high-resolution MAVRIC models to analogous but lower resolution Monte Carlo N-Particle Transport (MCNP) code models to estimate charged particle effects on models originally developed to study neutral particles only. This capability was initially developed to predict doses to various components of dynamic radioisotope power systems for spaceflight applications, but the tools and methods can be easily applied to other problem geometries and expanded to other space radiation environments.

79 ASTRONOMY AND ASTROPHYSICS↗