Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel application”

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,099 records · Page 61

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↗

Bayesian model updating with finite element vs surrogate models: Application to a miter gate structural system

Bayesian finite element (FE) model updating using direct model evaluations of large-scale high-fidelity FE models is extremely computationally expensive. Surrogate models can be used as fast emulators of FE models to accelerate the model calibration process. The physics/mechanics-based FE models are still the underpinning behind the surrogate models. Here, this paper evaluates the loss in accuracy and the gain in computational time while performing Bayesian model updating by using surrogate model evaluations compared to using direct FE model evaluations. This evaluation is crucial before entirely relying on surrogate models in model updating for structural health monitoring (SHM) and damage prognosis (DP) purposes. This paper also demonstrates Bayesian updating and surrogate model construction of large-scale high-fidelity FE models of infrastructure systems. In this regard, the miter gate structural system is considered as the testbed structure. Three predominant damage modes (loss of contact between gate and wall, loss of thickness due to corrosion, and loss of tension in the diagonal rods) are considered for model updating purposes. Bayesian model updating is performed using direct FE evaluations by leveraging parallel computing. Two types of surrogates, namely polynomial chaos expansion (PCE) and Gaussian process regression (GPR), are developed for the miter gate. Model updating is performed again using the trained surrogate models, and the updating results are compared with their counterparts obtained using the direct FE evaluation results. The posterior distribution of the FE model parameters obtained using the trained surrogates are sufficiently accurate with respect to the posterior obtained utilizing the direct FE evaluations. In addition, an approximate 4-fold decrease in the computational time was observed when using surrogate model evaluations instead of direct FE evaluations for model updating.

42 ENGINEERING↗

Optimizing High-Throughput Inference on Graph Neural Networks at Shared Computing Facilities with the NVIDIA Triton Inference Server

Abstract With machine learning applications now spanning a variety of computational tasks, multi-user shared computing facilities are devoting a rapidly increasing proportion of their resources to such algorithms. Graph neural networks (GNNs), for example, have provided astounding improvements in extracting complex signatures from data and are now widely used in a variety of applications, such as particle jet classification in high energy physics (HEP). However, GNNs also come with an enormous computational penalty that requires the use of GPUs to maintain reasonable throughput. At shared computing facilities, such as those used by physicists at Fermi National Accelerator Laboratory (Fermilab), methodical resource allocation and high throughput at the many-user scale are key to ensuring that resources are being used as efficiently as possible. These facilities, however, primarily provide CPU-only nodes, which proves detrimental to time-to-insight and computational throughput for workflows that include machine learning inference. In this work, we describe how a shared computing facility can use the NVIDIA Triton Inference Server to optimize its resource allocation and computing structure, recovering high throughput while scaling out to multiple users by massively parallelizing their machine learning inference. To demonstrate the effectiveness of this system in a realistic multi-user environment, we use the Fermilab Elastic Analysis Facility augmented with the Triton Inference Server to provide scalable and high-throughput access to a HEP-specific GNN and report on the outcome.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dual-pace transient heat conduction in vertically aligned carbon nanotube arrays induced by structure separation

Vertically aligned arrays of carbon nanotubes (CNTs) are attractive for a wide range of macroscopic applications which can exploit the remarkable properties of individual nanotubes. In this work, an abnormal behavior of CNT bundles in heat conduction is discovered under the transient electro-thermal characterization. The measured voltage change over the sample shows a dual-pace thermal response (DTR), which could not be fitted using a single thermal diffusivity heat transfer model. Instead, two thermal diffusivities are determined from the DTR phenomenon. Three rounds of cryogenic experiments are conducted to investigate the physics behind DTR phenomenon. It starts to emerge when the temperature is reduced to a certain level. After two rounds of cryogenic experiments, the DTR phenomenon becomes permanent from 295 K to 10 K. The nano-scale structure separation induced by the increased thermal strain leads to macroscale structure separation, which results in two parallel heat conduction paths responsible for the DTR phenomenon. By building a new parallel heat transfer model taking both the transient and steady-state electrical and thermal response into consideration, the area ratio of separated CNTs is determined. Furthermore, this result uncovers the existence of coiled morphology and nano-structure evolution of CNTs under cryogenic state and its effect on thermal and electrical transport, especially on their transient behaviors. Knowledge of the thermal transport in these CNTs arrays is important considering the fact that transient thermal response strongly affects their mechanical, optical, and electrical behaviors.

42 ENGINEERING↗

Efficient First-Order Algorithms for Large-Scale, Non-Smooth Maximum Entropy Models with Application to Wildfire Science

Maximum entropy (MaxEnt) models are a class of statistical models that use the maximum entropy principle to estimate probability distributions from data. Due to the size of modern data sets, MaxEnt models need efficient optimization algorithms to scale well for big data applications. State-of-the-art algorithms for MaxEnt models, however, were not originally designed to handle big data sets; these algorithms either rely on technical devices that may yield unreliable numerical results, scale poorly, or require smoothness assumptions that many practical MaxEnt models lack. In this paper, we present novel optimization algorithms that overcome the shortcomings of state-of-the-art algorithms for training large-scale, non-smooth MaxEnt models. Our proposed first-order algorithms leverage the Kullback–Leibler divergence to train large-scale and non-smooth MaxEnt models efficiently. For MaxEnt models with discrete probability distribution of n elements built from samples, each containing m features, the stepsize parameter estimation and iterations in our algorithms scale on the order of O(mn) operations and can be trivially parallelized. Moreover, the strong ℓ1 convexity of the Kullback–Leibler divergence allows for larger stepsize parameters, thereby speeding up the convergence rate of our algorithms. To illustrate the efficiency of our novel algorithms, we consider the problem of estimating probabilities of fire occurrences as a function of ecological features in the Western US MTBS-Interagency wildfire data set. Our numerical results show that our algorithms outperform the state of the art by one order of magnitude and yield results that agree with physical models of wildfire occurrence and previous statistical analyses of wildfire drivers.

Physics↗

Role of surface diffusion in formation of unique reactivity for graphite oxidation: Time-resolved measurements in a pulsed diffusion reactor

Quantification of oxidation kinetics is essential to develop graphitic materials for diverse applications: from refractories found in gas-cooled nuclear reactors to catalysts needed for chemical manufacturing. In this work, using well-defined highly oriented pyrolytic graphite, low-pressure isotopic transient experiments combined with controlled annealing periods, we resolve the role of surface diffusion and quantify oxidation kinetics with nanomole-precision. We observe an unexpected increase in reactivity following annealing which is explained by the role of surface diffusion increasing the probability for trapping mobile oxygen at more reactive edge sites. Here, the locus of adsorption and spillover to the basal plane is distinct from the trapping location creating a more active oxygen species. Isotopic products reflect the population dynamics of oxygen added at the edge and surface diffusion that relocates basal plane oxygen to more reactive edge sites. Since this process proceeds in parallel with direct oxidation reactions, it is not likely to be observed using steady-state or conventional ‘bulk’ characterization techniques. Our unique time-resolved non-equilibrium measurement in a well-defined transport regime, enables observation of three distinct behaviors: short-term deactivation due to the balance of rates in oxygen supply/product formation, reactivity increases due to surface diffusion and longer-term reactivity increase with oxygen accumulation.

36 MATERIALS SCIENCE↗

MEMS RF accelerators for nuclear energy and advanced manufacturing

Energetic ions are widely used to develop radiation hard structural materials and fuels for advanced nuclear energy. Ion beams are also used to generate neutrons for nuclear materials testing, as well as in a series of high impact, high value adding manufacturing processes (ranging from doping of semiconductors to hardening of materials against wear and to increase bio-compatibility of materials). But to date, high power beams of high energy ions are simply too expensive because they are delivered from large single beam accelerators. We have recently demonstrated that MEMS technology enables massively parallel, low cost batch fabrication of ion beams. We propose to scale intense ion accelerators based on MEMS (micro-electro mechanical systems) to high beam power (>10 kW). With our approach, hundreds to thousands of ion beamlets will be densely packed on silicon wafers. Ions are injected and accelerated across gaps formed in stacks of wafers, leading to uniquely high current densities for intense ion accelerators with variable kinetic energy (0.1 to 10 MeV) and variable ion species. Our team brings together experts in high power ion accelerators from Lawrence Berkeley National Laboratory and leading experts in MEMS technology from Cornell University. We will deliver a disruptive technology of low cost, flexible and massively scalable ion accelerators that will enable the rapid development of radiation hard nuclear materials for advanced nuclear energy and enable new applications in manufacturing.

43 PARTICLE ACCELERATORS↗

Physicochemical control of solvation and molecular assembly of charged amphiphilic oligomers at air-aqueous interfaces

Hypothesis : Understanding the rules that control the assembly of nanostructured soft materials at interfaces is central to many applications. We hypothesize that electrolytes can be used to alter the hydration shell of amphiphilic oligomers at the air-aqueous interface of Langmuir films, thereby providing a means to control the formation of emergent nanostructures. Experiments : Three representative salts – (NaF, NaCl, NaSCN) were studied for mediating the self-assembly of oligodimethylsiloxane methylimidazolium (ODMS-MIM + ) amphiphiles in Langmuir films. The effects of the different salts on the nanostructure assembly of these films were probed using vibrational sum frequency generation (SFG) spectroscopy and Langmuir trough techniques. Experimental data were supported by atomistic molecular dynamic simulations. Findings : Langmuir trough surface pressure – area isotherms suggested a surprising effect on oligomer assembly, whereby the presence of anions affects the stability of the interfacial layer irrespective of their surface propensities. In contrast, SFG results implied a strong anion effect that parallels the surface activity of anions. These seemingly contradictory trends are explained by anion driven tail dehydration resulting in increasingly heterogeneous systems with entangled ODMS tails and appreciable anion penetration into the complex interfacial layer comprised of headgroups, tails, and interfacial water molecules. These findings provide physical and chemical insight for tuning a wide range of interfacial assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Slow Electron Holes in the Earth's Magnetosheath

Here, we present a statistical analysis of electrostatic solitary waves observed aboard Magnetospheric Multiscale spacecraft in the Earth's magnetosheath. Applying single-spacecraft interferometry to several hundred solitary waves collected in about 2-minute interval, we show that almost all of them have the electrostatic potential of positive polarity and propagate quasi-parallel to the local magnetic field with plasma frame velocities of the order of 100 km/s. The solitary waves have typical parallel half-widths from 10 to 100 m that is between 1 and 10 Debye lengths and typical amplitudes of the electrostatic potential from 10 to 200 mV that is between 0.01% and 1% of local electron temperature. The solitary waves are associated with quasi-Maxwellian ion velocity distribution functions, and their plasma frame velocities are comparable with ion thermal speed and well below electron thermal speed. We argue that the solitary waves of positive polarity are slow electron holes and estimate the time scale of their acceleration, which occurs due to interaction with ions, to be of the order of one second. The observation of slow electron holes indicates that their lifetime was shorter than the acceleration time scale. We argue that multi-spacecraft interferometry applied previously to these solitary waves is not applicable because of their too-short spatial scales. The source of the slow electron holes and the role in electron-ion energy exchange remain to be established.

79 ASTRONOMY AND ASTROPHYSICS↗

First-order crosstalk mitigation in parallel quantum gates driven with multi-photon transitions

Here, we demonstrate an order of magnitude reduction in the sensitivity to optical crosstalk for neighboring trapped-ion qubits during simultaneous single-qubit gates driven with individual addressing beams. Gates are implemented via two-photon Raman transitions, where crosstalk is mitigated by offsetting the drive frequencies for each qubit to avoid first-order crosstalk effects from inter-beam two-photon resonance. The technique is simple to implement, and we find that phase-dependent crosstalk due to optical interference is reduced on the most impacted neighbor from a maximal fractional rotation error of 0.185(4) without crosstalk mitigation to ≤0.006 with the mitigation strategy. Furthermore, we characterize first-order crosstalk in the two-qubit gate and avoid the resulting rotation errors for the arbitrary-axis Mølmer–Sørensen gate via a phase-agnostic composite gate. Finally, we demonstrate holistic system performance by constructing a composite CNOT gate using the improved single-qubit gates and phase-agnostic two-qubit gate. This work is done on the Quantum Scientific Computing Open User Testbed; however, our methods are widely applicable for individual addressing Raman gates and impose no significant overhead, enabling immediate improvement for quantum processors that incorporate this technique.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of pre-training large language models on leadership-class supercomputers

Large language models (LLMs) have arisen rapidly to the center stage of artificial intelligence as the foundation models applicable to many downstream learning tasks. However, how to effectively build, train, and serve such models for many high-stake and first-principle-based scientific use cases are both of great interests and of great challenges. Moreover, pre-training LLMs with billions or even trillions of parameters can be prohibitively expensive not just for academic institutions, but also for well-funded industrial and government labs. Furthermore, the energy cost and the environmental impact of developing LLMs must be kept in mind. Here, in this work, we conduct a first-of-its-kind performance analysis to understand the time and energy cost of pre-training LLMs on the Department of Energy (DOE)’s leadership-class supercomputers. Employing state-of-the-art distributed training techniques, we evaluate the computational performance of various parallelization approaches at scale for a range of model sizes, and establish a projection model for the cost of full training. Our findings provide baseline results, best practices, and heuristics for pre-training such large models that should be valuable to HPC community at large. We also offer insights and optimization strategies for using the first exascale computing system, Frontier, to train models of the size of GPT-3 and beyond.

97 MATHEMATICS AND COMPUTING↗

Breakdown and Discharge Characteristics of High Repetition Frequency Nanosecond Pulsed Air Dielectric Barrier Discharge

This study explores the breakdown characteristics and discharge modes in a parallel plate air dielectric barrier discharge system using repetitive nanosecond pulses, with pulse repetition frequencies (PRFs) from 0.1 to 100 kHz. It examines how pulse parameters—PRF, pulse number and gas pressure—affect the memory effect, leveraging current and voltage measurements, fast imaging and optical emission spectroscopy. The findings show that higher PRFs lead to a reduction in breakdown voltage well below the streamer breakdown threshold. We argue that this effect may be attributed to the cumulative buildup of metastable species and negative ions in the discharge gap which could sustain free electrons in‐between the voltage pulses. Despite observing increased energy deposition with higher PRFs, the impact on filament formation was minimal, highlighting a strong dependence of discharge morphology on the accumulation of plasma‐produced species. This research studies provides valuable insights for controlling discharge regimes in applications such as plasma‐assisted combustion, surface treatment and air treatment by clarifying the interactions between discharge mechanisms at different PRFs and pressures.

42 ENGINEERING↗

Mechanisms for the Generation of Complex Fracture Networks: Observations From Slant Core, Analog Models, and Outcrop

We use observations of hydraulic fractures in core, outcrop attributes of natural hydraulic fractures, and analogue models, to address how hydraulic fracture networks evolve. A slant core from the Wolfcamp Formation—an unconventional shale hydrocarbon reservoir in the Permian Basin of West Texas—collected within 18 m and 30 m of two hydraulically stimulated horizontal wells, provided an opportunity to examine hydraulic fractures directly. In approximately 183 m of core, 309 calcite-sealed natural opening-mode fractures and 375 hydraulic fractures were identified. Many hydraulic fractures in the core show complex morphology, including twist-hackle segmentation, diversion, and bifurcation; these structures most commonly develop at lithological bed boundaries and mechanical heterogeneities such as natural fractures and concretions. An outcrop of bed-parallel pavements in the Cretaceous Boquillas Formation in West Texas contains opening-mode fractures that likely formed by natural hydraulic fracturing. Fracture traces provide evidence of twist-hackle segmentation, and are typically associated with bed boundaries and preexisting bed-parallel stylolites. A laboratory study of hydraulic fracturing of 33 synthetic blocks of gypsum and hydrostone revealed fracture steps, diversions, twist hackles, and multiple overlapping fractures together with information on fracture growth directions. These complexities in the fracture network were dominantly nucleated at inclusions used to simulate pre-existing fractures, and as a result of mechanical heterogeneity introduced by the wellbore and perforations. Collectively, our results show that complex fracture networks are produced in hydraulic fracturing of self-sourced reservoir strata. Mechanical stratigraphic boundaries and other heterogeneities are likely to enhance fracture network complexity through the processes of segmentation, diversion, and bifurcation. These processes create multiple fracture strands, resulting in an increased number of hydraulic fractures over those initiated, thereby increasing total fracture surface area. Our study provides insight into hydraulic fracture network propagation, and has applications for evaluation, completion, production, and fracture modeling of unconventional reservoirs.

54 ENVIRONMENTAL SCIENCES↗

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗

High performance power handling in the absence of an H-mode edge in negative triangularity DIII-D plasmas

Abstract Experiments performed during strongly-shaped high-power diverted negative triangularity (NT) experiments in DIII-D achieved detached divertor conditions and a transient-free edge, showcasing the potential for application of NT to a core-edge integrated reactor-like scenario and providing the first characterization of the parametric dependencies for detachment onset. Detached divertor conditions will be required in future devices to mitigate divertor heat fluxes. Access to dissipative divertor conditions was investigated via an increase in upstream density. Detachment onset at the outer strike point was achieved with H-mode level confinement H 98 − y 2 ∼ 1 and reactor-relevant normalized pressures β N ∼ 2 . Confinement degradation was observed with deeper detachment, associated with the loss of an electron temperature pedestal. Differences in geometry, radial transport, impact of cross field drifts are discussed to explain differences in access to detachment in NT discharges. Higher normalized densities, with respect to equivalent discharges in positive triangularity, were necessary to achieve detachment, partially explained by the shorter parallel connection length to the targets. The effect of cross-field particle drifts ( E × B , B × ∇ B ) on access to detachment was demonstrated by the lower upstream density needed to access detachment with ion B × ∇ B drift directed outside of the active divertor (Greenwald fraction f Gw ∼ 0.9–1.0 vs f Gw ∼ 1.3). The upstream density at detachment onset was observed to increase linearly with plasma current with ion B × ∇ B drift into the divertor, consistent with the observed narrowing of the scrape-off layer heat flux width λ q . Edge fluid simulations capture separatrix densities needed to achieve detachment in NT plasma and their dependence on drift direction. The ability to reproduce detachment dynamics in NT plasma increases the confidence in future design studies for NT divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Leveraging gradient weighted class activation mapping to improve classification effectiveness: Case study in transportation infrastructure characterization

Roadway “corners†are common for pedestrian use, whether designated with markings or not. Different types of markings have been deployed, ranging from simple parallel lines to more complex designs. Understanding the impact of different types of crosswalks is important for public safety. In this work we explore methods to improve the logging of marked crosswalk types. We used the Roadway Information Database from the Second Strategic Highway Research Project and used active learning methods with transfer learning to identify the crosswalk types (marked or unmarked). Upon completion we found our classifiers were unable to perform above roughly 94% correct classifications. To improve their efficacy, we separated the crosswalks into their “fine grained†types and used Gradient-Weighted Class Activation Mapping to isolate and study the features that classified the crosswalks. We compared this with sampled manually marked crosswalks and present findings. We believe this use case can represent a process to improve the active learning method for some visual machine learning applications.

Karnowski, Thomas↗

Direct Air Reactive Capture and Conversion for Utility-Scale Energy Storage (Final Report)

This final report for FEW0277 summarizes the work performed over the project performance period of October 2021 – March 2025. This project was funded under the “Reactive Capture and Conversion R&D” lab call released in FY2021. The goal of the project was to develop dual-function materials and process for capturing CO 2 from the atmosphere and converting it into CH 4 . The work was organized into four parallel tracks in 1) direct air capture materials synthesis and characterization, 2) catalysts for CO 2 conversion, 3) mechanistic investigations via ab initio simulations, and 4) process modeling, technoeconomic analysis, and lifecycle assessment. The project was split into two budget periods. The first budget period focused on development of amine-based materials, due to their known performance for CO 2 direct air capture and their potential to act synergistically with metal catalysts to enable a low-temperature methanation pathway. The second budget period focused on development of alkali-based materials and a simulated-moving-bed process for high conversion catalytic reduction of captured CO 2 to CH 4 . All project milestones were completed during the project performance period and are summarized in this report. Our work resulted in publication of eight peer-reviewed manuscripts, one patent application, and numerous presentations given at domestic and international conferences and invited academic department seminars.

03 NATURAL GAS↗