SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2022

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↗

DOE OSTI2022

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↗

DOE OSTI2019

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↗

DOE OSTI2020

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↗

DOE OSTI2026

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↗

DOE OSTI2018

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↗

DOE OSTI2023

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↗

DOE OSTI2016

OpenMP Parallelization and Optimization of Graph-Based Machine Learning Algorithms

In this paper, we investigate the OpenMP parallelization and optimization of two novel data classification algorithms. The new algorithms are based on graph and PDE solution techniques and provide significant accuracy and performance advantages over traditional data classification algorithms in serial mode. The methods leverage the Nystrom extension to calculate eigenvalue/eigenvectors of the graph Laplacian and this is a self-contained module that can be used in conjunction with other graph-Laplacian based methods such as spectral clustering. We use performance tools to collect the hotspots and memory access of the serial codes and use OpenMP as the parallelization language to parallelize the most time-consuming parts. Where possible, we also use library routines. We then optimize the OpenMP implementations and detail the performance on traditional supercomputer nodes (in our case a Cray XC30), and test the optimization steps on emerging testbed systems based on Intel’s Knights Corner and Landing processors. We show both performance improvement and strong scaling behavior. Finally, a large number of optimization techniques and analyses are necessary before the algorithm reaches almost ideal scaling.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

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↗

DOE OSTI2025

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↗

DOE OSTI2015

Simulation of multipactor on the rectangular grooved dielectric surface

Multipactor discharge on the rectangular grooved dielectric surface is simulated self-consistently by using a two-and-a-half dimensional (2.5 D) electrostatic particle-in-cell (PIC) code. Compared with the electromagnetic PIC code, the former can give much more accurate solution for the space charge field caused by the multipactor electrons and the deposited surface charge. According to the rectangular groove width and height, the multipactor can be divided into four models, the spatial distributions of the multipactor electrons and the space charge fields are presented for these models. It shows that the rectangular groove in different models gives very different suppression effect on the multipactor, effective and efficient suppression on the multipactor can only be reached with a proper groove size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Alpha particles diffusion due to charge changes

Alpha particles diffusion due to charge changes in a magnetized plasma is studied. Analytical calculations and numerical simulations are employed to show that this process can be very important in the pedestal-edge-SOL regions. This is the first study that presents clear evidence of the importance of atomic processes on the diffusion of alpha particles. A simple 1D model that includes inelastic collisions with plasma species, “cold” neutrals, and partially ionized species was employed. The code, which follows the exact particle orbits and includes the effect of inelastic collisions via a Monte Carlo type random process, runs on a graphic processor unit (GPU). The analytical and numerical results show excellent agreement when a uniform background (plasma and cold species) is assumed. The simulations also show that the gradients in the density of the plasma and cold species, which are large and opposite in the edge region, produce an inward flux of alpha particles. Calculations of the alpha particles flux reaching the walls or divertor plates should include these processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

New approach for precise computation of Lyman-α forest power spectrum with hydrodynamical simulations

Current experiments are providing measurements of the flux power spectrum from the Lyman-α forests observed in quasar spectra with unprecedented accuracy. Their interpretation in terms of cosmological constraints requires specific simulations of at least equivalent precision. In this paper, we present a suite of cosmological N-body simulations with cold dark matter and baryons, specifically aiming at modeling the low-density regions of the inter-galactic medium as probed by the Lyman-α forests at high redshift. The simulations were run using the GADGET-3 code and were designed to match the requirements imposed by the quality of the current SDSS-III/BOSS or forthcoming SDSS-IV/eBOSS data. They are made using either 2 × 768{sup 3} ≅ 1 billion or 2 × 192{sup 3} ≅ 14 million particles, spanning volumes ranging from (25 Mpc h{sup −1}){sup 3} for high-resolution simulations to (100 Mpc h{sup −1}){sup 3} for large-volume ones. Using a splicing technique, the resolution is further enhanced to reach the equivalent of simulations with 2 × 3072{sup 3} ≅ 58 billion particles in a (100 Mpc h{sup −1}){sup 3} box size, i.e. a mean mass per gas particle of 1.2 × 10{sup 5}M{sub ⊙} h{sup −1}. We show that the resulting power spectrum is accurate at the 2% level over the full range from a few Mpc to several tens of Mpc. We explore the effect on the one-dimensional transmitted-flux power spectrum of four cosmological parameters (n{sub s}, σ{sub 8}, Ω{sub m} and H{sub 0}) and two astrophysical parameters (T{sub 0} and γ) that are related to the heating rate of the intergalactic medium. By varying the input parameters around a central model chosen to be in agreement with the latest Planck results, we built a grid of simulations that allows the study of the impact on the flux power spectrum of these six relevant parameters. We improve upon previous studies by not only measuring the effect of each parameter individually, but also probing the impact of the simultaneous variation of each pair of parameters. We thus provide a full second-order expansion, including cross-terms, around our central model. We check the validity of the second-order expansion with independent simulations obtained either with different cosmological parameters or different seeds. Finally, a comparison to the one-dimensional Lyman-α forest power spectrum obtained with BOSS by [1] shows an excellent agreement.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

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↗

DOE OSTI2022

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↗

DOE OSTI2014

Irradiation instability at the inner edges of accretion disks

An instability can potentially operate in highly irradiated disks where the disk sharply transitions from being radially transparent to opaque (the 'transition region'). Such conditions may exist at the inner edges of transitional disks around T Tauri stars and accretion disks around active galactic nuclei. We derive the criterion for this instability, which we term the 'irradiation instability', or IRI. We also present the linear growth rate as a function of β, the ratio between radiation force and gravity, and c{sub s}, the sound speed of the disk, obtained using two methods: a semi-analytic analysis of the linearized equations and a numerical simulation using the GPU-accelerated hydrodynamical code PEnGUIn. In particular, we find that IRI occurs at β ∼ 0.1 if the transition region extends as wide as ∼0.05r, and at higher β values if it is wider. This threshold value applies to c{sub s} ranging from 3% of the Keplerian orbital speed to 5%, and becomes higher if c{sub s} is lower. Furthermore, in the nonlinear evolution of the instability, disks with a large β and small c{sub s} exhibit 'clumping', extreme local surface density enhancements that can reach over 10 times the initial disk surface density.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

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↗

DOE OSTI2015

SU-E-T-350: Effective Point of Measurement and Total Perturbation Correction P for Parallel-Plate Ion Chambers in High-Energy Photon Beams

Purpose: This paper aims to determine the effective point of measurement and the total perturbation correction p of parallel-plate chambers for clinical photon dosimetry. Methods: The effective point of measurement (EPOM) was calculated using the EGSnrc Monte Carlo code system with the EGSnrc user code egs- chamber. Depth dose curves of the ionization chambers were calculated in a water phantom for several high energy photon spectra (4, 6, 10, 15, 18 MV-X). Different normalization criterions (normalization to the maximum of the depth dose curve and normalization to the value in 10 cm depth) have been applied. The EPOM was determined by shifting the normalized depth dose curve of a small water voxel against the depth ionization curve until the disagreement (calculated by the root mean square deviation) reaches a minimum. In addition, the total perturbation correction p was calculated by the ratio of the dose to water and the product of the dose determined in the chamber and the water to air stopping power ratio. Results: The EPOM varied slightly depending on the chosen normalization criterion. For all chambers the necessary shift of the EPOM decreased linearly with increasing beam quality specifier TPR{sub 20/10}. For the Roos and NACP chamber, the results were positive suggesting that the chambers need to be shifted towards the focus. For the Markus chamber, the required shift was negative and for the Advanced Markus chamber partly negative and partly positive. The total perturbation correction p was almost independent of the depth. Only for regions below 1 cm the perturbation correction deviated significantly from unity. Conclusion: In the present study, the effective point of measurement and the total perturbation correction p was determined for four parallel-plate ionization chambers and five clinical relevant photon spectra. Applying the calculated EPOM, the residual perturbation correction p was mostly depth independent.

60 APPLIED LIFE SCIENCES↗