Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallel Applications”

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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Collaborative Research: A Model of Partially Ionized Plasma Flows with Kinetic Treatment of Neutral Atoms and Nonthermal Ions

Interactions of flows of partially ionized, magnetized plasma are frequently accompanied by the presence of both thermal and non-thermal (pickup) ion components. Such interactions cannot be modeled using traditional MHD equations and require more advanced approaches to treat them. If a nonthermal component of ions is formed due to charge exchange and collisions between the thermal (core) ions and neutrals, it experiences the action of magnetic field, its distribution function is isotropized, and it soon acquires the velocity of the ambient plasma without being thermodynamically equilibrated. This situation, e. g., takes place in the outer heliosphere –- the part of interstellar space beyond the solar system whose properties are determined by the solar wind interaction with the local interstellar medium. This is also possible in laboratory, at million degrees and above, when plasma is conducting electricity far too well, which makes Ohmic heating ineffective. To attain the target temperatures one needs additional heating eventually playing a dominant role. Among such sources is a so-called neutral particle beam heating. This is a wide-spread technique (Joint European Torus and International Thermonuclear Experimental Reactor experiments) based on the injection of powerful beams of neutral atoms into ohmically preheated plasma. In this project we have investigated the energy and density separation between the thermal and nonthermal components in the solar wind and interstellar plasmas. A new model has been developed in which we solve the ideal MHD equations for mixture of all ions and the kinetic Boltzmann equation to describe the transport of neutral atoms. As a separate capability, we can treat the flow of neutral atoms in a multi-component fashion, where neutral atoms born in each thermodynamically distinct region are governed by the Euler gas dynamic equations. We also describe the behavior of pickup ions either kinetically, using the Fokker--Planck equation, or as a separate fluid. Our numerical simulations have demonstrated that pickup ions play a major role in the interaction of the solar wind and (partially ionized) interstellar medium plasmas. Our teams have investigated the stability of the surface (the heliopause) that separates the solar wind from the local interstellar medium, the transport of galactic cosmic rays, the properties of the heliotail flow, and modifications to the bow wave in front of the heliopause due to charge exchange between the neutral H atoms born in the solar wind and interstellar ions. Modeling results have been validated against observational data, such as obtained by the Interstellar Boundary Explorer (IBEX), and made it possible to shed light on the structure of energetic neutral atom maps created by this spacecraft.. We have also demonstrated that charge-exchange modulated heliosphere is a source of anisotropy of the multi-TeV cosmic ray flux observed in a number of Earth-bound air shower experiments. Newly developed codes are implemented within a Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS), a publicly available code being developed by our team for over 12 years. MS-FLUKSS scales well up to 160,000 computing cores and has been ported on major supercomputers in the country. Efficient parallelization and data choreography in the continuum simulation modules are provided by Chombo, an adaptive mesh refinement framework managed by Phillip Colella's team at LBNL. We have implemented in-house, hybrid (MPI+OpenMP) parallelization of the kinetic modules that solve the Boltzmann equation with a Monte Carlo method. Currently, the kinetic modules are being rewritten to take advantage of the modern CPU-GPU supercomputer architecture. The scope of the project allowed us to enhance plasma research and education in such broad, multidisciplinary field as physics of partially ionized plasma and its application to space physics and fusion science. Besides the impact on the modeling of complex physical systems, our approach to computational resource management for complex codes utilizing multiple algorithm technologies appears to be a major advance on current approaches. The development of sophisticated resource management will be essential for all future modeling efforts that incorporate a diversity of scales and physical processes. Our effort provided leadership in promoting computational science and plasma physics within the UAH and FIT campuses and, through the training of a broad spectrum of scientists and engineers, fostering new technologies across the country.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Modeling & Analysis of a Small Hydropower Plant and Battery Energy Storage System Connected as a Microgrid

Puget Sound Energy (PSE) is a utility with headquarters in Bellevue, WA that provides electric and gas service to customers in many regions of western Washington State. PSE is undertaking the implementation of a battery energy storage system pilot project in Glacier, WA, a small town located in the North Cascade Mountains near Mt. Baker Ski Area. Glacier experiences power outages that are both somewhat frequent and lengthy (roughly 2.8 outages per year with a duration of 7.5 hours average). The primary cause of the outages is faults on a long 55 kV transmission line that serves the town of Glacier from the Kendall substation. This approximately 10-mile long line is shown in the upper left of Figure 1 below. The Kendall-Glacier transmission line extends through an area of heavy forest with few reasonable options for increasing the line’s reliability. The substation in Glacier (also shown on figure below) serves over 1000 customers that are mostly residential but includes over 50 small businesses. The town has approximately 250 year-round residents, and a peak in visitors during the winter ski season when storms are also more likely to cause power outages. Figure 1: Proposed Microgrid in Glacier, WA Given the reliability challenges in Glacier, PSE is assessing the feasibility of a 2MW/4.4MWh battery energy storage system (BESS) that would allow for improved reliability. The BESS would also allow for the ability to investigate other applications for utility-scale energy 5 storage. This energy storage system is being considered for operation in “planned islanding” and “microgrid” mode. “Planned islanding” indicates use of the BESS to provide power to part or potentially all of the load in the town of Glacier during an outage. “Microgrid” indicates utilization (during an outage) of the BESS in parallel with Nooksack Falls Hydroelectric Power Plant, a small (nameplate 1.5 MW/2000 HP) hydro power plant less than 10 miles from the town of Glacier. Nooksack is a run-of river small hydro generation plant with a gross head of 226 ft. The site has a 5-foot diameter penstock which enters the powerhouse and goes into a pelton turbine with six runners on one shaft. It was originally constructed by Stone & Webster in 1906, which makes it the second oldest operating power generation facility in Western Washington. It was placed on the National Register of Historic places in December of 1988. The plant operated until 1997 when a fire destroyed the generator. The project continued operations in 2003 after the replacement of the generator and other updates.

13 HYDRO ENERGY↗
DOE OSTI2024

Investigating Cu-Site Doped Cu–Sb–S Nanoparticles Using Photoelectron and Electron Paramagnetic Resonance Spectroscopy

Tetrahedrite (Cu 12 Sb 4 S 13 ) and famatinite (Cu 3 SbS 4 ) are good candidates for green energy applications because they possess promising thermoelectric and photovoltaic properties as well as contain earth-abundant and nontoxic constituents. Herein, X-ray photoelectron spectroscopy (XPS), ultraviolet photoelectron spectroscopy (UPS), and electron paramagnetic resonance spectroscopy (EPR) methods examined inherent electronic properties and interatomic magnetic interactions of Cu-site doped tetrahedrite and famatinite nanomaterials. An energy-efficient modified polyol method was utilized for the synthesis of tetrahedrite and famatinite nanoparticles doped on the Cu-site with Zn, Fe, Ni, Mn, and Co. This is the first parallel study of tetrahedrite and famatinite nanomaterials with XPS, UPS, and EPR methods alongside a systematic analysis of dopant-dependent effects on the electronic structure and magnetic interactions for each material. XPS showed that the Cu and Sb species in tetrahedrite and famatinite possess different oxidation states, while UPS characterization reveals larger dopant-dependent shifts in the work function for tetrahedrite nanoparticles (4.21 to 4.79 eV) than for famatinite nanoparticles (4.57 to 4.77 eV). Finally, all famatinite nanoparticles display an EPR signal, indicating trace amounts of paramagnetic Cu(II) present below the detection limit of XPS. For tetrahedrite, EPR signatures were observed only for the Zn-doped and Mn-doped nanoparticles, suggesting signal broadening from Cu–Cu spin exchange or spin–lattice relaxation. This study demonstrates the complementary nature of XPS and EPR techniques for studying the oxidation states of metals in solid-state nanomaterials. Comparing the electronic and magnetic properties of tetrahedrite and famatinite while studying the impact of dopant incorporation will guide future endeavors in designing sustainable, high-performance materials for renewable energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Scalability of OpenFOAM Density-Based Solver with Runge–Kutta Temporal Discretization Scheme

Compressible density-based solvers are widely used in OpenFOAM, and the parallel scalability of these solvers is crucial for large-scale simulations. In this paper, we report our experiences with the scalability of OpenFOAM’s native rhoCentralFoam solver, and by making a small number of modifications to it, we show the degree to which the scalability of the solver can be improved. The main modification made is to replace the first-order accurate Euler scheme in rhoCentralFoam with a third-order accurate, four-stage Runge-Kutta or RK4 scheme for the time integration. The scaling test we used is the transonic flow over the ONERA M6 wing. This is a common validation test for compressible flows solvers in aerospace and other engineering applications. Numerical experiments show that our modified solver, referred to as rhoCentralRK4Foam, for the same spatial discretization, achieves as much as a 123.2% improvement in scalability over the rhoCentralFoam solver. As expected, the better time resolution of the Runge–Kutta scheme makes it more suitable for unsteady problems such as the Taylor–Green vortex decay where the new solver showed a 50% decrease in the overall time-to-solution compared to rhoCentralFoam to get to the final solution with the same numerical accuracy. Finally, the improved scalability can be traced to the improvement of the computation to communication ratio obtained by substituting the RK4 scheme in place of the Euler scheme. All numerical tests were conducted on a Cray XC40 parallel system, Theta, at Argonne National Laboratory.

Li, Sibo↗
DOE OSTI2020

Validation of a Generalized Formulation for Load Sharing Behavior in Epicyclic Gears for Wind Turbines

In an ideal epicyclic gear set every parallel gear path transmits the same amount of torque. However, it is well known that certain manufacturing variations result in unequal load sharing between the parallel gear paths. Previous works have developed and validated a general closed-form analytical model of this phenomenon that describes the load sharing characteristics of epicyclic gear sets from three to six planets at any torque level. More recently, this analytical model has been reformulated to include the effects of gravity, carrier bearing clearance, and external applied moments, all of which are relevant to most gearboxes and their mounting configuration in horizontal-axis wind turbines. In this paper, the reformulated model is compared to load measurements collected from two similar wind turbine gearboxes with three-planet epicyclic gear sets. The resulting load sharing values are also compared to the mesh load factor requirements in the AGMA 6006 and IEC 61400-4 wind turbine gearbox design standards. Load sharing factors as high as 1.16 at extreme rotor moments and 1.08 with no moment were measured for the gearbox with cylindrical roller bearings, but the load sharing factor remained below 1.10 for the gearbox with tapered roller bearings The results show that in the wind turbine application, load sharing is not equal - even for three-planet systems with a floating central member because of the effects of gravity, rotor moments, and the resulting relative motions among the epicyclic gear components within the carrier bearing clearance.

closed-form analytical model↗
DOE OSTI2023

Surface Effects on Anisotropic Photoluminescence in One‐Dimensional Organic Metal Halide Hybrids

1D organic metal halide hybrids (OMHHs) exhibit strongly anisotropic optical properties, highly efficient light emission, and large Stokes shift, holding promise for novel photodetection and lighting applications. However, the fundamental mechanisms governing their unique optical properties and in particular the impacts of surface effects are not understood. Herein, 1D C 4 N 2 H 14 PbBr 4 by polarization‐dependent time‐averaged and time‐resolved photoluminescence (TRPL) spectroscopy, as a function of photoexcitation energy, is investigated. Surprisingly, it is found that the emission under photoexcitation polarized parallel to the 1D metal halide chains can be either stronger or weaker than that under perpendicular polarization, depending on the excitation energy. The excitation‐energy‐dependent anisotropic emission is attributed to fast surface recombination, supported by first‐principles calculations of optical absorption in this material. The fast surface recombination is directly confirmed by TRPL measurements, when the excitation is polarized parallel to the chains. The comprehensive studies provide a more complete picture for a deeper understanding of the optical anisotropy in 1D OMHHs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Priority-BF: A Task Manager for Priority-Based Scheduling

The increasing demand for computational resources, particularly in High-Performance Computing environments, necessitates to rethink how we handle job scheduling strategies. This work addresses the challenge of managing concurrent jobs with differing priorities on overloaded parallel systems, where strict QoS constraints are often difficult for users to define. Our solution relies on a qualitative description of priorities and pulls from two key approaches: the Easy-BF algorithm and the Conservative Backfilling algorithms. This solution improves the response time for high-priority jobs by 50% without affecting the overall system utilization. We show its applicability in several critical scenarios such as High-Performance Computing (HPC) resource management and in-situ computing.

Gainaru, Ana [ORNL]↗
DOE OSTI2018

Extension of Landau-fluid closure to weakly collisional plasma regime

Both collisionless and collisional Landau fluid (LF) closure have been developed and implemented in two-fluid plasma simulation framework BOUT++, taking full advantage of the fast non-Fourier method (Dimits, 2014). The spectral range of good fit can be conveniently extended to larger domain by adding more Lorentzians. The scaling for upper limit of wavenumber resolved in collisionless LF operator and the fitting coefficients for the collisional LF operator, the later is extended to weakly collisional regime (ν$^{*}_{e}$ ~ 0.01), are presented in this paper. The flux limited expression q$^{FL}_{∥e}$ is reduced to free streaming expression q$^{FS}_{∥e}$ in collisionless limit and gives classical Spitzer–Härm result q$^{SH}_{∥e}$ in collisional limit as we expected. The collisional LF closure recovers the classical Spitzer–Härm result q$^{SH}_{∥e}$ in collisional limit and is same as the collisionless LF closure in collisionless limit. The test cases also obviously reflect the nonlocal effects from LF closures and responding to different boundary conditions in open surface region. Our report will make it more applicable and possible to include the Landau damping kinetic effect in the fluid simulation of tokamak plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Crystallographic Texture, Structure, and Stress Transmission in Nugget Sandstone Examined With X‐Ray Tomography and Diffraction Microscopy

Subsurface processes in sandstones are controlled by porosity, permeability, and deformation mechanisms, all of which are controlled by a complex interplay of crystallographic rock texture, structure, and micromechanics. Texture, structure, and micromechanics have historically been studied using optical and electron microscopy of thin-sections. Here, we employed a new combination of in situ X-ray tomography and ray diffraction microscopy to study crystallographic texture, structure, and grain stresses in 3D. We examined these features in a sample of Nugget sandstone, a sandstone constituting hydrocarbon reservoirs across the American West. Our aims are threefold. First, we demonstrate the utility of X-ray diffraction microscopy probes for revealing texture, structure, and stress transmission in 3D. Second, we apply these techniques to Nugget sandstone and discuss findings in the context of prior work. Third, we study grain stress tensor evolution during mechanical compression to examine whether their heterogeneity and orientation evolution reflect that of inter-particle forces in granular materials. Our results show: (a) larger grains featured higher intra-granular misorientations, possibly from an increased prevalence of cements; (b) pores closed parallel to the loading direction and opened normal to loading; (c) grain stresses featured heterogeneity and orientations similar to inter-particle forces in non-cohesive granular materials; (d) grains featured compressive stresses in the loading direction and tensile stresses orthogonal to the loading direction, the latter resisting sample dilation and grain separation. Our work demonstrates the first known application of multi-modal X-ray tomography and diffraction microscopy to sandstone, providing new 3D insight into the nature of quartz cement and stress evolution.

high energy diffraction microscopy↗
DOE OSTI2026

High-rate heavy-ion tracking using MWPCs and PPACs with an Anti-Discharge Unit

This work reports on the development of two robust, heavy-ion beam tracking concepts operating at low pressure (< 15 Torr) for high rate applications (> 200 kHz). The first concept consists of a Multi-Wire Proportional Counter (MWPC) with a central anode consisting of 12 μm Au-plated Tungsten wires spaced 1 mm from each other. The anode grid is sandwiched between two segmented cathodes aligned orthogonally in the two dimensions for (x,y) particle localization. The second detector is a Parallel-Plate Avalanche Counter (PPAC) that uses the same readout geometry as the MWPC, but replaces the wire anode with a 150 nm silver layer deposited on both sides of a thin (< 1 mg/cm 2 ) polypropylene foil. Additionally, the bias circuitry for the PPAC anode central foil is equipped with an Anti-Discharge Unit (ADU) to prevent transitions from proportional operation to streamer formation, thereby avoiding damaging discharges. The localization capability of both detectors was tested with a low-rate alpha-particle source (241-Am). A position resolution of < 1 mm (FWHM) was achieved under stable, high-gas-gain (> 1000) operating conditions. Their performance in terms of detection efficiency as a function of the isotope charge (Z) was determined by irradiating the detectors with a cocktail beam (Z ≤ 15) with energy of ∼ 100 MeV/u. Full detection efficiency is maintained for all available fragments under optimal operational conditions (i.e., voltage bias). Full detection efficiency was achieved at rates above 200 kHz by irradiating the detectors with a 1 cm diameter 238 U beam at an energy of 143 MeV/u.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Umbrella sampling: a powerful method to sample tails of distributions

We present the umbrella sampling (US) technique and show that it can be used to sample extremely low-probability areas of the posterior distribution that may be required in statistical analyses of data. In this approach sampling of the target likelihood is split into sampling of multiple biased likelihoods confined within individual umbrella windows. We show that the US algorithm is efficient and highly parallel and that it can be easily used with other existing Markov Chain Monte Carlo (MCMC) samplers. The method allows the user to capitalize on their intuition and define umbrella windows and increase sampling accuracy along specific directions in the parameter space. Alternatively, one can define umbrella windows using an approach similar to parallel tempering. We provide a public code that implements US as a standalone PYTHON package. We present a number of tests illustrating the power of the US method in sampling low-probability areas of the posterior and show that this ability allows a considerably more robust sampling of multimodal distributions compared to standard sampling methods. We also present an application of the method in a real world example of deriving cosmological constraints using the supernova type Ia data. We show that US can sample the posterior accurately down to the approximate to 15 sigma credible region in the Omega(m) - Omega(Lambda) plane, while for the same computational effort the affine-invariant MCMC sampling implemented in the emcee code samples the posterior reliably only to approximate to 3 sigma.

cosmology: cosmological parameters↗
DOE OSTI2025

Entanglement-enhanced ac magnetometry in the presence of Markovian noise

Entanglement is a resource to improve the sensitivity of quantum sensors. In an ideal case, using an entangled state as a probe to detect target fields, we can beat the standard quantum limit by which all classical sensors are bounded. However, since entanglement is fragile against decoherence, it is unclear whether entanglement-enhanced metrology is useful in a noisy environment. Its benefit is indeed limited when estimating the amplitude of dc magnetic fields under the effect of parallel Markovian decoherence, where the noise operator is parallel to the target field. In this paper, on the contrary, we show an advantage to using an entanglement over the classical strategy under the effect of parallel Markovian decoherence when we try to detect ac magnetic fields. We consider a scenario to induce a Rabi oscillation of the qubits with the target ac magnetic fields. Although we can, in principle, estimate the amplitude of the ac magnetic fields from the Rabi oscillation, the signal becomes weak if the qubit frequency is significantly detuned from the frequency of the ac magnetic field. We show that, by using the Greenberger-Horne-Zeilinger (GHZ) states, we can significantly enhance the signal of the detuned Rabi oscillation even under the effect of parallel Markovian decoherence. Further, our method is based on the fact that the interaction time between the GHZ states and ac magnetic fields scales as 1/L to mitigate the decoherence effect, where L is the number of qubits, which contributes to improving the bandwidth of the detectable frequencies of the ac magnetic fields. Our results pave the way for new applications of entanglement-enhanced ac magnetometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Discorpy : algorithms and software for camera calibration and correction

Camera or lens-based detector calibration is essential for spatial accuracy in applications like dimensional tomography, optical metrology, and computer vision. Many methods and software exist yet there is still a lack of approaches that achieve both high accuracy and robustness while being easy to use and capable of handling a wide range of distortions. Radial lens distortion is common in high-resolution X-ray detector optics used in parallel-beam tomography at synchrotrons. Achieving sub-pixel accuracy requires calibrating with an optical target image. Although methods for characterizing radial distortion are well established, acquired images often also include perspective distortion and optical center offset. Here, we present our approaches to individually characterize and correct both types of distortion using a single calibration image, implemented in the Discorpy software.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Numbering-Up of Microscale Devices for Megawatt-Scale Supercritical Carbon Dioxide Concentrating Solar Power Receivers

Concentrated solar power (CSP) plants have the potential to reduce the consumption of nonrenewable resources and greenhouse gas emissions in electricity production. In CSP systems, a field of heliostats focuses solar radiation on a central receiver, and energy is then transferred to a thermal power plant at high temperature. However, maximum receiver surface fluxes are low (30–100 W cm -2 ) with high thermal losses, which has contributed to the limited market penetration of CSP systems. Recently, small (~4 cm 2 ), laminated micro pin-fin devices have shown potential to achieve concentrated surface fluxes over 100 W cm -2 using supercritical CO 2 as the working fluid. The present study explores the feasibility of using these microscale unit cells as building blocks for a megawatt-scale (250 MW thermal) open solar receiver through a numbering-up approach, where multiple microscale unit cell devices are connected in parallel. A multiscale model of the full-scale central receiver is developed. The model consists of interconnected unit cell and module level (i.e., multiple unit cells in parallel) submodels which predict local performance of the central receiver. Each full-scale receiver consists of 3000 micro pin-fin unit cells divided into 250 modules. The performance of three different full-scale receivers is simulated under representative operating conditions. The results show that the microscale unit cells have the potential to be numbered up to megawatt applications while providing high heat flux and thermal efficiency. At the design incident flux and surface emissivity, a global receiver efficiency of approximately 90% when heating sCO 2 from 550 °C to 650 °C at an average incident flux of 110 W cm -2 can be achieved.

14 SOLAR ENERGY↗
DOE OSTI2016

SU-F-J-147: Magnetic Field Dose Response Considerations for a Linac Monitor Chamber

Purpose: The impact of magnetic fields on the readings of a linac monitor chamber have not yet been investigated. Herein we examine the total dose response as well as any deviations in the beam parameters of flatness and symmetry when a Varian monitor chamber is irradiated within an applied magnetic field. This work has direct application to the development of Linac-MR systems worldwide. Methods: A Varian monitor chamber was modeled in the Monte Carlo code PENELOPE and irradiated in the presence of a magnetic field with a phase space generated from a model of a Linac-MR prototype system. The magnetic field strength was stepped from 0 to 3.0T in both parallel and perpendicular directions with respect to the normal surface of the phase space. Dose to each of the four regions in the monitor chamber were scored separately for every magnetic field adaptation to evaluate the effect of the magnetic field on flatness and symmetry. Results: When the magnetic field is perpendicular to the phase space normal we see a change in dose response with a maximal deviation (10–25% depending on the chamber region) near 0.75T. In the direction of electron deflection we expectedly see opposite responses in chamber regions leading to a measured asymmetry. With a magnetic field parallel to the phase space normal we see no measured asymmetries, however there is a monotonic rise in dose response leveling off at about +12% near 2.5T. Conclusion: Attention must be given to correct for the strength and direction of the magnetic field at the location of the linac monitor chamber in hybrid Linac-MR devices. Elsewise the dose sampled by these chambers may not represent the actual dose expected at isocentre; additionally there may be a need to correct for the symmetry of the beam recorded by the monitor chamber. Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Collaborative Research: A Model of Partially Ionized Plasma Flows with Kinetic Treatment of Neutral Atoms and Nonthermal Ions

Interactions of flows of partially ionized, magnetized plasma are frequently accompanied by the presence of both thermal and non-thermal (pickup) ion components. Such interactions cannot be modeled using traditional MHD equations and require more advanced approaches to treat them. If a nonthermal component of ions is formed due to charge exchange and collisions between the thermal (core) ions and neutrals, it experiences the action of magnetic field, its distribution function is isotropized, and it soon acquires the velocity of the ambient plasma without being thermodynamically equilibrated. This situation, e. g., takes place in the outer heliosphere - the part of interstellar space beyond the solar system whose properties are determined by the solar wind interaction with the local interstellar medium. This is also possible in laboratory, at million degrees and above, when plasma is conducting electricity far too well, which makes Ohmic heating ineffective. To attain the target temperatures one needs additional heating eventually playing a dominant role. Among such sources is a so-called neutral particle beam heating. This is a wide-spread technique (Joint European Torus and International Thermonuclear Experimental Reactor experiments) based on the injection of powerful beams of neutral atoms into ohmically preheated plasma. In this project we have investigated the energy and density separation between the thermal and nonthermal components in the solar wind and interstellar plasmas. A new model has been developed in which we solve the ideal MHD equations for mixture of all ions and the kinetic Boltzmann equation to describe the transport of neutral atoms. As a separate capability, we can treat the flow of neutral atoms in a multi-component fashion, where neutral atoms born in each thermodynamically distinct regions are governed by the Euler gas dynamic equations. We also describe the behavior of pickup ions either kinetically, using the Fokker–Planck equation, or as a separate fluid. Our numerical simulations have demonstrated that pickup ions play a major role in the interaction of the solar wind and (partially ionized) interstellar medium plasmas. Our teams have investigated the stability of the surface (the heliopause) that separates the solar wind from the local interstellar medium, the transport of galactic cosmic rays, the properties of the heliotail flow, and modifications to the bow wave in front of the heliopause due to charge exchange between the neutral H atoms born in the solar wind and interstellar ions. Modeling results have been validated against observational data, such as obtained by the Interstellar Boundary Explorer (IBEX), and made it possible to shed light on the structure of energetic neutral atom maps created by this spacecraft.. We have also demonstrated that charge-exchange modulated heliosphere is a source of anisotropy of the multi-TeV cosmic ray flux observed in a number of Earth-bound air shower experiments. Newly developed codes are implemented within a Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS), a publicly available code being developed by our team for over 12 years. MS-FLUKSS scales well up to 160,000 computing cores and has been ported on major supercomputers in the country. Efficient parallelization and data choreography in the continuum simulation modules are provided by Chombo, an adaptive mesh refinement framework managed by Phillip Colella’s team at LBNL. We have implemented in-house, hybrid (MPI+OpenMP) parallelization of the kinetic modules that solve the Boltzmann equation with a Monte Carlo method. Currently, the kinetic modules are being rewritten to take advantage of the modern CPU-GPU supercomputer architecture. The scope of the project allowed us to enhance plasma research and education in such broad, multidis- ciplinary field as physics of partially ionized plasma and its application to space physics and fusion science. Besides the impact on the modeling of complex physical systems, our approach to computational resource management for complex codes utilizing multiple algorithm technologies appears to be a major advance on current approaches. The development of sophisticated resource management will be essential for all future modeling efforts that incorporate a diversity of scales and physical processes. Our effort provided leadership in promoting computational science and plasma physics within the UAH and FIT campuses and, through the training of a broad spectrum of scientists and engineers, foster new technologies across the country.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Novel Geometric Operations for Linear Programming

This report summarizes the work performed under the project "Linear Programming in Strongly Polynomial Time." Linear programming (LP) is a classic combinatorial optimization problem heavily used directly and as an enabling subroutine in integer programming (IP). Specifically IP is the same as LP except that some solution variables must take integer values (e.g. to represent yes/no decisions). Together LP and IP have many applications in resource allocation including general logistics, and infrastructure design and vulnerability analysis. The project was motivated by the PI's recent success developing methods to efficiently sample Voronoi vertices (essentially finding nearest neighbors in high-dimensional point sets) in arbitrary dimension. His method seems applicable to exploring the high-dimensional convex feasible space of an LP problem. Although the project did not provably find a strongly-polynomial algorithm, it explored multiple algorithm classes. The new medial simplex algorithms may still lead to solvers with improved provable complexity. We describe medial simplex algorithms and some relevant structural/complexity results. We also designed a novel parallel LP algorithm based on our geometric insights and implemented it in the Spoke-LP code. A major part of the computational step is many independent vector dot products. Our parallel algorithm distributes the problem constraints across processors. Current commercial and high-quality free LP solvers require all problem details to fit onto a single processor or multicore. Our new algorithm might enable the solution of problems too large for any current LP solvers. We describe our new algorithm, give preliminary proof-of-concept experiments, and describe a new generator for arbitrarily large LP instances.

97 MATHEMATICS AND COMPUTING↗