SEARCH · Engineering Papers

Results for “distributed and parallel processing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2024

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

DOE OSTI2023

HBMax: Optimizing Memory Efficiency for Parallel Influence Maximization on Multicore Architectures

The goal of influence maximization is to select k most-influential vertices or seeds in a network, where influence is defined by a given diffusion process. The problem has a number of important applications such as viral marketing, information spread, and epidemic control. Although computing optimal seed set is NP-Hard, due to the submodular nature of the problem efficient approximation algorithms exist. However, even state-of-the-art parallel implementations are limited by a sampling step that incurs large memory footprints. This in turn limits the problem size reach and approximation quality. In this work, we study the memory footprint of the sampling process collecting reverse reachability information in the IMM algorithm over large real-world social networks. We present an adaptive and memory-efficient optimization approach for a state-of-the-art multi-threaded parallel influence maximization algorithm. Our approach,HuffMax, uses a portion of the reverse reachable (RR) sets collected by the algorithm to learn the characteristics of the graph. Then, it compresses the intermediate reverse reachability information with Huffman coding, and queries directly on the compressed data to preserve the memory savings obtained through compression. We also propose an efficient sampling strategy based on the distribution of RR sets, which can further reduce the computation time for typical social networks with long-tail distributions. Considering a NUMA architecture, we scale up our solution on 128-core CPUs and reduce the memory footprint by up to 45.7% with negligible time overhead (or even faster) and without perceivable loss of accuracy.

Chen, Xinyu↗

DOE OSTI2020

10 MWE CDCL Large Pilot Plant – pre-FEED Study

Coal Direct Chemical Looping (CDCL) is an advanced oxy-combustion technology that has the potential to substantially reduce the energy penalty and the cost of electricity (COE) for coal-fired power generation with CO₂ capture. The Babcock & Wilcox Company (B&W) and The Ohio State University (OSU) have been collaborating on the development of an iron oxide oxygen-carrier based chemical looping technology for clean power generation with inherent carbon capture. In this process, coal is dried and pulverized prior to being transported into a moving-bed reducer. In the reducer, coal reacts with the oxygen-carrier particles, forming combustion byproducts, predominantly CO₂ and H₂O, while reducing the iron oxide oxidation state from Fe₂O₃ to a mixture of FeO and Fe. The reduced state particles are then transported to a combustor reactor and re-oxidized with air. Following the oxidation, the oxygen-carrier particles are regenerated, and a large amount of heat is released for steam production. The produced steam is sent to a turbine for electricity generation. Meanwhile, the CO₂–rich stream leaving the reducer is cooled, cleaned, and compressed for subsequent pipeline transportation and sequestration. By combining air separation and fuel conversion into a single system, the CDCL technology enables the intensification of oxy-combustion processes by eliminating the energy and cost intensive cryogenic air separation unit and thereby results in higher overall plant efficiencies and lower COE’s. The use of a moving-bed reducer results in high conversions of volatile hydrocarbons and high CO₂ purity, which reduces the cost of downstream CO₂ purification for sequestration or utilization. The Babcock & Wilcox Company in collaboration with The Ohio State University, Johnson & Matthey, The Electric Power Research Institute, and Dover Light & Power completed a Preliminary Front-End Engineering and Design (Pre-FEED) study of a 10 MWe coal-direct chemical looping (CDCL) pilot plant. The planned system is a modular 10 MWe CDCL large pilot facility consisting in 4 modules of 2.5 MWe each, working in parallel and to be hosted within the current structure at the City of Dover’s Municipal Power Plant. The CDCL system can achieve auto-thermal operation and includes a sub-critical steam cycle for power generation. The pilot system was designed to demonstrate full commercial operation at a reduced scale. The coal distribution per plan area is a representative slice of larger commercial arrangements. The system includes a CO₂ recycle system, but it does not include a compression system. The large pilot includes all the environmental control equipment and oxygen carrier and ash handling systems. As part of the project, the Team performed laboratory testing and carried out multiple pilot test campaigns to obtain design and performance information at the 250 kWth CDCL pilot facility at the Babcock & Wilcox Company’s Research Center. The Team demonstrated sustained operations at designed coal inputs, high coal conversion, high CO₂ purity, heat generation on the combustor, low carbon carryover between reactors and low particle attrition. Emissions generated in the reducer reactor were identified as SO₂ and NO x . The commercial manufacturing cost of oxygen carrier particle was evaluated by JM. A particle manufacturing report was generated and submitted to the DOE. Based on the results from the pilot tests and the pre-FEED design efforts, a techno-economic analysis was performed. The study shows that the CDCL process is a promising carbon-friendly technology capable of producing electricity with high efficiency. The estimated COE of the supercritical CDCL plant is $83.3 / MWh, which meets DOE’s target of less than 30% increase in COE when compared to a supercritical PC plant without CO₂ capture. This is the lowest among the existing carbon capture technologies (post-combustion and oxy-combustion) for fossil fuel power plant. CDCL is evaluated to be the most promising technology for carbon capture from the economic aspect.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Optimizing High Performance Markov Clustering for Pre-Exascale Architectures

HipMCL is a high-performance distributed memory implementation of the popular Markov Cluster Algorithm (MCL) and can cluster large-scale networks within hours using a few thousand CPU-equipped nodes. It relies on sparse matrix computations and heavily makes use of the sparse matrix-sparse matrix multiplication kernel (SpGEMM). The existing parallel algorithms in HipMCL are not scalable to Exascale architectures, both due to their communication costs dominating the runtime at large concurrencies and also due to their inability to take advantage of accelerators that are increasingly popular. In this work, we systematically remove scalability and performance bottlenecks of HipMCL. We enable GPUs by performing the expensive expansion phase of the MCL algorithm on GPU. Additionally, we propose a CPU-GPU joint distributed SpGEMM algorithm called pipelined Sparse SUMMA and integrate a probabilistic memory requirement estimator that is fast and accurate. Furthermore, we develop a new merging algorithm for the incremental processing of partial results produced by the GPUs, which improves the overlap efficiency and the peak memory usage. We also integrate a recent and faster algorithm for performing SpGEMM on CPUs. We validate our new algorithms and optimizations with extensive evaluations. With the enabling of the GPUs and integration of new algorithms, HipMCL is up to 12.4x faster, being able to cluster a network with 70 million proteins and 68 billion connections just under 15 minutes using 1024 nodes of ORNL's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

Design consideration in constructing high performance embedded Knowledge-Based Systems (KBS)

As the hardware trends for artificial intelligence (AI) involve more and more complexity, the process of optimizing the computer system design for a particular problem will also increase in complexity. Space applications of knowledge based systems (KBS) will often require an ability to perform both numerically intensive vector computations and real time symbolic computations. Although parallel machines can theoretically achieve the speeds necessary for most of these problems, if the application itself is not highly parallel, the machine's power cannot be utilized. A scheme is presented which will provide the computer systems engineer with a tool for analyzing machines with various configurations of array, symbolic, scaler, and multiprocessors. High speed networks and interconnections make customized, distributed, intelligent systems feasible for the application of AI in space. The method presented can be used to optimize such AI system configurations and to make comparisons between existing computer systems. It is an open question whether or not, for a given mission requirement, a suitable computer system design can be constructed for any amount of money.

Dalton, Shelly D.↗

NASA NTRS1996

Highly Efficient Amplifier for Ka-Band Communications

An amplifier developed under a Small Business Innovation Research (SBIR) contract will have applications for both satellite and terrestrial communications. This power amplifier uses an innovative series bias arrangement of active devices to achieve over 40-percent efficiency at Ka-band frequencies with an output power of 0.66 W. The amplifier is fabricated on a 2.0- by 3.8-square millimeter chip through the use of Monolithic Microwave Integrated Circuit (MMIC) technology, and it uses state-of-the-art, Pseudomorphic High-Electron-Mobility Transistor (PHEMT) devices. Although the performance of the MMIC chip depends on these high-performance devices, the real innovations here are a unique series bias scheme, which results in a high-voltage chip supply, and careful design of the on-chip planar output stage combiner. This design concept has ramifications beyond the chip itself because it opens up the possibility of operation directly from a satellite power bus (usually 28 V) without a dc-dc converter. This will dramatically increase the overall system efficiency. Conventional microwave power amplifier designs utilize many devices all connected in parallel from the bias supply. This results in a low-bias voltage, typically 5 V, and a high bias current. With this configuration, substantial I(sup 2) R losses (current squared times resistance) may arise in the system bias-distribution network. By placing the devices in a series bias configuration, the total current is reduced, leading to reduced distribution losses. Careful design of the on-chip planar output stage power combiner is also important in minimizing losses. Using these concepts, a two-stage amplifier was designed for operation at 33 GHz and fabricated in a standard MMIC foundry process with 0.20-m PHEMT devices. Using a 20-V bias supply, the amplifier achieved efficiencies of over 40 percent with an output power of 0.66 W and a 16-dB gain over a 2-GHz bandwidth centered at 33 GHz. With a 28-V bias, a power level of 1.1 W was achieved with a 12-dB gain and a 36-percent efficiency. This represents the best reported combination of power and efficiency at this frequency. In addition to delivering excellent power and gain, this Ka-band MMIC power amplifier has an efficiency that is 10 percent greater than existing designs. The unique design offers an excellent match for spacecraft applications since the amplifier supply voltage is closely matched to the typical value of spacecraft bus voltage. These amplifiers may be used alone in applications of 1 W or less, or several may be combined or used in an array to produce moderate power, Ka-band transmitters with minimal power combining and less thermal stress owing to the combination of excellent efficiency and power output. The higher voltage operation of this design may also save mass and power because the dc-dc power converter is replaced with a simpler voltage regulator.

Source record↗

NASA NTRS2020

Assessment of Holographic Microscopy for Quantifying Marine Particle Size and Concentration

Holographic microscopy has emerged as a tool for in situ imaging of microscopic organisms and other particles in the marine environment: appealing because of relatively larger sampling volumes and simpler optical configurations compared to similar imaging systems. However, its quantitative capabilities have remained uncertain, in part because hologram reconstruction and image recognition have required manual operation. Here we assess the quantitative skill of our automated hologram processing pipeline, the CCV Pipeline, to evaluate the size and concentration measurements of environmental and cultured assemblages of marine plankton particles, and microspheres. Over one-million particles, ranging from 10-200 microns equivalent spherical diameter (ESD), imaged by the 4-Deep HoloSea digital inline holographic microscope (DIHM) are analyzed. These measurements were collected in parallel with FlowCam, Imaging FlowCytobot, and manual microscope identification. Once corrections for particle location and nonuniform illumination were developed and applied, the DIHM showed an underestimate in ESD of about 3-10%, but successfully reproduced the size-spectral-slope from environmental samples, and the size distribution of monocultures and microspheres. DIHM concentrations (order 1-1000 particles~ml-1) showed linear agreement (r^2=0.73) with the other instruments, but individual comparisons at times had large uncertainty.

holograpy↗

DOE OSTI2026

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗

NASA NTRS2020

Assessment of holographic microscopy for quantifying marine particle size and concentration

Holographic microscopy has emerged as a tool for in situ imaging of microscopic organisms and other particles in the marine environment: appealing because of the relatively larger sampling volume and simpler optical configuration compared to other imaging systems. However, its quantitative capabilities have so far remained uncertain, in part because hologram reconstruction and image recognition have required manual operation. Here, we assess the quantitative skill of our automated hologram processing pipeline (CCV Pipeline), to evaluate the size and concentration measurements of environmental and cultured assemblages of marine plankton particles, and microspheres. Over 1 million particles, ranging from 10 to 200 μm in equivalent spherical diameter, imaged by the 4‐Deep HoloSea digital inline holographic microscope (DIHM) are analyzed. These measurements were collected in parallel with a FlowCam (FC), Imaging FlowCytobot (IFCB), and manual microscope identification. Once corrections for particle location and nonuniform illumination were developed and applied, the DIHM showed an underestimate in ESD of about 3% to 10%, but successfully reproduced the size spectral slope from environmental samples, and the size distribution of cultures (Dunaliella tertiolecta, Heterosigma akashiwo, and Prorocentrum micans) and microspheres. DIHM concentrations (order 1 to 1000 particles ml(exp −1)) showed a linear agreement (r(exp 2) = 0.73) with the other instruments, but individual comparisons at times had large uncertainty. Overall, we found the DIHM and the CCV Pipeline required extensive manual correction, but once corrected, provided concentration and size estimates comparable to the other imaging systems assessed in this study. Holographic cameras are mechanically simple, autonomous, can operate at very high pressures, and provide a larger sampling volume than comparable lens‐based tools. Thus, we anticipate that these characterization efforts will be rewarded with novel discovery in new oceanic environments.

Noah L. Walcutt↗

NASA NTRS2003

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗

NASA NTRS1982

On the structure of the Io Torus

Since the plasma flow near Io is reduced, neutrals originating in charge exchange are not energetic enough to leave the Jovian system and are therefore distributed over an extensive region, as indicated by the sodium cloud. New ions subsequently created in the distributed neutral atomic cloud as a result of charge exchange or electron impact ionization are picked up by the corotating magnetic field, and the radial current driven by the pickup process cannot close in the Io torus, so that it must instead be connected to the planetary ionosphere by field-aligned currents. These field-aligned currents will flow away from the equator at the outer edge of the neutral cloud, and towards it at the inner edge. It is found that the Jovian ionospheric photoelectrons also cannot supply the current flowing away from the equator, so that torus ions accelerated by a parallel electric field could be involved. The parallel potential drop is large enough to push the torus into the Jovian atmosphere, explaining both the sharp, discontinuous change of flux tube content and ion temperature at L equals 5.6, as well as the generation of Auroral-type hiss at that point.

Goertz, C. K.↗

NASA NTRS1985

Heating of thermal helium in the equatorial magnetosphere - A simulation study

The heating of He(+) ions observed by the GEOS satellite at L of about seven starting from an initial anisotropic proton distribution is studied. A simulation of electromagnetic ion cyclotron waves in a plasma consisting of hot anisotropic protons, a dominant thermal hydrogen plasma, and a minority species of helium ions is utilized. It is shown that the thermal ions are heated to 100 times their initial temperature due to the growth of large-amplitude ion cyclotron waves. He(+) ions are heated mainly in the perpendicular direction, and they are heated more than cold H(+) ions. The He(+) ions are heated in a two-step process. First, they are set into oscillations by the growing wave until some of the ions reach a parallel velocity on the order of the resonant velocity; then, strong heating occurs.

Omura, Y.↗

NASA NTRS2002

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 km or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed-shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗

NASA NTRS2002

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 kin or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed- shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗

NASA NTRS2017

Electron Scattering by High-Frequency Whistler Waves at Earth's Bow Shock

Electrons are accelerated to non-thermal energies at shocks in space and astrophysical environments. While different mechanisms of electron acceleration have been proposed, it remains unclear how non-thermal electrons are produced out of the thermal plasma pool. Here, we report in situ evidence of pitch-angle scattering of non-thermal electrons by whistler waves at Earths bow shock. On 2015 November 4, the Magnetospheric Multiscale (MMS) mission crossed the bow shock with an Alfvn Mach number is approximately 11 and a shock angle of approximately 84deg. In the ramp and overshoot regions, MMS revealed bursty enhancements of non-thermal (0.52 keV) electron flux, correlated with high-frequency (0.2 - 0.4 Omega(sub ce), where Omega(sub ce) is the cyclotron frequency) parallel-propagating whistler waves. The electron velocity distribution (measured at 30 ms cadence) showed an enhanced gradient of phase-space density at and around the region where the electron velocity component parallel to the magnetic field matched the resonant energy inferred from the wave frequency range. The flux of 0.5 keV electrons (measured at 1ms cadence) showed fluctuations with the same frequency. These features indicate that non-thermal electrons were pitch-angle scattered by cyclotron resonance with the high-frequency whistler waves. However, the precise role of the pitch-angle scattering by the higher-frequency whistler waves and possible nonlinear effects in the electron acceleration process remains unclear.

Oka, M.↗

DOE OSTI2023

AI ATAC 1: An Evaluation of Prominent Commercial Malware Detectors

This work presents an evaluation of six prominent commercial endpoint malware detectors, a network malware detector, and a file-conviction algorithm from a cyber technology vendor. The evaluation was administered as the first of the Artificial I ntelligence Applications t o Autonomous Cybersecurity (AI ATAC) prize challenges, funded by / completed in service of the US Navy. The experiment employed 100K files (50/50% benign/malicious) with a stratified distribution of file types, including ~1K zero-day program executables (increasing experiment size two orders of magnitude over previous work). We present an evaluation process of delivering a file to a fresh virtual machine donning the detection technology, waiting 90s to allow static detection, then executing the file and waiting another period for dynamic detection; this allows greater fidelity in the observational data than previous experiments, in particular, resource and time-to-detection statistics. To execute all 800K trials (100K files × 8 tools), a software framework is designed to choreograph the experiment into an automated, time-synced, and reproducible workflow with substantial parallelization. Software with base classes for this framework are provided. A cost-benefit model was configured to integrate the tools’ detection statistics into a comparable quantity by simulating costs of use. This provides a ranking methodology for cyber competitions and a lens for reasoning about the varied statistical results. The results provide insights on state of commercial malware detection.

Bridges, Robert↗

DOE OSTI2024

Laser-assisted radiative recombination in a cold hydrogen plasma

Abstract We study the process of laser-assisted radiative recombination of an electron with a proton in a cold hydrogen plasma employing the semiclassical Kramers’ approach which involves calculation of classical trajectories in combined laser and Coulomb fields and the use of the correspondence principle. Due to the Coulomb focusing effect, recombination is the most effective when the initial electron momentum is parallel to the laser polarization. Orders of magnitude enhancement of the cross section, as compared to the laser-free case, is observed in this case. With increasing angle between the electron momentum and polarization, the recombination cross section drops. However, even after averaging over Maxwellian velocity distribution we obtain a substantial enhancement of the recombination rate constant, as compared to the zero-field case. For the field intensities in the range 30–350 MW cm −2 , the enhancement occurs in the region of the radiation wavelength from 5 to 20 µ m and for the plasma temperature from 20 to 300 K.

Fabrikant, I. I. (ORCID:0000000293849454)↗

NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗