Search NASA⌕ Search

SEARCH · Search NASA

Results for “Beam code development”

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 163 records · Page 9

A Comparison of Model Calculation and Measurement of Absorbed Dose for Proton Irradiation

With the increase in the amount of time spent EVA that is necessary to complete the construction and subsequent maintenance of ISS, it will become increasingly important for ground support personnel to accurately characterize the radiation exposures incurred by EVA crewmembers. Since exposure measurements cannot be taken within the organs of interest, it is necessary to estimate these exposures by calculation. To validate the methods and tools used to develop these estimates, it is necessary to model experiments performed in a controlled environment. This work is such an effort. A human phantom was outfitted with detector equipment and then placed in American EMU and Orlan-M EVA space suits. The suited phantom was irradiated at the LLUPTF with proton beams of known energies. Absorbed dose measurements were made by the spaceflight operational dosimetrist from JSC at multiple sites in the skin, eye, brain, stomach, and small intestine locations in the phantom. These exposures are then modeled using the BRYNTRN radiation transport code developed at the NASA Langley Research Center, and the CAM (computerized anatomical male) human geometry model of Billings and Yucker. Comparisons of absorbed dose calculations with measurements show excellent agreement. This suggests that there is reason to be confident in the ability of both the transport code and the human body model to estimate proton exposure in ground-based laboratory experiments.

Zapp, N.↗

Single-Crystal X-Ray Spectropolarimeter

X-ray spectropolarimetry is a diagnostics technique suitable for measuring intense electron beams and strong magnetic fields in hot, dense plasmas. To increase the accuracy achievable, we investigated the possibility of measuring both orthogonal components of linear polarization simultaneously using a single crystal. We demonstrated experimentally for the first time that crystals with cubic symmetry can be used as polarizing beam splitters. We also found that the polarization-splitting planes in cubic crystals can be accessed in two crystal orientations, a property that can be used to decrease the clutter in polarized spectra. In parallel, we developed a ray-tracing code and are working to add the crystal reflectivity to instrument modeling codes, to help us interpret the spectra and to design spectroscopic and spectropolarimetric instruments. This work is motivated by the need to perform measurements of important quantities in difficult experimental conditions, such as the axial magnetic field at stagnation in magnetized liner inertial fusion loads.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A theoretical investigation of laser-sustained plasma thruster

A numerical code has been successfully developed for the investigation of thruster performance using a laser-sustained hydrogen plasma as the propellant. The plasma was sustained using a 10.6-micron CO2 laser beam focused at different positions within the thruster. The physical model assumed that plasma is in thermodynamical equilibrium (LTE), and geometric ray tracing was adopted to describe the laser beam. The steady-state, axisymmetric, Navier-Stokes equations coupled with the laser power absorption process have been solved numerically. A pressure based Navier-Stokes numerical solver using body-fitted coordinates was used to calculate the laser-supported rocket flow which includes both subsonic and supersonic flow regions. From the limited parametric study, which did not try to optimize the rocket performance, it was found that better performance was obtained when the laser beam was focused closer to the rocket throat.

Jeng, San-Mou↗

Nonintrusive spectroscopic techniques for supersonic/hypersonic aerodynamics and combustion diagnostics

This paper presents an overview of the primary nonintrusive diagnostic techniques being developed by the NASA Langley Research Center to address the validation needs of Computational Fluid Dynamic (CFD) codes. The techniques include absorption in the UV and IR, Laser Induced Fluorescence, electron beam fluorescence, and a number of scattering techniques including Rayleigh, spontaneous Raman, and several coherent Raman spectroscopies. Most of the techniques are highly specialized, require complex data interpretation, and can satisfy only a few of the CFD needs. For these reasons, the evolving trend in flowfield diagnostics appears to favor a mode in which the diagnostic researcher, the experimental aerodynamicist, and the CFD community jointly define experiments based on the aeronautical requirements and on available diagnostic techniques.

Exton, R. J.↗

Development of a Benchmark Example for Delamination Fatigue Growth Prediction

The development of a benchmark example for cyclic delamination growth prediction is presented and demonstrated for a commercial code. The example is based on a finite element model of a Double Cantilever Beam (DCB) specimen, which is independent of the analysis software used and allows the assessment of the delamination growth prediction capabilities in commercial finite element codes. First, the benchmark result was created for the specimen. Second, starting from an initially straight front, the delamination was allowed to grow under cyclic loading in a finite element model of a commercial code. The number of cycles to delamination onset and the number of cycles during stable delamination growth for each growth increment were obtained from the analysis. In general, good agreement between the results obtained from the growth analysis and the benchmark results could be achieved by selecting the appropriate input parameters. Overall, the results are encouraging but further assessment for mixed-mode delamination is required

Krueger, Ronald↗

Development of a Benchmark Example for Delamination Fatigue Growth Prediction

The development of a benchmark example for cyclic delamination growth prediction is presented and demonstrated for a commercial code. The example is based on a finite element model of a Double Cantilever Beam (DCB) specimen, which is independent of the analysis software used and allows the assessment of the delamination growth prediction capabilities in commercial finite element codes. First, the benchmark result was created for the specimen. Second, starting from an initially straight front, the delamination was allowed to grow under cyclic loading in a finite element model of a commercial code. The number of cycles to delamination onset and the number of cycles during stable delamination growth for each growth increment were obtained from the analysis. In general, good agreement between the results obtained from the growth analysis and the benchmark results could be achieved by selecting the appropriate input parameters. Overall, the results are encouraging but further assessment for mixed-mode delamination is required.

Krueger, Ronald↗

Parametric raytracing modeling for NSTX-U scenario development with high harmonic fast waves and neutral beam injection

High harmonic fast waves (HHFW) are a versatile heating and current drive tool for scenario development. Extensive modeling scans were performed to find optimal parameters for different uses of HHFW in National Spherical Tokamak Experiment (NSTX-U). Scans of plasma density, temperature, magnetic field, and antenna phasing were performed both with and without neutral beam injection. For speed of calculation, the ray-tracing code GENRAY coupled to the quasilinear Fokker–Planck code CQL3D was used. CQL3D allows for a more accurate description of the fast ion population, as well as for quasilinear effects such as HHFW-induced modifications of the distribution function. Best current drive results are obtained at elevated electron temperatures and with the lowest k φ phasing. Adding neutral beams however typically strongly reduces the HHFW current drive efficiency at the low density cases due to HHFW absorption on beam ions. Results of this parametric study will feed into scenario development and predict-first whole-shot modeling of NSTX-U discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma Photonics for the Control of Kilojoule Lasers

The main goal of this project was to develop plasma photonic devices (plasma lenses or ellipsoidal plasma mirrors (EPMs)) for controlling the focus of high-energy (kJ-class) laser systems. In order to defocus a kJclass laser beam, a plasma optic must produce a very specific plasma density profile. For a plasma lens, this profile must be peaked on axis and fall off parabolicly with radius. We developed four plasma-lens target designs that were predicted to yield the ideal plasma profile required for defocusing kJ-class laser beams. Particle-in-cell (PIC) simulations using the code OSIRIS were used to explore the design. The resulting plasma lens must also be underdense, which sets an upper limit on the solid density of the plasma-lens targets to ~ 1 mg/cm 3 or less. This low of a density is highly challenging to fabricate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASTRAN DMAP Fuzzy Structures Analysis: Summary of Research

The main proposed tasks of Cooperative Agreement NCC1-382 were: (1) developing MSC/NASTRAN DMAP language scripts to implement the Soize fuzzy structures approach for modeling the dynamics of complex structures; (2) benchmarking the results of the new code to those for a cantilevered beam in the literature; and (3) testing and validating the new code by comparing the fuzzy structures results to NASA Langley experimental and conventional finite element results for two model test structures representative of aircraft fuselage sidewall construction: (A) a small aluminum test panel (SLP, single longeron panel) with a single longitudinal stringer attached with bolts; and (B) a 47 by 72 inch flat aluminum fuselage panel (AFP, aluminum fuselage panel) including six longitudinal stringers and four frame stiffeners attached with rivets.

Sparrow, Victor W.↗

Numerical Simulation of Two-grid Ion Optics Using a 3D Code

A three-dimensional ion optics code has been developed under NASA's Project Prometheus to model two grid ion optics systems. The code computes the flow of positive ions from the discharge chamber through the ion optics and into the beam downstream of the thruster. The rate at which beam ions interact with background neutral gas to form charge exchange ions is also computed. Charge exchange ion trajectories are computed to determine where they strike the ion optics grid surfaces and to determine the extent of sputter erosion they cause. The code has been used to compute predictions of the erosion pattern and wear rate on the NSTAR ion optics system; the code predicts the shape of the eroded pattern but overestimates the initial wear rate by about 50%. An example of use of the code to estimate the NEXIS thruster accelerator grid life is also presented.

ion optics↗

High-Efficiency Inductive Output Tubes Using a Third Harmonic Drive on the Grid

In this article, we discuss the use of a third harmonic component to the drive voltage on the grid of an inductive output tube (IOT). High-efficiency IOTs are typically characterized by efficiencies up to 70%–75%. However, the achievement of efficiencies greater than 80% would substantially reduce the operating costs of next-generation accelerators. In order to achieve this goal, we consider the addition of a third harmonic component to the drive signal on the grid of an IOT. The use of a third harmonic drive component in IOT guns has been considered in order to apply such a gun as the injector of radio frequency linear accelerators (RF linacs). Furthermore, we consider that the IOT will be used to provide the rf power to drive RF linacs and apply the third harmonic with the intention of increasing the efficiency of the RF output of the IOT. We consider a model IOT with a 700-MHz resonant cavity and using an annular beam with a voltage of 30 kV, an average current of 6.67 A yielding a perveance of about 1.3 μP. We simulate this IOT using the NEMESIS simulation code which has been successfully validated by comparison with the K5H90W-2 IOT developed by Communications and Power Industries. It is found that the effect of the third harmonic on the efficiency is greatest when the phase of the third harmonic is shifted by radians with respect to the fundamental drive signal and with third harmonic powers greater than about 50% that of the fundamental drive power. For the present example, we show that efficiencies approaching 86% are possible by this means.

43 PARTICLE ACCELERATORS↗

Electron beam irradiated foil expansion simulated by a coupled HYDRA/MCNP6 application

Simulating the hydrodynamic evolution of a temporally pulsed electron beam with multiple metal foils requires coupling the nearly monoenergetic, relativistic electron beam to the foil target. This deposition may be effectively calculated by Monte Carlo transport. Once this energy is deposited on a computational mesh representing the foil target, highly developed radiation-hydrodynamic codes, such as HYDRA, can describe the subsequent evolution. Currently, HYDRA does not have an internal Monte Carlo electron transport capability. To partially rectify this deficiency, a PYTHON-based scheme to link MCNP6 to HYDRA has been developed for quadrilateral meshes. Results from a seven-foil tantalum target configuration are presented to demonstrate this capability.

42 ENGINEERING↗

End-to-end differentiable digital twin for the IOTA/FAST facility

As the design complexity of modern accelerators grows, there is more interest in using controllable-fidelity simulations that have fast execution time or can yield additional insights about accelerator state. One notable example of additional information are gradients of physical observables with respect to design parameters produced by differentiable simulations. The IOTA/FAST facility has recently begun a program to implement and experimentally validate an end-to-end digital twin to serve as a virtual accelerator test stand, allowing for rapid prototyping of new software and experiments with minimal beam time costs. In this contribution we will discuss our plans and progress. Specifically, we will cover the selection and benchmarking of both physics and ML codes for linac and ring simulation, the development of generic interfaces between surrogate and physics-based sections, and presenting the control interface as either a deterministic event loop or a fully asynchronous EPICS soft input/output controller. We will also discuss challenges in model calibration and uncertainty quantification, as well as future plans to implement larger proton accelerators like PIPII and Booster.

Kuklev, N. [Fermilab]↗

Modeling betatron radiation using particle-in-cell codes for plasma wakefield accelerator diagnostics

The analysis of plasma wakefield acceleration experimental measurements, particularly in the characterization of photons emitted through the betatron radiation mechanism, requires the development of accurate numerical models. These computational models are crucial for supporting modern instrumentation designed to measure the single-shot, double-differential angular-energy radiation spectra resulting from interactions between beams and plasmas. Motivated by the needs of such applications, this paper presents detailed numerical models of betatron radiation generated in beam-plasma acceleration experiments. These models are based on the integration of the Liénard-Wiechert (LW) potentials, applied to computed particle trajectories. The particle trajectories are generated using three distinct methods: first, by tracking particles through idealized fields in the blowout regime of PWFA; second, by obtaining trajectories using the fast quasistatic particle-in-cell (PIC) code quickpic; and third, obtaining trajectories from the fully self-consistent PIC code osiris. To ensure the accuracy and reliability of these models, the paper includes various benchmark tests using analytical expressions, as well as employing the PIC code epoch, which takes an alternative approach by using a Monte Carlo quantum electrodynamics (QED)-based radiation model. Additionally, the paper presents simulations of the expected experimental betatron radiation spectra, taking into account parameters relevant to PWFA and plasma photocathode experiments at the SLAC FACET-II facility.

Yadav, M. [University of California, Los Angeles, ↗

3D simulations of inertial confinement fusion implosions part 1: inline modeling of polarized cross beam energy transfer and subsequent drive anomalies on OMEGA and NIF

Inertial confinement fusion experiments are sensitive to cross-beam energy transfer (CBET), a laser-plasma instability that redistributes laser energy in the coronal plasma through self-generated ion acoustic wave (IAW) gratings. The detailed CBET coupling depends on the polarization state of the crossing wavefields. CBET itself can also scramble the beam polarizations by inducing ellipticity through the IAW grating, and rotating the seed polarization toward that of the pump. We develop a ray-based model that describes the polarized CBET coupling and that is compatible with the framework of 3D inline radiative hydrodynamics simulations. The model is implemented in the ASTER/IFRIIT code and verified against an academic test case and an offline polarized CBET post-processor. It is then applied to the detailed configuration of the distributed polarization rotator system on OMEGA, where results highlight how polarized CBET induces significant low modes in the collisional absorption source term. Finally, the modeling is applied to a simple indirect-drive configuration, comparing CBET calculations with 96 unpolarized or polarized beams with 24 unpolarized quads. It is shown that these cases produce similar power amplification per cone of beams grouped with similar polar angles. However, the 96 beam geometry itself is found to reduce azimuthal variations in quad power after the interaction and favors beams with larger polar angles within the cones, an effect that is amplified by the polarized CBET. Application of the model to inline calculations of OMEGA implosions are presented in a companion paper.

3D inline modeling↗

Tutorial: The equations of electron emission and their evaluation

Electron emission and transport through and over potential barriers is an essential process requiring modeling and simulation to meet the design needs and characterization of an exceedingly broad range of technologically important devices and processes. The simulation and description of thermal, field, and photoemission, and the related concerns of space–charge affected electron flow, often make use of specialized formulations developed in the early days of quantum mechanics. Advancements in the utilization of electron sources and particularly the simulation of devices and applications using advanced particle-in-cell and trajectory methods for beam optics codes create a strong need for a pedagogical account of the emission models to ensure correct numerical evaluation of their equations. This Tutorial starts from simple phenomenological accounts and progressively builds to comprehensive models emphasizing straightforward and often rapid calculation. It recommends formulations to supplant the canonical Richardson–Laue–Dushman (thermal), Fowler–Nordheim (field), Fowler–DuBridge (photo), and Baroody (secondary) equations and provides a useful formulation of space–charge affected flow commonly described by the Child–Langmuir relation that takes into account cathode dependence on surface field.

Physics↗

Doppler lidar signal and turbulence study

Comparison of the second moments of the Doppler lidar signal with aircraft and tower measured parameters is being carried out. Lidar binary data tapes were successfully converted to ASCII Code on the VAX 11/780. These data were used to develop the computer programs for analyzing data from the Marshall Space Flight Center field test. Raw lidar amplitude along the first 50 forward and backward beams of Run No. 2, respectively was plotted. Plotting techniques for the same beams except with the amplitude thresholded and range corrected were developed. Plotting routines for the corresponding lidar width of the first 50 forward and backward beams were also established. The relationship between raw lidar amplitude and lidar width was examined. The lidar width is roughly constant for lidar amplitudes less than 120 dB. A field test with the NASA/MSFC ground based Doppler lidar, the instrumented NASA B-57B gust gradient aircraft, and the NASA/MSFC eight tower array was carried out. The data tape for the lidar was received and read. The aircraft data and tower data are being digitized and converted to engineering units. Velocities computed sequentially along each of the lidar beams beginning at 16:40:00, May 12, 1983 were plotted for Run No. 1.

Frost, W.↗