SEARCH · Engineering Papers

Results for “Monte Carlo Code”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2013

Obtaining identical results with double precision global accuracy on different numbers of processors in parallel particle Monte Carlo simulations

We describe and compare different approaches for achieving numerical reproducibility in photon Monte Carlo simulations. Reproducibility is desirable for code verification, testing, and debugging. Parallelism creates a unique problem for achieving reproducibility in Monte Carlo simulations because it changes the order in which values are summed. This is a numerical problem because double precision arithmetic is not associative. Parallel Monte Carlo, both domain replicated and decomposed simulations, will run their particles in a different order during different runs of the same simulation because the non-reproducibility of communication between processors. In addition, runs of the same simulation using different domain decompositions will also result in particles being simulated in a different order. In [1], a way of eliminating non-associative accumulations using integer tallies was described. This approach successfully achieves reproducibility at the cost of lost accuracy by rounding double precision numbers to fewer significant digits. This integer approach, and other extended and reduced precision reproducibility techniques, are described and compared in this work. Increased precision alone is not enough to ensure reproducibility of photon Monte Carlo simulations. Non-arbitrary precision approaches require a varying degree of rounding to achieve reproducibility. For the problems investigated in this work double precision global accuracy was achievable by using 100 bits of precision or greater on all unordered sums which where subsequently rounded to double precision at the end of every time-step.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2012

Development of array-type prompt gamma measurement system for in vivo range verification in proton therapy

Purpose: In vivo range verification is one of the most important parts of proton therapy to fully utilize its benefits delivering high radiation dose to tumor, while sparing the normal tissue with the so-called Bragg peak. Currently, however, range verification method is not used in clinics. The purpose of the present study is to optimize and evaluate the configuration of an array-type prompt gamma measurement system on determining distal dose edge for in vivo range verification of proton therapy. Methods: To effectively measure the prompt gammas against the background gammas, the Monte Carlo simulations with the MCNPX code were employed in optimizing the configuration of the measurement system, and the Monte Carlo method was also used to understand the effect of the background gammas, mainly neutron capture gammas, in the measured gamma distribution. To reduce the effect of the background gammas, the optimized energy window of 4-10 MeV in measuring the prompt gammas was employed. A parameterized source was used to maximize computation speed in the optimization study. A simplified test measurement system, using only one detector moving from one measurement location to the next, was constructed and applied to therapeutic proton beams of 80-220 MeV. For accurate determination of the distal dose edge, the sigmoidal curve-fitting method was applied to the measured distributions of the prompt gammas, and then, the location of the half-value between the maximum and minimum value in the curve-fitting was determined as the distal dose edge and compared with the beam range assessed by the proton dose distribution. Results: The parameterized source term employed in optimization process improved the calculation speed by up to {approx}300 times. The optimization study indicates that an array-type measurement system with 3, 2, 2, and 150 mm for scintillator thickness, slit width, septal thickness, and slit length, respectively, can effectively measure the prompt gamma distributions minimizing the contribution of background gammas. The present results show that a few hundred counts of prompt gammas can be easily obtained by measuring 10 s at each measurement location for proton beams of {approx}4 nA. The distal dose edges determined by the prompt gamma distribution are 5.45, 14.73, and 27.74 cm for proton beams of 5.17 (80 MeV), 14.99 (150 MeV), and 27.38 (220 MeV) cm, respectively. Conclusions: The results show that the array-type measurement system can measure prompt gamma distributions from a therapeutic proton beam within a short measurement time, and that the distal dose edge can be determined within a few millimeters of error without using any sophisticated analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2019

About the Impact of the Unresolved Resonance Region in Monte Carlo Simulations of Sodium Fast Reactors

In the last few years, and within the framework of different European projects, KENO-VI code from SCALE system has been employed to perform detailed continuous-energy Mon-te Carlo transport calculations for advanced fast reactors. The core characterization of both the sodium-cooled ASTRID and the lead-cooled ALFRED reactors was performed during the FP7 cross-cutting ESNII+ project; more recently, core calculations for the sodium-cooled Superphénix reactor and the improved European Sodium Fast Reactor design were performed within the HORIZON2020 ESFR-SMART project. In all cases, the effective multiplication factor predicted by KENO-VI was systematically higher (around 400-500 pcm) than the values computed by MCNP and Serpent Monte Carlo codes, using the same nu-clear data library.In order to provide insight into the origin of the observed discrepancies, a simplified 2D MOX-fueled SFR pin-cell benchmark has been launched. The multiplication factor, as well as 1-group and VITAMINJ 175-group cross-sections computed by KENO-VI, Serpent and MCNP codes employing ENDF/B-VII.1 data library, have been compared.Significant differences between KENO-VI and the other codes have been found in the un-resolved resonance regions of 239Pu and 241Pu capture and production cross sections, while negligible differences appeared outside those energy ranges. On the other hand, calculations without using probability tables have shown very good agreement. Quantita-tive comparison is presented and analyzed, along with a discussion of the impact of the probability-table treatment in the three codes for MOX-fueled systems with typical SFR spectrum.

Jimenez carrascosa, A↗

DOE OSTI2018

An angular biasing method using arbitrary convex polyhedra for Monte Carlo radiation transport calculations

This paper presents a new method for performing angular biasing in Monte Carlo radiation transport codes using arbitrary convex polyhedra to define regions of interest toward which to project particles (DXTRAN regions). The method is derived and is implemented using axis-aligned right parallelepipeds (AARPPs) and arbitrary convex polyhedra. Attention is also paid to possible numerical complications and areas for future refinement. A series of test problems are executed with void, purely absorbing, purely scattering, and 50% absorbing/50% scattering materials. For all test problems tally results using AARPP and polyhedral DXTRAN regions agree with analog and/or spherical DXTRAN results within statistical uncertainties. In cases with significant scattering the figure of merit (FOM) using AARPP or polyhedral DXTRAN regions is lower than with spherical regions despite the ability to closely fit the tally region. This is because spherical DXTRAN processing is computationally less expensive than AARPP or polyhedral DXTRAN processing. Thus, it is recommended that the speed of spherical regions be considered versus the ability to closely fit the tally region with an AARPP or arbitrary polyhedral region. It is also recommended that short calculations be made prior to final calculations to compare the FOM for the various DXTRAN geometries because of the influence of the scattering behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Criticality concerns of a group actinide co-crystallization separations approach to used nuclear fuel recycling

We report the co-crystallization of hexavalent actinides with uranyl nitrate hexahydrate has been established as a viable approach to achieving group actinide separation in used nuclear fuel recycle while avoiding the difficulties associated with the solvent-extraction process altogether. A preliminary neutronics safety analysis has been performed to alleviate criticality concerns associated with this group hexavalent actinide co-crystallization approach. The basic uranyl nitrate hexahydrate material compositions used are representative of used fuel found at commercial pressurized water reactors all around the U.S. This material is subsequently doped with various percentages of either typical reactor grade plutonium or with the total transuranic vector, replacing the uranium atoms in the crystal lattice structure. Infinite and effective neutron multiplication factors were calculated using the state-of-the-art Monte Carlo neutron transport code MCNP in a variety of plausible material configurations, assuming an acceptable safety limit for criticality to be an effective neutron multiplication factor less than 0.95. The results indicate that the crystal structures form a very subcritical intrinsic geometry even with six water molecules being incorporated into the crystal structure. In conclusion, the results of this research suggest that significant opportunity exists for the safe transfer of the group actinide separation approach from a laboratory to an operational scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

Measurement of leakage neutron spectra for zirconium with D-T neutrons and validation of evaluated nuclear data

Integral neutronics experiments have been a powerful tool and method for examining the reliability of evaluated nuclear data. In the present paper, the accuracy of evaluated nuclear data for Zirconium has been validated by comparing measured leakage neutron spectra with calculated ones. Leakage neutron spectra from the irradiation of D-T neutrons on Zirconium samples were experimentally measured at 60° and 120° by using a time-of-flight method. Theoretical calculations are carried out by Monte Carlo neutron transport code MCNP-4C with evaluated nuclear data of the ENDF/B-VII.0, ENDF/B-VIII.0, JENDL-4.0 and CENDL-3.1 libraries. From the comparisons, it is found that the calculations with JENDL-4.0 gives better agreements with the experiments at 60° while a large discrepancies are observed in the 4–13.5 MeV energy range at 120°. However, the calculated spectra using the data from ENDF/B-VII.0, ENDF/B-VIII.0 and CENDL-3.1 libraries showed some discrepancies with the measured ones, especially caused by the difference of the elastic and inelastic angular distributions and the secondary neutron energy distributions in the continuum inelastic scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗

DOE OSTI2013

Excited State Electronic Properties of Sodium Iodide and Cesium Iodide

We compute from first principles the dielectric function, loss function, lifetime and scattering rate of quasiparticles due to electronic losses, and secondary particle spectrum due to plasmon decay in two scintillating alkali halides, sodium iodide and cesium iodide. Particular emphasis is placed on quasiparticles within several multiples of the band gap from the band edges. A theory for the decay spectra of plasmons and other electronic excitations in crystals is presented. Applications to Monte Carlo radiation transport codes are discussed.

Campbell, Luke W↗

DOE OSTI2014

Reconstruction of NSTX midplane neutral density profiles from visible imaging data

In this report the experimental determination of neutral densities in tokamak plasmas from line radiation is only accurate in the narrow region in which both the excitation rate and neutral density are significant. We describe an alternative procedure using the DEGAS 2 Monte Carlo neutral transport code to invert light emission data obtained from a tangentially viewing camera, yielding absolute radial profiles of deuterium atoms and molecules at midplane. That the neutral source in these simulations can be adequately characterized as a uniform flux at the vacuum vessel wall is demonstrated by the similarity of the shapes of the simulated and observed brightness profiles. A second test is obtained by comparing the resulting neutral pressures at the vessel walls with data from midplane micro-ion gauges. We also show that the simulated camera image is insensitive to variations in the spatial distribution of the neutral source.

(PSI-21) NSTX DEGAS↗

DOE OSTI2026

Surrogate modeling of Monte Carlo radiation transport with convolutional neural networks for shielding optimization

Here, we present a machine learning (ML)-based surrogate model using convolutional neural networks (CNN) designed to emulate the attenuation of neutron fields as they pass through various shielding materials. This model can compute the outgoing neutron flux almost instantaneously and achieves reasonable accuracy compared to traditional Monte Carlo (MC)-based codes, which are computationally intensive. This emulator alleviates the complexity of neutron radiation transport through shielding materials by reducing the dimensionality and enables shielding optimization for a known radiation environment. This optimization process, which would have taken an unrealistic timeline due to several complex radiation transport simulations, can now be achieved in minutes, thus increasing computational capabilities in radiation shielding assessment. We demonstrate the applications of this emulator in computing effective dose rates and optimizing shielding solutions for a heavy-ion accelerator facility, such as the Facility for Rare Isotope Beams, where secondary neutrons produced via beam interactions dominate the radiation environment.

accelerator shielding↗

DOE OSTI2020

Assessment of CASL VERA for BWR analysis and application to accident tolerant SiC/SiC channel box

Application of the Virtual Environment for Reactor Applications (VERA) to BWR analysis is assessed in this paper by comparing results to those calculated using other widely-used modeling tools, namely the U.S. Nuclear Regulatory Commission’s PARCS/PATHS and the Serpent Monte Carlo particle transport code. Additionally, VERA is used to calculate 3-D temperature and fast neutron flux distributions in silicon carbide (SiC) fiber-reinforced, SiC matrix composite (SiC/SiC) BWR channel boxes, which are being studied as an Accident Tolerant Fuel core structural material concept. We found that the code-to-code comparisons were favorable, and the SiC/SiC channel box evaluation demonstrates the many advanced modeling features VERA provides while also highlighting the non-uniformity in fast neutron flux distributions that can play a role in potential SiC/SiC channel box deformation. Traditional BWR analysis tools do not have the calculation fidelity necessary for coupled assessment of flux and temperature gradients in a SiC/SiC channel box. VERA is a state-of-the-art modeling environment that was developed to increase the safety and economic competitiveness of nuclear power through improved modeling accuracy. While VERA has already been deployed in the nuclear industry for PWR applications, the current study is a vital initial step in the extensive development, validation, and verification that VERA must go through to be useful for BWR applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

SN 2021adxl: A luminous nearby interacting supernova in an extremely low-metallicity environment

SN 2021adxl is a slowly evolving, luminous, Type IIn supernova with asymmetric emission line profiles, similar to the well-studied SN 2010jl. We present extensive optical, near-ultraviolet, and near-infrared photometry and spectroscopy covering ∼1.5 years post discovery. SN 2021adxl occurred in an unusual environment, atop a vigorously star-forming region that is offset from its host galaxy core. The appearance of Ly α and O II , as well as the compact core, would classify the host of SN 2021adxl as a “Blueberry” galaxy, analogous to higher redshift, low-metallicity, star-forming dwarf “Green Pea” galaxies. Using several abundance indicators, we find a metallicity of the explosion environment of only ∼0.1 Z ⊙ , the lowest reported metallicity for a Type IIn SN environment. SN 2021adxl reaches a peak magnitude of M r ≈ −20.2 mag and since discovery, SN 2021adxl has faded by only ∼4 magnitudes in the r band with a cumulative radiated energy of ∼1.5 × 10 50 erg over 18 months. SN 2021adxl shows strong signs of interaction with a complex circumstellar medium, seen by the detection of X-rays, revealed by the detection of coronal emission lines, and through multi-component hydrogen and helium profiles. In order to further understand this interaction, we model the H α profile using a Monte Carlo electron scattering code. The blueshifted high-velocity component is consistent with emission from a radially thin spherical shell resulting in the broad emission components due to electron scattering. Using the velocity evolution of this emitting shell, we find that the SN ejecta collide with circumstellar material of at least ∼5 M ⊙ assuming a steady-state mass-loss rate of ∼4 − 6 × 10 −3 M ⊙ yr −1 for the first ∼200 days of evolution. SN 2021adxl was last observed to be slowly declining at ∼0.01 mag d −1 , and if this trend continues, SN 2021adxl will remain observable after its current solar conjunction. Continuing the observations of SN 2021adxl may reveal signatures of dust formation or an infrared excess, similar to that seen for SN 2010jl.

Brennan, S. J. (ORCID:0000000313256235)↗

DOE OSTI2022

Details of the neutral energy distribution and ionization source using spectrally resolved Balmer-alpha measurements on DIII-D

Spectrally resolved passive Balmer- α (D- α, H- α) measurements from the DIII-D 16 channel edge main-ion charge exchange recombination system confirm the presence of higher energy neutrals (“thermal” neutrals) in addition to the cold neutrals that recycle off the walls in the edge region of DIII-D plasmas. Charge exchange between thermal ions and edge neutrals transfers energy and momentum between the populations giving rise to thermal neutrals with energies approximating the ions in the pedestal region. Multiple charge exchange events in succession allow an electron to effectively take a random walk, transferring from ion to ion, providing a pathway of increasing energy and velocity, permitting a neutral to get deeper into the plasma before a final ionization event that contributes to the ion and electron particle fueling. Spectrally resolved measurements provide information about the density and velocity distribution of these neutrals, which has been historically valuable for validating Monte Carlo neutral models, which include the multi stage charge exchange dynamics. Here, in this study, a multi-channel set of such measurements is used to specifically isolate the details of the thermal neutrals that are responsible for fueling inside the pedestal top. Being able to separate the thermal from the cold emission overcomes several challenges associated with optical filter-based neutral density measurements. The neutral dynamics, deeper fueling by the thermal neutrals, and spectral measurement are modeled with the FIDASIM Monte Carlo collisional radiative code, which also produces synthetic spectra with a shape that is in close agreement with the measurements. By scaling the number of neutrals in the simulation to match the intensity of the thermal emission, we show it is possible to obtain local neutral densities and ionization source rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Development of Two-Step Method for Fast Chloride MSR Neutronics

This work presents neutronics models of a small and large fast-spectrum molten chloride-salt reactor. The models are similar to designs being pursued by industry, and they may serve as generic preconceptual and simplified neutronics models that provide information for decision making in licensing-related areas. Here, the two models were created using Serpent, a Monte Carlo neutron transport code, and Moltres, a neutron diffusion core simulator tool. Specifically, this study focused on exploring the applicability of diffusion theory to fast molten salt reactor (MSR) models, the capabilities of an open-source, MSR-oriented simulation tool (Moltres), and optimal energy-group structures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Monte Carlo Analysis of the Performance of the ITER Diagnostic Residual Gas Analyzer

The composition of exhausted gas is a key parameter in long-pulse plasma fusion experiments, and its evolution shall be monitored at timescales relevant to plasma dynamics and plasma-wall interactions. A diagnostic residual gas analyzer (DRGA) is a multisensor instrument particularly suited to these studies, and ITER will adopt DRGAs in the equatorial and in the divertor tokamak regions. In this work, we have revisited the design of the ITER divertor DRGA through simple vacuum analytical considerations supported by simulations conducted with Molflow+, a test particle Monte Carlo (TPMC) simulation code commonly used in the particle accelerator community. Here, starting with recommendations on the manufacturing of the vacuum piping of the DRGA, this work is followed by a complete vacuum characterization of the diagnostic vacuum setup (pressure profiles at base pressure and during sampling, orifice diameter, and length optimization), and finally, the in-vessel residence time of the most important gas species is simulated. These studies have allowed us to give insights into some experimental results recently found on the prototype DRGA installed in the Wendelstein W7-X stellarator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Inferring Neutron Star Merger Ejecta Morphology with Kilonovae

In this study we incorporate a new grid of kilonova simulations produced by the Monte Carlo radiative transfer code SuperNu in an inference pipeline for astrophysical transients and evaluate their performance. These simulations contain four different two-component ejecta morphology classes. We analyze follow-up observational strategies by Vera Rubin Observatory in optical and James Webb Space Telescope (JWST) in mid-infrared (MIR). Our analysis suggests that, within these strategies, it is possible to discriminate between the four different morphologies only when late-time JWST observations in MIR are available. We conclude that follow-ups by the new Vera Rubin Observatory alone are not sufficient to determine ejecta morphology. Additionally, we make comparisons between surrogate models based on radiative transfer simulation grids by SuperNu and POSSIS, by analyzing the historic kilonova AT2017gfo that accompanied the gravitational wave event GW170817. We show that both SuperNu and POSSIS models provide similar fits to photometric observations but their qualitative interpretations differ.

79 ASTRONOMY AND ASTROPHYSICS↗