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 397 records · Page 22

Developments on two lithium vapor-box linear test-stand experiments

The lithium vapor-box divertor is a possible fusion power exhaust solution.It uses condensation pumping to create a gradient of vapor density in a divertor slot; this should allow a stable detachment front without active feedback.As initial explorations of the concept, two test stands which take the form of three connected cylindrical stainless steel boxes are being developed: one without plasma at PPPL, to test models of lithium evaporation and flow; and one for the linear plasma device Magnum-PSI (at DIFFER in Eindhoven, The Netherlands) to test the ability of a lithium vapor cloud to induce volumetric detachment and redistribute the plasma power.The first experiment uses boxes with diameters of 6 cm, joined by apertures with diameters of 2.2 cm. Up to 1 g of Li is placed in one box, which is heated to up to 600 degrees C. The Li evaporates, then flows to and condenses in the two other, cooler boxes over several minutes. The quantity of Li transported is assessed by weighing the boxes before and after the heating cycle, and is compared to the quantity predicted to flow for the box at its measured temperature using a Direct Simulation Monte Carlo code, SPARTA. With good experimental conditions, the two values agree to within 15%.The experiment on Magnum-PSI is in the conceptual design stage.The design is assessed by simulations using the code B2.5-Eunomia.They show that when the hydrogen-ion plasma beam, with n_e = 4e20 per cubic meter, T_e = 1.5 eV, and r = 1 cm, is passed through a 16 cm long, 12 Pa, 625 degree C Li vapor cloud, the plasma heat flux and pressure on the target are significantly reduced compared to the case without Li.With the Li present, the plasma is cooled by excitation of Li neutrals followed by radiation until it volumetrically recombines, lowering the heat flux from 3.7 MW/m^2 to 0.13 MW/m^2, and the pressure is reduced by 93%, largely by collisions of hydrogen ions with neutral Li.

divertor↗

Development of deterministic transport methods for low energy neutrons for shielding in space

Transport of low energy neutrons associated with the galactic cosmic ray cascade is analyzed in this dissertation. A benchmark quality analytical algorithm is demonstrated for use with BRYNTRN, a computer program written by the High Energy Physics Division of NASA Langley Research Center, which is used to design and analyze shielding against the radiation created by the cascade. BRYNTRN uses numerical methods to solve the integral transport equations for baryons with the straight-ahead approximation, and numerical and empirical methods to generate the interaction probabilities. The straight-ahead approximation is adequate for charged particles, but not for neutrons. As NASA Langley improves BRYNTRN to include low energy neutrons, a benchmark quality solution is needed for comparison. The neutron transport algorithm demonstrated in this dissertation uses the closed-form Green's function solution to the galactic cosmic ray cascade transport equations to generate a source of neutrons. A basis function expansion for finite heterogeneous and semi-infinite homogeneous slabs with multiple energy groups and isotropic scattering is used to generate neutron fluxes resulting from the cascade. This method, called the FN method, is used to solve the neutral particle linear Boltzmann transport equation. As a demonstration of the algorithm coded in the programs MGSLAB and MGSEMI, neutron and ion fluxes are shown for a beam of fluorine ions at 1000 MeV per nucleon incident on semi-infinite and finite aluminum slabs. Also, to demonstrate that the shielding effectiveness against the radiation from the galactic cosmic ray cascade is not directly proportional to shield thickness, a graph of transmitted total neutron scalar flux versus slab thickness is shown. A simple model based on the nuclear liquid drop assumption is used to generate cross sections for the galactic cosmic ray cascade. The ENDF/B V database is used to generate the total and scattering cross sections for neutrons in aluminum. As an external verification, the results from MGSLAB and MGSEMI were compared to ANISN/PC, a routinely used neutron transport code, showing excellent agreement. In an application to an aluminum shield, the FN method seems to generate reasonable results.

Ganapol, Barry↗

The lightcraft project

Rensselaer Polytechnic Institute has been developing a transatmospheric 'Lightcraft' technology which uses beamed laser energy to propel advanced shuttle craft to orbit. In the past several years, Rensselaer students have analyzed the unique combined-cycle Lightcraft engine, designed a small unmanned Lightcraft Technology Demonstrator, and conceptualized larger manned Lightcraft - to name just a few of the interrelated design projects. The 1990-91 class carried out preliminary and detailed design efforts for a one-person 'Mercury' Lightcraft, using computer-aided design and finite-element structural modeling techniques. In addition, they began construction of a 2.6 m-diameter, full-scale engineering prototype mockup. The mockup will be equipped with three robotic legs that 'kneel' for passenger entry and exit. More importantly, the articulated tripod gear is crucial for accurately pointing at, and tracking the laser relay mirrors, a maneuver that must be performed just prior to liftoff. Also accomplished were further design improvements on a 6-inch-diameter Lightcraft model (for testing in RPI's hypersonic tunnel), and new laser propulsion experiments. The resultant experimental data will be used to calibrate Computational Fluid Dynamic (CFD) codes and analytical laser propulsion models that can simulate vehicle/engine flight conditions along a transatmospheric boost trajectory. These efforts will enable the prediction of distributed aerodynamic and thruster loads over the entire full-scale spacecraft.

Messitt, Don G.↗

A Mixed-Mode I/II Fracture Criterion and Its Application in Crack Growth Predictions

A crack tip opening displacement (CTOD)-based, mixed mode fracture criterion is developed for predicting the onset and direction of crack growth. The criterion postulates that crack growth occurs in either the Mode I or Mode II direction, depending on whether the maximum in either the opening or the shear component of CTOD, measured at a specified distance behind the crack tip, attains a critical value. For crack growth direction prediction, the proposed CTOD criterion is shown to be equivalent to seven commonly used crack growth criteria under linearly elastic and asymptotic conditions. Under elastic-plastic conditions the CTOD criterion's prediction of the dependence of the crack growth direction on the crack-up mode mixity is in excellent agreement with the Arcan test results. Furthermore, the CTOD criterion correctly predicts the existence of a crack growth transition from mode I to mode II as the mode mixity approaches the mode II loading condition. The proposed CTOD criterion has been implemented in finite element crack growth simulation codes Z1P2DL and FRANC2DL to predict the crack growth paths in (a) a modified Arcan test specimen and fixture made of AL 2024-T34 and (b) a double cantilever beam (DCB) specimen made of AL 7050. A series of crack growth simulations have been carried out for the crack growth tests in the Arcan and DCB specimens and the results further demonstrate the applicability of the mixed mode CTOD fracture criterion crack growth predictions and residual strength analyses for airframe materials.

Sutton, Michael A.↗

Time-Averaged Velocity, Temperature and Density Surveys of Supersonic Free Jets

A spectrally resolved molecular Rayleigh scattering technique was used to simultaneously measure axial component of velocity U, static temperature T, and density p in unheated free jets at Mach numbers M = 0.6,0.95, 1.4 and 1.8. The latter two conditions were achieved using contoured convergent-divergent nozzles. A narrow line-width continuous wave laser was passed through the jet plumes and molecular scattered light from a small region on the beam was collected and analyzed using a Fabry-Perot interferometer. The optical spectrum analysis air density at the probe volume was determined by monitoring the intensity variation of the scattered light using photo-multiplier tubes. The Fabry-Perot interferometer was operated in the imaging mode, whereby the fringe formed at the image plane was captured by a cooled CCD camera. Special attention was given to remove dust particles from the plume and to provide adequate vibration isolation to the optical components. The velocity profiles from various operating conditions were compared with that measured by a Pitot tube. An excellent comparison within 5m's demonstrated the maturity of the technique. Temperature was measured least accurately, within 10K, while density was measured within 1% uncertainty. The survey data consisted of centerline variations and radial profiles of time-averaged U, T and p. The static temperature and density values were used to determine static pressure variations inside the jet. The data provided a comparative study of jet growth rates with increasing Mach number. The current work is part of a data-base development project for Computational Fluid Dynamics and Aeroacoustics codes that endeavor to predict noise characteristics of high speed jets. A limited amount of far field noise spectra from the same jets are also presented. Finally, a direct experimental validation was obtained for the Crocco-Busemann equation which is commonly used to predict temperature and density profiles from known velocity profiles. Data presented in this paper are available in ASCII format upon request.

Panda, Jayanta↗

General Tool for Evaluating High-Contrast Coronagraphic Telescope Performance Error Budgets

The Coronagraph Performance Error Budget (CPEB) tool automates many of the key steps required to evaluate the scattered starlight contrast in the dark hole of a space-based coronagraph. The tool uses a Code V prescription of the optical train, and uses MATLAB programs to call ray-trace code that generates linear beam-walk and aberration sensitivity matrices for motions of the optical elements and line-of-sight pointing, with and without controlled fine-steering mirrors (FSMs). The sensitivity matrices are imported by macros into Excel 2007, where the error budget is evaluated. The user specifies the particular optics of interest, and chooses the quality of each optic from a predefined set of PSDs. The spreadsheet creates a nominal set of thermal and jitter motions, and combines that with the sensitivity matrices to generate an error budget for the system. CPEB also contains a combination of form and ActiveX controls with Visual Basic for Applications code to allow for user interaction in which the user can perform trade studies such as changing engineering requirements, and identifying and isolating stringent requirements. It contains summary tables and graphics that can be instantly used for reporting results in view graphs. The entire process to obtain a coronagraphic telescope performance error budget has been automated into three stages: conversion of optical prescription from Zemax or Code V to MACOS (in-house optical modeling and analysis tool), a linear models process, and an error budget tool process. The first process was improved by developing a MATLAB package based on the Class Constructor Method with a number of user-defined functions that allow the user to modify the MACOS optical prescription. The second process was modified by creating a MATLAB package that contains user-defined functions that automate the process. The user interfaces with the process by utilizing an initialization file where the user defines the parameters of the linear model computations. Other than this, the process is fully automated. The third process was developed based on the Terrestrial Planet Finder coronagraph Error Budget Tool, but was fully automated by using VBA code, form, and ActiveX controls.

Marchen, Luis F.↗

SpaRibs Geometry Parameterization for Wings with Multiple Sections using Single Design

The SpaRibs topology of an aircraft wing has a significant effect on its structural behavior and stability as well as the flutter performance. The development of additive manufacturing techniques like Electron Beam Free Form Fabrication (EBF3) has made it feasible to manufacture aircraft wings with curvilinear spars, ribs (SpaRibs) and stiffeners. In this article a new global-local optimization framework for wing with multiple sections using curvilinear SpaRibs is described. A single design space is used to parameterize the SpaRibs geometry. This method has been implemented using MSC-PATRAN to create a broad range of SpaRibs topologies using limited number of parameters. It ensures C0 and C1 continuities in SpaRibs geometry at the junction of two wing sections with airfoil thickness gradient discontinuity as well as mesh continuity between all structural components. This method is advantageous in complex multi-disciplinary optimization due to its potential to reduce the number of design variables. For the global-local optimization the local panels are generated by an algorithm which is totally based on a set algebra on the connectivity matrix data. The great advantage of this method is that it is completely independent of the coordinates of the nodes of the finite element model. It is also independent of the order in which the elements are distributed in the FEM. The code is verified by optimizing of the CRM Baseline model at trim condition at Mach number equal to 0.85 for five different angle of attack (-2deg, 0deg,2deg,4deg and 6deg). The final weight of the wing is 19,090.61 lb. This value is comparable to that obtained by Qiang et al. 6 (19,269 lb).

wing sections↗

NASA Tech Briefs, July 2013

Dielectrophoresis-Based Particle Sensor Using Nanoelectrode Arrays; Multi-Dimensional Damage Detection for Surfaces and Structures; ULTRA: Underwater Localization for Transit and Reconnaissance Autonomy; Autonomous Cryogenic Leak Detector for Improving Launch Site Operations; Submillimeter Planetary Atmospheric Chemistry Exploration Sounder; Method for Reduction of Silver Biocide Plating on Metal Surfaces; Silicon Micromachined Microlens Array for THz Antennas; Forward-Looking IED Detector Ground Penetrating Radar; Fully Printed, Flexible, Phased Array Antenna for Lunar Surface Communication, Battery Charge Equalizer with Transformer Array; An Efficient, Highly Flexible Multi-Channel Digital Downconverter Architecture; Dimmable Electronic Ballast for a Gas Discharge Lamp; Conductive Carbon Nanotube Inks for Use with Desktop Inkjet Printing Technology; Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates; High-Performance, Low-Temperature-Operating, Long-Lifetime Aerospace Lubricants; Carbon Nanotube Microarrays Grown on Nanoflake Substrates; Differential Muon Tomography to Continuously Monitor Changes in the Composition of Subsurface Fluids; Microgravity Drill and Anchor System; 20 Granular Media-Based Tunable Passive Vibration Suppressor; 21 Miga Aero Actuator and 2D Machined Mechanical Binary Latch; Micro-XRF for In Situ Geological Exploration of Other Planets; Hydrogen-Enhanced Lunar Oxygen Extraction and Storage Using Only Solar Power; Uplift of Ionospheric Oxygen Ions During Extreme Magnetic Storms; Miniaturized, High-Speed, Modulated X-Ray Source; Hollow-Fiber Spacesuit Water Membrane Evaporator 25 High-Power Single-Mode 2.65-micrometers InGaAsSb/AlInGaAsSb Diode Lasers; Optical Device for Converting a Laser Beam Into Two Co-aligned but Oppositely Directed Beams; A Hybrid Fiber/Solid-State Regenerative Amplifier with Tunable Pulse Widths for Satellite Laser Ranging; X-Ray Diffractive Optics; SynGenics Optimization System (SynOptSys); 29 CFD Script for Rapid TPS Damage Assessment; radEq Add-On Module for CFD Solver Loci-CHEM; Science Opportunity Analyzer (SOA) Version 8; 30 Autonomous Byte Stream Randomizer; Distributed Engine Control Empirical/Analytical Verification Tools; Dynamic Server-Based KML Code Generator Method for Level-of-Detail Traversal of Geospatial Data; Automated Planning of Science Products Based on Nadir Overflights and Alerts for Onboard and Ground Processing; Linked Autonomous Interplanetary Satellite Orbit Navigation; Risk-Constrained Dynamic Programming for Optimal Mars Entry, Descent, and Landing; Scheduling Operations for Massive Heterogeneous Clusters; Deepak Condenser Model (DeCoM); Flight Software Math Library; Recirculating 1-K-Pot for Pulse-Tube Cryostats; 35 Method for Processing Lunar Regolith Using Microwaves; Wells for In Situ Extraction of Volatiles from Regolith (WIEVR); and Estimating the Backup Reaction Wheel Orientation Using Reaction Wheel Spin Rates Flight Telemetry from a Spacecraft.

Source record↗

Free-Space Optical Communication for CubeSats in Low Lunar Orbit (LLO)

The NASA ARTEMIS Program will include LunaNet, a highly extensible, open architecture, lunar communications and navigation network. A constellation of CubeSats in Low Lunar Orbit (LLO), 100 km, could form an optical communications and navigation network as part of LunaNet, with terminals on the lunar surface, including mobile ones such as with astronauts and rovers. The proposed CubeSat nodes should provide data relay and navigational aid services. The proposed effort herein is to develop a fine pointing capability for laser beam pointing to augment body pointing by CubeSats. Body pointing was used by Aerospace Corporation for the CubeSats in LEO in NASA’s Optical Communications and Sensors Demonstration (OCSD) program [1]. Previously, this fine pointing capability was computer simulated for the OCSD program [2,3]. With fine pointing, the spot size on the Earth was reduced by a factor of eight with a reduction in laser output power by a factor of sixty-four, thereby mitigating the thermal load challenge on the CubeSats. The same reductions in spot size and laser output power can be achieved for CubeSats in LLO. A new method is described for optical data transmissions from satellites, which uses laser arrays for laser beam pointing. It combines a lens system and an array of vertical-cavity, surface-emitting lasers and photodetectors, an VCSEL/Photodetector Array, (both mature technologies), in a novel way. This system is applied to CubeSats in low lunar orbit, (LLO), which use body pointing. Also, It may be able to replace current architectures which use dynamical systems, (i.e., moving parts) to point the laser, and which may also use vibration isolation platforms. The computer simulations used the optics code, OpticStudio, from Zemax, LLC, which has the capabilities to model the laser source and diffraction effects from wave optics. These capabilities make it possible to model laser beam propagation over long space communication distances.

Goorjian, P. M.↗

Hall2De Simulations of a Magnetic Nozzle

The two-dimensional, fluid-based Hall thruster code Hall2De is adapted to simulate a low temperature propulsive magnetic nozzle. Electron cyclotron resonance (ECR) thrusters developed at the Office National d’Etudes et de Recherches Aérospatiales (ONERA) in France and the University of Michigan are investigated. For each device, the simulated domain is defined from the nozzle throat to a location 10 thruster radii downstream. A numerical mesh is defined based on the measured magnetic field, and the governing fluid equations for singlycharge ions, electrons, and neutrals are solved. The boundary conditions for the throat of the magnetic nozzle are based on measured values of electron temperature and ion velocity. Multiple boundary conditions for the electron temperature in the downstream regions of the domain are considered. The simulated ion velocity and potential agrees well with experiment, while the simulated divergence efficiencies and thrusts are 25% lower and 45-55% higher than experiment. When a Neumann boundary condition for electron temperature is applied to the far-plume boundaries, the nozzle was seen to be purely isothermal. This suggests that the thermal conductivity along field lines in real devices is much lower than classical simulations predict, and anomalous resistivity may be needed to reproduce the temperature distribution observed in experiment. The larger beam divergence observed in our simulation may also be related to the lack of anomalous resistivity.

Jorns, Benjamin A.↗

Fundamental Mixing and Combustion Experiments for Propelled Hypersonic Flight

The first experiment is a study of a coaxial jet discharging into stagnant laboratory air, with center jet of a mixture of 5% oxygen and 95% helium by volume and coflow jet of air. The exit flow pressure of both center-jet and coflow nozzles is 1 atmosphere. The presence of oxygen in the center jet is to allow the use of an oxygen flow-tagging technique (RELIEF4) to obtain non-intrusive velocity measurements. Both jets are nominally Mach 1.8, but, because of the greater speed of sound, the center jet velocity is more than twice that of the coflow. The mixing layer which forms between the center jet and the coflow near the nozzle exit is compressible, with a calculated convective Mach number of approximately 0.7. This geometry has several advantages: The streamwise development of the flow is generally dominated by turbulent stresses (rather than pressure forces), and thus calculations are sensitive to turbulence modeling. It includes features present in supersonic combustors, including a compressible mixing layer near the nozzle exit and a light-gas/air plume downstream. Since it is a free jet, it provides easy access for both optical instrumentation and probes. Since it is axisymmetric, it requires fewer experimental measurements to fully characterize, and calculations can be performed with more modest computer resources. However, weak shock waves formed at the nozzle exit strengthen and turn normal as they approach the axis, complicating the flow. Care is thus taken in the design of the facility to provide as near as possible to 1-D flow at the exit of both center and coflow nozzles, and to minimize the strength of waves generated at the nozzle exit. Results from this experiment are compared to CFD solutions obtained by VULCAN, a previously developed code used in engine analysis. The second experiment is a study of a supersonic combustor consisting of a diverging duct with single downstream-angled wall injector. Thus, the geometry is relatively simple and large regions of subsonic recirculating flow are avoided. The nominal entrance Mach number is 2 and the enthalpy of the test gas (hot air "simulant") is nominally that of Mach 7 flight. It was believed, on the basis of calculations performed that this would produce mixing-limited flow, that is to say, one for which chemical reaction to equilibrium proceeds at a much greater rate than mixing. It later proved that this was not the case. The primary experimental technique employed is coherent anti-Stokes Raman spectroscopy, known by its acronym CARS. The species probed is molecular nitrogen and the quantity measured is temperature. Intrusive probes, such as Pitot, total temperature, hot-wire, etc., are not used due to access difficulty and high heat flux in the combustor, and because they may alter the flow. CARS has several advantages over other optical methods. It is a relatively mature and well-understood technique. Signal levels are relatively high and the signal is in the form of a coherent (laser) beam that can be collected through small windows. Incoherent (non-CARS) interferences are rejected by spatial filtering.

Diskin, G. S.↗

Analytic model of electron transport through and over non-linear barriers

Tunneling barriers are an essential component of electron sources, sensors, detectors, and vacuum nanoelectronics, and a pivotal factor in their performance, but the barriers themselves routinely depart from the analytic models used to model their behavior. Here, a new formalism is developed to analytically and accurately model emission through and over barriers associated with depletion layers, nanotip barriers, and MIM structures. The transmission probability for depletion layer and metal-insulator-metal (MIM) and metaloxide-semiconductor (MOS) barriers is accurately modeled as the electron energy exceeds the barrier height using approaches designed for rapid implementation demanded by simulation codes, and extensible to general barriers. The models supersede conventional thermal and field models in depletion and MIM/MOS barrier studies. Thermal-field methods are used to treat the transmission probability, and shape factor methods to treat the tunneling factor. Analytic formulae for current density are obtained. The methods ease device simulation and characterization of current-voltage relations for emerging technologically interesting barriers with better accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Build-Up Approach to Updating the Mock Quiet Spike Beam Model

When a new aircraft is designed or a modification is done to an existing aircraft, the aeroelastic properties of the aircraft should be examined to ensure the aircraft is flight worthy. Evaluating the aeroelastic properties of a new or modified aircraft can include performing a variety of analyses, such as modal and flutter analyses. In order to produce accurate results from these analyses, it is imperative to work with finite element models (FEM) that have been validated by or correlated to ground vibration test (GVT) data, Updating an analytical model using measured data is a challenge in the area of structural dynamics. The analytical model update process encompasses a series of optimizations that match analytical frequencies and mode shapes to the measured modal characteristics of structure. In the past, the method used to update a model to test data was "trial and error." This is an inefficient method - running a modal analysis, comparing the analytical results to the GVT data, manually modifying one or more structural parameters (mass, CG, inertia, area, etc.), rerunning the analysis, and comparing the new analytical modal characteristics to the GVT modal data. If the match is close enough (close enough defined by analyst's updating requirements), then the updating process is completed. If the match does not meet updating-requirements, then the parameters are changed again and the process is repeated. Clearly, this manual optimization process is highly inefficient for large FEM's and/or a large number of structural parameters. NASA Dryden Flight Research Center (DFRC) has developed, in-house, a Mode Matching Code that automates the above-mentioned optimization process, DFRC's in-house Mode Matching Code reads mode shapes and frequencies acquired from GVT to create the target model. It also reads the current analytical model, as we11 as the design variables and their upper and lower limits. It performs a modal analysis on this model and modifies it to create an updated model that has similar mode shapes and frequencies as those of the target model. The Mode Matching Code output frequencies and modal assurance criteria (MAC) values that allow for the quantified comparison of the updated model versus the target model. A recent application of this code is the F453 supersonic flight testing platform, NASA DFRC possesses a modified F-15B that is used as a test bed aircraft for supersonic flight experiments. Traditionally, the finite element model of the test article is generated. A GVT is done on the test article ta validate and update its FEM. This FEM is then mated to the F-15B model, which was correlated to GVT data in fall of 2004, A GVT is conducted with the test article mated to the aircraft, and this mated F-15B/ test article FEM is correlated to this final GVT.

Herrera, Claudia Y.↗

Understanding extraction limits of plasma cathodes with experiment and simulation

The project is focused on computations for enhanced ionization near the exit orifice of a plasma cathode and how the cathode extraction electrode geometry and stray magnetic field structure influence available extractable current. The computational effort will take place at the Princeton Collaborative Research Laboratory (PCRF). The specific modeling platform to be utilized from PCRF is the electrostatic direct implicit particle in cell (EDIPIC) code. EDIPIC will be used to study ionization processes in the electron extraction region of the cathode. The control of electron energy is a key aspect to minimizing power losses in the plasma. Exiting electrons ionize gas leaving the orifice and in the process produce the dense plasma from which large amounts of current can be extracted. Project activities include studying processes that determine the energy distribution of the exiting electrons thus providing insight into how to optimize the source. The model will be validated and developed using experimental measurements as warranted. The familiarity and expertise of the team at PPPL with this type of plasma and the challenges posed greatly enhances the likelihood of gaining insight into the extraction sheath under the duration of the proposed activity. The effort is expected to shed a great deal of insight into electron extraction processes and pathways for optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗

Shielding Analysis of a Preclinical Bremsstrahlung X-ray FLASH Radiotherapy System within a Clinical Radiation Therapy Vault

A preclinical radiotherapy system producing FLASH dose rates with 12 MV bremsstrahlung x rays is being developed at Stanford University and SLAC National Accelerator Laboratory. Because of the high expected workload of 6,800 Gy w –1 at the isocenter, an efficient shielding methodology is needed to protect operators and the public while the preclinical system is operated in a radiation therapy vault designed for 6 MV x rays. Here, in this study, an analysis is performed to assess the shielding of the local treatment head and radiation vault using the Monte Carlo code FLUKA and the empirical methodology given in the National Council on Radiation Protection and Measurements Report 151. Two different treatment head shielding designs were created to compare single-layer and multilayer shielding methodologies using high-Z and low-Z materials. The multilayered shielding methodology produced designs with a 17% reduction in neutron fluence leaking from the treatment head compared to the single layered design of the same size, resulting in a decreased effective dose to operators and the public. The conservative assumptions used in the empirical methods can lead to over-shielding when treatment heads use polyethylene or multilayered shielding. High-Z/Low-Z multilayered shielding optimized via Monte Carlo is shown to be effective in the case of treatment head shielding and provide more effective shielding design for external beam radiotherapy systems that use 12 MV bremsstrahlung photons. Modifications to empirical methods used in the assessment of MV radiotherapy systems may be warranted to capture the effects of polyethylene in treatment head shielding.

61 RADIATION PROTECTION AND DOSIMETRY↗

Development of a reduced model for energetic particle transport by sawteeth in tokamaks

We report that the sawtooth instability is known for inducing transport and loss of energetic particles (EPs), and for generating seed magnetic islands that can trigger tearing modes. Both effects degrade the overall plasma performance. Several theories and numerical models have been previously developed to quantify the expected EP transport caused by sawteeth, with various degrees of sophistication to differentiate the response of EPs at different energies and on different orbits (e.g. passing vs. trapped), although the analysis is frequently limited to a single time slice during a tokamak discharge. This work describes the development and initial benchmark of a framework that enables a reduced model for EP transport by sawteeth retaining the full EP phase-space information. The model, implemented in the ORBIT hamiltonian particle-following code, can be used either as a standalone post-processor taking input data from codes such as TRANSP, or as a pre-processor to compute transport coefficients that can be fed back to TRANSP for time-dependent simulations including the effects of sawteeth on EPs. The advantage of the latter approach is that the evolution of the EP distribution can be simulated quantitatively for sawtoothing discharges, thus enabling a more accurate modeling of sources, sinks and overall transport properties of EP and thermal plasma species for comprehensive physics studies that require detailed information of the fast-ion distribution function and its evolution over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗