Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow analysis”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Development of Methods for In-line Monitoring of Yield Stress During the Transfer of Radioactive Waste - 20435

Radioactive waste slurries at the Savannah River Site (SRS) and Hanford often behave as non-Newtonian fluids with a yield stress. The yield stress increases the size of equipment needed to transport or mix these slurries and can increase the risk of process upsets. This property is often measured by collection of samples and transportation to a laboratory, which can increase the risk of personnel exposure and even change the rheology. In addition, slurry rheology has been observed to change with time. Real-time in-line monitoring (RTIM) of rheological properties of slurries can have substantial benefits to critical processes which are currently in operation or will be part of future plans in Department of Energy (DOE) complex processes. The need for an appropriate RTIM technology that can potentially minimize or eliminate numerous sampling operations was recently identified. This research at Florida International University focuses on addressing the need for a reliable RTIM technology at US Department of Energy sites. The objective of this work is to investigate and develop appropriate and reliable technologies that can perform accurate and easy measurements of yield stress with minimal changes to the flow conditions at US Department of Energy sites. Use of ultrasound in longitudinal and shear modes was considered as a method for monitoring of yield stress variations with no disturbance to the flow. Measurement of shear wave velocity (in shear mode) and monitoring of sound velocity variations with start of flow from rest (longitudinal mode) are the underlying principles. Measurement of pressure loss and liquid rise (due to static differential pressure) during flow of slurries were two other methods that were considered and investigated. In addition, analysis of stability depth of a penetrometer in slurry medium (penetrometer method) revealed a potential suitability. This paper presents efforts associated with investigation of different methods mentioned earlier for kaolin-water mixtures at different concentrations. Theoretical analysis showed that all methods under investigation were capable of monitoring yield stress in a desirable range of 0 to 70 pascal. A bench scale test setup was constructed for investigation of pressure loss and liquid rise methods. Simulants at different concentrations were created for initial testing campaigns. For ultrasonic testing, effects of yield stress on sound propagation speed in longitudinal and shear modes were investigated for simulants placed inside PVC pipes. Parameters such as simulant concentration and pipe diameter were varied to find optimal values. It was found that spacing (or medium thickness) of 1/8 inch or less was necessary to observe echoes with sufficient strength. In addition, a significant effect of yield stress on sound propagation velocity was observed. For the liquid rise method, kaolin-water simulant at 1.17 specific gravity was pumped through a 1/4 inch pipe and liquid height was measured in a vertical riser branching off a three-way connection. Static gauge pressure was varied within [0.3 to 1.1] psi range by increase of pumping pressure and flow. Results showed consistent yield stress data within 10% deviation from values reported by a rheometer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Impact of climate change on water availability and its propagation through the Western U.S. power grid

Climate impact assessment for power systems has tended to focus on the effects of water availability on generation without explicitly considering how these impacts propagate through the interconnected grid to affect generation in other locations. We investigate the impact of projected water availability on the operations of the contemporary Western U.S. power grid, focusing specifically on regional interdependencies. Hydrologic simulations derived from three GCMs (CCSM4, INMCM4 and GFDL-CM3) and two radiative scenarios (RCP4.5 and RCP8.5) are used to force a large-scale, distributed water management model (MOSART-WM), which translates water availability into power generation constraints at hydropower plants and water-dependent thermo-electric plants. Power system dynamics are evaluated using the production cost model PLEXOS. We find that the inter-regional connections across the Western U.S. electricity infrastructure are key in managing drastic variations in regional water availability. Net generation in the Northwest and Northern California is substantially sensitive to local water availability. Generation from the Desert Southwest emerges as a critical resource to compensate for variations in water availability in those regions. Climate change impacts on Northwestern water availability therefore drive future changes in other regions’ generation and regional power flows. The large inter-annual hydrological variability of Northern California exacerbates grid stress arising from long-term changes in water availability in the Northwest. While the regional power flow directions seem insensitive to long-term variations in water availability, the analysis highlights the need to consider other compounding regional factors such as changes in Southern California net load or changes in regional fuel prices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Effect of a liquid-type temperature rise inhibitor on cement hydration

A liquid-type temperature rise inhibitor (L-TRI) based on modified sorbitol was developed to solve the thermal cracking issue by modulating the heat release from cement hydration during the early age. Both the exothermic process and phase evolution of cement pastes blended with/without L-TRI were investigated by the combination of isothermal calorimetry, in-situ XRD, ICP-OES and SEM. Also, the adsorption of the L-TRI in the cement paste was identified by total organic carbon analysis. The testing results suggested that the L-TRI affected the heat flow curve in two ways. It prolonged the induction period as a retarding effect. It also reduced the heat flow during the acceleration period and decreased the maximum heat flow as a depressing effect. The retarding effect was mainly caused by the complexation of unabsorbed molecules in the pore solution, and the depressing effect was caused by inhibiting the growth of C-S-H.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.

Anisotropic diffusion↗
DOE OSTI2024

Implementation of a Drift Flux Model into SAM with Development of a Verification and Validation Test Suite for Modeling of Noncondensable Gas Mixtures

The advanced thermal-hydraulic system code, System Analysis Module (SAM), was originally developed for the modeling of single-phase flow in advanced reactors. It has since been expanded to include a four-equation drift flux model for the modeling of two-phase flows containing a noncondensable gas. The model was expanded to support the modeling of molten salt reactor (MSR) designs in which the fuel is directly dissolved in the circulating coolant. These designs have shown that circulating gas bubbles can play an important role in the management of fission products and the operational behavior of the reactor. A drift flux model was implemented to more accurately capture the localized behavior of the void in the core and its impact on the mass transfer of fission products. A thorough assessment of the new model was performed by developing a verification and validation test suite. Verification problems were designed to test all major terms in the new governing equations. The new model converged to the correct solution at the expected order of accuracy for all verification cases. The validation cases included a wide range of flow and void conditions in different pipe geometries. Although higher void experiments show a slight underprediction of void by the drift flux model, experiments that aim to reproduce Molten Salt Reactor Experiment (MSRE) experimental conditions show good agreement with the model. The gas transport model was activated for a SAM model of the MSRE to demonstrate that it can be used in a more complex model. Finally, this gas transport model will be used along with an interfacial area transport equation being implemented in SAM for the prediction of mass transport behavior in MSR conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Hyphae move matter and microbes to mineral microsites: Integrating the hyphosphere into conceptual models of soil organic matter stabilization

Associations between soil minerals and microbially derived organic matter (often referred to as mineral-associated organic matter or MAOM) form a large pool of slowly cycling carbon (C). Here, the rhizosphere, soil immediately adjacent to roots, is thought to control the spatial extent of MAOM formation because it is the dominant entry point of new C inputs to soil. However, emphasis on the rhizosphere implicitly assumes that microbial redistribution of C into bulk (non-rhizosphere) soils is minimal. We question this assumption, arguing that because of extensive fungal exploration and rapid hyphal turnover, fungal redistribution of soil C from the rhizosphere to bulk soil minerals is common, and encourages MAOM formation. First, we summarize published estimates of fungal hyphal length density and turnover rates and demonstrate that fungal C inputs are high throughout the rhizosphere–bulk soil continuum. Second, because colonization of hyphal surfaces is a common dispersal mechanism for soil bacteria, we argue that hyphal exploration allows for the non-random colonization of mineral surfaces by hyphae-associated taxa. Third, these bacterial communities and their fungal hosts determine the chemical form of organic matter deposited on colonized mineral surfaces. Collectively, our analysis demonstrates that omission of the hyphosphere from conceptual models of soil C flow overlooks key mechanisms for MAOM formation in bulk soils. Moving forward, there is a clear need for spatially explicit, quantitative research characterizing the environmental drivers of hyphal exploration and hyphosphere community composition across systems, as these are important controls over the rate and organic chemistry of C deposited on minerals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Thermal Hydraulic Experimental Test Article - Fiscal Year 2023 (Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is a facility that is used to develop sodium components and instrumentation as well as acquire experimental data for validation of reactor thermal hydraulic and safety analysis codes. The facility simulates nominal conditions as well as protected/unprotected loss of flow accidents in a sodium-cooled fast reactor (SFR). High fidelity distributed temperature profiles of the developed flow field may be acquired with Rayleigh backscatter based optical fiber temperature sensors. The facility was designed in partnership with systems code experts to tailor the experiment to ensure the most relevant and highest quality data for code validation. THETA is comprised of a traditional primary coolant and secondary coolant system. The primary system is submerged in the pool of sodium and consists of a pump, electrically heated core, intermediate heat exchanger, and connected piping and thermal barriers (redan). The secondary system, located outside of the sodium pool, consists of a pump, sodium to air heat exchanger, and connected piping and valves. To date a test matrix has been completed utilizing the primary system of THETA. These tests, along with computational fluid dynamics and systems code models, determined the heat transfer across the core barrel and intermediate heat exchanger outlet was too great to effectively represent scaled thermal hydraulic phenomena of a liquid metal cooled reactor. Therefore, a significant effort was made to remove the primary system from the METL 28” test vessel #4 and clean the residual sodium from the primary system to facilitate upgrades. Thermal insulation was then incorporated in the core barrel and intermediate heat exchanger outlets. The primary system was then replaced, and a series of tests were performed to assess the performance of the thermal insulation. With primary system testing and upgrades complete, the secondary system could then be brought online. The tube side of the shell-and-tube intermediate heat exchanger was installed onto the primary system flange to begin installing the secondary system. The support structure for the secondary system was then erected on the METL mezzanine alongside the THETA primary system to facilitate installation of the secondary system components (sodium-to-air heat exchanger, flowmeter and pump). The piping and expansion tank were welded into the secondary system. Non-destructive examination of the secondary system welds was completed in order to satisfy ASME B31.3 pipe code for class M process fluids. The heating system and insulation were then added to prepare the system to be filled with sodium. The ancillary electrical equipment was installed which included the pump control box, blower VFD, pipe heater control system, etc. The secondary system will be filled, and a test matrix will be completed in early FY2024.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

DOE Final Report on Virginia Tech’s contribution to “Tokamak Disruption Simulation”

This work was performed by Virginia Tech in collaboration with multiple institutions led by Los Alamos National Laboratory as a part of the Tokamak Disruption Simulation SciDAC (Scientific Discovery through Advanced Computing) project supported jointly by the Department of Energy Office of Science and Ad- vanced Scientific Computing Research. This report summarizes Virginia Tech’s contributions to the Sci- DAC project. Virginia Tech researchers (presently University of Washington researchers) focused on the fundamental role of plasma-material interaction on transport, which could then have macroscopic effects on simulations of tokamak disruptions. The plasma sheath, which regulates plasma particle and energy fluxes to the wall, is an essential component in the study of plasma-material interaction (PMI). Understanding sheath theory by accounting for finite sheath thickness and transport in the vicinity of the sheath entrance can sig- nificantly modify typical assumptions that are made in the Bohm speed analysis, where the Bohm speed provides the lower bound of the plasma exit flow speed. Our work provides a modified Bohm speed formu- lation that accounts for the critical role of transport and for applications that are away from the asymptotic limits that are typically assumed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

MRCI Subtask 2.4/2.5: Regional/Subregional Analysis and Risk Assessment Final Technical Summary Report

The objective of the Midwest Regional Carbon Initiative (MRCI) project is to implement a collaborative Regional Initiative (RI) to accelerate the deployment of carbon capture, and storage (CCS) in the Midwest-Northeastern quadrant of the United States. This report is a Technical Summary report describing work performed on Tasks 2.4 (Conducting Regional/Subregional Analysis) and 2.5 (Assessing and Managing Risk) during the MRCI project. In Task 2.4, detailed numerical reservoir simulation models were developed for selected carbon storage (CS) systems identified under Task 2.1 (Battelle, 2021a) in the MRCI study area. The objective of Task 2.4 is to demonstrate a dynamic modeling methodology for evaluating the suitability of the (selected) CS systems in the MRCI region for hosting a commercial-scale storage project. In this study, CO2 injectivity was evaluated for different CS systems with an annual injection rate of 1 million metric tonnes (MMT) of CO2 considered as the minimum requirement for a commercial scale project. The objective of Task 2.5 is to assess key risks associated with storing CO2 in the different CS systems across the MRCI region and to demonstrate a method(s) for assessing these risks that may be used by developers of future CO2 storage projects in the region. The risk analysis was limited to evaluating two types of leakage risks (i.e., wellbore leakage, flow across unfractured caprock) at three modeled sites considered in Task 2.4.

MRCI,Report,Summary,Technical Challenges,dynamic m↗
DOE OSTI2025

Iodine Capture Studies of Copper- and Bismuth-Based Sorbents

The release of radioiodine, one of several radionuclides of concern when recycling used nuclear fuel (UNF), is an important consideration in the fuel cycle. In this study, two sorbent materials, Cu 0 -polyacrylonitrile (PAN) and Bi 0 -PAN, were tested as solid sorbent candidates for iodine capture. Experiments using a thin bed of sorbent material were exposed to vaporized iodine for over 300 hours (~2 weeks) under varied conditions in a dynamic flow environment. The overall performance was monitored in real-time using thermogravimetric analysis, and the materials were subsequently characterized for surface and bulk analysis using scanning electron microscopy – energy-dispersive spectroscopy (SEM-EDS) and powder x-ray diffraction (pXRD), respectively. Iodine (in the form of I 2 ) is expected to be released primarily in the dissolver off-gas (DOG) stream; therefore, this study demonstrates the effects of elemental iodine (I 2 ), water vapor (H 2 O), and nitrogen dioxide (NO 2 ). When exposed to ‘ideal’ conditions (in which I 2 is carried by dry air), Cu 0 -PAN and Bi 0 -PAN behave differently, with TGA analysis indicating that Cu 0 sorption performance is higher than that of Bi 0 , as evidenced by a larger mass change. Under ‘harsh’ conditions—such as a gaseous feed containing I 2 , H 2 O, and NO 2 vapors carried by air,—iodine capture performance for both Cu 0 - and Bi 0 -PAN are affected. SEM-EDS and pXRD analysis of these materials is discussed herein.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

My vehicle is a data mine

In this talk we explore how analysis of vehicle data provides information of individual vehicle behaviors, information of other vehicles in the flow of traffic, and insights into the behavior of drivers. Over the last two decades traditional passenger vehicles have been transformed from integrated two-port electrical nodes to cyber-physical systems of communicating computational nodes whose individual state and control variables are shared on a standard controller area network (CAN) bus. As driver assistance systems have crept into vehicles as safety features, driver behaviors can be observed through analysis of the data streams on the CAN bus as these new nodes communicate with one another. The properties of these data streams, as well as architectures and approaches to gather the data, are important to consider when drawing conclusions on the relevance of the data in making decisions at varying levels of the information hierarchy. We will demonstrate several technical challenges associated with these data collection processes, as well as preliminary results that demonstrate application relevance of the data to behavior, traffic, and systems domains.

42 ENGINEERING↗
DOE OSTI2026

Moments of parton distribution functions of the pion from lattice QCD using gradient flow

We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.

lattice QCD↗
DOE OSTI2025

Processing behavior evolution of recycled polypropylene: An integrated experimental and Computer-Aided engineering simulation study

Polypropylene (PP) comprises 21% of global plastics production and 18% of plastics waste, yet less than 1% of solid-waste PP is recycled in the United States (U.S.), representing significant environmental and economic challenges. Mechanical recycling, the most prevalent recycling method, subject's materials to thermomechanical stresses, which typically degrade polymer properties, affecting the quality of polymer products. This study replicates the impact of mechanical recycling through multiple extrusion cycles to examine the effects on PP's processing behavior. Dynamic scanning calorimetry (DSC) measurements showed stable melting behavior across all processing conditions, while crystallization analysis exhibited consistent shifts in kinetic parameters. Rheological characterization demonstrated progressive viscosity reductions through successive cycles, particularly pronounced at elevated reprocessing temperatures. Here, the integration of this experimental data into injection molding simulations showed that recycled PP maintains viable processing characteristics. Our findings establish quantitative correlations between processing history and material behavior, enabling optimization of processing parameters directly rather than relying on trial-and-error approaches. While these results reflect idealized recycling conditions with minimal contamination, they provide a framework for understanding fundamental property evolution during mechanical recycling.

42 ENGINEERING↗
DOE OSTI2024

Once-Through Steam Generator Model Analysis Using Python and Advanced Optimization Tools (Summer Internship Report)

This study focuses on the parametric analysis of design parameters for a once-through steam generator (OTSG) model, using python and advanced optimization tools to facilitate applications such as the flowing autoclave steam generator (FASG) test cases. Building on previous research involving another OTSG with a different design, this project aims to enhance our understanding of how steam generators (SGs) behave and how their outputs are influenced by changes in design. The reason for this design change is to allow for more precise modeling and optimization of SG performance, to provide a comparative analysis between the two designs, and to set up the model for integration with the FASG test case. The OTSG python-model is a mathematical representation (including fluid flow and heat transfer equations/models/correlations) of a steam-generating unit in a pressurized water reactor-type small modular reactor system. Design studies involve changing the model’s input design parameters to observe the resulting effects on the output of the system. By using advanced optimization tools, such as the Risk Analysis Virtual Environment (RAVEN) developed at Idaho National Laboratory, detailed design parametric studies and model optimization were performed. Six input parameters—pressure, temperature, and mass flow rate for the inlet of the primary-side (hot fluid) and secondary-side (cold fluid), respectively, of the SG—were randomly perturbed via RAVEN’s Monte Carlo Sampler module, using uniform distributions (i.e., ±1%, ±5% and ±10% relative changes) for 600 samples. The analysis provides valuable insights into SG optimization and can be used for sensor placement optimization to effectively monitor and obtain experimental data in other tests.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Investigation of Anchor Link Magnetohydrodynamic Effects in the Toroidally Symmetric Lead-Lithium (TSLL) Blanket Concept

This paper presents a comprehensive study on the magnetohydrodynamic (MHD) flow in a slotted channel with a cylindrical anchor link, which is an imperative component of the toroidally symmetric lead-lithium (TSLL) blanket concept. Following the validation of the MHD solver implemented in COMSOL and a proposed subtraction approach to compute the anchor link pressure drop, the effects of computational domain size on the pressure drop and velocity distribution are examined. The results show that the pressure drop associated with the anchor link and the maximum flow velocity follow an asymptotic trend as the domain width increases, with wall-induced pressure drop being more dominant than that of the anchor link. The velocity distribution analysis revealed the formation of an internal boundary layer extending along the magnetic field direction, which is a unique feature of the investigated MHD flow. An estimation of the total MHD pressure drop associated with the array of anchor links under the TSLL blanket conditions suggests ~0.3 MPa, which is significantly lower than the recommended maximum allowable blanket pressure drop of 2 MPa. In conclusion, this work offers valuable insights into the anchor link–associated MHD phenomena and serves as a foundation for further development of the TSLL blanket concept for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Role of deep levels and barrier height lowering in current-flow mechanism in 150 μ m thick epitaxial n-type 4H–SiC Schottky barrier radiation detectors

Schottky barrier detectors (SBDs) require larger surface areas than conventional electronics to increase the detection efficiency although such SBDs manifest large diode ideality factors due to inhomogeneous areal distribution of surface barrier height (SBH). Inhomogeneous SBH distributions lead to various current flow mechanisms in SBDs, which need to be identified to optimize detector performance. In this Letter, we identify the current flow mechanism in large area Schottky barrier diodes for radiation detection fabricated on 150 μm thick n-4H–SiC epitaxial layers. The analysis of temperature-dependent forward current–voltage (I–V–T) characteristics of SBDs revealed two linear regions in current–voltage curves up to 450 K, one corresponding to the current flow through a low barrier patch, while the other corresponds to that of average barrier distribution. Applying a SBH distribution model to the reverse I–V–T characteristics, an activation energy of 0.76 eV for the current flow over the Schottky barrier was calculated. Additionally, the activation energy did not directly correspond to any of the defect levels observed from the deep level transient spectroscopy (DLTS). Above 450 K, a Schottky type barrier lowering suggested a current flow through a low barrier patch of ≈ 0.8 eV. The absence of any SBH lowering below 450 K indicated that the current corresponded to a neutrally charged trap level at ≈ 0.6 eV below the conduction band edge, which was consistent with DLTS measurements revealing the presence of an electron trap level Z1/2 at 0.59 eV below the conduction band edge.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

MPLBM-UT: Multiphase LBM library for permeable media analysis

MPLBM-UT is a specialized lattice-Boltzmann library that makes running single- and two-phase flow simulations in porous media accessible to everyone. We provide a suite of tools to pre-process computational domains for simulation, to set up custom boundary conditions, to run simulations, to post-process simulation outputs, and to visualize simulation results and data. All of these tools are easily accessible to users through the mplbm_utils Python package included in and automatically installed with MPLBM-UT. The high-performance, highly parallel library Palabos is used as the solver backend. MPLBM-UT is easily deployed in a variety of systems, from laptops to supercomputer clusters. MPLBM-UT also features multiple examples and benchmark templates that allow for fast prototyping of different porous media problems. We also provide an interface for reading in different file types and downloading domains from the Digital Rocks Portal to perform simulations.

97 MATHEMATICS AND COMPUTING↗