Engineering PapersSearch

SEARCH · Engineering Papers

Results for “eagles”

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 37 records · Page 2

Beyond the Surface: Non-Invasive Low-Field NMR Analysis of Microbially-Induced Calcium Carbonate Precipitation in Shale Fractures

Microbially-induced calcium carbonate precipitation (MICP) is a biological process in which microbially-produced urease enzymes convert urea and calcium into solid calcium carbonate (CaCO 3 ) deposits. MICP has been demonstrated to reduce permeability in shale fractures under elevated pressures, raising the possibility of applying this technology to enhance shale reservoir storage safety. For this and other applications to become a reality, non-invasive tools are needed to determine how effectively MICP seals shale fractures at subsurface temperatures. In this study, two different MICP strategies were tested on 2.54 cm diameter and 5.08 cm long shale cores with a single fracture at 60 °C. Flow-through, pulsed-flow MICP-treatment was repeatedly applied to Marcellus shale fractures with and without sand (“proppant”) until reaching approximately four orders of magnitude reduction in apparent permeability, while a single application of polymer-based “immersion” MICP-treatment was applied to an Eagle Ford shale fracture with proppant. Low-field nuclear magnetic resonance (LF-NMR) and X-Ray computed microtomography (micro-CT) techniques were used to assess the degree of biomineralization. With the flow-through approach, these tools revealed that while CaCO 3 precipitation occurred throughout the fracture, there was preferential precipitation around proppant. Without proppant, the same approach led to premature sealing at the inlet side of the core. In contrast, immersion MICP-treatment sealed off the fracture edges and showed less mineral precipitation overall. This study highlights the use of LF-NMR relaxometry in characterizing fracture sealing and can help guide NMR logging tools in subsurface remediation efforts.

MICP

On the optimal sizing of power take-off systems for wave energy converters

The power take-off (PTO) system in a wave energy converter (WEC) is the means by which the energy in an ocean wave is converted into useful energy. There has been a general lack of technology convergence in the wave energy industry and PTO design has been largely device-specific with the design process determined by the WEC developer. An optimized PTO design can significantly affect the efficiency, reliability, performance, and overall cost and viability of a device. There remains a need in the marine energy industry for best practices and efficient design processes to be identified as it could greatly benefit developers and researchers in designing an optimal WEC PTO. This research utilizes the open-source Matlab-Simulink based software, WEC-Sim, which was developed by the National Renewable Energy Laboratory, and the Eagle, a high-performance computing system, to model and optimize two WECs of different archetypes and then investigates if any trends exist that can be exploited for greater efficiency in optimal PTO design. The results from this research indicate that there may be an optimal sizing for WEC PTOs for both power and PTO force rating and that substantial reductions can be made in the cost of a device without significant loss in the energy produced. Furthermore this optimal sizing may be independent of WEC type and deployment location. With an optimal, or near optimal, approach to WEC PTO sizing demonstrated, a methodology is proposed to address the challenge of nameplate ratings in the wave energy industry.

13 HYDRO ENERGY

Seasonality and Declining Intensity of Methane Emissions from the Permian and Nearby US Oil and Gas Basins

We quantify weekly methane emissions and trends from oil and gas production in the US Permian Basin for 2019-2023, and in nearby basins for 2022-2023, by analytical inversion of Tropospheric Monitoring Instrument (TROPOMI) satellite observations with the Integrated Methane Inversion (IMI) at 25 km resolution. Permian oil and gas emissions averaged 4.0 ± 1.1 Tg a-1 over 2019-2023, with large seasonal variation but little interannual variability. Methane intensity fell from 5.2 to 3.2% as production surged. Intensity in the New Mexico Permian fell from 4.5 to 2.1%, approaching the state's 2026 target of <2%. Emissions were on average 50 ± 10% higher in winter than summer, which we corroborate with Permian Basin Tower Network measurements, Insight M aircraft data, and GHGSat satellite observations. This seasonality may be driven in part by higher winter emissions from liquid storage tanks due to decreased separator efficiency in cold conditions. Similar but weaker seasonality along with decreasing emissions and intensities is found in weekly inversions for the Anadarko, Barnett, Eagle Ford, and Haynesville basins in 2022-2023. Our work suggests that better weatherization of oil and gas facilities could significantly reduce methane emissions.

Permian

Experimental Characterization of Hydrogen Diffusion in Shale Rocks for Geologic Storage Applications

As global energy systems undergo a transition to cleaner alternatives, geologic hydrogen storage has emerged as a promising solution for large-scale energy storage. A critical factor in determining the feasibility of this approach is the effectiveness of caprock formations, such as shale, in preventing hydrogen migration. This study investigates the diffusion behavior of hydrogen through shale to assess its suitability as a caprock for geologic hydrogen storage. Using a novel double-seal core holder design and a through-diffusion apparatus, hydrogen diffusion was measured through shale rock from the Eagle Ford and Wolfcamp Formations under dry conditions. These measurements were complemented by microstructural and mineralogical analyses using low-pressure nitrogen adsorption and X-ray diffraction. The effective diffusion coefficient of hydrogen in these shale caprocks ranged from 2.51 × 10 –8 to 9.85 × 10 –8 m 2 /s. Notably, we observed that the diffusion behavior was more related to the pore network structure and could not be attributed to differences in the total pore volume between shale types alone. Here, to further understand the role of pore network complexity, a fractal pore model was developed to correlate tortuosity with the fractal dimension of the pore structure (a measure of pore network complexity). The proposed model closely matched tortuosity values obtained from diffusion experiments, outperforming existing theoretical tortuosity–porosity correlations. These findings provide key quantitative parameters needed to assess the feasibility of geologic hydrogen storage as well as insights that can be applied to hydrogen storage in a range of geologic formations.

08 HYDROGEN

A single diiron enzyme catalyses the oxidative rearrangement of tryptophan to indole nitrile

Abstract Nitriles are uncommon in nature and are typically constructed from oximes through the oxidative decarboxylation of amino acid substrates or from the derivatization of carboxylic acids. Here we report a third nitrile biosynthesis strategy featuring the cyanobacterial nitrile synthase AetD. During the biosynthesis of the eagle-killing neurotoxin, aetokthonotoxin, AetD transforms the 2-aminopropionate portion of 5,7-dibromo-l-tryptophan to a nitrile. Employing a combination of structural, biochemical and biophysical techniques, we characterized AetD as a non-haem diiron enzyme that belongs to the emerging haem-oxygenase-like dimetal oxidase superfamily. High-resolution crystal structures of AetD together with the identification of catalytically relevant products provide mechanistic insights into how AetD affords this unique transformation, which we propose proceeds via an aziridine intermediate. Our work presents a unique template for nitrile biogenesis and portrays a substrate binding and metallocofactor assembly mechanism that may be shared among other haem-oxygenase-like dimetal oxidase enzymes.

Chemistry

The S-PLUS Fornax Project (S+FP): A first 12-band glimpse of the Fornax galaxy cluster

ABSTRACT The Fornax galaxy cluster is the richest nearby (D ∼ 20 Mpc) galaxy association in the southern sky. As such, it provides a wealth of opportunities to elucidate on the processes where environment holds a key role in transforming galaxies. Although it has been the focus of many studies, Fornax has never been explored with contiguous homogeneous wide-field imaging in 12 photometric narrow and broad bands like those provided by the Southern Photometric Local Universe Survey (S-PLUS). In this paper, we present the S-PLUS Fornax Project (S+FP) that aims to comprehensively analyse the galaxy content of the Fornax cluster using S-PLUS. Our data set consists of 106 S-PLUS wide-field frames (FoV∼1.4 × 1.4 deg2) observed in five Sloan Digital Sky Survey-like ugriz broad bands and seven narrow bands covering specific spectroscopic features like [O ii], Ca ii H+K, Hδ, G band, Mg b triplet, Hα, and the Ca ii triplet. Based on S-PLUS specific automated photometry, aimed at correctly detecting Fornax galaxies and globular clusters in S-PLUS images, our data set provides the community with catalogues containing homogeneous 12-band photometry for ∼3 × 106 resolved and unresolved objects within a region extending over ∼208 deg2 (∼5 Rvir in RA) around Fornax’ central galaxy, NGC 1399. We further explore the eagle and IllustrisTNG cosmological simulations to identify 45 Fornax-like clusters and generate mock images on all 12 S-PLUS bands of these structures down to galaxies with M⋆ ≥ 108 M⊙. The S+FP data set we put forward in this first paper of a series will enable a variety of studies some of which are briefly presented.

Astronomy & Astrophysics

Evaluating quenching in cosmological simulations of galaxy formation with spectral covariance in the optical window

ABSTRACT Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of the observables. Here, we follow an alternative approach based on the spectral covariance in a model-independent way. We build upon previous work by Sharbaf et al. that studied the covariance of high-quality SDSS (Sloan Digital Sky Survey) continuum-subtracted spectra in a relatively narrow range of velocity dispersion ($\sigma \in [100,150]$ km s$^{-1}$). Here, the same analysis is applied to synthetic data from the eagle and IllustrisTNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral covariance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age–metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with active galactic nucleus (AGN) feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid prescriptions, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasize that this methodology only relies on the spectral covariance to assess whether simulations provide a true representation of galaxy formation.

Sharbaf, Z. (ORCID:0009000450545946)

Estimating Critical Customer Outages Resulting from Extreme Hurricanes

US power outage data has been collected by organizations such as Oak Ridge National Laboratory (ORNL) through Environment for Analysis Geo-Located Energy Infrastructure (EAGLE-I: freely available) and poweroutage.us (commercial data: available to purchase). However, these sources do not provide information specific to outages of critical customers. Critical customers include entities, facilities, and individuals whose continuous access to electricity is essential for public safety, emergency response, disaster recovery, the well-being of vulnerable populations, public safety and order, and public utilities such as natural gas, communications, water and sanitation. Identification and geolocation of critical customers is crucial for understanding and addressing the effects of power outages on essential services and ensuring that necessary measures are taken to maintain their operations during power disruptions. This work is a first step towards estimating the occurrences of critical customer outages and developing a critical customer power outage data repository. This work estimates outage incidents of critical customers through spatiotemporal mapping of power outage data, weather data, building data, and critical infrastructure network data. Our results show that critical customer effects vary across different counties. We provide appropriate mathematical explanations and simplifications to define and systematize the proposed approach.

Bhusal, Narayan [ORNL] (ORCID:0000000222752145)

MAPLE v.1.0

SAND2025-00659O MAPLE is a software tool that uses epigenomic data to predict gene expression. MAPLE uses a set of epigenomic modifications to determine the effect on gene expression in a specific subset of species. The algorithm can be trained on additional species and epigenomic modifications, enhancing its predictive capabilities. EAGLE employs a hybrid neural network architecture, featuring a convolutional front-end and a multi-head attention layer, to process pre-processed signal data as input. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Davis IV, Warren

Data-driven Community-centered Resilient Assessment and Planning Toolkit for Nexus of Energy and Water (DCRAPT-NEW)

Urban areas, including Detroit and Pittsburgh, have suffered significant dual outages of the electrical and water infrastructure in the past decade due, in part, to the increasing number of extreme weather events. With increasing temperatures and rainfall intensity, these regions need to prepare for increasing extreme events through community-based energy and water resilience analysis, planning, and enhancement. This project developed a suite of open-source, open-access, community-centered, data-driven assessment and distributed energy resource (DER) and planning tools for energy and water resilience enhancement in urban areas. Through establishing a multi-level community awareness and engagement mechanism and a comprehensive collection of power outage and flooding data, an innovative group of community energy and water resilience assessment and planning tools have been developed for a wide range of users with differing and variable sets of data available to them. The developed tools include (1) DOE EAGLE-I data-driven, deep-learning assisted resilience assessment and DER planning tools at the county level with socioeconomic factors incorporated; (2) Utility annual power outage data-driven tools for long term resilience assessment and DER planning and 15-min power outage data-driven tools for short term resilience assessment and planning; (3) Detailed engineering tools for energy and water systems resilience assessment and planning when the system topology and component fragility curves are available; (4) Alternative Resiliency Metric Calculation that extracts and separates outage and restoration processes; and (5) Co-optimization tools that evaluate the resilience of the power and sewage system and allow users to conduct joint planning with energy and wastewater systems. The developed tools provide planners, decision-makers, and stakeholders with powerful capabilities to systematically evaluate system/community resilience and optimal and actionable guidance for enhancing resilience while prioritizing DER investments. The tools have been used and validated in Detroit and Pittsburgh and can be used in other areas of the nation. In addition, this project will (1) advance the knowledge and applications of machine-learning methods in analyzing and fusing different layers of information and generating meaningful data points such as generating rare weather events; (2) significantly improve the energy and water resilience of the identified communities in Detroit and Pittsburgh and prepare for more frequent and severe weather conditions; (3) help communities assess extreme weather event impacts and address short-term and long-term resilience-related issues The developed tools have been made public via GitHub and demonstrated to community stakeholders and utility companies via the two annual workshops and numerous community engagement meetings. The project outcomes are also disseminated through publications in various journals and conference proceedings, and presentations at top conferences.

13 HYDRO ENERGY

Microwave-Assisted Dehydroaromatization of Flare Gas: Reactor Modeling, Plant-Wide Simulation and Economic Feasibility Analysis

Flaring is widely practiced in the oil, gas, and petrochemical sectors to ensure safety during upsets and maintenance but emits large amounts of GHGs, causing energy and economic losses. In the U.S., about one-third of Bakken gas (~250 MMSCFD) and ~100 MMSCFD from Eagle Ford are flared. Recovering this gas is essential for sustainability. Existing recovery methods—compression and reinjection (EOR), conversion to NGL, LNG/CNG, GTL, and GTW—are often limited by flowrate, composition, and variability, especially in unconventional wells. This study develops a microwave-assisted dehydroaromatization (DHA) process to convert flare gas into benzene, toluene, ethylene, and naphthalene. A laboratory reactor model is scaled up into a modular plant-wide system. Techno-economic (TEA) and life-cycle (LCA) analyses evaluate performance and sustainability, with sensitivity studies on plant capacity, electricity cost, and catalyst price confirming strong economic potential.

dehydroaromatization

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

Christensen, Matthew [Pacific Northwest National L

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

54 ENVIRONMENTAL SCIENCES

On a Critical Acceleration Scale of Dark Matter in ΛCDM and Dynamical Dark Energy

Abstract Universal acceleration a 0 emerges in various empirical laws, yet its fundamental nature remains unclear. Using Illustris and Virgo N -body simulations, we focus on the velocity and acceleration fluctuations in collisionless dark matter involving long-range gravity. For comparison, in the kinetic theory of gases, molecules undergo random elastic collisions involving short-range interactions, where only velocity fluctuations are relevant. Hierarchical structure formation proceeds through the merging of smaller halos to form larger halos, which facilitates a continuous energy cascade from small to large halos at a constant rate ε u ≈ −10 −7 m 2 s −3 . Velocity fluctuations involve a critical velocity u c ∝ (1 + z ) −3/4 . Acceleration fluctuations involve a critical acceleration a c ∝ (1 + z ) 3/4 . Two critical quantities are related by the rate of energy cascade ε u ≈ − a c u c /[2(3 π ) 2 ], where factor 3 π is from the angle of incidence during merging. With critical velocity u c on the order of 300 km s −1 at z = 0, the critical acceleration is determined to be a c 0 ≡ a c ( z = 0) ≈ 10 −10 m s −2 , suggesting a c might explain the universal acceleration a 0 ≈ 10 −10 m s −2 in the empirical Tully–Fisher relation or modified Newtonian dynamics. The redshift evolution a c ∝ (1 + z ) 3/4 is in good agreement with Magneticum and EAGLE simulations and in reasonable agreement with limited observations. This suggests a larger a 0 at a higher redshift such that galaxies of fixed mass rotate faster at a higher redshift. Note that for dark energy (DE) density ρ DE 0 ≈ a c 0 2 / G = 1 0 − 10 J m −3 , we postulate an entropic origin of the DE from acceleration fluctuations of dark matter, analogous to the gas pressure from velocity fluctuations. This leads to a dynamical DE coupled to the structure evolution involving a relatively constant DE density followed by a slow weakening phase, suggesting possible deviations from the standard ΛCDM paradigm.

N-body simulations

Correlating Power Outage Spread with Infrastructure Interdependencies During Hurricanes

Power outages caused by extreme weather events, such as hurricanes, can significantly disrupt essential services and delay recovery efforts, underscoring the importance of enhancing our infrastructure's resilience. This study investigates the spread of power outages during hurricanes by analyzing the correlation between the network of critical infrastructure and outage propagation. We leveraged datasets from Hurricanemapping.com, the North American Energy Resilience Model Interdependency Analysis (NAERM-IA), and historical power outage data from the Oak Ridge National Laboratory (ORNL)'s EAGLE-I system. Our analysis reveals a consistent positive correlation between the extent of critical infrastructure components accessible within a certain number of steps (k-hop distance) from initial impact areas and the occurrence of power outages in broader regions. This insight suggests that understanding the interconnectedness among critical infrastructure elements is key to identifying areas indirectly affected by extreme weather events.

Bose, Avishek

Dark matter decay signals in cosmic filaments

Cosmic filaments form the backbone of the cosmic web, yet their properties and evolution remain uncertain. Using the EAGLE hydrodynamical simulations, we investigate the dark matter density profiles in filaments and their implications for dark matter decay signals. We show that GeV-scale dark matter particles decaying into electron-positron pairs can produce detectable radio synchrotron emission. By leveraging stacked filament radio data, we place stringent constraints on the dark matter decay lifetime, improving existing limits by up to two orders of magnitude for strong filamentary magnetic fields.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Validation of the NLR Pumped Storage Hydropower Cost Model

The National Laboratory of the Rockies (NLR) first released its pumped storage hydropower (PSH) cost model in 2023 as the most detailed bottom-up PSH cost model available to the public. It is available both as a spreadsheet and an interactive web tool, enabling users with a variety of PSH interests to transparently characterize costs of alternative PSH sites and designs. The PSH cost model cannot replace detailed site-level studies and design, but it is important to validate it against other industry PSH cost estimates. The initial model methodology report validated the cost model for a single proposed site, the Eagle Mountain Project in California. This slide deck documents an expanded validation exercise using cost data from six other sites: Goldendale (Washington), Seminoe (Wyoming), Gordon Butte (Montana), Swan Lake (Oregon), White Pine (Oregon), and Lewis Ridge (Kentucky). It compares itemized costs from Federal Energy Regulatory Commission (FERC) applications and other reported costs with NLR PSH cost model outputs after customizing inputs for each site. The validation exercise finds that the NLR model's conservative indirect cost assumptions often drive overall cost overestimation, with direct cost comparisons typically agreeing more closely. All cost model estimates are well within an Association for the Advancement of Cost Engineering (AACE) Class 5 estimation range (-50% to +100%), with five within the AACE Class 4 range (-30% to +50%) and four being within 15%. This result is considered reasonable performance for a parametric model applied at a preliminary design stage.

13 HYDRO ENERGY

Advanced Computing Annual Report 2024

In fiscal year (FY) 2024, the National Renewable Energy Laboratory (NREL) took a major leap forward with the completed full buildout of Kestrel, the Office of Energy Efficiency and Renewable Energy's newest high-performance computing (HPC) system. Kestrel is already supporting science across the portfolio, bringing roughly 44 petaflops of computing power, which is more than five times the capacity of our previous supercomputer, Eagle. By delivering greater GPU capacity, Kestrel enables faster progress in artificial intelligence (AI) and opens new avenues in energy research - from defining long-term planning scenarios to accommodate a growing power system to material discovery to improving energy efficiency in photovoltaics (PV). Across the portfolio, research is being accelerated by Kestrel's impressive power. During FY24, 427 projects and more than 700 researchers used NREL's HPC, supporting the U.S. Department of Energy's Office of Energy Efficiency and Renewable Energy across 13 funding areas. Through these collaborations, researchers produced more than 450 technical outputs, including 195 articles in peer-reviewed publications, pushing the boundaries of science and engineering. This year's report features new sections spotlighting the expanding roles of Artificial Intelligence and Accelerated Computing. We also introduce an early career section to celebrate the accomplishments of our up-and-coming researchers, whose pioneering work is shaping the future of energy. We hope you enjoy the new insights and discoveries highlighted in these pages.

97 MATHEMATICS AND COMPUTING