Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 1,441 records · Page 80

Hyperparameter Setting for a Marked Multidimensional Hawkes Process with Dissimilar Decays

We provide further details for using Lim, et al.'s marked multidimensional Hawkes processes with dissimilar decays. We first describe what makes these different from other Hawkes processes, then describe each model hyperparameter and how to initialize it informed by the input data and any prior biases. We derived tighter bounds than Lim, et al. for faster convergence of rejection sampling. The resulting hyperparameters and bounds have been helpful against both synthetic and real-world datasets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nek5000: improvements in the available RANS models, meshing, tutorials, and training

This year, the Nuclear Energy Advanced Modeling Simulation program (NEAMS) thermal-hydraulics report for Nek5000 NRC- and verification and validation (V&V)-driven development focuses on following areas of code application and improvement. First we have continued improvements of RANS modeling capabilities in Nek5000 including improved k-tau model focusing mostly on wallfunction initial implementation with spectral element method (SEM) and initiating investigation of an alternative approach XSEM that greatly reduces discretization errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Summary of UMI Chemical Work

Our initial work for the Urban Modeling Integration (UMI) project was to validate the QUIC model dense gas dispersion of chlorine advection up a ramp. The validation task was to use QUIC to model an experiment done in a wind tunnel. Naturally the QUIC model uses a slumping cylindrical model where inside the cylinder the gas is treated as dense but outside ideal. In trying to validate using the wind tunnel experiment, due to the scale of the wind tunnel experiment and ramp, The model formed a very clear, discrete, concentration barrier at the edge of the cylinder that likely doesn’t line up well with the real dense gas advection/diffusion. This effect is much less noticeable at a larger scale such as an instantaneous release of a 25 ton tank.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Design Method for Low-Pressure Venturi Nozzles

The purpose of this work is to provide empirical design models for low-pressure, subsonic Venturi nozzles. Experimentally validated simulations were used to determine the effect of nozzle geometry and operating conditions on the suction ratio (ratio of suction mass flow rate to motive mass flow rate) of low-pressure, subsonic Venturi nozzles, over a wide range of geometries and operating conditions, through a parametric study. The results of the parametric study were used to develop seven empirical models, each with a different range of applicability or calculating a different indicator of nozzle performance (i.e., suction ratio, momentum ratio, or dynamic pressure ratio), of the Venturi nozzles using a constrained multi-variable global optimization method. Of the seven empirical models, the best models were found to be those for low- (less than one) and high-suction ratios (greater than one), with mean absolute percentage errors of 5% and 18%, respectively. These empirical models provide a design tool for subsonic, low-pressure Venturi nozzles that is more than an order of magnitude more accurate than a governing equation approach or conventional flow head calculations. These newly-developed empirical models can be applied for initial nozzle design when precise suction ratios are required.

O’Hern, Hannah↗

The entire virial radius of the fossil cluster RXJ 1159 + 5531. II. Dark matter and baryon fraction

In this second paper on the entire virial region of the relaxed fossil cluster RXJ 1159+5531, we present a hydrostatic analysis of the azimuthally averaged hot intracluster medium (ICM) using the results of Su et al. For a model consisting of ICM, stellar mass from the central galaxy (BCG), and an NFW dark matter (DM) halo, we obtain a good description of the projected radial profiles of ICM emissivity and temperature that yield precise constraints on the total mass profile. The BCG stellar mass component is clearly detected with a K-band stellar mass-to-light ratio, M{sub ⋆}/L{sub K}=0.61±0.11 M{sub ⊙}/L{sub ⊙}, consistent with stellar population synthesis models for a Milky Way initial mass function. We obtain a halo concentration, c{sub 200}=8.4±1.0, and virial mass, M{sub 200}=(7.9±0.6)×10{sup 13} M{sub ⊙}. For its mass, the inferred concentration is larger than most relaxed halos produced in cosmological simulations with Planck parameters, consistent with RXJ 1159+5531 forming earlier than the general halo population. The baryon fraction at r {sub 200}, f{sub b,200}=0.134±0.007, is slightly below the Planck value (0.155) for the universe. However, when we take into account the additional stellar baryons associated with non-central galaxies and the uncertain intracluster light (ICL), f{sub b,200} increases by ≈0.015, consistent with the cosmic value and therefore no significant baryon loss from the system. The total mass profile is nearly a power law over a large radial range (∼0.2–10 R {sub e}), where the corresponding density slope α obeys the α−R{sub e} scaling relation for massive early-type galaxies. Performing our analysis in the context of MOND still requires a large DM fraction (85.0%±2.5% at r = 100 kpc) similar to that obtained using the standard Newtonian approach. The detection of a plausible stellar BCG mass component distinct from the NFW DM halo in the total gravitational potential suggests that ∼10{sup 14} M{sub ⊙} represents the mass scale above which dissipation is unimportant in the formation of the central regions of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Galaxy-scale gravitational lens candidates from the hyper suprime-cam imaging survey and the galaxy and mass assembly spectroscopic survey

We present a list of galaxy-scale lens candidates including a highly probable interacting galaxy-scale lens in the Hyper Suprime-Cam (HSC) imaging survey. We combine HSC imaging with the blended-spectra catalog from the Galaxy And Mass Assembly (GAMA) survey to identify lens candidates, and use lens mass modeling to confirm the candidates. There are 45 matches between the HSC S14A-0b imaging data release and the GAMA catalog. We separate lens and lensed arcs using color information, and exclude those candidates with small image separations (<1.″0, estimated with the lens/source redshifts from the GAMA survey) that are not easily resolved with ground-based imaging. After excluding these, we find 10 probable lens systems. There is one system with an interacting galaxy pair, HSC J084928+000949, that has a valid mass model. We predict the total mass enclosed by the Einstein radius of ∼0.″72 (∼1.65 kpc) for this new expected lens system to be ∼10{sup 10.59} M{sub ⊙}. Using the photometry in the grizy bands of the HSC survey and stellar population synthesis modeling with a Salpeter stellar initial mass function, we estimate the stellar mass within the Einstein radius to be ∼10{sup 10.46} M{sub ⊙}. We thus find a dark matter mass fraction within the Einstein radius of ∼25%. Further spectroscopy or high-resolution imaging would allow confirmation of the nature of these lens candidates. The particular system with the interacting galaxy pair, if confirmed, would provide an opportunity to study the interplay between dark matter and stars as galaxies build up through hierarchical mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

A Census of Star Formation in the Outer Galaxy. II. The GLIMPSE360 Field

We have conducted a study of star formation in the outer Galaxy from 65°< l < 265° in the region observed by the GLIMPSE360 program. This Spitzer warm mission program mapped the plane of the outer Milky Way with IRAC at 3.6 and 4.5 μm. We combine the IRAC, Wide-field Infrared Survey Explorer (WISE), and Two Micron All Sky Survey catalogs and our previous results from another outer Galaxy survey and identify a total of 47,338 young stellar objects (YSOs) across the field spanning >180° in Galactic longitude. Using the DBSCAN method on the combined catalog, we identify 618 clusters or aggregations of YSOs having five or more members. We identify 10,476 class I, 29,604 class II, and 7325 anemic class II/class III YSOs. The ratio of YSOs identified as members of clusters was 25,528/47,338, or 54%. We found that 100 of the clusters identified have previously measured distances in the WISE H ii survey. We used these distances in our spectral energy distribution (SED) fitting of the YSOs in these clusters, of which 96 had YSOs with <3σ fits. We used the derived masses from the SED model fits to estimate the initial mass function (IMF) in the inner and outer Galaxy clusters; dividing the clusters by galactocentric distances, the slopes were Γ = 1.87 ± 0.31 above 3 M {sub ⊙} for R {sub Gal} < 11.5 kpc and Γ = 1.15 ± 0.24 above 3 M {sub ⊙} for R {sub Gal} > 11.5 kpc. The slope of the combined IMF was found to be Γ = 1.92 ± 0.42 above 3 M {sub ⊙}. These values are consistent with each other within the uncertainties and with literature values in the inner Galaxy high-mass star formation regions. The slopes are likely also consistent with a universal Salpeter IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

THE BURSTY STAR FORMATION HISTORIES OF LOW-MASS GALAXIES AT 0.4 < z < 1 REVEALED BY STAR FORMATION RATES MEASURED FROM H β AND FUV

We investigate the burstiness of star formation histories (SFHs) of galaxies at 0.4 < z < 1 by using the ratio of star formation rates (SFRs) measured from H β and FUV (1500 Å) (H β -to-FUV ratio). Our sample contains 164 galaxies down to stellar mass ( M {sub *}) of 10{sup 8.5} M {sub ⊙} in the CANDELS GOODS-N region, where Team Keck Redshift Survey Keck/DEIMOS spectroscopy and Hubble Space Telescope /WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the ratio of H β - and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H β -to-FUV ratio of our sample increases with M {sub *} and SFR. The median ratio is ∼0.7 at M {sub *} ∼ 10{sup 8.5} M {sub ⊙} (or SFR ∼ 0.5 M {sub ⊙} yr{sup −1}) and increases to ∼1 at M {sub *} ∼ 10{sup 10} M {sub ⊙} (or SFR ∼ 10 M {sub ⊙} yr{sup −1}). At M {sub *} < 10{sup 9.5} M {sub ⊙}, our median H β -to-FUV ratio is lower than that of local galaxies at the same M {sub *}, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation for our results, and the importance of the burstiness increases as M {sub *} decreases. Due to sample selection effects, our H β -to-FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

79 ASTRONOMY AND ASTROPHYSICS↗

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying Human Behavior and Decision Errors in Security Screening Operations

The potential for human errors in conducting security related screening operations can lead to inadvertent and adverse decision outcomes. This paper overviews an initial mathematical framework designed to model and quantify the various human factors and decision outcomes that may occur in conducting security screening operations, such as U.S. port of entry radiological and nuclear (rad/nuc) security screening. This framework is based on the Human Error Assessment and Reduction (HEART) technique. As applied here, the framework incorporates a set of rules for human engagement, including a prescribed concept of operations (CONOPS) and deviations that may occur from this established CONOPS due to inadvertent personnel decision errors. We also review some of the various factors that may adversely influence such decisions by security screening personnel. Some of these factors include current workload, environmental conditions, training, and various other intangible factors. Using the HEART methodology, we translate each of these factors into error producing conditions, their effects on error, an assessed proportion of effects, and finally the overall probability of human error at each stage of the screening process. We then include a small scale example to demonstrate the methodology and results based on an assumed set of input conditions at a notional port of entry.

Brigantic, Robert T.↗

Sensitivities to neutrino non-standard interactions from the disappearance of muon (anti)neutrinos at the NOvA experiment

The experimental confirmation of the phenomenon of neutrino oscillations provided evidence for a non-zero value for the mass of neutrinos, a property not initially predicted by the Standard Model of Physics. Investigations of the existence of the so called Physics Beyond the Standard Model is a subject of growing interest and effort in the field of particle physics. Interactions of neutrinos with matter are well described by the Standard Model, however in the presence of new Physics, additional phenomena regarding neutrino properties would be possible, including non-standard interactions of neutrinos with matter. The NOvA experiment, at Fermilab, has been investigating neutrino oscillations through studies of ν µ produced at the Fermilab accelerator facilities, and also comprise a suitable environment for investigation of additional phenomena taking part in the neutrino oscillation framework. We show in this thesis the results of the sensitivity studies to neutrino non standard interactions in the NOvA experiment through investigation of the disappearance of muon (anti)neutrinos during their evolution in the NOvA 810 km baseline. We investigate the effect of the NSI parameters, |ε µτ | and δ µτ , on the determination of the standard oscillation parameters sin 2 (θ 23 ) and Δ$m^2_{32}$, as well as the relation between each parameter, for both neutrino mass hierarchies. This study makes use of an exposure of 9.48 × 10 20 POT for neutrino beam (FHC), and 12.33 × 10 20 POT for antineutrino beam (RHC). No deviations from the standard oscillation scenario is found at the 90% CL limit. The joint neutrino and antineutrino fit yields a limit of |ε µτ | ≤ 0.086 (0.085) for a fixed value of δ µτ = 0, and |ε µτ | ≤ 1.054 (1.053) for the case where both NSI parameters are fitted, for the normal (inverted) neutrino mass hierarchy, at the 90% CL. In the current sensitivity, all values of δ µτ are found to be allowed. The 90% CL interval for the standard oscillation parameters are Δ$m^2_{32}$ ϵ [2.379, 2.580] and sin 2 (θ 23 ) ϵ [0.420, 0.602] for the normal hierarchy scenario, and Δ$m^2_{32}$ ϵ [-2.639, - 2.438] and sin 2 (θ 23 ) ϵ [0.418, 0.602] for the inverted hierarchy scenario.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

MELCOR Model Development of the General Atomics Fast Modular Reactor

General Atomics Electromagnetic Systems (GA-EMS) is currently developing a fast modular reactor (FMR) capable of providing safe, carbon-free electricity by 2035. As part of the development process, researchers at the University of Wisconsin-Madison (UWM) use MELCOR 2.2 to develop a simulation model of the proposed design to perform accident analyses. This paper summarizes the efforts of researchers at UWM in the creation of the FMR model. This model was built based on initial FMR design parameters. The design parameters will change as the design evolves and is refined.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Active Learning Meets Foundation Models: Fast Remote Sensing Data Annotation for Object Detection

Object detection in remote sensing demands extensive, high-quality annotations—a process that is both labor-intensive and time-consuming. In this work, we introduce a real-time active learning and semi-automated labeling framework that leverages foundation models to streamline dataset annotation for object detection in remote sensing imagery. For example, by integrating a Segment Anything Model (SAM), our approach generates mask-based bounding boxes that serve as the basis for dual sampling: (a) uncertainty estimation to pinpoint challenging samples, and (b) diversity assessment to ensure broad data coverage. Furthermore, our Dynamic Box Switching Module (DBS) addresses the well-known cold start problem for object detection models by replacing its suboptimal initial predictions with SAM-derived masks, thereby enhancing early-stage localization accuracy. Extensive evaluations on multiple remote sensing datasets plus a real-world user study, demonstrate that our framework not only reduces annotation effort, but also significantly boosts detection performance compared to traditional active learning sampling methods. The code for training and the user interface will be made available.

Burges, Marvin [ORNL] (ORCID:0000000312690769)↗

Transport Upscaling under Flow Heterogeneity and Matrix-Diffusion in Three-Dimensional Discrete Fracture Networks

For this work, we investigate the combined effects of network scale flow variability and retention due to matrix-diffusion on the scaling behavior of transport through fractured media. Two of the principal mechanisms controlling the transport of solutes through fractured low-permeability media are broad distributions of flow velocities and retention times in the solid matrix. We study the relative impact of these two processes under different initial conditions using a set of three-dimensional discrete fracture network simulations. We use these simulations to develop and calibrate an upscaled continuous time random walk (CTRW) approach for advective transport based on an Ornstein-Uhlenbeck model for the particle velocities that accounts for the fracture-matrix coupling using a compound Poisson process. This CTRW model can be conditioned on the initial solute distribution and allows to observe late-time scaling behavior at distances beyond what is feasible to observe using high-fidelity direct numerical simulations. We determine that the initial distribution of particles leads to marked differences in the persistent long-term scale behavior in the solute travel time distributions, even those undergoing retention due to matrix diffusion through implementation and analysis of the model.

54 ENVIRONMENTAL SCIENCES↗

Behavior of Bilayer Leaflets in Asymmetric Model Membranes: Atomistic Simulation Studies

Spatial organization within lipid bilayers is an important feature for a range of biological processes. Leaflet compositional asymmetry and lateral lipid organization are just two of the ways in which membrane structure appears to be more complex than initially postulated by the fluid mosaic model. This raises the question of how the phase behavior in one bilayer leaflet may affect the apposing leaflet and how one begins to construct asymmetric model systems to investigate these interleaflet interactions. In this paper, we report on all-atom molecular dynamics simulations (a total of 4.1 μs) of symmetric and asymmetric bilayer systems composed of liquid-ordered (Lo) or liquid-disordered (Ld) leaflets, based on the nanodomain-forming POPC/DSPC/cholesterol system. We begin by analyzing an asymmetric bilayer with leaflets derived from simulations of symmetric Lo and Ld bilayers. In this system, we observe that the properties of the Lo and Ld leaflets are similar to those of the Lo and Ld leaflets in corresponding symmetric systems. However, it is not obvious that mixing the equilibrium structures of their symmetric counterparts is the most appropriate way to construct asymmetric bilayers nor that these structures will manifest interleaflet couplings that lead to domain registry/antiregistry. We therefore constructed and simulated four additional asymmetric bilayer systems by systematically adding or removing lipids in the Ld leaflet to mimic potential density fluctuations. We find that the number of lipids in the Ld leaflet affects its own properties, as well as those of the apposing Lo leaflet. Collectively, the simulations reveal the presence of weak acyl chain interdigitation across bilayer leaflets, suggesting that interdigitation alone does not contribute significantly to the interleaflet coupling in nonphase-separated bilayers of this chemical composition. Finally, however, the properties of both leaflets appear to be sensitive to changes in in-plane lipid packing, possibly providing a mechanism for interleaflet coupling by modulating local density and/or curvature fluctuations.

59 BASIC BIOLOGICAL SCIENCES↗

Cloud, Aerosol, and Complex Terrain Interactions (CACTI) Science Plan

General circulation models and downscaled regional models exhibit persistent biases in deep convective initiation location and timing, cloud top height, stratiform area and precipitation fraction, and anvil coverage. Despite important impacts on the distribution of atmospheric heating, moistening, and momentum, nearly all climate models fail to represent convective organization, while system evolution is not represented at all. Improving representation of convective systems in models requires characterization of their predictability as a function of environmental conditions, and this characterization depends on observing many cases of convective initiation, non-initiation, organization, and non-organization. The Cloud, Aerosol, and Complex Terrain Interactions (CACTI) experiment in the Sierras de Córdoba mountain range of north-central Argentina is designed to improve understanding of cloud life cycle and organization in relation to environmental conditions so that cumulus, microphysics, and aerosol parameterizations in multiscale models can be improved. The Sierras de Córdoba range has a high frequency of orographic boundary-layer clouds, many reaching congestus depths, many initiating into deep convection, and some organizing into mesoscale systems uniquely observable from a single fixed site. Some systems even grow upscale to become among the deepest, largest, and longest-lived in the world. These systems likely contribute to an observed regional trend of increasing extreme rainfall, and poor prediction of them likely contributes to a warm, dry bias in climate models downstream of the Sierras de Córdoba range in a key agricultural region. Many environmental factors influence the convective life cycle in this region, including orographic, low-level jet, and frontal circulations, surface fluxes, synoptic vertical motions influenced by the Andes, cloud detrainment, and aerosol properties. Local and long-range transport of smoke resulting from biomass burning as well as blowing dust are common in the austral spring, while changes in land surface properties as the wet season progresses impact surface fluxes and boundary-layer evolution on daily and seasonal time scales that feed back to cloud and rainfall generation. This range of environmental conditions and cloud properties coupled with a high frequency of events makes this an ideal location for improving our understanding of cloud-environment interactions. The following primary science questions will be addressed through coordinated first U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1), mobile C-band Scanning ARM Precipitation Radar (C-SAPR2), ARM Aerial Facility (AAF) Gulfstream-1 (G-1), and guest instrument observations: How are the properties and life cycles of orographically generated cumulus humulis, mediocris, and congestus clouds affected by environmental kinematics, thermodynamics, aerosols, and surface properties? How do these cloud types alter these environmental conditions? How do environmental kinematics, thermodynamics, and aerosols impact deep convective initiation, upscale growth, and mesoscale organization? How are soil moisture, surface fluxes, and aerosol properties altered by deep convective precipitation events and seasonal accumulation of precipitation? This multi-faceted experiment involves a long-term 7-month extended operational period (EOP, 1 October, 2018–30 April, 2019) as well as a 1.5-month intensive operational period (IOP, 30 October–13 December) that will include G-1 observations coinciding with the international multi-agency Remote Sensing of Electrification, Lightning, and Meso-scale/Micro-scale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign.

Convection, Andes, Argentina, cloud life cycle, cl↗

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Extremely Bright GRB 160625B with Multiple Emission Episodes: Evidence for Long-term Ejecta Evolution

GRB 160625B is an extremely bright GRB with three distinct emission episodes. By analyzing its data observed with the Gamma-Ray Burst Monitor (GBM) and Large Area Telescope (LAT) on board the Fermi mission, we find that a multicolor blackbody (mBB) model can be used to fit very well the spectra of the initial short episode (Episode I) within the hypothesis of photosphere emission of a fireball model. The time-resolved spectra of its main episode (Episode II), which was detected with both GBM and LAT after a long quiescent stage (∼180 s) following the initial episode, can be fitted with a model comprising an mBB component plus a cutoff power-law (CPL) component. This GRB was detected again in the GBM and LAT bands with a long extended emission (Episode III) after a quiescent period of ∼300 s. The spectrum of Episode III is adequately fitted with CPL plus single power-law models, and no mBB component is required. These features may imply that the emission of the three episodes are dominated by distinct physics processes, i.e., Episode I is possible from the cocoon emission surrounding the relativistic jet, Episode II may be from photosphere emission and internal shock of the relativistic jet, and Episode III is contributed by internal and external shocks of the relativistic jet. On the other hand, both X-ray and optical afterglows are consistent with the standard external shocks model.

79 ASTRONOMY AND ASTROPHYSICS↗