Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Modeling”

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 811 records · Page 45

Symmetry tuning of a near one-dimensional 2-shock platform for code validation at the National Ignition Facility

Here, we introduce a new quasi 1-D implosion experimental platform at the National Ignition Facility designed to validate physics models as well as to study various Inertial Confinement Fusion aspects such as implosion symmetry, convergence, hydrodynamic instabilities, and shock timing. The platform has been developed to maintain shell sphericity throughout the compression phase and produce a round hot core at stagnation. This platform utilizes a 2-shock 1 MJ pulse with 340 TW peak power in a near-vacuum AuHohlraum and a CH ablator capsule uniformly doped with 1% Si. We also performed several inflight radiography, symmetry capsule, and shock timing experiments in order to tune the symmetry of the capsule to near round throughout several epochs of the implosion. Finally, adjusting the relative powers of the inner and outer cones of beams has allowed us to control the drive at the poles and equator of the capsule, thus providing the mechanism to achieve a spherical capsule convergence. Details and results of the tuning experiments are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Modeling of Physical Vapor Transport in Contactless Crystal Growth Geometry

Growth from the vapor under conditions of limited contact with the walls of the growth ampoule is beneficial for the quality of the growing crystal due to reduced stress and contamination which may be caused by interactions with the growth container. The technique may be of a particular interest for studies on crystal growth under microgravity conditions: elimination of some factors affecting the crystal quality may make interpretation of space-conducted processes more conclusive and meaningful. For that reason, and as a part of our continuing studies on 'contactless' growth technique, we have developed a computational model of crystal growth process in such system. The theoretical model was built, and simulations were performed using the commercial computational fluid dynamics code, (CFD) ACE. The code uses an implicit finite volume formulation with a gray discrete ordinate method radiation model which accounts for the diffuse absorption and reflection of radiation throughout the furnace. The three-dimensional model computes the heat transfer through the crystal, quartz, and gas both inside and outside the ampoule, and mass transport from the source to the crystal and the sink. The heat transport mechanisms by conduction, natural convection, and radiation, and mass transport by diffusion and convection are modeled simultaneously and include the heat of the phase transition at the solid-vapor interfaces. As the thermal boundary condition, temperature profile along the walls of the furnace is used. For different thermal profiles and furnace and ampoule dimensions, the crystal growth rate and development of the crystal-vapor and source-vapor interfaces (change of the interface shape and location with time) are obtained. Super/under-saturation in the ampoule is determined and critical factors determining the 'contactless' growth conditions are identified and discussed. The relative importance of the ampoule dimensions and geometry, the furnace dimensions and its temperature, and the properties of the grown material are analyzed. The results of the simulations are compared with related experimental results on growth of CdTe, CdZnTe, ZnTe, PbTe, and PbSnTe crystals by this technique.

Palosz, W.↗

Effect of dislocation pile-up on size-dependent yield strength in finite single-crystal micro-samples

Recent research has explained that the steeply increasing yield strength in metals depends on decreasing sample size. In this work, we derive a statistical physical model of the yield strength of finite single-crystal micro-pillars that depends on single-ended dislocation pile-up inside the micro-pillars. We show that this size effect can be explained almost completely by considering the stochastic lengths of the dislocation source and the dislocation pile-up length in the single-crystal micro-pillars. The Hall–Petch-type relation holds even in a microscale single-crystal, which is characterized by its dislocation source lengths. Our quantitative conclusions suggest that the number of dislocation sources and pile-ups are significant factors for the size effect. They also indicate that starvation of dislocation sources is another reason for the size effect. Moreover, we investigated the explicit relationship between the stacking fault energy and the dislocation “pile-up” effect inside the sample: materials with low stacking fault energy exhibit an obvious dislocation pile-up effect. Our proposed physical model predicts a sample strength that agrees well with experimental data, and our model can give a more precise prediction than the current single arm source model, especially for materials with low stacking fault energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Public-Private-Academic Partnership to Advance Solar Power Forecasting

Executive Summary NOAA is making major contributions to the solar forecasting project in three areas. First, it is improving its forecasts of solar irradiance, clouds, and aerosols in its numerical weather prediction models. Second, it is providing advanced satellite products for DOE's FOA awardees to use in their forecast systems. Third, it is using high-quality ground-based measurements from SURFRAD and ISIS stations to verify and validate forecast model output. This reports covers results from all three areas for the period May 1, 2014 - April 30, 2015. Modeling In its modeling effort, NOAA continues work to improve the skill of solar forecasts from the Earth System Research Lab (ESRL) research versions of the 13-km Rapid Refresh (RAP) and the 3-km High-Resolution Rapid Refresh (HRRR) models, which are in turn transitioned into operations at the National Centers for Environmental Prediction (NCEP). A major milestone was achieved in September 2014 with the initial operational implementation of the HRRR at NCEP. In the ESRL research versions of the models, testing and development, in both real-time runs and retrospective experiments, is guided by an extensive in-house verification system. Early in the SFIP project, we developed the capability to verify our model forecasts against the high-quality surface radiation measurements from the SURFRAD and ISIS networks. This highlighted some shortcomings with the RAP and HRRR forecasts of incoming shortwave radiation. Most of our effort during Phase 1 of SFIP was focused on addressing these problems with a variety of model system improvements. The RAP and HRRR models during the warm season of 2014 had a noticeable warm and dry bias in near-surface conditions over most of the central and eastern United States, and our new SURFRAD/ISIS verification revealed that there was also a large excess of incoming global horizontal irradiance in the models. We hypothesized that a lack of cloud cover (particularly low-level cloud cover) in the models was resulting in too much heating of the land surface. This, in turn, caused unrealistically strong surface heat fluxes and turbulent mixing in the planetary boundary layer (PBL), which further reduced the already deficient cloud cover. We addressed these issues with a combination of data assimilation system modifications and model physics improvements. Many of our data assimilation changes were made with a view towards improving the near-term representation of clouds and precipitation. One of these changes involved better accounting for regions of weak reflectivity in the RAP cloud / hydrometeor assimilation system, in order to improve the representation of light precipitation in the RAP initial conditions and provide more realistic initial cloud cover. Additional modifications more accurately accounted for radar beam blockage and data gaps (particularly in the western United States), which improves shorter lead times forecasts of clouds and precipitation. We have also tested the assimilation of new data sources within the RAP and the HRRR, including radar radial velocity data and surface mesonet observations. Within the HRRR, we have tested the cycling of the 3-km land surface fields to allow a higher-resolution treatment of land surface processes. In terms of model physics development for SFIP, we have implemented a shallow cumulus scheme within the RAP, and have made numerous improvements to the Mellor-Yamada-Nakanishi-Niino (MYNN) PBL scheme to address insufficient low-level cloud cover in the models. We have conducted tests incorporating the radiation effects of (parameterized) boundary-layer clouds within the modified MYNN PBL scheme (independent of the convective schemes). The Grell-Freitas-Olson shallow cumulus scheme has also been tested within the 3-km HRRR. Finally, we have also modified the RUC land surface model (LSM) treatment of the vegetation wilting point, reducing it to increase evapotranspiration and increase cloud cover in the boundary layer. All of these changes work in tandem to significantly improve the model forecasts of cloud cover, incoming shortwave radiation, and near-surface temperature and moisture. Satellite The role of NOAA/NESDIS in the Solar Forecasting Improvement Project is to provide Advanced Satellite Products (ASPs) for the two forecasting teams at NCAR and IBM. The ASPs are cloud, surface, and atmosphere products derived from geostationary satellite imagery at the highest possible spatial and temporal resolution - such quantities as cloud mask, cloud probability, cloud transmission, cloud top height, cloud top temperature, cloud effective particle size, etc. Ancillary data, such as elevation and numerical weather prediction fields are provided in the files at the same resolution as well. There are at this time 147 different variables in the ASP output, including quality flags and processing information. The main goals for Year 1 of the project were to implement an Advanced Satellite Products system for the use of the IBM and NCAR teams, begin validation, and make any needed changes based on feedback from the teams. ASP files are being produced every GOES Imager acquisition, which occur on a 15-30 minute schedule. Processing is done on a dedicated computer, with a turn-around time of 8-21 minutes from image acquisition to results available on ftp. Several helpful visualizations of the data are also created for users on web pages. Users have been provided with a document titled "User's Guide for 1km Cloud Products Derived from GOES Imager Data using CLAVR-x", which discusses the basics of the source imagery, the process by which it is turned into Advanced Satellite Products, and considerations users should make when using the data. Validation of selected variables from the older 4km version of the products was also included. Future work will concentrate on validation of the 1km products and improving the turn-around time, product variety, and product quality as needed. Ground Observations In the ground-based measurement effort, NOAA's main objectives are to provide high quality radiation products for validation and verification of short-term to day-ahead solar forecasts. More specifically for the three year project, our goals include (1) Maintaining and providing data from our 7 SURFRAD and 7 ISIS; (2) Update ISIS radiation measurements from 3 min to 1 min data: (3) Purchase and install new pyrheliometers for direct solar irradiance measurements at the 7 SURFRAD sites; (4) Building, testing, and deploying two mobile SURFRAD stations at two utility plants in collaboration with DOE sponsored partners, and includes ongoing maintenance and processing of the data at the mobile sites; (5) Upgrading the data acquisition and communications at 7 SURFRAD sites and 7 ISIS sites; (6) Providing radiation data at the 7 SURFRAD sites in near real-time; (7) Develop and provide aerosol optical depth and cloud images and cloud fraction at our two mobile sites; (8) Provide data recovery rates each year; (9) Provide temporally and spatially averaged radiation products for comparison to HRRR and RAP solar forecasts and advanced satellite products; (10) Provide a data-set for analysis of conversion of direct and diffuse to sloped surfaces; (11) and as time permits develop and provide spectral solar irradiance, cloud optical depth and spectral albedo from the mobile sites. Milestones this year include working with the DOE sponsored teams to find locations to deploy two mobile SURFRAD stations. One existing unit was deployed at a 30MW PV facility in the San Luis Valley in collaboration with Xcel and the NCAR team in August, 2014. The second unit was built and tested at our facilities in Boulder, CO and deployed near Green Mountain Power's Education Center in Rutland, VT in collaboration with Green Mountain Power and the IBM Team in October, 2014. Data processing was implemented and the radiation data from these two mobile sites have been made available on our ftp server in near real-time. We also are providing images and cloud fraction from the TSI cameras for these two mobile sites on our ftp site. Another milestone was upgrading our data acquisition and communication systems at 7 SURFRAD and 7 ISIS sites. We accelerated our schedule for these upgrades to provide timely radiation products. These upgrades allow more reliable and near-real time radiation data delivery to the DOE sponsored teams to meet their goals. Lastly, we changed the data rate at the ISIS sites from 3 min to 1 min.

14 SOLAR ENERGY↗

Experimental Investigation on Metallic Droplet Behavior in Molten BOF Slag

In order to better understand the metallic droplet behavior during a slag treatment process, a physical modeling based on the similarity principle was performed in a transparent scaled-down vessel at room temperature. Paraffin oil, 20 wt pct copper sulfate solution, and compressed air were used to simulate the molten slag, metallic droplet, and carrier gas, respectively. The droplets injected into paraffin oil during the experiment were captured by a high speed camera and were analyzed by Image Pro Plus software to obtain the droplet size distribution. The critical droplet size in the physical modeling and slag treatment process is quantitatively correlated. The results show that droplet breakage phenomenon is dominant over its coalescence in the current industrial practice, and droplet breakage is enhanced with increasing gas flow rate and/or lance depth. No significant effect of the nozzle configuration was found on the droplet breakage and coalescence. The droplet size distribution varies with the lance position. Gas flow rate and lance depth are the most important factors for droplet breakage, the extent of which can be reduced through a proper selection of the operational conditions. A linear relationship between the droplet size and the input energy flux is obtained.

36 MATERIALS SCIENCE↗

Multiscale high-order/low-order (HOLO) algorithms and applications

We review the state of the art in the formulation, implementation, and performance of so-called high-order/low-order (HOLO) algorithms for challenging multiscale problems. HOLO algorithms attempt to couple one or several high-complexity physical models (the high-order model, HO) with low-complexity ones (the low-order model, LO). The primary goal of HOLO algorithms is to achieve nonlinear convergence between HO and LO components while minimizing memory footprint and managing the computational complexity in a practical manner. Key to the HOLO approach is the use of the LO representations to address temporal stiffness, effectively accelerating the convergence of the HO/LO coupled system. The HOLO approach is broadly underpinned by the concept of nonlinear elimination, which enables segregation of the HO and LO components in ways that can effectively use heterogeneous architectures. The accuracy and efficiency benefits of HOLO algorithms are demonstrated with specific applications to radiation transport, gas dynamics, plasmas (both Eulerian and Lagrangian formulations), and ocean modeling. Across this broad application spectrum, HOLO algorithms achieve significant accuracy improvements at a fraction of the cost compared to conventional approaches. It follows that HOLO algorithms hold significant potential for high-fidelity system scale multiscale simulations leveraging exascale computing.

97 MATHEMATICS AND COMPUTING↗

Scaling Containment on a Large Centrifuge

Containment science concerns the trapping or leakage of radioactive cavity gases from underground nuclear explosions (UNEs). The most physically relevant data for validation of containment science comes from legacy UNE testing or large field-scale chemical explosions. Field-scale tests can inform on the degree of residual stress imparted to host geology and possible formation of an accompanying gas-containing “stress cage”, for example, as a function of chemical explosion yield. However, field testing can be expensive, difficult to conduct, and challenging to instrument with limitations on data coverage (i.e., the number of boreholes; difficult coring through damage zones). Small-scale laboratory tests are typically simpler to conduct and more thoroughly characterize, but may not include the right scale of containment processes for stress-cage or chimney formation. This report presents the theory and fundamentals of large centrifuge physical modeling, including a history of explosive or similar testing relevant to containment science investigations on centrifuges. Subscale models with embedded chemical explosives in the enhanced gravity of a large centrifuge can represent hundreds of meters of depth and large explosions not otherwise attainable in the laboratory as based on scaling of length, energy, and other processes by the g-factor of the centrifuge (i.e., the number of times larger the centrifugal force is than the gravitational force on the Earth’s surface). The centrifuge uniquely joins the simplicity and exhaustiveness of laboratory-scale characterization with the physics of field-scale processes. This report is part of Sandia National Laboratories’ (SNL’s) effort to prepare for physical modeling with approximately one-meter tall geologic models of stress cage and/or chimney formation on its 29-foot [8.84 m] radius centrifuge with a load capacity of 1.6 million 𝑔-pounds [726 𝑔-ton]. The scaling relationships herein will inform model design decisions and performance requirements for in-flight sensors.

42 ENGINEERING↗

Three-dimensional simulation of H-mode plasmas with localized divertor impurity injection on Alcator C-Mod using the edge transport code EMC3-EIRENE

Experiments in Alcator C-Mod to assess the level of toroidal asymmetry in divertor conditions resulting from poloidally and toroidally localized extrinsic impurity gas seeding show a weak toroidal peaking (∼1.1) in divertor electron temperatures for high-power enhanced D-alpha H-mode plasmas. This is in contrast to similar experiments in Ohmically heated L-mode plasmas, which showed a clear toroidal modulation in the divertor electron temperature. Modeling of these experiments using the 3D edge transport code EMC3-EIRENE [Y. Feng et al., J. Nucl. Mater. 241, 930 (1997)] qualitatively reproduces these trends, and indicates that the different response in the simulations is due to the ionization location of the injected nitrogen. Low electron temperatures in the private flux region (PFR) in L-mode result in a PFR plasma that is nearly transparent to neutral nitrogen, while in H-mode the impurities are ionized in close proximity to the injection location, with this latter case yielding a largely axisymmetric radiation pattern in the scrape-off-layer. The consequences for the ITER gas injection system are discussed. Quantitative agreement with the experiment is lacking in some areas, suggesting potential areas for improving the physics model in EMC3-EIRENE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vectorized Monte Carlo methods for reactor lattice analysis

Some of the new computational methods and equivalent mathematical representations of physics models used in the MCV code, a vectorized continuous-enery Monte Carlo code for use on the CYBER-205 computer are discussed. While the principal application of MCV is the neutronics analysis of repeating reactor lattices, the new methods used in MCV should be generally useful for vectorizing Monte Carlo for other applications. For background, a brief overview of the vector processing features of the CYBER-205 is included, followed by a discussion of the fundamentals of Monte Carlo vectorization. The physics models used in the MCV vectorized Monte Carlo code are then summarized. The new methods used in scattering analysis are presented along with details of several key, highly specialized computational routines. Finally, speedups relative to CDC-7600 scalar Monte Carlo are discussed.

Brown, F. B.↗

Laser assisted crystallization of ferromagnetic amorphous ribbons: A multimodal characterization and thermal model study

This paper focuses on laser-based de-vitrification of amorphous soft magnetic Fe-Si-B ribbons and its consequent influence on the magnetic properties. Laser processing resulted in a finer scale of crystallites due to rapid heating and cooling during laser annealing compared to conventional furnace annealing process. A significant increase in saturation magnetization is observed for laser-annealed ribbons compared to both as-received and furnace annealed samples coupled with an increase in coercivity compared to the as received samples. The combined effect of thermal histories and stresses developed during laser annealing results in the formation of nano-crystalline phase along the laser track. The phase evolution is studied by micro-XRD and TEM analysis. Solute partitioning and compositional variation within the phases are obtained by Local Electrode Atom probe analysis. The evolution of microstructure is rationalized using a Finite Element based heat transfer multi-physics model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using Bayesian Networks for Sensitivity Analysis of Complex Biogeochemical Models

Sensitivity analysis is a vital tool in numerical modeling to identify important parameters and processes that contribute to the overall uncertainty in model outputs. We developed a new sensitivity analysis method to quantify the relative importance of uncertain model processes that contain multiple uncertain parameters. The method is based on the concepts of Bayesian networks (BNs) to account for complex hierarchical uncertainty structure of a model system. We derived a new set of sensitivity indices using the methodology of variance-based global sensitivity analysis with the Bayesian inference. The framework is capable of representing the detailed uncertainty information of a complex model system using BNs and affords flexible grouping of different uncertain inputs given their characteristics and dependency structures. We have implemented the method on a real-world biogeochemical model at the groundwater-surface water interface within the Hanford Site’s 300 Area. The uncertainty sources of the model were first grouped into forcing scenario and three different processes based on our understanding of the complex system. The sensitivity analysis results indicate that both the reactive transport and groundwater flow processes are important sources of uncertainty for carbon-consumption predictions. Within the groundwater flow process, the structure of geological formations is more important than the permeability heterogeneity within a given geological formation. Our new sensitivity analysis framework based on BNs offers substantial flexibility for investigating the importance of combinations of interacting uncertainty sources in a hierarchical order, and it is expected to be applicable to a wide range of multi-physics models for complex systems.

Dai, Heng↗

Physically Based Modeling of PEMFC Cathode Catalyst Layers: Effective Microstructure and Ionomer Structure–Property Relationship Impacts

Abstract This work presents a pseudo-two-dimensional proton exchange membrane fuel cell (PEMFC) model incorporating Nafion ionomer structure–property relationships in the cathode catalyst layer (CL) to capture and explain losses at low Pt loading. Structural data from neutron reflectometry and thin film Nafion conductivity measurements predict variations in the oxygen diffusion coefficient and ionic conductivity with changing CL ionomer thickness and Pt loading. By including these structure–property relationships, predicted polarization curves agree closely with previously published experimental data from cells with Pt loadings between 0.025 and 0.2 mg/cm2. Results demonstrate that structure–property relationships based on physically measurable ionomer and CL properties provide a feasible interpretation of PEMFC CL phenomena for a range of Pt loadings and help explain previously unaccounted-for losses at low Pt. Results also show that simulations must account for surface species coverage variations in order to properly capture the kinetic losses. Finally, results suggest that an increase in ionomer thickness surrounding the C/Pt surfaces may lead to improved cell performance due to improved ionic conductivity.

Electrochemistry↗

Online Bayesian State Estimation for Real-Time Monitoring of Growth Kinetics in Thin Film Synthesis

Rapid validation of newly predicted materials through autonomous synthesis requires real-time adaptive control methods that exploit physics knowledge, a capability that is lacking in most systems. Here, in this study, we demonstrate an approach to enable real-time control of thin film synthesis by combining in situ optical diagnostics with a Bayesian state estimation method. We developed a physical model for film growth and applied the direct filter (DF) method for real-time estimation of nucleation and growth rates during pulsed laser deposition (PLD). We validated the approach using simulated and experimental reflectivity data for WSe 2 growth and ultimately deployed the algorithm on an autonomous PLD system during the growth of 1T'-MoTe 2 . The DF robustly estimates growth parameters in real time at early stages of growth, down to 15% monolayer area coverage. This fusion of in situ diagnostics, data assimilation, and physical modeling opens new opportunities in adaptive control of synthesis trajectories toward desired material states.

36 MATERIALS SCIENCE↗

Using onset times from frequent seismic surveys to understand fluid flow at the Peace River Field, Canada

With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.

58 GEOSCIENCES↗

Event generators for high-energy physics experiments

We provide an overview of the status of Monte-Carlo event generators for high-energy particle physics. Guided by the experimental needs and requirements, we highlight areas of active development, and opportunities for future improvements. Particular emphasis is given to physics models and algorithms that are employed across a variety of experiments. These common themes in event generator development lead to a more comprehensive understanding of physics at the highest energies and intensities, and allow models to be tested against a wealth of data that have been accumulated over the past decades. A cohesive approach to event generator development will allow these models to be further improved and systematic uncertainties to be reduced, directly contributing to future experimental success. Event generators are part of a much larger ecosystem of computational tools. They typically involve a number of unknown model parameters that must be tuned to experimental data, while maintaining the integrity of the underlying physics models. Making both these data, and the analyses with which they have been obtained accessible to future users is an essential aspect of open science and data preservation. It ensures the consistency of physics models across a variety of experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A unique, ring-like radio source with quadrilateral structure detected with machine learning

ABSTRACT We report the discovery of a unique object in the MeerKAT Galaxy Cluster Legacy Survey (MGCLS) using the machine learning anomaly detection framework astronomaly. This strange, ring-like source is 30′ from the MGCLS field centred on Abell 209, and is not readily explained by simple physical models. With an assumed host galaxy at redshift 0.55, the luminosity (1025 W Hz−1) is comparable to powerful radio galaxies. The source consists of a ring of emission 175 kpc across, quadrilateral enhanced brightness regions bearing resemblance to radio jets, two ‘ears’ separated by 368 kpc, and a diffuse envelope. All of the structures appear spectrally steep, ranging from −1.0 to −1.5. The ring has high polarization (25 per cent) except on the bright patches (<10 per cent). We compare this source to the Odd Radio Circles recently discovered in ASKAP data and discuss several possible physical models, including a termination shock from starburst activity, an end-on radio galaxy, and a supermassive black hole merger event. No simple model can easily explain the observed structure of the source. This work, as well as other recent discoveries, demonstrates the power of unsupervised machine learning in mining large data sets for scientifically interesting sources.

Astronomy & Astrophysics↗

Verification of Reactor Activation Modelling with Physical Characterization Data and the Impact on Long Term Assessments - 20303

Whiteshell Laboratories, located in Manitoba, Canada, provided research facilities for the Canadian nuclear industry since the early 1960's. The facility was centered on Whiteshell Reactor 1 (WR-1), an organically cooled, heavy water moderated nuclear reactor. WR-1 was safely shut down in 1985, defueled and drained, and has been in a safe Storage with Surveillance (Deferred Decommissioning) state until present day. In support of an Environmental Assessment for proposed in-situ disposal of WR-1, detailed characterization information regarding the reactor core components (calandria, calandria tubes, and pressure tubes) was necessary to provide confidence that the long-term impacts of the facility can be assessed adequately. A simplified estimate of the total activation of core components was completed in 1992 using conservative assumptions, expected reactor configuration and computer models built with WIMS-CRNL, ONEDANT and ORIGEN-S. Due to the simplification, and age of the work, there was considerable doubt that the modelled inventory estimate of the core components was sufficient to support long-term performance modeling of the planned in-situ decommissioning. Additional data would be required to provide this confidence. Several options were considered for characterizing the WR-1 core components, but it was decided to collect only a limited number of samples in an effort to validate the 1992 simplified estimate as bounding. The limited number of samples reduced risks to workers, both radiological and non-radiological, and significantly reduced the costs required to confidently define and bound the radiological inventory of activation products in the core. The results of the sampling were compared to the simplified activation estimation. There was good agreement between the modelled and sampled results, with all major contributors to total activity showing similar relative quantities. The samples were generally lower than predicted by the models, with few exceptions. Notable exceptions included Nb-94, which were higher in the sample results than the model. This was determined to be caused by the use of Niobium containing steel alloys in the stainless steel pressure tube that were not accounted for in the models. The comparison of results provided confidence that the original inventory estimates produced through the simplified model are conservative and bounding, and therefore suitable for use in the long-term assessment of in situ disposal of the WR-1 reactor. This paper provides details of the simplified computer model results, their comparison to results of the sampling, and the changes adopted in the assessment of the in-situ decommissioning of WR-1 as a result. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effects of thermal fluctuations on thermal inflation

The mechanism of thermal inflation, a relatively short period of accelerated expansion after primordial inflation, is a desirable ingredient for a certain class of particle physics models if they are not to be in contention with the cosmology of the early Universe. Though thermal inflation is most simply described in terms of a thermal effective potential, a thermal environment also gives rise to thermal fluctuations that must be taken into account. We numerically study the effects of these thermal fluctuations using lattice simulations. We conclude that though they do not ruin the thermal inflation scenario, the phase transition at the end of thermal inflation proceeds through phase mixing and is therefore not accompanied by the formations of bubbles nor appreciable amplitude of gravitational waves.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗