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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Numerical simulations of the ablative Rayleigh-Taylor instability in planar inertial-confinement-fusion targets using the FastRad3D code

The ablative Rayleigh-Taylor (RT) instability is a central issue in the performance of laser-accelerated inertial-confinement-fusion targets. Historically, the accurate numerical simulation of this instability has been a challenging task for many radiation hydrodynamics codes, particularly when it comes to capturing the ablatively stabilized region of the linear dispersion spectrum and modeling ab initio perturbations. Here, we present recent results from two-dimensional numerical simulations of the ablative RT instability in planar laser-ablated foils that were performed using the Eulerian code FastRad3D. Our study considers polystyrene, (cryogenic) deuterium-tritium, and beryllium target materials, quarter- and third-micron laser light, and low and high laser intensities. An initial single-mode surface perturbation is modeled in our simulations as a small modulation to the target mass density and the ablative RT growth-rate is calculated from the time history of areal-mass variations once the target reaches a steady-state acceleration. By performing a sequence of such simulations with different perturbation wavelengths, we generate a discrete dispersion spectrum for each of our examples and find that in all cases the linear RT growth-rate γ is well described by an expression of the form γ=α [kg/(1+ϵ kLm)]1/2−βkVa, where k is the perturbation wavenumber, g is the acceleration of the target, Lm is the minimum density scale-length, Va is the ablation velocity, and ϵ is either one or zero. The dimensionless coefficients α and β in the above formula depend on the particular target and laser parameters and are determined from two-dimensional simulation results through the use of a nonlinear curve-fitting procedure. While our findings are generally consistent with those of Betti et al. (Phys. Plasmas 5, 1446 (1998)), the ablative RT growth-rates predicted in this investigation are somewhat smaller than the values previously reported for the same target and laser parameters. It is speculated that differences in the equation-of-state and opacity models are largely responsible for the discrepancy. Resolution of this issue awaits the development of better experimental diagnostics capable of measuring small-wavelength (5–20 μm) perturbation growth due to the ablative RT instability in the linear regime.

Bates, J. W. (ORCID:0000000188087240)↗
DOE OSTI2013

A maximum likelihood approach to estimating correlation functions

We define a maximum likelihood (ML for short) estimator for the correlation function, ξ, that uses the same pair counting observables (D, R, DD, DR, RR) as the standard Landy and Szalay (LS for short) estimator. The ML estimator outperforms the LS estimator in that it results in smaller measurement errors at any fixed random point density. Put another way, the ML estimator can reach the same precision as the LS estimator with a significantly smaller random point catalog. Moreover, these gains are achieved without significantly increasing the computational requirements for estimating ξ. We quantify the relative improvement of the ML estimator over the LS estimator and discuss the regimes under which these improvements are most significant. We present a short guide on how to implement the ML estimator and emphasize that the code alterations required to switch from an LS to an ML estimator are minimal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Aligning thermal and current quenches with a high density low-Z injection

The conventional approach for thermal quench (TQ) mitigation in a tokamak disruption is through a high- Z impurity injection that radiates away the plasma’s thermal energy before it reaches the wall. The downside is a robust Ohmic-to-runaway current conversion due to the radiatively clamped low post-thermal-quench electron temperature. An alternative approach is to deploy a low- Z (either deuterium or hydrogen) injection that aims to slow down the TQ, and ideally aligns it with the current quench (CQ). This approach has been investigated here via 3D MHD simulations using the PIXIE3D code. By boosting the hydrogen density, a fusion-grade plasma is dilutionally cooled at approximately the original pressure. Energy loss to the wall is controlled by a Bohm outflow condition at the boundary where the magnetic field intercepts a thin plasma sheath at the wall, in addition to Bremsstrahlung bulk losses. Robust MHD instabilities proceed as usual, while the collisionality of the plasma has been greatly increased and parallel transport is now in the Braginskii regime. The main conclusion of this study is that the decreased transport loss along open field lines due to a sufficient low- Z injection slows down the TQ rate to the order of 20 ms, aligned with the CQ timescale for a 15 MA ITER plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Two-photon double ionization with finite pulses: Application of the virtual sequential model to helium

As a step toward the full ab initio description of two-photon double ionization processes, we present a finite-pulse version of the virtual-sequential model for polyelectronic atoms. The model relies on the ab initio description of the single ionization scattering states of both the neutral and ionized target system. As a proof of principle and a benchmark, the model is applied to the helium atom using the newstock atomic photoionization code. The results of angularly integrated observables, which are in excellent agreement with existing time-dependent Schrödinger equation (TDSE) simulations, show how the model is able to capture the role of electron correlation in the nonsequential regime, and the influence of autoionizing states in the sequential regime, at a comparatively modest computational cost. The model also reproduces the two-particle interference with ultrashort pulses, which is within reach of current experimental technologies. Furthermore, the model shows the modulation of the joint energy distribution in the vicinity of autoionizing states, which can be probed with extreme-ultraviolet pulses of duration much longer than the characteristic lifetime of the resonance. Furthermore, the formalism discussed here applies also to polyelectronic atoms and molecules, thus opening a window on nonsequential and sequential double ionization in these more complex systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Monte Carlo Driven MDI Optimization at a Muon Collider

A Muon Collider represents a very interesting possibility for a future machine to explore the energy frontier in particle physics. However, to reach the needed luminosity, beam intensities of the order of 10⁹–10¹² muons per bunch are needed. In this context, the Beam-Induced Background must be taken into account for its effects on magnets and detector. Several mitigation strategies can however be conceived. In this view, it is of crucial importance to develop a flexible tool that allows to easily reconstruct the machine geometry in a Monte Carlo code, allowing to simulate in detail the interaction of muon decay products in the machine, while being able to change the machine optics itself to find the best configuration. In this contribution, a possible approach to such a purpose is presented, based on FLUKA for the Monte Carlo simulation and on LineBuilder for the geometry reconstruction. Results based on the 1.5 TeV machine optics developed by the MAP collaboration are discussed, as well as a first approach to possible mitigation strategies.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Uniaxial Compression of Cellular Materials at a 10 -1 s -1 Strain Rate Simultaneously with Synchrotron X-ray Computed Tomographic Imaging

The topic is presented as a series of slides. Motivation for the work included the following: X-ray tomography is a fantastic technique for characterizing a material’s starting structure as well as for non-destructive, in situ experiments to investigate material response; 3D X-ray tomography is needed to fully characterize the morphology of cellular materials; and synchrotron micro-CT can capture 3D images without pausing experiment. Among the conclusions reached are these: High-rate radiographic and tomographic imaging (0.25 s 3D frame rate) using synchrotron CT can capture full 3D images of hyper-elastic materials at a 10 -2 strain rate; dynamic true in situ uniaxial loading can be accurately captured; the three stages of compression can be imaged: bending, buckling, and breaking; implementation of linear modeling is completed; meshes have been imported into LANL modeling codes--testing and validation is underway and direct comparison and validation between in situ data and modeled mechanical response is possible.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quantifying Moss Response to Metal Contaminant Exposure Using Laser-Induced Fluorescence

Tracing sources of contamination, including potentially toxic elements (PTEs), has historically been achieved through sampling and analysis of soil or biota, which are labor-intensive, costly, and destructive methods. Thus, availability of a non-destructive in situ remote sensing method for monitoring metals deposited in biota is of great interest. Laser-induced fluorescence (LIF) is an emerging spectroscopic and imaging technique that documents changes in molecular energy level in plants as a biological response to metal contamination. For a proof-of-concept study and preliminary experiment, moss was selected for experimentation due to its long history of use in tracing atmospheric deposition of PTEs. Consecutive treatments of copper chloride (CuCl2) were administered to three moss samples, simulating wet deposition every 48 h over 10 days until reaching cumulative Cu concentrations of 2.690 to 8.075 μmol/cm2. While these Cu amounts are above environmentally relevant concentrations, they allowed the best conditions for testing and fine tuning of the imaging and data processing protocols presented in this paper. Moss fluorescence was induced using both 532 nm green and 355 nm UV lasers. A CMOS camera captured images of the LIF response, and red–green–blue (RGB) decimal code values were extracted for each pixel in the images, and pixel densities of color channels from treated and untreated moss samples were compared. Results show a shift towards lower color decimal codes corresponding to increased Cu concentration. We developed and contrasted multiple quantitative analyses of color distributions and demonstrated that LIF shows great promise for remote sensing of Cu accumulation in moss at μmol/cm2 levels. Though currently, the method would be limited to highly toxic sites, it illustrates the possibility and provides a framework for development of higher-sensitivity methods to detect nmol/cm2 that are viable for urban contamination level monitoring.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2014

SU-E-J-53: Dosimetric Evaluation at Volumetric Modulated Arc Therapy for Treatment of Prostate Cancer Using Single Or Double Arcs

Purpose: Evaluate and compare retrospective prostate treatment plan using Volumetric Modulated Arc Therapy (RapidArc™ - Varian) technique with single or double arcs at COI Group. Methods: Ten patients with present prostate and seminal vesicle neoplasia were replanned as a target treatment volume and a prescribed dose of 78 Gy. A baseline planning, using single arc, was developed for each case reaching for the best result on PTV, in order to minimize the dose on organs at risk (OAR). Maintaining the same optimization objectives used on baseline plan, two copies for optimizing single and double arcs, have been developed. The plans were performed with 10 MV photon beam energy on Eclipse software, version 11.0, making use of Trilogy linear accelerator with Millenium HD120 multileaf collimator. Comparisons on PTV have been performed, such as: maximum, minimum and mean dose, gradient dose, as well as the quantity of monitor units, treatment time and homogeneity and conformity index. OARs constrains dose have been evaluated, comparing both optimizations. Results: Regarding PTV coverage, the difference of the minimum, maximum and mean dose were 1.28%, 0.7% and 0.2% respectively higher for single arc. When analyzed the index of homogeneity found a difference of 0.99% higher when compared with double arcs. However homogeneity index was 0.97% lower on average by using single arc. The doses on the OARs, in both cases, were in compliance to the recommended limits RTOG 0415. With the use of single arc, the quantity of monitor units was 10,1% lower, as well as the Beam-On time, 41,78%, when comparing double arcs, respectively. Conclusion: Concerning the optimization of patients with present prostate and seminal vesicle neoplasia, the use of single arc reaches similar objectives, when compared to double arcs, in order to decrease the treatment time and the quantity of monitor units.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

DaMaSCUS: the impact of underground scatterings on direct detection of light dark matter

Conventional dark matter direct detection experiments set stringent constraints on dark matter by looking for elastic scattering events between dark matter particles and nuclei in underground detectors. However these constraints weaken significantly in the sub-GeV mass region, simply because light dark matter does not have enough energy to trigger detectors regardless of the dark matter-nucleon scattering cross section. Even if future experiments lower their energy thresholds, they will still be blind to parameter space where dark matter particles interact with nuclei strongly enough that they lose enough energy and become unable to cause a signal above the experimental threshold by the time they reach the underground detector. Therefore in case dark matter is in the sub-GeV region and strongly interacting, possible underground scatterings of dark matter with terrestrial nuclei must be taken into account because they affect significantly the recoil spectra and event rates, regardless of whether the experiment probes DM via DM-nucleus or DM-electron interaction. To quantify this effect we present the publicly available Dark Matter Simulation Code for Underground Scatterings (DaMaSCUS), a Monte Carlo simulator of DM trajectories through the Earth taking underground scatterings into account. Our simulation allows the precise calculation of the density and velocity distribution of dark matter at any detector of given depth and location on Earth. The simulation can also provide the accurate recoil spectrum in underground detectors as well as the phase and amplitude of the diurnal modulation caused by this shadowing effect of the Earth, ultimately relating the modulations expected in different detectors, which is important to decisively conclude if a diurnal modulation is due to dark matter or an irrelevant background.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Enhancement of ECCD by the current condensation effect for stabilizing large magnetic islands caused by neoclassical tearing modes in tokamak plasmas

The radio frequency current condensation effect reported in Reiman and Fisch (2018 Phys. Rev. Lett. 121 225001) is modeled in the nonlinear resistive magnetohydrodynamic code. A series of numerical investigations have been performed to investigate the enhancement of electron cyclotron current drive (ECCD) by the current condensation effect during the control of neoclassical tearing mode (NTM) in tokamak plasmas. In the numerical model, both the parallel transport and the perpendicular transport of electron temperature are considered. The EC driven current and driven perturbed electron temperature can nonlinearly evolve within the given magnetic configuration and eventually reach saturation states. The input power threshold of ECCD and the fold bifurcation phenomenon are numerically verified via nonlinear simulations. The numerical results show good agreements with the analytical results. Moreover, spatial distributions of EC current for the two solutions at different condensed level are displayed. The control effectiveness of ECCD for large NTM islands has been evaluated while considering the current condensation effect. While taking into account current condensation effect, for a sufficiently large input power, a larger island can be more effectively stabilized than a smaller one, which suggests a reassessment of the previous idea that the ECCD should always be turned on as early as possible. The potential physics mechanism behind the ECCD control have all been discussed in detail. Furthermore, the condensation effect is found to have favorable effects on the radial misalignment of ECCD. In the consideration of the situation for extremely localized control needs, a highly peaked heating profile is adopted to verify that the fold bifurcation phenomenon still exists and the current condensation effect can still take effect in this extreme condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Divertor detachment characterization in negative triangularity discharges in DIII-D via 2D edge-plasma transport modeling

Edge fluid modeling of the first divertor-plasma detachment experiments in negative triangularity (NT) discharges in the DIII-D tokamak is presented using the 2D multifluid edge transport code UEDGE, including cross-field particle drifts. Experiments were performed where the lower single-null magnetic equilibrium had a strong NT (δ≈−0.5), that is, where the magnetic X-point is at a larger major radius than the core magnetic axis. Divertor-plasma detachment was induced by increasing the core plasma density in DIII-D via intrinsic gas puffing. Here density scans are performed with UEDGE to reach a detached plasma and to quantitatively recover the experimental rollover of the ion saturation current on the outer divertor target plate. The simulations cover experiments with both signs of the toroidal magnetic field, B T , where the ion magnetic Grad-B drifts are directed into (forward B T ) and out of (reverse B T ) the divertor region. Consistent with experiments with neutral beam power injection, the NT simulations reproduce: 40% higher density is needed to reach detachment onset with forward B T compared with reverse B T , and the absence of deep detachment is found with reverse B T . Similarly, comparison between Ohmic discharges in NT and positive triangularity (PT) shaping confirms that a substantially higher density is needed to achieve detachment in NT than in PT, with NT requiring an line-average density of at least the Greenwald density or higher. Simulation results suggest that higher densities are needed to reach detachment in negative compared to PT because these discharges have a shorter midplane-to-target distance along the total magnetic field B, a shorter outer divertor poloidal leg length (0.06 m vs 0.2 m), and reduced radial transport [near-scrape-off layer (SOL) D ⊥ /χ ⊥ =0.3/0.5 vs D ⊥ /χ ⊥ =1.0/1.0, all in (m 2 /s)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantifying deformation during Zry-4 burst testing: a comparison of BISON and a combined in-situ digital image correlation and infrared thermography method

The development and verification of nuclear fuel-cladding performance codes require extensive testing to establish empirical correlations, necessitating more accelerated testing methods that can replace large validation studies. Here, in the present work, a technique was developed and implemented to experimentally quantify the relationship between temperature, cladding internal pressure, and strain in-situ during a burst test, with the specific aim of generating data for code validation and model refinement. Digital image correlation was used to measure hoop, axial, and radial strain, and infrared thermography to quantify axial temperature gradients. Several key experimental modifications to traditional LOCA burst testing were necessary, but are shown to effectively have little impact on burst temperatures. The time/temperature dependent pressure and time/spatially dependent thermal gradient data were used with BISON to simulate cladding burst and strain rates. Although the measured DIC strains reach appreciable amounts at slightly lower temperatures than BISON simulated strains, there is good agreement between the measured and simulated strains in terms of magnitude and rates leading up to burst. However, in the few seconds prior to and during burst, the BISON simulated strains are significantly lower than the measured strains, indicating the potential for model improvement. The combined experiment and simulation technique may be applied to any new cladding concept to accelerate fuel qualification.

36 MATERIALS SCIENCE↗
DOE OSTI2016

High current DC negative ion source for cyclotron

A filament driven multi-cusp negative ion source has been developed for proton cyclotrons in medical applications. In Cs-free operation, continuous H{sup −} beam of 10 mA and D{sup −} beam of 3.3 mA were obtained stably at an arc-discharge power of 3 kW and 2.4 kW, respectively. In Cs-seeded operation, H{sup −} beam current reached 22 mA at a lower arc power of 2.6 kW with less co-extracted electron current. The optimum gas flow rate, which gives the highest H{sup −} current, was 15 sccm in the Cs-free operation, while it decreased to 4 sccm in the Cs-seeded operation. The relationship between H{sup −} production and the design/operating parameters has been also investigated by a numerical study with KEIO-MARC code, which gives a reasonable explanation to the experimental results of the H{sup −} current dependence on the arc power.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Tunnel ionization within a one-dimensional, undriven plasma sheath

In high density, high temperature plasmas, the plasma sheath that develops can result in extremely high electric fields, on the order of tens to hundreds of V/nm. Under the right conditions, these electric fields can reach magnitudes that can increase the probability of electron tunneling ionization to occur, resulting in one or more electron-ion pairs. The presence of tunneling ionization can then modify the development of the plasma sheath, as well as properties such as the ion and electron densities and plasma potential. The tunnel ionization process for hydrogen atoms is demonstrated, in this work, as implemented in a Sandia National Laboratories, particle-in-cell code Aleph. Results are presented for the application of the tunnel ionization process to a one-dimensional, undriven plasma sheath. Additional results for cases that consider warm ions and neutrals, the inclusion of electron–neutral collisions, and the injection of neutral particles, as well as the application to various plasma devices, will be discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

SU-E-T-439: An Improved Formula of Scatter-To-Primary Ratio for Photon Dose Calculation

Purpose: Scatter-to-primary ratio (SPR) is an important dosimetric quantity that describes the contribution from the scatter photons in an external photon beam. The purpose of this study is to develop an improved analytical formula to describe SPR as a function of circular field size (r) and depth (d) using Monte Carlo (MC) simulation. Methods: MC simulation was performed for Mohan photon spectra (Co-60, 4, 6, 10, 15, 23 MV) using EGSNRC code. Point-spread scatter dose kernels in water are generated. The scatter-to-primary ratio (SPR) is also calculated using MC simulation as a function of field size for circular field size with radius r and depth d. The doses from forward scatter and backscatter photons are calculated using a convolution of the point-spread scatter dose kernel and by accounting for scatter photons contributing to dose before (z' d) reaching the depth of interest, d, where z' is the location of scatter photons, respectively. The depth dependence of the ratio of the forward scatter and backscatter doses is determined as a function of depth and field size. Results: We are able to improve the existing 3-parameter (a, w, d0) empirical formula for SPR by introducing depth dependence for one of the parameter d0, which becomes 0 for deeper depths. The depth dependence of d0 can be directly calculated as a ratio of backscatter-to-forward scatter doses for otherwise the same field and depth. With the improved empirical formula, we can fit SPR for all megavoltage photon beams to within 2%. Existing 3-parameter formula cannot fit SPR data for Co-60 to better than 3.1%. Conclusion: An improved empirical formula is developed to fit SPR for all megavoltage photon energies to within 2%.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

EMPIRE-PIC Code Verification of a Cold Diode

This report presents the code verification of EMPIRE-PIC to the analytic solution to a cold diode which was first derived by Jaffe. The cold diode was simulated using EMPIRE-PIC and the error norms were computed based on the Jaffe solution. The diode geometry is one-dimensional and uses the EMPIRE electrostatic field solver. After a transient start-up phase as the electrons first cross the anode-cathode gap, the simulations reach an equilibrium where the electric potential and electric field are approximately steady. The expected spatial order of convergence for potential, electric field and particle velocity are observed.

42 ENGINEERING↗
DOE OSTI2016

DUST CAPTURE AND LONG-LIVED DENSITY ENHANCEMENTS TRIGGERED BY VORTICES IN 2D PROTOPLANETARY DISKS

We study dust capture by vortices and its long-term consequences in global two-fluid inviscid disk simulations using a new polar grid code RoSSBi. We perform the longest integrations so far, several hundred disk orbits, at the highest resolution attainable in global disk simulations with dust, namely, 2048 × 4096 grid points. We vary a wide range of dust parameters, most notably the initial dust-to-gas ratio ϵ varies in the range of 10{sup −4}–10{sup −2}. Irrespective of the value of ϵ , we find rapid concentration of the dust inside vortices, reaching dust-to-gas ratios of the order of unity inside the vortex. We present an analytical model that describes this dust capture process very well, finding consistent results for all dust parameters. A vortex streaming instability develops, which invariably causes vortex destruction. After vortex dissipation large-scale dust rings encompassing a disk annulus form in most cases, which sustain very high dust concentration, approaching ratios of the order of unity; they persist as long as the duration of the simulations. They are sustained by a streaming instability, which manifests itself in high-density dust clumps at various scales. When vortices are particularly long-lived, rings do not form but dust clumps inside vortices can survive a long time and would likely undergo collapse by gravitational instability. Rings encompass almost an Earth mass of solid material, while even larger masses of dust do accumulate inside vortices in the earlier stage. We argue that rapid planetesimal formation would occur in the dust clumps inside the vortices as well as in the post-vortex rings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Integrated Decision Methodology for Final Disposal of Radioactive Sources

Radioactive sources are used throughout the world, and their beneficial uses are many. Sources are at the core of many applications in different branches of industry, medicine, agriculture and research. Sources are present in a very wide range of equipment that is used for, inter alia, cancer treatment, killing bacteria in food, sterilizing medical supplies, measuring instruments such as gauges used to measure soil moisture and soil density and many other components, irradiating seeds for enhancing food production, protecting buildings from lightning strikes, mapping underground sources of water, prospecting for oil and gas reserves, measuring density of soil for construction projects, or even detecting smoke. At their “cradle” and when in use, radioactive sources are usually properly managed and controlled. It is when a source has reached the end of its useful life when they are at a higher risk of being under inadequate control, poorly managed, or even becoming orphan. Sources that have reached the end of their useful life must be carefully and safely managed. The “grave” of disused sources needs a proper disposal process to prevent them from posing a potential threat to people or the environment. To protect the public and the environment from the potential hazards of ionizing radiation, and to prevent disused sources from becoming orphan, a “cradle-to-grave” control of radioactive sources is essential, as promoted by the Code of Conduct on the Safety and Security of Radioactive Sources. This approach requires a national policy and strategy, an adequate legal and regulatory framework, and adequate resources and infrastructure. Many IAEAMember States have expressed the need for effective solutions to this challenge. There is also a gap related to the life cycle costs of potential disposal practices. This paper proposes a multi-attribute decision making methodology to compare disposal options using both economic and non-economic criteria leading to a preferred option for disposal and a strategy for final disposal of disused radioactive sources that is specific and optimized for each IAEA Member State.

Beatty, Randy↗