Search NASA⌕ Search

SEARCH · Search NASA

Results for “Injection testing”

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 1,117 records · Page 62

A learning-augmented approach for AC optimal power flow

Because of the high nonlinearity of AC optimal power flow (OPF), numerous efforts have been made in recent decades to find efficient methods. Machine learning (ML) has proven to significantly reduce the computational costs in many real-world problems. Thus, this paper develops a learning-augmented method for solving AC OPF, which integrates both power network equations and ML to yield near-optimal solutions. More specifically, ML models are developed to first predict bus voltage magnitudes and angles. Then, physics-based network equations are employed to calculate the power injection at different buses. Three ML algorithms, i.e., random forest, multi-target decision tree, and extreme learning machine, are explored and compared. To evaluate the efficiency of the proposed learning-augmented AC OPF solver, the MATPOWER Interior Point Solver is adopted as a baseline. Case studies on both 500-bus and 4918-bus test networks show that the proposed learning-augmented method has reduced the computational time by 15–100 times depending on the network size with a minimal loss in optimality.

42 ENGINEERING↗

The Dynamics of Megafire Smoke Plumes in Climate Models: Why a Converged Solution Matters for Physical Interpretations

As the climate system warms, megafires have become more frequent with devastating effects. A byproduct of these events is the creation of smoke plumes that can rise into the stratosphere and spread across the globe where they reside for many months. To gain a deeper understanding of the plume dynamics, global climate simulations of a megafire were performed at a wide range of grid spacings from 2.0° down to 7 km, including a 7 km nonhydrostatic experiment. The analysis focuses on how the resolved dynamics affects the specification of the plume characteristics such as injection height and black carbon (BC) mass. Prior studies initialize the smoke plume at one or a few grid points and this is shown here to produce severely dissipative dynamics. In order to validate such simulations with observations, enhancements of the plume characteristics to offset the dissipation is necessary. Using a numerically converged simulation, sensitivity tests show that to approximate the observed stratospheric lifetime, a reduction in BC fraction by 50% is necessary for external mixtures. The vorticity dynamics of the plume is also analyzed with a Lagrangian budget to understand the mechanisms responsible for the evolution of a collocated anticyclonic vortex. The results can be distilled down into a simple conceptual model. As the plume rises, the air diverges at the top of the updraft where the largest concentrations of smoke are found. This divergence induces a dilution of the background cyclonic absolute vorticity producing an anticyclonic vortex. Vortex decay occurs from opposite arguments.

54 ENVIRONMENTAL SCIENCES↗

Integrated plasma facing component calorimetry for measurement of shot integrated deposited energy in the NSTX-U

Here, the upgrade to the National Spherical Torus eXperiment (NSTX-U) increases the injected neutral beam power up to 12 MW and the plasma current up to I p = 2 MA for plasma durations up to 5 s. The graphite plasma facing components have been re-designed to handle greater heat and energy fluxes than were seen in NSTX using a castellated design. We present the experimental testing and validation of a castellated graphite target, similar to the prototype tile design, instrumented with thermocouples at various depths in the castellation. During testing, incident heat flux is provided by a programmed electron beam system and surface temperatures are measured via infrared thermography directly viewing the target surface. It was found that the thermocouple response scaled linearly with the measured surface temperature rise regardless of thermocouple depth in the castellation. A sensitivity of 14.3 °C/kJ of deposited energy was found when treating individual castellations as a semi-infinite solid.

47 OTHER INSTRUMENTATION↗

Southeast Regional Carbon Sequestration Partnership (SECARB): Analyze Data Collected, Revise Surface Monitoring Plan to Maximize Benefit and Continue Data Collection

Objectives of the shallow groundwater study at the Cranfield site since 2008, as a part of the SECARB Phase III early test, were to 1) assess potential impacts of CO 2 leakage on groundwater quality; and 2) test and develop monitoring technologies for CO 2 leakage detection. The study did not intend to seek CO2 leakage signals via monitoring groundwater chemistry in the Cranfield shallow aquifer. Instead, quarterly or semi-annually groundwater sampling from the shallow aquifer was used to characterize spatial and temporal variations in geochemical parameters because CO 2 injection into the deep reservoir started in 2008.

01 COAL, LIGNITE, AND PEAT↗

Using Controlled-Source Electromagnetic (CSEM) for CO 2 Storage Monitoring in the North Dakota CarbonSAFE Project

Conference presentation for International Meeting for Applied Geophysics & Energy (IMAGE), September 26 – October 1, 2021, Denver, CO. Geophysical methods are key for characterizing the geologic formations to store CO 2 and monitor the injected CO 2 over time to ensure containment. In the integrated multimeasurement geophysical approach considered for this project, it is expected that the controlled-source electromagnetic (CSEM) method is a strong contributor to mapping the CO 2 movement. A feasibility study of the CSEM method, including 1D and 3D modeling and a field noise test, was conducted to determine its effectiveness in monitoring CO 2 in the Broom Creek and Deadwood Formations. The study results demonstrate that the CSEM method can be used for CO 2 storage monitoring in the study area.

20 FOSSIL-FUELED POWER PLANTS↗

Method and apparatus to produce high specific impulse and moderate thrust from a fusion-powered rocket engine

A system and method for producing and controlling high thrust and desirable specific impulse from a continuous fusion reaction is disclosed. The resultant relatively small rocket engine will have lower cost to develop, test, and operate that the prior art, allowing spacecraft missions throughout the planetary system and beyond. The rocket engine method and system includes a reactor chamber and a heating system for heating a stable plasma to produce fusion reactions in the stable plasma. Magnets produce a magnetic field that confines the stable plasma. A fuel injection system and a propellant injection system are included. The propellant injection system injects cold propellant into a gas box at one end of the reactor chamber, where the propellant is ionized into a plasma. The propellant and fusion products are directed out of the reactor chamber through a magnetic nozzle and are detached from the magnetic field lines producing thrust.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The FNAL Booster 2nd Harmonic RF Cavity

A second harmonic RF cavity which uses perpendicularly biased garnet for frequency tuning is currently being constructed for use in the Fermilab Booster. The cavity will operate at twice the fundamental RF frequency, from ~76 - 106 MHz, and will be turned on only during injection, and transition or extraction. Its main purpose is to reduce beam loss as required by Fermilab’s Proton Improvement Plan (PIP). After three years of optimization and study, the cavity design has been finalized and all constituent parts have been received. We discuss the design aspects of the cavity and its associated systems, component testing, and status of the cavity construction.

43 PARTICLE ACCELERATORS↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

A novel framework for hosting capacity analysis with spatio-temporal probabilistic voltage sensitivity analysis

Smart grids are envisioned to accommodate high penetration of distributed photovoltaic (PV) generation, which may cause adverse grid impacts in terms of voltage violations. Therefore, PV Hosting capacity is being used as a planning tool to determine the maximum PV installation capacity that causes the first voltage violation and above which would require infrastructure upgrades. Additionally, traditional methods of Hosting capacity analysis are scenario based and computationally complex as they rely on iterative load flow algorithms that require investigating a large number of scenarios for accurate assessment of PV impacts. Therefore, this paper presents a computationally efficient analytical approach to compute the probability distribution of voltage change due to random behavior of randomly located multiple distributed PVs. The proposed approach is based on Spatio-temporal probabilistic voltage sensitivity analysis that exploits both spatial and temporal uncertainties associated with PV injections. Thereafter, the derived distribution is used to quantify voltage violations for various PV penetration levels and subsequently determine the hosting capacity of the system without the need to examine large number of scenarios. Results of the proposed framework are validated via conventional load flow based simulation approach on the IEEE 37 and IEEE 123 node test systems.

42 ENGINEERING↗

Characterization of the Triassic Newark Basin of New York and New Jersey for geologic storage of carbon dioxide

The Newark Basin is a Triassic-aged rift basin underlying densely populated, industrialized sections of New York, New Jersey and Pennsylvania. The Basin is an elongate half-graben encompassing an area of more than 7,510 square-kilometers (2,900 square-miles), and could represent a key storage component for commercial scale management of carbon dioxide emissions via geologic sequestration. The project team first acquired published reports, surface and subsurface maps, and seismic data, which formed the basis for a three-dimensional model framework for the northern end of the Basin incorporating stratigraphic, hydrologic, and water quality data. Field investigations included drilling, coring, and logging of two stratigraphic test borings in Clarkstown, NY (Exit 14 Tandem Lot Well No. 1), drilled to a depth of 2,099 meters (6,885 feet); and Palisades, NY (Lamont-Doherty Earth Observatory Test Well No. 4) drilled to a depth of 549 meters (1,802 feet). Two two-dimensional seismic reflection data lines arrayed perpendicularly were acquired by Schlumberger/WesternGeco to help characterize the structure and stratigraphy and as part of pre-drilling field screening activities for the deep stratigraphic borehole. A total of 47 meters (155 feet) of continuous whole core was recovered from the Tandem Lot boring from depths of 1,393 meters (4,570 feet) to 1,486 meters (4,877 feet). Twenty-five horizontal rotary cores were collected in mudstones and sandstones in the surface casing hole and fifty-two cores were taken in various lithologies in the deep borehole. Rotary core plugs were analyzed by Weatherford Laboratories for routine and advanced testing. Rotary core plug trim end thin sections were evaluated by the New York State Museum for mineralogical analysis and porosity estimation. Using core samples, Lawrence Berkley National Laboratory designed and completed laboratory experiments and numerical modeling analyses to characterize the dissolution and reaction of carbon dioxide with formation brine and minerals, and resulting effects on injection rate, pressure, effective storage volume, and carbon dioxide migration within a prospective sandstone reservoir. $$Three potential porous and permeable sandstone units were identified in the Passaic Formation at the New York State Thruway Exit 14 location. Potential Flow Unit 1, at a depth of 643 meters (2,110 feet) to 751 meters (2,465 feet); Potential Flow Unit 2 at a depth of 853 meters (2,798 feet) to 1,000 meters (3,280 feet); and Potential Flow Unit 3, at a depth of 1,114 meters (3,655 feet) to 1,294 meters (4,250 feet). Reactive transport simulations of interactions between carbon dioxide, brine and formation minerals were carried out to evaluate changes in formation water chemistry, mineral precipitation and dissolution reactions, and any potential resulting effects on formation permeability. The experimental and modeling analyses suggest that mineral precipitation and dissolution reactions (within the target formation) are not expected to lead to significant changes to the underground hydrologic system over time frames (~30 years) typically relevant for carbon dioxide injection operations. Key findings of this basin characterization study include an estimate of carbon dioxide storage capacity in the Newark Basin. Assuming an average porosity of twelve percent and an aquifer volume of 6.1E+12 meters3, calculated ranges of likely storage capacity range from 1.9 – 20.2 gigatonnes under high temperature (low carbon dioxide density) conditions; and 2.9 – 30.2 gigatonnes under low temperature (low carbon dioxide density) conditions. Intra-basin faulting, geometry of the Palisades Sill, and the presence of altered meta-sediments above and below the Sill, increase potential compartmentalization within the basin. A structural/stratigraphic trap type may occur where porous/permeable sediments are cross-cut by the Palisades Sill. Potential injection intervals are present within the Stockton Formation of the Newark Basin. Additional porous/permeable intervals may be present within sandstones of the Passaic Formation, increasing projected storage capacity. Deeper wedges of strata are likely present in the deeper portions of the basin in southern New York and into northern New Jersey. Abundant mudstones are present in the Passaic, Lockatong, and Stockton Formations. These intervals have the requisite petrophysical properties to form effective primary and secondary containment intervals to industrial-scale sequestration of carbon dioxide in the Newark Basin. Hydro-thermally altered meta-sediments in the region immediately surrounding the top and base of the Palisades Sill is devoid of porosity/permeability and forms an additional effective lateral/vertical sealing cap rock.

58 GEOSCIENCES↗

Benchmark Exercise Report for Experimental Study of Bubble Scrubbing in Water Coolant Pool

Mechanistic assessments of radionuclide release during postulated accidents are expected to be included in advanced reactor license applications. The mechanistic source term (MST) provides an opportunity for vendors to realistically evaluate the radiological consequences of an incident, and may aid in justifying reduced emergency planning zones and plant sites. However, the development of MSTs for advanced nuclear reactors is challenging because there are numerous phenomena that can affect the transport and retention of radionuclides. As part of a trial MST assessment for a metal-fueled, pool-type sodium cooled fast reactor (SFR), led by Argonne National Laboratory, a simplified radionuclide transport code (SRT code) was developed, which includes models to estimate the quantity of fission product aerosols scrubbed in the sodium pool during postulated accident scenarios. In a pool-type SFR, when fission products are released into the coolant pool due to failure of fuel pins, most of the radionuclides are scrubbed by the coolant pool, but some have the potential to migrate to the cover gas region through entrainment within gas bubbles. The SRT code contains a model that evaluates this scrubbing behavior and calculates the fraction of fission product aerosols that reach the cover gas. Due to a lack of available validation data for sodium pool scrubbing, the U.S. Department of Energy funded an experiment at the University of Wisconsin-Madison to measure aerosol scrubbing by injecting air bubbles containing aerosol into a coolant pool. Prior to performing an experiment with liquid sodium, a water loop experiment was performed. Their experiment evaluated the effect of changing the aerosol size, aerosol density, aerosol concentration, bubble size, and pool depth on the aerosol scrubbing efficiency of the pool. In this benchmark experiment, the base tests were conducted by repeated tests of isolated bubbles. Afterwards, more prototypic tests with bubble swarms were performed to evaluate the interactions between the bubbles. The bubble swarm test was able to confirm that a larger amount of aerosol scrubbing occurred than the single bubble test. It was also confirmed that as the bubble size, aerosol density, and pool height increase, the extent of pool scrubbing also increases and does not change with the aerosol concentration. In addition, since the degree of scrubbing is the lowest at aerosol sizes between 0.01 and 1 μm, that is, the largest amount of aerosol is emitted, it was confirmed that the analysis of this size in MST is the most important. This benchmark experiment informs the direction of future sodium experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Economics of Big Area Addtiive Manufacturing

Case studies on the economics of Additive Manufacturing (AM) suggest that processing time is the dominant cost in manufacturing. Most additive processes have similar performance metrics: small part sizes, low production rates and expensive feedstocks. Big Area Additive Manufacturing is based on transitioning polymer extrusion technology from a wire to a pellet feedstock. Utilizing pellets significantly increases deposition speed and lowers material cost by utilizing low cost injection molding feedstock. The use of carbon fiber reinforced polymers eliminates the need for a heated chamber, significantly reducing machine power requirements and size constraints. We hypothesize that the increase in productivity coupled with decrease in feedstock and energy costs will enable AM to become more competitive with conventional manufacturing processes for many applications. As a test case, we compare the cost of using traditional fused deposition modeling (FDM) with BAAM for additively manufacturing composite tooling.

additive manufacturing↗

Deep Controlled Source Electro-Magnetic Sensing: A Cost Effective, Long-Term Tool for Sequestration Monitoring

The proposed system was designed to operate as a permanent, autonomous monitoring and data collection system that can provide much higher temporal data density than can be achieved economically with 3-Dimensional (3D) seismic surveys. It can operate over broad areas for long periods of time providing full 3D data sets on a monthly basis at a very low cost. By borrowing techniques commonly used in marine CSEM, structural information from background seismic surveys can be incorporated into the CSEM modeling to provide high resolution models of CO 2 progression within reservoirs. The system uses borehole-based vertical-electric-dipole sources placed at reservoir depths in the formation. The electric and magnetic fields induced by this source are received on the surface using an array of stations. The project was conducted in three phases. Phase I demonstrated the feasibility of the system to collect static/reference data at the Ketzin CO 2 storage pilot site in Germany. In Phase I, numerical modeling was used to determine the optimal configurations and requirements for sensor sensitivity and data accuracy. Based on the model results, existing hardware and software were modified. The CSEM system was then field tested at the Ketzin site. The data were imaged and the results were compared with independent studies of the reservoir and overburden geo-electrical characteristics. Phase II demonstrated the ability to provide sensitive, cost-effective measurement of changes in reservoir properties and changes in the overlying formations using a second round of measurements at the Ketzin site. A prototype autonomous recording system was developed and tested as a subset of the measurement points. Phase III of the project quantified the advantages (and disadvantages) of the fully autonomous data collection subsystems by comparing them with repeated measurements made with mobile stations. The Phase III also provided an additional time point in measuring post-closure changes in and above the Ketzin CO 2 injection site. Validation of the CSEM field survey results were completed using a method such as other alternative and available field site data.

58 GEOSCIENCES↗

Alfven eigenmodes and fast ion transport in negative triangularity DIII-D plasmas

The first energetic particle experiments in negative triangularity tokamak plasmas have been carried out on DIII-D. Alfvén eigenmode (AE) activity and associated fast ion transport comparable to that in positive triangularity is observed during the current ramp portion of all plasmas with early beam heating indicating negative triangularity does not confer a special advantage with respect to AE induced transport. In these discharges, a range of mode activity driven by the sub-Alfvénic (Vbeam/VA < 0.5) 80 kV neutral beams is found including beta induced Alfvén acoustic eigenmodes, beta induced Alfvén eigenmodes, reversed shear Alfvén eigenmodes (RSAEs) and toroidicity induced Alfvén eigenmodes (TAEs). Mode intermittency and possibly chirping appears to be more common than comparable positive triangularity and/or oval discharges but overall, the unstable spectra and mode amplitudes observed on magnetics, CO2 interferometry, and electron cyclotron emission in a set of matched positive and negative triangularity cases at moderate beam power is similar. Large levels of Alfvén eigenmode induced fast ion transport are found in both positive and negative triangularity with up to 70% central fast ion pressure deficits relative to classical predictions early during the current ramp phase for discharges with 3 MW injected 80 kV neutral beam power. The deficit in both cases is reduced toward zero as the current penetrates and eventually reaches classical levels during current flattop with qmin ≈ 1. Fast ion transport in negative triangularity plasmas measured using a beam modulation technique show very similar levels to those measured in oval plasmas over the range of tested beam powers (Pbeam ≈ 2–7 MW). This similarity is found despite quite different unstable spectra in the highest beam power case where a mixture of coherent and quasi-coherent modes in the TAE frequency range are observed in negative triangularity and a spectrum of more typical narrowband TAEs and RSAEs is observed for the oval plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predictive Engineering Tools for Injection-Molded Long-Carbon-Fiber Thermoplastic Composites (FY 2015 Third Quarterly Report)

During the third quarter of FY 2015, the following technical progress has been made toward project milestones: 1) Magna oversaw the tool build and prepared the molding plan for the complex part of Phase II. 2) PlastiComp hosted a visit by Magna and Toyota on April 23rd to finalize the molding scope and schedule. The plan for molding trials including selection of molding parameters for both LFT and D-LFT for the U-shape complex part was established. 3) Toyota shipped the U-shape complex part tool to Magna on May 28th, 2015. 4) Plasticomp provided 30wt% LCF/PP and 30wt% LCF/PA66 compounded pellets to Magna for molding the complex part. 5) Magna performed preliminary molding trials on June 2nd, 2015 to validate wall thickness, fill profile, tool temperature and shot size requirements for the complex part. 6) Magna performed the first complex part run on June 16th and 17th, 2015 at Magna’s Composite Centre of Excellence in Concord, ON, Canada. Dale Roland of Plasticomp, and Umesh Gandhi of Toyota also attended the molding. 7) Magna discussed and finalized the plan with PNNL and the team for cutting samples from molded parts at selected locations for fiber orientation and length measurements. 8) Magna provided the computer-aided design (CAD) files of the complex parts with and without ribs to PNNL and Autodesk to build the corresponding ASMI models for injection molding simulations. Magna also provided the actual parameters used. 9) Plasticomp’s provided knowledge and experience of molding LCF materials essential to the successful molding of the parts including optimization of fill speed, tool temperatures, and plasticizing conditions for the 30wt% LCF/PP and 30wt% LCF/PA66 materials in both rib and non-rib versions. 10) Magna molded additional parts for evaluation of mechanical property testing including torsional stiffness on June 29th and 30th, 2015 at Magna’s Composite Center of Excellence. 11) Toyota began preparation for the torsion test of the specimens. Preparation of a computer-aided engineering (CAE) model to predict the performance is in progress. 12) Autodesk fixed an error in the implementation of the proper orthogonal decomposition (POD) calculation of fiber length that had caused the ASMI solution to crash and provided an updated build of ASMI containing the fix. 13)Autodesk reviewed and provided feedback for the complex part molding and measurement locations. 14) Autodesk provided support to set up the workflow for ASMI-ABAQUS® analysis, and provided a fix and workaround for a bug in the ASMI-ABAQUS® output command. 15) Autodesk helped build ASMI analysis models for the complex parts with and without ribs. 16) Autodesk worked on improving the orientation prediction accuracy in the shearing layer for 3D meshes based on comparison to measured data of the plaque moldings. 17) PNNL installed a new ASMI version received from Autodesk and performed comparative analyses to assess mid-plane versus 3D fiber length predictions using the full fiber length model and the reduced-order model (ROM) using POD. 18) PNNL presented the project scope, accomplishments, significant results and future plans to DOE and the USCAR Materials Tech Team on June 3rd, 2015. 19) PNNL discussed the cutting of samples from molded parts and finalized a plan with Magna and the team suggesting the sample size, locations and number of samples per location. 20) PNNL and Autodesk built ASMI models for the complex parts with and without ribs, and preliminary analyses of the part with ribs were conducted using the actual molding parameters received from Magna. 21) PNNL worked on a procedure to extract fiber orientation and length results from a 3D ASMI analysis to a 3D ABAQUS model. This procedure is essential to import ASMI fiber orientation and length to a 3D ABAQUS model of the part allowing future part structural analysis for weight reduction study.

36 MATERIALS SCIENCE↗

Operator splitting method for simulation of dynamic flows in natural gas pipeline networks

Here, we develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for mass flow and pressure on a metric graph, where turbulent losses of momentum are modeled by phenomenological Darcy-Weisbach friction. Mass flow balance is maintained through the boundary conditions at the network nodes, where natural gas is injected or withdrawn from the system. Gas flow through the network is controlled by compressors boosting pressure at the inlet of the adjoint pipe. Our operator splitting numerical scheme is unconditionally stable and it is second order accurate in space and time. The scheme is explicit, and it is formulated to work with general networks with loops. We test the scheme over range of regimes and network configurations, also comparing its performance with performance of two other state of the art implicit schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Beam simulation studies of ECR beam extraction and low energy beam transport for FRIB

To meet the beam power requirements of 400 kW at the fragmentation target for facility for Rare Isotope Beams (FRIB), simultaneous acceleration of two-charge states should be used for heavier ions. These intense multi-charged ion beams will be produced by a 28 GHz electron cyclotron resonance (ECR) ion source at a high voltage of 35 kV. After extraction, the ion beam will be pre-accelerated to 12 keV/u with a 50 kV platform, transported down to an achromatic charge state selection (CSS) system followed by a vertical transport line, and then injected into a radio frequency quadrupole accelerator. The TRACK code developed at ANL is used to perform the simulations of the ECR beam extraction and low energy beam transport for FRIB. In this study, we include the magnetic field of ECR ion source into simulations. Different initial beam conditions as well as different space charge neutralization levels are tested for the ECR beamline. The beam loss in CSS system and the corresponding protective measures are discussed. The detailed results about the beam dynamic simulation and beam loss in CSS system will be presented in this paper.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Alfvén eigenmodes and fast ion transport in negative triangularity DIII-D plasmas

The first energetic particle experiments in negative triangularity tokamak plasmas have been carried out on DIII-D. Alfvén eigenmode (AE) activity and associated fast ion transport comparable to that in positive triangularity is observed during the current ramp portion of all plasmas with early beam heating indicating negative triangularity does not confer a special advantage with respect to AE induced transport. In these discharges, a range of mode activity driven by the sub-Alfvénic ( V beam / V a < 0.5) 80 kV neutral beams is found including beta induced Alfvén acoustic eigenmodes, beta induced Alfvén eigenmodes, reversed shear Alfvén eigenmodes (RSAEs) and toroidicity induced Alfvén eigenmodes (TAEs). Mode intermittency and possibly chirping appears to be more common than comparable positive triangularity and/or oval discharges but overall, the unstable spectra and mode amplitudes observed on magnetics, CO2 interferometry, and electron cyclotron emission in a set of matched positive and negative triangularity cases at moderate beam power is similar. Large levels of Alfvén eigenmode induced fast ion transport are found in both positive and negative triangularity with up to 70% central fast ion pressure deficits relative to classical predictions early during the current ramp phase for discharges with 3 MW injected 80 kV neutral beam power. The deficit in both cases is reduced toward zero as the current penetrates and eventually reaches classical levels during current flattop with q min ≈ 1. Fast ion transport in negative triangularity plasmas measured using a beam modulation technique show very similar levels to those measured in oval plasmas over the range of tested beam powers ( P beam ≈ 2–7 MW). Furthermore, this similarity is found despite quite different unstable spectra in the highest beam power case where a mixture of coherent and quasi-coherent modes in the TAE frequency range are observed in negative triangularity and a spectrum of more typical narrowband TAEs and RSAEs is observed for the oval plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗