SEARCH · Search NASA

Results for “Simulated flux procedure”

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.

149 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2004

Role of Turbulent Prandtl Number on Heat Flux at Hypersonic Mach Number

Present simulation of turbulent flows involving shock wave/boundary layer interaction invariably overestimates heat flux by almost a factor of two. One possible reason for such a performance is a result of the fact that the turbulence models employed make use of Morkovin's hypothesis. This hypothesis is valid for non-hypersonic Mach numbers and moderate rates of heat transfer. At hypersonic Mach numbers, high rates of heat transfer exist in regions where shock wave/boundary layer interactions are important. As a result, one should not expect traditional turbulence models to yield accurate results. The goal of this investigation is to explore the role of a variable Prandtl number formulation in predicting heat flux in flows dominated by strong shock wave/boundary layer interactions. The intended applications involve external flows in the absence of combustion such as those encountered in supersonic inlets. This can be achieved by adding equations for the temperature variance and its dissipation rate. Such equations can be derived from the exact Navier-Stokes equations. Traditionally, modeled equations are based on the low speed energy equation where the pressure gradient term and the term responsible for energy dissipation are ignored. It is clear that such assumptions are not valid for hypersonic flows. The approach used here is based on the procedure used in deriving the k-zeta model, in which the exact equations that governed k, the variance of velocity, and zeta, the variance of vorticity, were derived and modeled. For the variable turbulent Prandtl number, the exact equations that govern the temperature variance and its dissipation rate are derived and modeled term by term. The resulting set of equations are free of damping and wall functions and are coordinate-system independent. Moreover, modeled correlations are tensorially consistent and invariant under Galilean transformation. The final set of equations will be given in the paper.

Xiao, X.↗

DOE OSTI2025

Deep potential molecular dynamics simulations of ion-enhanced etching of silicon by atomic chlorine

The continued development of plasma-assisted processing techniques requires a fundamental understanding of plasma-surface interactions. Molecular dynamics (MD) simulations have been employed to complement experimental studies and better understand the properties of such systems. Recently, machine learning (ML) methods have enabled the development of ab initio-based interatomic potentials, which can be generalized to complex combinations of multiple atom types. In this work, we use ML potentials developed using the Deep Potential Molecular Dynamics (DeepMD) framework to provide a model of ion-enhanced etching of Si by Cl atoms. We demonstrate the importance of proper selection of the training data set to the accuracy of the DeepMD model and compare our results to MD results using empirical potentials, as well as to experimental measurements. Exposure of undoped Si at 300 K to thermal Cl atoms yields a steady-state Cl coverage of 1.25 monolayers, which is slightly lower than the value obtained in previous experimental studies. Predictions of Si etch yields by simultaneous Cl atom and Ar + ion impacts as a function of ion energy, neutral to ion flux ratio, and angle of incidence of the ions are in reasonably good agreement with classical MD results and experimental measurements. Finally, etch yields and SiCl x mixed layer thicknesses during simultaneous bombardment of the Si(100) surface by Cl atoms and Cl + ions are in good agreement with experimental data. In conclusion, the present work is a necessary condition for the extension of the DeepMD procedure to more complex systems of interest in plasma-surface interactions.

Artificial neural networks↗

DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Revisions of ORNL 188 W Process Based on Nonradiological Experiments

Tungsten-188 is in widespread use in the 188 W(t 1/2 = 69 d)/ 188 Re(t 1/2 =16.9 h) biomedical generator. Oak Ridge National Laboratory (ORNL) has been providing this product to the world since 1999. At ORNL, 188 W is produced via irradiation in ORNL’s High Flux Isotope Reactor (HFIR). Enriched 186 W targets in the form of sintered metallic pellets or rings achieve a compact loading in the irradiation vessel, providing a high yield per unit target. The enrichment of the target is >90% 186 W and this isotope undergoes double neutron capture to produce the desired 188W product. While 188 W is produced by neutron bombardment, 191 Os(t 1/2 = 15.4 d) is simultaneously produced as a byproduct requiring separation from 188 W by post irradiation treatment. In the current processing pathway, the irradiated W metal pellets or rings are first converted into an oxide form of WO 3 by heating the irradiated W metal target at 750°C in a quartz reaction vessel inside a vertical furnace under a constant flow of air. During heating, W metal reacts with oxygen in the air to produce WO 3 , which is soluble in 6 M NaOH for later purification process. This oxidation process also converts 188 Os (the decay daughter of 188 W) and 191 Os (the irradiation produced byproduct) into OsO 4 , a highly volatile and toxic gas. The gaseous effluents driven from the quartz reaction vessel are passed through a scrubbing array to remove OsO 4 before the air is discharged from the process. This heterogeneous oxidation method simultaneously achieves two goals: 1). converting metal target to a soluble oxide form and 2). separating volatile OsO 4 away from the solid WO 3 product, although harmful 191 OsO 4 is unfavorable but being taken care of by the down road scrubbing array. In the past twenty years the existing OsO 4 scrubbing array served well in preventing OsO 4 from being released into the environment, until July of 2020 when a minor amount of 191 Os was found to have been released into the environment which resulted in a standby of building 4501 from July 2020 to February 2021.Since October 20 of 2020 a team of researchers from groups of RSTD and NEFD were organized to work for the following tasks: (1) Fully understand the root causes of the July incident of 191 Os release and create a redesigned 191 Os scrubbing array to eliminate the possibility of 191 Os release in future 188 W processes; (2) In Phase-1 tests, select a correct air (O 2 ) flow rate in 188 W process, that provides sufficient oxygen for the oxidation of tungsten metal, while still allowing sufficient resonance time for 191 Os removal by the scrubbers; (3) In Phase-2 tests, confirm the efficiency of the redesigned scrubbing array to absorb excessive amount of non-rad OsO 4 (generated from Os metal powder) at the selected air flow; (4) In Phase-3 tests, simulate the hot cell 188 W process by heating pressed non-rad W-Os pellets in a quartz vessel of new design and the new array under operation parameters selected in Phase-1 and -2; (5) Complete a new operation procedure for 188 W process in hot cell and a TM report as a summary of the Os Mitigation Project. This TM report summarizes the investigation on root causes from technical aspects in July 191 Os release and the correspondent improvement towards the redesigned 191 Os scrubbing system. The report presents the efforts made on selection of optimal operational parameter with the new scrubbing array and the quartz vessel of new design, based on experiment data obtained in Phase-1, -2 and -3 tests. Further renovations for 188 W process and additional improvement of W target treatment will also be discussed.

07 ISOTOPE AND RADIATION SOURCES↗