Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large scale power system”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Safe Deep Reinforcement Learning for Active Distribution System Model Predictive Control with EVs and DERs

The temporal and spatial mismatch between PV generation and electric vehicle (EV) charging and discharging may cause voltage violations in active distribution networks. Despite the widespread use of deep reinforcement learning (DRL) in power system optimization and control, it lacks guarantees on constraint satisfaction during both training and deployment. This paper proposes a Lagrangian-based safe DRL approach for model predictive control (MPC) of active distribution systems with large-scale integration of PVs, EVs, and energy storage systems (ESSs). A Transformer-LSTM time-series model is proposed to forecast EV charging demand, which is then formulated as a constraint to ensure charging requirements are met. Using this prediction, a Lagrangian-based safe soft actor-critic (SAC) framework is developed for real-time control in a three-phase unbalanced distribution system, enforcing voltage safety constraints while optimizing the cumulative net reward. By integrating the forecasting model with multi-period constraints, the proposed framework jointly coordinates PV systems, EV charging and discharging, and ESS scheduling within the MPC horizon. Numerical experiments on a modified IEEE 123-bus system with real-world data show that, under a high PV penetration scenario, the proposed method increases the net reward by 30.74% and reduces average voltage violations from 0.0011 p.u. to 0.0002 p.u. compared with standard SAC. Compared with the optimal power flow (OPF) approach, it achieves similar voltage security while yielding lower line losses. It also maintains real-time control capability, reducing operation latency to 53.21 ms per 15-minute control interval. The proposed method remains effective under varying PV/EV penetrations and load conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

2020 Hydrogen and Fuel Cell Technologies Market Report

This report features summaries of relevant hydrogen and fuel cell activity in 2020, consistent with prior annual reports across applications and sectors. However, 2020 was no ordinary year. Despite pandemic hardships, the hydrogen and fuel cell industry not only stepped up to assist with humanitarian and healthcare efforts but saw continued momentum in several growing market sectors with numerous partnerships, initiatives, and investments. Using hydrogen and fuel cell technologies to provide an environmental and economic pathway to clean, reliable power was a key feature of world economic recovery and energy infrastructure plans. Today, hydrogen energy and fuel cell systems are being developed and deployed in a range of stationary power, portable power, and transportation applications around the world. Throughout 2020, there were significant achievements in each sector, building on years of research and development (R&D), field testing, large-scale demonstration projects, and real-world customer experience. The heavy-duty truck market in particular continued to make substantial gains, with new players entering the space and a coalition of companies pledging to produce and deploy more than 100,000 fuel cell trucks in the next two decades. To provide hydrogen for trucks and other end users, larger and larger electrolyzer projects were announced throughout the year. Collectively, there were more than 430 megawatts (MW) of electrolyzer systems either deployed, ordered, or proposed around the world. This 2020 report also provides special “Spotlights” on the collective industry response to the COVID-19 pandemic as well as international commitments and hydrogen strategies, providing overviews of significant activity during and after the initial onset of COVID-19 around the world.

08 HYDROGEN↗
DOE OSTI2024

Wet chemistry-synthesized Fe/mixed ferrite soft magnetic composites for high-frequency power conversion

State-of-the-art soft magnetic (SM) alloy systems such as electrical steels and permalloys exhibit large eddy current losses at high-frequency (kHz range and above), limiting their application in fast switching devices. Eddy current losses can be reduced for metallic alloys through a decrease in characteristic length scale or embedding them in an insulation layer. This work proposes the development of fine-scale core materials consisting of metallic SM nano/microparticles with magnetic, inorganic insulation layers, and their synthesis using a wet chemistry-based, scalable method for high frequency and high-power applications. More specifically, a magnetic ferrite coating (Ni 0.5-x Mn x Zn 0.5 Fe 2 O 4 , x = 0.0–0.5) was applied via deposition of ferrite powder produced using a “wet chemistry-based” co-precipitation method; the ferrite was distributed through the Fe micropowder via either ball milling or ultrasonic mixing. Powder cores were prepared by compaction of the synthesized composites. Relative magnetic permeability and core loss were measured at excitation frequencies to 200 kHz. A core loss of 127 kW/m 3 was measured at excitation frequency of 100 kHz and magnetization of 0.02 T. This value improves on the 199.3 kW/m 3 reported in literature at identical excitation conditions for a compact formed from a composite comprising Fe microparticles coated with low-permeability Ni 0.5 Zn 0.5 Fe 2 O 4 . Ultrasonic mixing resulted in slightly lower core loss than ball milling, possibly because ball milling causes loss-increasing deformation of the Fe micropowder. XRD and SEM were used to observe composite composition and core cross-section microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Flexibility Options: A Proposed Product for Managing Imbalance Risk

The presence of variable renewable energy resources with uncertain outputs in day-ahead electricity markets results in additional balancing needs in real-time. Addressing those needs cost-effectively and reliably within a competitive market with unbundled products is challenging as both the demand for and the availability of flexibility depends on day-ahead energy schedules. Existing approaches for reserve procurement usually rely either on oversimplified demand curves that do not consider how system conditions that particular day affect the value of flexibility, or on bilateral trading of hedging instruments that are not co-optimized with day-ahead schedules. This article proposes a new product, ‘Flexibility Options', to address these two limitations. The demand for this product is endogenously determined in the day-ahead market and it is met cost-effectively by considering real-time supply curves for product providers, which are co-optimized with the energy supply. As we illustrate with numerical examples and mathematical analysis, the product addresses the hedging needs of participants with imbalances cost-effectively, provides a less intermittent revenue stream for participants with flexible outputs, promotes value-driven pricing of flexibility, and ensures that the system operator is revenue-neutral. This article provides a comprehensive design that can be further tested and applied in large-scale systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

To What Extent Will Decarbonization Deepen the Conversation Between Industry and the Grid?

Decarbonization - the transition away from un-mitigated fossil fuel combustion throughout the economy - requires big changes from both power and process systems. On the power system side, those changes are expected to include large increases in variable generation, e.g., from wind and solar, which has near-zero marginal costs and at large shares can produce infrequent but consequential energy droughts. On the process systems side, industries are investigating their options for direct and indirect electrification, the latter exemplified by replacing fossil fuel inputs with zero-carbon, energy-carrying chemicals like hydrogen and ammonia produced via electrochemical processes. The economic features of these changes within the larger context of power and process systems suggest that their realization could be accompanied by a paradigm shift in how industrial facilities interact with the grid. For example, the dominant type of demand participation in power markets could change from today's focus on load reductions at peak times to a new focus on shifting electricity use, enabled in part by large-scale product storage, to take advantage of renewable energy that would otherwise be curtailed and to avoid consumption during high-price energy droughts. This talk will describe these and other possible design and operational approaches from grid and industrial economic perspectives, culminating in an enumeration of open problems that lie at the interface of today and tomorrow's power and process systems.

co-design↗
DOE OSTI2025

Coherency-Constrained Spectral Clustering for Power Network Reduction

This paper presents a methodology for reducing the complexity of large-scale power network models using spectral clustering, aggregation of electrical components, and cost function approximation. Two approaches are explored using unconstrained and constrained spectral clustering to determine areas for effective system reduction. Once the system areas are determined, both loads and generators by type are aggregated, and their new cost function is approximated through polynomial curve-fitting or statistical methods. The performance of reduced networks is evaluated in terms of their ability to follow the true daily cost of the original system over a 24-hour period considering a set of several days. Two test systems are taken as test beds. Application of the methodology to a modified version of the IEEE 39-bus system reduces it from 17 generators to a 4-bus system and 9 generators with about 93% of accuracy. Similarly, the IEEE 118-bus system is reduced from 19 generators to a 3-bus system with three aggregated units achieving over 99% of accuracy. These findings address scalability challenges and enhance accuracy for high and mid-loading level conditions, and by aggregating thermal units with similar cost functions.

42 ENGINEERING↗
DOE OSTI2022

High Fidelity CFD Simulations Supporting the KP-FHR

Kairos Power, LLC, is developing its version of the Fluoride-cooled High-temperature Reactor, the KP-FHR. The design uses a pebble bed core with fluoride salt as a coolant. The pebbles used in the KP-FHR have a diameter of 4 cm, with a shell fuel region where TRISO particles are embedded. A Pebble bed core design is adopted by several Gen IV reactors, They boast many benefits, such as fuel integrity, highly efficient heat transfer, and passive safety. However, it is challenging to accurately predict temperature and flow inside a pebble bed. Traditional approaches use the porous media model, which regards the pebble bed as a continuous medium, but with different temperature fields representing different levels, such as the fluid temperature, pebble surface temperature, and pebble center temperature. Empirical heat transfer correlations are adopted to calculate the heat transfer coefficient between different phases. However, empirical correlations are usually validated with experimental data, which usually lacks detail inside the pebble bed. The available experimental data is also generally at a high Reynolds number, which falls outside of the conditions of KP-FHR. Explicit computational fluid dynamics (CFD) simulations of randomly packed pebble beds have only become feasible recently. This is thanks to the rapid development of computational power and scalable algorithms. In this work, we used the Spectral Element Method (SEM) CFD code NekRS to simulate the randomly packed pebble bed in a cylindrical container. NekRS, which is the GPU variant of Nek5000, but refactored to utilize the computational power of GPUs using the OCCA library to run on hybrid architecture high performance computing systems. It was initially developed with the libParamunal library, but truncated and tuned for large-scale turbulence simulation. As a result, the SEM reaches higher precision with the same degrees of freedom by using a high-order Lagrange polynomial basis distributed on Gauss-Lobatto-Legendre quadrature inside each element, compared to lower-order methods, such the Finite Volume Method and Finite Element Method. The report is divided into five parts. We start with a general discussion of the pebble bed reactor, along with a specific investigation into the KP-FHR. The second part presents the numerical methodology. In the third part, we study a modular pebble bed with 1741 pebbles in a container of 7 pebble-diameter radius. Beyond LES simulations done by NekRS, we also leveraged the thermal radiation model in OpenFOAM to study heat transfer under no-forced-flow scenarios. Then, in the fourth part we simulated a pebble bed similar to the size of the Hermes Test Reactor. The total number of pebbles is in these simulations is 34,374. The container radius is 14 pebble-diameters. Finally, the report concludes in part five, with a discussion of future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Resilient buildings for fire-adapted landscapes: EE and flexible loads integrated with solar and storage microgrids

Energy efficiency (EE) and flexible loads can be part of a resilient buildings packagewhen they are combined with solar and storage in a clean energy microgrid to reduce the carbonfootprint of buildings and enable resilience to extreme events. Recent large wildfires and an emerging understanding of the fire-adapted landscapes in the American West underscore the urgency of work towards scaling and commercializing these systems. For example, a recent power shutoff in Northern California (October 9-12, 2019) resulted in 738k customers disconnected at the peak of the outage and is emblematic of challenges to come. Our paper reports on insights from a clean energy microgrid deployment pilot that integrates a 50 kWAC PV array, a 109 kW / 174 kWh battery system, switchgear to safely isolate from the regional power system, and communicating controllers for HVAC and refrigeration. The project will be commissioned in May 2020 and is sited at a critical infrastructure site in rural Northern California – in this case a gasoline station with convenience store. Our experience and results shed light on capabilities of integrated microgrids to provide value to customers during blue sky conditions and resilience during black sky days with high fire risk, and what opportunities and barriers exist for scaling these integrated microgrid systems in the near term. We use a simulation model to estimate how EE and flexibility can extend the run time of solar and storage, improving the reliability of power at critical sites.

Alstone, Peter↗
DOE OSTI2021

Distributional Deep Reinforcement Learning-Based Emergency Frequency Control

Emergency frequency control is one of the most critical approaches to maintain power system stability after major disturbances. With the increasing number of grid-connected renewable energy sources, existing model-based methods of frequency control are facing up with challenges of computational speed and scalability for large-scale systems. In this paper, the emergency frequency control problem is formulated as a Markov Decision Process (MDP) and solved through a novel Distributional Deep Reinforcement Learning (DDRL) method, namely the distributional soft actor critic (DSAC) method. Compared with other RL methods that only estimate the mean value, the proposed DSAC model estimates the distribution of value function over returns. This advancement can lead to more insights and knowledge for the agent, with the benefit of a much faster and more stable learning process, and the improved frequency control performance. Here, the simulation results on IEEE 39-bus and IEEE 118-bus systems demonstrate the effectiveness and robustness of proposed models, as well as the advantage compared to other state-of-the-art DRL algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

How do North American weather regimes drive wind energy at the sub-seasonal to seasonal timescales?

Abstract There has been an increasing need for forecasting power generation at the subseasonal to seasonal (S2S) timescales to support the operation, management, and planning of the wind-energy system. At the S2S timescales, atmospheric variability is largely related to recurrent and persistent weather patterns, referred to as weather regimes (WRs). In this study, we identify four WRs that influence wind resources over North America using a universal two-stage procedure approach. These WRs are responsible for large-scale wind and power production anomalies over the CONUS at the S2S timescales. The WR-based reconstruction explains up to 40% of the monthly variance of power production over the western United States, and the explanatory power of WRs generally increases with the increase of timescales. The identified relationship between WRs and power production reveals the potential and limitations of the regional WR-based wind resource assessment over different regions of the CONUS across multiple timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Scalable learning of potentials to predict time-dependent Hartree–Fock dynamics

We propose a framework to learn the time-dependent Hartree–Fock (TDHF) inter-electronic potential of a molecule from its electron density dynamics. Although the entire TDHF Hamiltonian, including the inter-electronic potential, can be computed from first principles, we use this problem as a testbed to develop strategies that can be applied to learn a priori unknown terms that arise in other methods/approaches to quantum dynamics, e.g., emerging problems such as learning exchange–correlation potentials for time-dependent density functional theory. We develop, train, and test three models of the TDHF inter-electronic potential, each parameterized by a four-index tensor of size up to 60 × 60 × 60 × 60. Two of the models preserve Hermitian symmetry, while one model preserves an eight-fold permutation symmetry that implies Hermitian symmetry. Across seven different molecular systems, we find that accounting for the deeper eight-fold symmetry leads to the best-performing model across three metrics: training efficiency, test set predictive power, and direct comparison of true and learned inter-electronic potentials. All three models, when trained on ensembles of field-free trajectories, generate accurate electron dynamics predictions even in a field-on regime that lies outside the training set. To enable our models to scale to large molecular systems, we derive expressions for Jacobian-vector products that enable iterative, matrix-free training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Dynamic Validation of CNN-Based Surrogate Models for Inverter-Based Resources in Open-Source Solvers

Traditionally, distribution system planning has focused on steady-state analyses, with limited consideration of dynamic behavior. However, as large or medium-scale inverter-based resources (IBRs), particularly grid-following (GFL) inverters in commercial or industry buildings, become more prevalent, understanding their dynamic impact is essential for grid planning and operation. This article presents an innovative deep-learning (DL)-approach using convolutional neural networks technique to model the GFL inverters. Developed from real grid-tied commercial IBR transient data, these dynamic DL models overcome proprietary constraints by requiring minimal knowledge of internal converter physics while maintaining high accuracy and flexibility. To demonstrate their applicability, the models were incorporated into GridLAB-D, an open-source, three-phase distribution analysis tool. This integration enables dynamic simulations of large-scale distribution networks with high IBR penetration stability analysis. Rigorous testing and validation, aligned with industry standards, confirmed the reliability and efficiency of this approach, paving the way for enhanced planning and operational assessments of modern power systems.

Deep-learning↗
DOE OSTI2021

Markets and Economic Requirements for Fission Batteries and Other Nuclear Systems

Fission Batteries (FBs) are nuclear reactors defined by five characteristics which enable large-scale deployment: cost competitive, standardized sizes for economic mass production, easy installation and removal, secure and safe unattended operation with high reliability. FBs are not defined by technology or power level. Technical and market considerations suggest that most FBs will produce 20 to 30 MWt. This proceedings reports on the outcomes of two workshops that were held in January 2021 to better define markets and economic challenges for FBs. Three major markets were identified. The largest market is the industrial and commercial heat market. There are about 4000 industrial users (excluding utilities) that require more than one megawatt of heat. The number of customers versus size of heat demand was determined. In a low-carbon world there is the potential for many additional customers—including expanded biofuels production and district heat. The second market is for non-grid electricity. This includes co-generation plants that produce heat and electricity for a single customer. The third market is the maritime market with ~100,000 ships worldwide. In the United States, natural gas is the low-cost energy option today and will remain so unless constraints or taxes impact its use. If restrictions on greenhouse gas emissions, the FB competition includes natural gas with carbon capture, biofuels, hydrogen and grid electricity. Natural gas with carbon capture is not economically viable on a small scale. Biofuels may be expensive but may be the economically preferred option for locations with small energy demands of a few megawatts. Hydrogen is a potential competitor with many of the characteristics of natural gas. Grid electricity is not a competitive source of heat. For FBs to be economically competitive, the price of delivered heat must be $20-50/MWh ($6-15/million BTU). The economically competitive range for non-grid electricity is estimated at $70-100/MWh. These electricity prices are competitive with the retail prices of electricity in many parts of the United States for the customer. FBs are not expected to be competitive selling wholesale electricity to the grid. To achieve the aforementioned cost targets for heat and electricity markets, FB designers must (1) maximize the power output within the constraints of a FB (e.g., truck transportability, passive decay heat removal), (2) drastically reduce the size of onsite staff, (3) adopt core designs with low fuel costs (enrichment and fabrication), and (4) develop a system design that is efficiently manufactured in factories. The business case depends upon more than being just a replacement for natural gas. The largest incentives for adoption of FBs is where they create new markets and new sources of revenue. An example is the paper and pulp industry that burns biomass wastes to provide heat and electricity to make paper. An external heat source could meet the demand for heat and electricity by the paper process and enable converting waste biomass into liquid biofuels rather than burning to provide heat. Other markets, such as data centers, are driven by special energy requirements such as extreme reliability. Most customers are not in the energy business but need heat and electricity to produce a product—a manufactured good, education, retail sales (shopping malls), marine transport or some other product. As a consequence, there will be large incentives to lease rather than own FBs. Leasing avoids the regulatory challenges that remain with the owner of the FB. Leasing creates large incentives for FP standardization of sizes and transportability to maintain the value of the FB at the end of the lease—similar to the leasing of jet engines and aircraft. The economic constraints combined with technical constraints suggest competitive FBs will likely have outputs exceeding 10 MWt. There appear to be little incentives for very long-lived reactor cores because such machines require much larger inventories of fuel. Maintenance requirements and the options to provide technology updates may favor shorter lifetimes (~5 years). The assessment is that there is the potential for FBs to be economically viable and play a major role in global decarbonization in three markets: heat, non-grid electricity and maritime applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Technical Impacts of Light-Duty and Heavy-Duty Transportation Electrification on a Coordinated Transmission and Distribution System

In this study, we propose a strategy to model the required spatiotemporal charging demand from light-duty (LD) and medium- and heavy-duty (MHD) electric vehicles (EVs) using actual transportation data by mapping the demand for the required EV charging to a realistic and coordinated distribution and transmission electric grid at the predicted times of the day to study their impact on the power system in a variety of load, weather, and EV penetration scenarios. This work is the first study that includes the actual weather data and transportation data with realistic and coordinated distribution and transmission grid data in a large industry-scale level study. The main goal of this study is to identify possible issues and required upgrades in the electric grid, caused by an increase in EV integration. The transmission case study is a large grid with 6717 buses over a Texas footprint, and the distribution grid is over Houston, a city in Texas, covering over three million customers. The resulting overloads and voltage violations experienced in the system are discussed, and required planning upgrades to avoid these issues are suggested.

AC optimal power flow (AC-OPF)↗
DOE OSTI2022

Peer-to-Peer Energy Trading under Network Constraints Based on Generalized Fast Dual Ascent

We report the wide deployment of renewable energy resources, combined with a more proactive demand-side management, is inducing a new paradigm in both power system operation and electricity market trading, which especially boosts the emergence of the peer-to-peer (P2P) market. A more flexible local market mechanism is highly desirable in response to fast changes in renewable power generation at the distribution network level. Moreover, large-scale implementation of P2P energy trading inevitably affects the secure and economic operation of the distribution network. This paper presents a new P2P electricity trading framework with distribution network security constraints considered using the generalized fast dual ascent method. First, an event-driven local P2P market framework is presented to facilitate short-term or immediate local energy transactions. Then, the sensitivity analysis of nodal voltage and network loss with respect to nodal power injections is used to evaluate the impacts of P2P transactions on the distribution network, which ensures the secure operation of the distribution system. Thereby, the external operational constraints are internalized, and the cost of P2P energy trading can be appropriately allocated in an endogenous way. Moreover, a generalized fast dual ascent method is employed to implement distributed market-clearing efficiently. Finally, numerical results indicate that the proposed model could guarantee secure operation of the distribution system with P2P energy trading, and the solution method enjoys good convergence performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Large-Signal Stability Improvement of Parallel Grid-Forming Inverter-Driven Black Start: Preprint

With the rapid increase of inverter-based resources in modern grids, advanced grid-forming (GFM) inverter capa- bilities, such as system restoration and operation under faults, are urgently needed to realize power electronics-dominant grids at scale. One such capability is inverter-driven black start using GFM inverters. This paper analyzes the ability of two recently proposed advanced GFM controls to help GFM inverters sustain system-wide, off-nominal conditions and remain synchronized until they can overcome the momentary overloading as more GFMs join the process without generator sequence coordination or communications and finally stabilize the grid. Through an extensive set of 1,200 full-order electromagnetic transient simu- lations, we evaluate the black-start process while employing the various GFM inverter controls. The results show that the GFM current limiter and primary control have a significant impact on the stability of the system during dynamic operating conditions and thus impact the success of inverter-driven system restoration.

current limit↗
DOE OSTI2022

Deep Learning Based Superconducting Radio-Frequency Cavity Fault Classification at Jefferson Laboratory

This work investigates the efficacy of deep learning (DL) for classifying C100 superconducting radio-frequency (SRF) cavity faults in the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab. CEBAF is a large, high-power continuous wave recirculating linac that utilizes 418 SRF cavities to accelerate electrons up to 12 GeV. Recent upgrades to CEBAF include installation of 11 new cryomodules (88 cavities) equipped with a low-level RF system that records RF time-series data from each cavity at the onset of an RF failure. Typically, subject matter experts (SME) analyze this data to determine the fault type and identify the cavity of origin. This information is subsequently utilized to identify failure trends and to implement corrective measures on the offending cavity. Manual inspection of large-scale, time-series data, generated by frequent system failures is tedious and time consuming, and thereby motivates the use of machine learning (ML) to automate the task. This study extends work on a previously developed system based on traditional ML methods (Tennant and Carpenter and Powers and Shabalina Solopova and Vidyaratne and Iftekharuddin, Phys. Rev. Accel. Beams, 2020, 23, 114601), and investigates the effectiveness of deep learning approaches. The transition to a DL model is driven by the goal of developing a system with sufficiently fast inference that it could be used to predict a fault event and take actionable information before the onset (on the order of a few hundred milliseconds). Because features are learned, rather than explicitly computed, DL offers a potential advantage over traditional ML. Specifically, two seminal DL architecture types are explored: deep recurrent neural networks (RNN) and deep convolutional neural networks (CNN). We provide a detailed analysis on the performance of individual models using an RF waveform dataset built from past operational runs of CEBAF. In particular, the performance of RNN models incorporating long short-term memory (LSTM) are analyzed along with the CNN performance. Furthermore, comparing these DL models with a state-of-the-art fault ML model shows that DL architectures obtain similar performance for cavity identification, do not perform quite as well for fault classification, but provide an advantage in inference speed.

97 MATHEMATICS AND COMPUTING↗