Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reach 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.

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

At least 415 records · Page 23

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counterpropagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scattering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. Finally, we compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Thermostructural Analysis of a Diboride Composite Leading Edge

In an effort to support the design of zirconium diboride composite leading edges for hypersonic vehicles, a finite element model (FEM) of a prototype leading edge was created and finite element analysis (FEA) was employed to assess its thermal and structural response to aerothermal boundary conditions. Unidirectional material properties for the structural components of the leading edge, a continuous fiber reinforced diboride composite, were computed with COSTAR. These properties agree well with those experimentally measured. To verify the analytical approach taken with COSMOS/M, an independent FEA of one of the leading edge assembly components was also done with COSTAR. Good agreement was obtained between the two codes. Both showed that a unidirectional lay-up had the best margin of safety for a simple loading case. Both located the maximum stress in the same region and ply. The magnitudes agreed within 4 percent. Trajectory based aerothermal heating was then applied to the leading edge assembly FEM created with COSMOS/M to determine steady state temperature response, displacement, stresses, and contact forces due to thermal expansion and thermal strains. Results show that the leading edge stagnation line temperature reached 4700 F. The maximum computed failure index for the laminated composite components peaks at 4.2, and is located at the bolt flange in layer 2 of the side bracket. The temperature gradient in the tip causes a compressive stress of 279 ksi along its width and substantial tensile stresses within its depth.

Kowalski, Tom↗

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counter- propagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scat- tering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. We compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhancement of the Earth Science and Remote Sensing Group's Website and Related Projects

The major problem addressed throughout the term was the need to update the group's current website, as it was outdated and required streamlining and modernization. The old Gateway to Astronaut Photography of the Earth website had multiple components, many of which involved searches through expansive databases. The amount of work required to update the website was large and due to a desired release date, assistance was needed to help build new pages and to transfer old information. Additionally, one of the tools listed on the website called Image Detective had been underutilized in the past. It was important to address why the public was not using the tool and how it could potentially become more of a resource for the team. In order to help with updating the website, it was necessary to first learn HTML. After assisting with small edits, I began creating new pages. I utilized the "view page source" and "developer" tools in the internet browser to observe how other websites created their features and to test changes without editing the code. I then edited the code to create an interactive feature on the new page. For the Image Detective Page I began an evaluation of the current page. I also asked my fellow interns and friends at my University to offer their input. I took all of the opinions into account and wrote up a document regarding my recommendations. The recommendations will be considered as I help to improve the Image Detective page for the updated website. In addition to the website, other projects included the need for additional, and updated image collections, along with various project requests. The image collections have been used by educators in the classroom and the impact crater collection was highly requested. The glaciers collection focused mostly on South American glaciers and needed to include more of the earth's many glaciers. The collections had not been updated or created due to the fact that related imagery had not been catalogued. The process of cataloging involves identifying the center point location of the image and feature identification. Other project needs included collecting night images of India in for publishing. Again, many of the images were not catalogued and the database was lacking in night time imagery for that region. The last project was to calculate the size of mega fans in South Africa. Calculating the fan sizes involved several steps. To expedite the study, calculations needed to be made after the base maps had been created. Using data files that included an outline of the mega fans on a topographic map, I opened the file in Photoshop, determined the number of pixels within the outlined area, created a one degree squared box, determined the pixels within the box, converted the pixels within the box to kilometers, and then calculated the fan size using this information. Overall, the internship has been a learning experience for me. I have learned how to use new programs and I developed new skills. These These skills can help me as I enter into the next phase of my career. Learning Photoshop and HTML in addition to coding in Dreamweaver are highly sought after skills that are used in a variety of fields. Additionally, the exposure to different aspects of the team and working with different people helped me to gain a broader set of skills and allowed me to work with people with different experiences. The various projects I have worked on this summer have directly benefitted the team whether it was completing projects they did not have the time to do, or by helping the team reach deadlines sooner. The new website will be the best place to see all of my work as it will include the newly designed pages and will feature my updates to collections.

Coffin, Ashley↗

Evolving Applications of Aerosol Sealant

Solving the issues facing the residential building industry today requires developing not only new technologies but also new, innovative applications of existing high-performance technologies. The U.S. Department of Energy's (DOE's) Building America Program is using this approach to develop a low-cost, easy solution for whole-home envelope air sealing. Aerosol sealant is a proven method for sealing ducts quickly and efficiently, especially in hard-to-reach places. This technology, originally developed at Lawrence Berkeley National Laboratory in 1994 with support from DOE, seals 70%-90% of duct leaks (LBNL n.d.), saving $600-$850 per year in energy costs for an average household (DOE 2011). A new, innovative application takes this technology out of the duct system and into the conditioned space. Recent Building America projects have validated a method for whole-house air sealing using the aerosol sealant technology. With tighter envelope requirements in new codes and above-code programs, adequately air sealing by hand is more labor-intensive and costly than ever. The aerosol sealing method is designed to eliminate the guesswork and automate this process. The continuous measurement during installation allows builders to dial in the exact air tightness they need to meet their performance goals. Initial field testing shows that this method regularly achieves air leakage levels less than 1 ACH50 as measured by a blower door test. In a current Building America project, the Center for Energy and Environment is performing iterative field tests in 20 homes to identify the optimal integration and sequencing of this method in production home building processes to achieve a lower installation cost than traditional air sealing. This method has the potential to improve quality assurance of envelope sealing and significantly reduce labor costs compared with traditional air sealing approaches.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Simulated Time Lags of Hinode/XRT and SDO/AIA Lightcurves as an Indication of Loop Heating Scenario

The precise nature of the heating mechanism (location, duration) in coronal loops is still a matter of enormous research. We present results from a one-dimensional (1D) hydrodynamic loop simulation of a coronal loop which was run using different parameters such as loops length (50, 200 and 500 light-seconds), maximum temperature reached (3 million degrees Kelvin and 10 million degrees Kelvin), and abundances. For each scenario the model outputs were used to calculate the corresponding lightcurves as seen by X-ray telescope/Be-thin filter and various Extreme Ultra Violet Atmospheric Imaging Assembly channels. The lag time between the peak of these lightcurves was computed using cross-correlation and plotted as a function of loop length. Additional results were computed using the zero-dimensional Enthalpy-Based Thermal Evolution of Loops (EBTEL) code in order to test the compatibility of the two codes and to investigate additional loop lengths. Initial results indicate that the long (greater than 5000 seconds) lags observed in the approximately 100 light-seconds loops of active regions can only be reproduced using photospheric abundances and much longer loop lengths. This result suggests that the observed time lags cannot be completely explained by impulsive heating.

nanoflares↗

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Higher-order space-charge stability in anisotropic beams: Vlasov-Poisson derivation, refined dispersion relations, and stability charts

The Hofmann stability chart is used to screen working points in space-charge-dominated linacs. We identify two errors in its published higher-order dispersion relations: missing $(1\mp2\hatη^2/α)$ factors in the third-order $S^4$ coupling residues, and a sign error in the stated isotropic reduction of the fourth-order relation. Both corrections follow from Hofmann's Vlasov-Poisson equations without fitted parameters. They reproduce coherent tune-shift coefficients in the author's later monograph that the printed forms miss by 24% and 127%. Mode-resolved figures from a published application agree with the corrected relations and reject the printed forms, indicating an inconsistency between the 1998 equations and the calculations underlying those tested figures. We quantify the effect on the non-oscillatory stability chart. Inside the adopted $S^2\le10$ comparison domain, printed and corrected forms disagree on 0.73-2.11% of cells, with no preferred direction. Among excluded cells, disagreement reaches 22%, and the printed relation over-predicts instability at every sampled anisotropy. This concentration may help explain why the errors persisted, although it does not establish their historical cause. For PIP-II, the corrected chart flags four of thirty-two evaluable periods, including one on a third-order odd branch missed by a second-order screen. This count covers non-oscillatory modes only and remains conditional on an unresolved factor-five disagreement between two codes on transverse emittance growth.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

Primary transcripts: From the discovery of RNA processing to current concepts of gene expression ‐ Review

The main purpose of this review is to recall for investigators - and in particular students -, some of the early data and concepts in molecular genetics and biology that are rarely cited in the current literature and are thus invariably overlooked. There is a growing tendency among editors and reviewers to consider that only data produced in the last 10–20 years or so are pertinent. However this is not the case. In exact science, sound data and lucid interpretation never become obsolete, and even if forgotten, will resurface sooner or later. In the field of gene expression, covered in the present review, recent post-genomic data have indeed confirmed many of the earlier results and concepts developed in the mid-seventies, well before the start of the recombinant DNA revolution. Human brains and even the most powerful computers, have difficulty in handling and making sense of the overwhelming flow of data generated by recent high-throughput technologies. This was easier when low throughput, more integrative methods based on biochemistry and microscopy dominated biological research. Nowadays, the need for organising concepts is ever more important, otherwise the mass of available data can generate only {sup b}uilding ruins{sup }- the bricks without an architect. Concepts such as pervasive transcription of genomes, large genomic domains, full domain transcripts (FDTs) up to 100 kb long, the prevalence of post-transcriptional events in regulating eukaryotic gene expression, and the 3D-genome architecture, were all developed and discussed before 1990, and are only now coming back into vogue. Thus, to review the impact of earlier concepts on later developments in the field, I will confront former and current data and ideas, including a discussion of old and new methods. Whenever useful, I shall first briefly report post-genomic developments before addressing former results and interpretations. Equally important, some of the terms often used sloppily in scientific discussions will be clearly defined. As a basis for the ensuing discussion, some of the issues and facts related to eukaryotic gene expression will first be introduced. In chapter 2 the evolution in perception of biology over the last 60 years and the impact of the recombinant DNA revolution will be considered. Then, in chapter 3 data and theory concerning the genome, gene expression and genetics will be reviewed. The experimental and theoretical definition of the gene will be discussed before considering the 3 different types of genetic information - the {sup T}riad{sup }- and the importance of post-transcriptional regulation of gene expression in the light of the recent finding that 90% of genomic DNA seems to be transcribed. Some previous attempts to provide a conceptual framework for these observations will be recalled, in particular the {sup C}ascade Regulation Hypothesis{sup }(CRH) developed in 1967–85, and the {sup G}ene and Genon{sup }concept proposed in 2007. A knowledge of the size of primary transcripts is of prime importance, both for experimental and theoretical reasons, since these molecules represent the primary units of the {sup R}NA genome{sup }on which most of the post-transcriptional regulation of gene expression occurs. In chapter 4, I will first discuss some current post-genomic topics before summarising the discovery of the high Mr-RNA transcripts, and the investigation of their processing spanning the last 50 years. Since even today, a consensus concerning the real form of primary transcripts in eukaryotic cells has not yet been reached, I will refer to the viral and specialized cellular models which helped early on to understand the mechanisms of RNA processing and differential splicing which operate in cells and tissues. As a well-studied example of expression and regulation of a specific cellular gene in relation to differentiation and pathology, I will discuss the early and recent work on expression of the globin genes in nucleated avian erythroblasts. An important concept is that the primary transcript not only embodies protein-coding information and regulation of its expression, but also the 3D-structure of the genomic DNA from which it was derived. The wealth of recent post-genomic data published in this field emphasises the importance of a fundamental principle of genome organisation and expression that has been overlooked for years even though it was already discussed in the 1970–80ties. These issues are addressed in chapter 5 which focuses on the involvement of the nuclear matrix and nuclear architecture in DNA and RNA biology. This section will make reference to the Unified Matrix Hypothesis (UMH), which was the first molecular model of the 3D organisation of DNA and RNA. The chapter on the {sup R}NA-genome and peripheral memories{sup }discusses experimental data on the ribonucleoprotein complexes containing pre-mRNA (pre-mRNPs) and mRNA (mRNPs) which are organised in nuclear and cytoplasmic spaces respectively. Finally, {sup O}utlook {sup }will enumerate currently unresolved questions in the field, and will propose some ideas that may encourage further investigation, and comprehension of available experimental data still in need of interpretation. In chapter 8, some propositions and paradigms basic to the authors own analysis are discussed. {sup I}n conclusion{sup }the raison d'être of this review is recalled and positioned within the overall framework of scientific endeavour.

60 APPLIED LIFE SCIENCES↗

Multiphysics Simulations of MSRE with NEAMS Thermal Hydraulics Tools

This report documents the benchmarks being developed and simulations performed using tools and codes developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing MSRE experimental data. In FY23, three main work scopes were investigated under the NEAMS MSR work package at ANL. The first scope investigated the Griffin-SAM coupling model for simulating the pump startup transient experiment of MSRE. The analyses start with a simple model (single-channel, single-lattice), gradually adding more details (multi-channel, full-core) into the model. The results show that the reactivity loss curve is very sensitive to the axial boundary conditions and the radial core discretization. The simple model can predict a similar reactivity trend as that of the more sophisticated model, which is likely due to error cancellation. Accurately modeling the axial boundary condition may further improve the reactivity trend but would require significant efforts to generate the mesh of the MSRE inlet and upper plenum. The core channel radial discretization for the Griffin-SAM coupled model also depends on the flow distribution. Given the complex geometry in the inlet plenum, the flow distribution needed to be calculated from CFD analysis, which was performed using the NekRS code. This analysis employed a MSRE CAD model developed by Copenhagen Atomics. The CAD model was disassembled to keep the inlet plenum region only, which was subsequently cleaned and modified so that the mesh generated is under the memory limit. The results are merged to a few radial regions to show that the flow rate is highest in the central region. This would be useful for future improvement of the Griffin-SAM coupling model of the MSRE core. The last task investigated is tritium transport modeling using the standalone SAM code. This task aimed to initiate the effort to demonstrate and validate the tritium transport model implemented in SAM. The preliminary investigation employed an MSRE model consisting of the primary loop. Three tritium transport pathways were examined including the retention in the graphite, the permeation through the HX tube wall, and the removal from the off-gas system. The results compare well with the MSRE data, but improvements are still needed on the initial conditions (i.e., the present state may not have reached equilibrium), the boundary conditions, the off-gas system modeling, and a better numerical strategy to reach the equilibrium state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fiber motion in highly confined flows of carbon fiber and non-Newtonian polymer [Advanced modeling of non-Newtonian polymer and carbon fibers]

Inks compounded of short carbon fibers suspended in polymer resin can be extruded to produce composite materials during additive manufacturing or 3D printing processes. The flow process induces anisotropic orientation of the fibers which is set into the matrix and significantly affects the mechanical and physical properties of the final product. Therefore, the flow of fiber suspensions needs to be understood in order to predict the orientation distribution of the fibers during such manufacturing processes. There is still a lack of knowledge for extrusion of the complex mixture of a non-Newtonian polymer containing a high-volume fraction of fibers with high fiber aspect ratio, where both inter-particle, fluid-particle and fiber-wall interactions are computationally evaluated. This paper presents predictive numerical simulations and experimental results of such confined flow in a concentrated regime. The code is based on Lagrange multiplier technique and resolves each particle and interaction between the fibers and surrounding fluid and nozzle walls. We investigate numerically how the fiber length impacts the fiber alignment during extrusion. We found that a fiber length above 67% of the nozzle's diameter induces dramatic change in the fiber flow, causing fibers to concentrate at the nozzle boundaries with very low concentration of fibers in the center of the nozzle. It was found that the best fiber alignment is reached for fiber lengths equal to 40–50% of the nozzle diameters. Numerical findings are supported by experimental results. In conclusion, this work improves understanding of fiber orientation during 3D printing and is an important milestone for the prediction of the complex mechanical properties of additively manufactured fiber composites.

36 MATERIALS SCIENCE↗

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗

Conceptual study on using Doppler backscattering to measure magnetic pitch angle in tokamak plasmas

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

beam tracing↗

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗

A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos

Quasi-N-body simulations, such as FastPM, provide a fast way to simulate cosmological structure formation, but have yet to adequately include the effects of massive neutrinos. In this work, we present a method to include neutrino particles in FastPM, enabling computation of the CDM and total matter power spectra to percent-level accuracy in the non-linear regime. The CDM-neutrino cross-power can also be computed at a sufficient accuracy to constrain cosmological observables. To avoid the shot noise that typically plagues neutrino particle simulations, we employ a quasi-random algorithm to sample the relevant Fermi-Dirac distribution when setting the initial neutrino thermal velocities. We additionally develop an effective distribution function to describe a set of non-degenerate neutrinos as a single particle to speed up non-degenerate simulations. The simulation is accurate for the full range of physical interest, M ν ≲ 0.6eV, and applicable to redshifts z ≲ 2. Such accuracy can be achieved by initializing particles with the two-fluid approximation transfer functions (using the REPS package). Convergence can be reached in ~ 25 steps, with a starting redshift of z=99. Probing progressively smaller scales only requires an increase in the number of CDM particles being simulated, while the number of neutrino particles can remain fixed at a value less than or similar to the number of CDM particles. In turn, the percentage increase in runtime-per-step due to neutrino particles is between ~ 5-20% for runs with 1024 3 CDM particles, and decreases as the number of CDM particles is increased. The code has been made publicly available, providing an invaluable resource to produce fast predictions for cosmological surveys and studying reconstruction.

79 ASTRONOMY AND ASTROPHYSICS↗

Reconstruction for Liquid Argon TPC Neutrino Detectors Using Parallel Architectures [Slides]

Neutrinos are particles that interact rarely, so identifying them requires large detectors which produce lots of data. Processing this data with the computing power available is becoming more difficult as the detectors increase in size to reach their physics goals. In liquid argon time projection chambers (TPCs) the charged particles from neutrino interactions produce ionization electrons which drift in an electric field towards a series of collection wires, and the signal on the wires is used to reconstruct the interaction. The MicroBooNE detector currently collecting data at Fermilab has 8000 wires, and planned future experiments like DUNE will have 100 times more, which means that the time required to reconstruct an event will scale accordingly. Modernization of liquid argon TPC reconstruction code, including vectorization, parallelization and code portability to GPUs, will help to mitigate these challenges. The liquid argon TPC hit finding algorithm within the L arSoft framework used across multiple experiments has been vectorized and parallelized. This increases the speed of the algorithm on the order of ten times within a stand alone version on Intel architectures. This new version has been incorporated back into LarSoft so that it can be generally used. These methods will also be applied to other low level reconstruction of the wire signals such as the deconvolution. The applications and performance of this modernized liquid argon TPC wire reconstruction will be presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

BISON microstructure-based pulverization criterion in high burnup structure

To improve the economics of commercial nuclear power production, utilities are seeking to increase the allowable burnup limit of UO$_2$ fuel. One of the main factors that contributes to the current burnup limit of 62 GWd/MTU in commercial light water reactors (LWRs) is the risk of fine fragmentation or pulverization during a loss of coolant accident (LOCA). Pulverization primarily occurs at high burnups, especially when the high burnup structure (HBS) has formed. To allow the industry to pursue increased burnup and develop mitigation strategies, it is essential to have improved capability to predict the onset of pulverization. However, the mechanism of pulverization is not well understood, and the existing predictive capabilities implemented in the BISON fuel performance code are empirical in nature. In this report, mesoscale simulations are used to improve understanding of the formation mechanism of the HBS and how it responds during a LOCA transient, and inform development of a BISON pulverization criterion. A phase-field model was used to simulate the evolution of bubble pressure as a result of HBS formation. The simulations showed that gas atoms diffuse from grain interiors to the new grain boundaries created during HBS formation, and diffuse rapidly along these grain boundaries to reach existing bubbles. This causes an increase in bubble pressure in existing bubbles, leading to bubble growth during steady-state operation. To simulate the response of HBS bubbles to a LOCA transient, a newly developed phase-field model was used; in agreement with preliminary results from FY20, bubble size did not change significantly during the duration of the transient. A phase-field fracture model was used to study fragmentation patterns in the HBS, including using input from the phase-field model as initial conditions. Phase-field fracture simulations were used to determine a pulverization criterion for BISON. A function for the critical pressure for pulverization to occur was fit to data from the phase-field fracture simulations, and this function was implemented as a material property in BISON. For comparison, an analytical criterion for pulverization was developed and implemented within the same material property.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evidence for Surprising Heavy Nitrogen Isotopic Enrichment in Comet 46P/Wirtanen's Hydrogen Cyanide

46P/Wirtanen is a Jupiter-family comet, probably originating from the solar system's Kuiper Belt, that now resides on a 5.4 yr elliptical orbit. During its 2018 apparition, comet 46P passed unusually close to the Earth (within 0.08 au), presenting an outstanding opportunity for close-up observations of its inner coma. Here we present observations of HCN, H13CN, and HC15N emission from 46P using the Atacama Compact Array. The data were analyzed using the SUBLIME non-LTE radiative transfer code to derive 12C/13C and 14N/15N ratios. The HCN/H13CN ratio is found to be consistent with a lack of significant 13C fractionation, whereas the HCN/HC15N ratio of 68 ± 27 (using our most conservative 1σ uncertainties), indicates a strong enhancement in 15N compared with the solar and terrestrial values. The observed 14N/15N ratio is also significantly lower than the values of ∼140 found in previous comets, implying a strong 15N enrichment in 46P's HCN. This indicates that the nitrogen in Jupiter-family comets could reach larger isotopic enrichments than previously thought, with implications for the diversity of 14N/15N ratios imprinted into icy bodies at the birth of the solar system.

Short period comets↗