Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Building performance simulation”

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 235 records · Page 13

Retrofitting Buildings with Solar-Reflective Roofs and Walls and its Impact on Peak Power Demand

Buildings are a major consumer of electricity in the United States and a significant portion of the consumption comes from heating, ventilation, and air-conditioning (HVAC) applications. Passive cooling strategies in the building envelope help to reduce the energy consumption for HVAC as well as peak electricity demand. Although being one of the most cost-effective passive cooling strategies, modern reflective roofing and reflective exterior wall technology is not well documented for its impact on peak demand. This study utilized whole building energy simulations on residential and commercial building prototype models to quantify the impact of cool roofs and cool exterior walls. The analysis was performed in three climate zones with varying insulation levels and solar reflectances for roofs and exterior walls. For both the residential and commercial buildings, the baseline building had a roof solar reflectance value of 0.10 and an exterior wall solar reflectance value of 0.25. The results from the simulations show that roofs and exterior walls with higher reflectance values increase cooling energy savings but can also increase heating energy consumption. The impact of changes in solar reflectances was greater in buildings with low roof/wall insulation levels compared to roofs/walls with higher insulation levels. A baseline for the simulations was set with the roof and exterior wall solar reflectances set at 0.1 and 0.25, respectively and simulations having varying roof and exterior wall thermal resistances were compared to the baselines.

14 SOLAR ENERGY↗

Zero Export Feeder Through Transactive Markets

This project is aimed at creating a transactive energy market to address the challenges faced by utility providers when increasing distributed energy resource (DER) adoption in their service area. One major challenge is mitigating export back to the grid during times of excess production. The transactive energy market operates at the distribution level and balances the supply and demand on the feeder, thus maintaining a zero-energy export at the primary feeder head. The market participants in this case are the residential customers on the feeder, who bid into the market. Building controls are then optimized based on the settled price. Market performance was demonstrated in this study by simulating different levels of DER penetration on a selected Pepco feeder. The feeder successfully achieved a zero export while providing cost-effective electricity to the participants, demonstrating that this market design can enable high DER penetration on existing feeders.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Dynamic Decarbonization through Autonomous Physics-Centric Deep Learning and Optimization of Building Operations (Abstract only)

This project directly addresses the primary goal of Area of Interest 2 in the CRADA call: to advance optimization-based integrated energy management systems in commercial and residential buildings. Pacific Northwest National Laboratory (PNNL) and its industry partner PassiveLogic aim to accomplish this by reaching three key objectives. First, to ensure a broad impact in the building controls industry, PNNL will extend its open-source library for predictive control synthesis by augmenting its capabilities with data-driven self-learning of building models and auto-calibration of predictive controllers. The effort will focus on building use cases selected in collaboration with PassiveLogic. The team will specifically address the development of methods for data-driven adaptation of building models, investigation of model architectures that best address specific building types, and automated synthesis of differentiable predictive controllers that optimize diverse objectives. Second, PNNL will collaborate with PassiveLogic to integrate the aforementioned methods with PasiveLogic’s advanced controls platform. The collaborative integration effort will inform the developments under the first objective by providing specific data on the attainable performance of model learning on resource-constrained edge computing platforms. This software integration effort will increase the technical maturity of the developed libraries by exploring the use of software integration tools and methods. Third, PNNL and PassiveLogic will work to improve the technology readiness of the developed predictive controllers by testing their performance in relevant test environments, such as high-fidelity simulation, hardware in the loop, and actual test buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Tank 48H Tetraphenylborate Mitigation: Simulant Studies using Sodium Permanganate

Tank 48H currently holds legacy material containing organic tetraphenylborate (TPB) compounds from the operation of the In-Tank Precipitation (ITP) process. TPB was added during the ITP process to precipitate the otherwise soluble cesium as insoluble cesium TPB (CsTPB), but excessive benzene generation from TPB degradation curtailed this treatment method. The contents of Tank 48H, which include approximately 26.000 kg of potassium TPB (KTPB) and trace CsTPB, are not compatible with the waste treatment facilities at the Savannah River Site (SRS) since the organic content and the associated flammability issues pose a challenge to the salt processing and sludge processing facilities within the liquid waste system. An in-tank process to remove (or decompose) TPB safely would be of great value. Previous testing at Savannah River National Laboratory (SRNL) demonstrated the destruction of glycolate via chemical oxidation using sodium permanganate with simulated and radioactive waste. Scoping tests were performed to study the destruction of TPB to determine if the contents of Tank 48H would be amenable to the same type of destruction. Partial destruction of TPB was observed in Tank 48H simulants under mild conditions (e g., pH 11, room temperature) with no definitive indication of benzene generation. To build upon the success of the scoping tests, an additional study was requested to provide a better understanding of the underlying chemistry for Tank 48H content destruction using sodium permanganate. Three experiments were performed with Tank 48H simulants at 40 °C to determine the efficacy of using sodium permanganate for TPB destruction. Three starting pH values were selected: 1) pH 11 for comparison with the previous work at room temperature, 2) pH 10 as the minimum pH recommended by the Corrosion Control Program (CCP) for in-tank processing, and 3) pH 8 to determine the effectiveness of TPB destruction at near neutral pH. While below the allowable pH for the CCP, the experiment at pH 8 was performed to study the TPB-Permanganate reaction under more extreme conditions and further verify the potential for out-of-tank processing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Two-Stage Reinforcement Learning Policy Search for Grid-Interactive Building Control

This paper develops an intelligent grid-interactive building controller, which optimizes building operation during both normal hours and demand response (DR) events. To avoid costly on-demand computation and to adapt to non-linear building models, the controller utilizes reinforcement learning (RL) and makes real-time decisions based on a near-optimal control policy. Learning such a policy typically amounts to solving a hard non-convex optimization problem. We propose to address this problem with a novel global-local policy search method. In the first stage, an RL algorithm based on zero-order gradient estimation is leveraged to search for the optimal policy globally, due to its scalability and the potential to escape some poor performing local optima. The obtained policy is then fine-tuned locally to bring the first-stage solution closer to that of the original unsmoothed problem. Experiments on a simulated five-zone commercial building demonstrate the advantages of the proposed method over existing learning approaches. They also show that the learned control policy outperforms a pragmatic linear model predictive controller (MPC) and approaches the performance of an oracle MPC in testing scenarios. Using a state-of-the-art advanced computing system, we demonstrate that the controller can be learned and deployed within hours of training.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Experimental characterization and potential energy savings of insulated cladding for U.S. residential buildings

This study evaluates a novel insulated cladding composed of foamed cement and vinyl siding that is designed as a retrofit product to enhance thermal performance in U.S. residential buildings. Experimental measurements showed that foamed cement exhibits a thermal conductivity of 0.0334–0.0365 W/m-K (R-value of 4.0–4.3 °F⋅ft 2 ⋅h/BTU), which is significantly higher than that of conventional fiber cement boards. Although mechanical testing confirmed the material’s brittleness and low strength, it also indicated its suitability for nonstructural insulation applications. Water vapor permeability testing demonstrated effective moisture resistance with the polymer coating. To evaluate the insulated cladding’s energy performance, we used the U.S. Department of Energy’s prototype single-family building model and simulated across 50 locations spanning 16 International Energy Conservation Code climate zones. The R-7 (R-value of 7 °F⋅ft 2 ⋅h/BTU per inch) retrofit cladding resulted in heating energy savings up to 275 therms (Fairbanks, AK) and cooling energy savings up to 1,547 kWh (Phoenix, AZ). The energy cost savings varied by region, but the highest annual cost savings ($\$$434/year) were observed in Santa Maria, CA, and the highest annual cost savings percentage (37 %) were noted in Monterey, CA. On a national scale, the insulated retrofit provided an average annual energy cost reduction of 16 %, highlighting the cladding system’s broad applicability and financial viability.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Benefits assessment of cool skin and ventilated cavity skin: Saving energy and mitigating heat and grid stress

This study assessed the energy-saving and climate-adaptive potential of cool skin and ventilated cavity skin facade technologies in Seoul's high-rise apartment buildings. We created weather scenarios for historical, mid-term future, and long-term future conditions using Coordinated Regional Downscaling EXperiment (CORDEX) method. Building energy simulations were conducted on a South Korean high-rise apartment model to evaluate their performance under different weather conditions. The results indicate that cool skin and ventilated cavity skin technologies can save cooling energy during summers but lead to heating energy penalties in winters. Ventilated cavity skin outperforms cool skin, offering better cooling energy savings and reduced heating penalties. Combining both technologies yields the highest overall energy savings, with 7 %, 9 %, and 10 % cooling energy savings for cool skin, ventilated cavity skin, and the combined package, respectively. However, cool skin increases heating energy consumption by 5 %, while ventilated cavity skin has minimal impact on heating energy. These envelope technologies also reduce peak electricity demand by at least 5 %, 8 %, and 9 %, respectively. They contribute to heat stress reduction, enhance resilience, and decrease extreme heat risks for occupants during power outages by at least 18 % under various weather conditions. Considering the prevalence of aging high-rise apartments in South Korea, adopting these envelope renovation strategies can effectively reduce cooling loads, enhance thermal comfort, and boost resilience under future climates, while avoiding costly reconstruction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Deep reinforcement learning control for co-optimizing energy consumption, thermal comfort, and indoor air quality in an office building

With the recent demand for decarbonization and energy efficiency, advanced HVAC control using Deep Reinforcement Learning (DRL) becomes a promising solution. Due to its flexible structures, DRL has been successful in energy reduction for many HVAC systems. However, only a few researches applied DRL agents to manage the entire central HVAC system and control multiple components in both the water loop and the air loop, owing to its complex system structures. Moreover, those researches have not extended their applications by incorporating the indoor air quality, especially both CO2 and PM2.5concentrations, on top of energy saving and thermal comfort, as achieving those objectives simultaneously can cause multiple control conflicts. What's more, DRL agents are usually trained on the simulation environment before deployment, so another challenge is to develop an accurate but relatively simple simulator. Therefore, we propose a DRL algorithm for a central HVAC system to co-optimize energy consumption, thermal comfort, indoor CO2 level, and indoor PM2.5 level in an office building. To train the controller, we also developed a hybrid simulator that decoupled the complex system into multiple simulation models, which are calibrated separately using laboratory test data. The hybrid simulator combined the dynamics of the HVAC system, the building envelope, as well as moisture, CO2, and particulate matter transfer. Three control algorithms (rule-based, MPC, and DRL) are developed, and their performances are evaluated on the hybrid simulator environment with a realistic scenario (i.e., with stochastic noises). The test results showed that, the DRL controller can save 21.4 % of energy compared to a rule-based controller, and has improved thermal comfort, reduced indoor CO2 concentration. The MPC controller showed an 18.6 % energy saving compared to the DRL controller, mainly due to savings from comfort and indoor air quality boundary violations caused by unmeasured disturbances, and it also highlights computational challenges in real-time control due to non-linear optimization. Finally, we provide the practical considerations for designing and implementing the DRL and MPC controllers based on their respective pros and cons.

Guo, Fangzhou↗

Floating Block Method for Quantum Monte Carlo Simulations

Quantum Monte Carlo simulations are powerful and versatile tools for the quantum many-body problem. In addition to the usual calculations of energies and eigenstate observables, quantum Monte Carlo simulations can in principle be used to build fast and accurate many-body emulators using eigenvector continuation or design time-dependent Hamiltonians for adiabatic quantum computing. Furthermore, these new applications require something that is missing from the published literature, an efficient quantum Monte Carlo scheme for computing the inner product of ground state eigenvectors corresponding to different Hamiltonians. In this work, we introduce an algorithm called the floating block method, which solves the problem by performing Euclidean time evolution with two different Hamiltonians and interleaving the corresponding time blocks. We use the floating block method and nuclear lattice simulations to build eigenvector continuation emulators for energies of 4 He, 8 Be, 12 C, and 16 O nuclei over a range of local and nonlocal interaction couplings. From the emulator data, we identify the quantum phase transition line from a Bose gas of alpha particles to a nuclear liquid.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Designing reinforcement learning algorithms for building HVAC control: From experimental observation to simulation comparisons

Advanced supervisory-level control with reinforcement learning (RL) is regarded as a promising solution for HVAC systems to minimize energy consumption while maintaining thermal comfort and indoor air quality. However, most RL applications were conducted in the simulation environment rather than real-world HVAC systems. This paper developed a value-based RL controller termed Deep Q-Network (DQN) for a typical central HVAC system and evaluated its performance in a building test facility. By comparing DQN with a rule-based controller, the study not only demonstrated the cases where DQN could properly maintain indoor comfort but also discussed possible reasons why DQN failed in some other situations. Recognizing the limitations of value-based RL algorithms from the experimental tests, a simulation study was conducted to compare DQN with an alternative RL approach, an actor–critic algorithm termed Deep Deterministic Policy Gradient (DDPG). In scenarios with a relatively large action space, DDPG outperformed DQN by requiring fewer computational resources and achieving better thermal comfort, lower energy consumption, and more stable control actions. The findings suggest that the ability of DDPG to handle continuous control variables more effectively allows for faster convergence in training and more precise control in practice, which enhances the overall efficiency and reliability of the HVAC system.

Guo, Fangzhou↗

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗

Adaptive Control for Residential HVAC Systems to Support Grid Services

In this paper, we propose an adaptive control strategy for residential buildings' heating ventilation and air conditioning (HVAC) systems to support grid services. The proposed control strategy captures the time-varying characteristics of the building and weather environments. It assigns to each building a local controller to maintain the indoor temperature within a desired comfort level, while the load aggregator allocates the assigned reference power provided by the distribution system operator to support a specific grid service, such as demand peak reduction, balancing supply and demand, and consuming the solar photovoltaic (PV) power locally. Simulation results are presented to verify the performance of the proposed control strategy in terms of maintaining the indoor temperatures of several buildings within the desired comfort bands while tracking the reference power profile. Also, a comparison with a recently developed model-free control strategy is presented.

Wu, Tumin↗

Influence of Phase-Change Materials and Operating Periods on Load Shifting for Thermal Storage Heat Pumps in European Climate Zones

As the reliance on heat pumps (HPs) for space conditioning increases worldwide, ex-tensive grid strains during utility peak hours are expected. Significant research efforts focus on investigating the use of integrated thermal energy storage (TES) systems to reduce peak energy demands. Conventionally, TES systems are often integrated to store and provide either cooling or heating, but not both. In this work, a single room-temperature phase-change material (PCM)-TES was integrated for dual-mode opera-tion with an R290 air-to-water 10.5 kW (3-Ton) HP. Two PCM melting temperatures (17°C, 22°C) were assessed to investigate their impact on system performance and cost reductions for Milan, Italy (4A); Barcelona, Spain (3A); and Oslo, Norway (6A). HP and HP-TES performance maps were generated using Modelica-based transient models and co-simulated with a prototype residential building using the Spawn of En-ergyPlusTM with simple rule-based controls. Our simulations indicate that PCM-17°C provided annual cost savings at higher peak-to-off-peak cost ratios, as the demand re-duction dominated the savings potential. The recharge energy costs and requirements became more dominant at lower cost ratios, making PCM-22°C more beneficial at lower cost ratios. Finally, despite Oslo's heating-dominated climate, the utility peak hours favor cooling demand reductions, making PCM-17°C more appealing. In all cases, the cost savings potential stabilized between 4 – 10% regardless of cost ratio, suggesting that system-level component and operation control optimization is required to maximize cost savings further.

25 ENERGY STORAGE↗

Cellulose nanofibrils and nanocrystals in confined flow: Single-particle dynamics to collective alignment revealed through scanning small-angle x-ray scattering and numerical simulations

Nanostructured materials made through flow-assisted assembly of proteinaceous or polymeric nanosized fibrillar building blocks are promising contenders for a family of high-performance biocompatible materials in a wide variety of applications. Optimization of these processes relies on improving our knowledge of the physical mechanisms from nano- to macroscale and especially understanding the alignment of elongated nanoparticles in flows. Here, we study the full projected orientation distributions of cellulose nanocrystals (CNCs) and nanofibrils (CNFs) in confined flow using scanning microbeam SAXS. For CNCs, we further compare with a simulated system of dilute Brownian ellipsoids, which agrees well at dilute concentrations. However, increasing CNC concentration to a semidilute regime results in locally arranged domains called tactoids, which aid in aligning the CNC at low shear rates, but limit alignment at higher rates. Similarly, shear alignment of CNF at semidilute conditions is also limited owing to probable bundle or flock formation of the highly entangled nanofibrils. Finally, this work provides a quantitative comparison of full projected orientation distributions of elongated nanoparticles in confined flow and provides an important stepping stone towards predicting and controlling processes to create nanostructured materials on an industrial scale.

36 MATERIALS SCIENCE↗

Comparative Analysis of HEATNETS for Geothermal Network Performance: Preprint

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer, rather than also comparing control systems responses. HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study shows that HEATNETS, while not intended to replace other modeling tools, can be a unique modeling tool for the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗

Comparative Analysis of HEATNETS for Geothermal Network Performance

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification and validation process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study focuses on a validation of HEATNETS using results from TRNSYS. HEATNETS is not intended to replace other modeling tools, but this work demonstrates, via a comparison with an industry standard code, that HEATNETS can be a unique, high-level and rapid modeling tool for estimating the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗

Subdiffusive High-Pressure Hydrogen Gas Dynamics in Elastomers

We report elastomeric rubber materials serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. High-pressure exposure and subsequent rapid decompression often lead to cavitation and stress-induced damage of the elastomer due to localization of the hydrogen gas. Here, we use all-atom classical molecular dynamics simulations to assess the impact of compositional variations on gas diffusion within the commonly used elastomer ethylene–propylene–diene monomer (EPDM). With the aim to build a predictive understanding of precursors to cavitation and to motivate material formulations that are less sensitive to hydrogen-induced failure, we perform systematic simulations of gas dynamics in EPDM as a function of temperature, gas concentration, and cross-link density. Our simulations reveal anomalous, subdiffusive hydrogen motion at pressure and intermediate times. We identify two groups of gas with different mobilities: one group exhibiting high mobility and one group exhibiting low mobility due to their motion being impeded by the polymer. With decreasing temperatures, the low-mobility group shows increased gas localization, the necessary precursor for cavitation damage in these materials. At lower temperatures, increasing cross-link density led to greater hydrogen gas mobility and a lower fraction of caged hydrogen, indicating that increasing cross-link density may reduce precursors to cavitation. Finally, we use a two-state kinetic model to determine the energetics associated with transitions between these two mobility states.

36 MATERIALS SCIENCE↗

Landmark-Warped Emulators for Models with Misaligned Functional Response

Many computer models output functional data, and in some cases, these functional data have similar, but misaligned, shape characteristics. In this paper, we introduce a general approach for building emulators for computer models that output misaligned functional data when key values in the functional response (landmarks) can be easily identified. This approach has two main parts: modeling the aligned (using the landmarks) functional data, and modeling the functions that map the misaligned data to the aligned space (warping functions). As the warping functions are required to be monotonic, we give special attention to modeling monotonic functional response data. We discuss how our approach can be easily applied for a variety of typical emulators, such as Gaussian processes, Bayesian multivariate adaptive regression splines, and Bayesian additive regression trees, and how sensitivity analysis can be performed. We demonstrate our approach by building emulators for two applications: (1) a high-energy-density physics computer model used to simulate inertial confinement fusion ignition experiments, where model outputs are highly misaligned, and (2) a multiphysics continuum hydrocode used to simulate high-velocity impact experiments, where model outputs are only slightly misaligned. In case (1) traditional methods cannot be applied, while in (2) they can be applied, but the proposed method performs significantly better.

97 MATHEMATICS AND COMPUTING↗