Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Baryonic imprints on DM haloes: the concentration–mass relation in the C amels simulations

ABSTRACT The physics of baryons in haloes, and their subsequent influence on the total matter phase space, has a rich phenomenology and must be well understood in order to pursue a vast set of questions in both cosmology and astrophysics. We use the Cosmology and Astrophysics with MachinE Learning Simulation (Camels) suite to quantify the impact of four different galaxy formation parameters/processes (as well as two cosmological parameters) on the concentration–mass relation, cvir−Mvir. We construct a simulation-informed non-linear model for concentration as a function of halo mass, redshift, and six cosmological/astrophysical parameters. This is done for two galaxy formation models, IllustrisTNG and Simba, using 1000 simulations of each. We extract the imprints of galaxy formation across a wide range in mass $M_{\rm vir}\in [10^{11}, 10^{14.5}] \, {\rm M}_\odot \, h^{-1}$ and in redshift z ∈ [0, 6] finding many strong mass- and redshift-dependent features. Comparisons between the IllustrisTNG and Simba results show the astrophysical model choices cause significant differences in the mass and redshift dependence of these baryon imprints. Finally, we use existing observational measurements of cvir−Mvir to provide rough limits on the four astrophysical parameters. Our non-linear model is made publicly available and can be used to include Camels-based baryon imprints in any halo model-based analysis.

Astronomy & Astrophysics↗
DOE OSTI2025

Litter Production and Foliar Nutrient Resorption in Pioneer and Non-Pioneer Species in a Selective Logging Experiment in the Central Amazon, BIONTE, ZF-2, Manaus, 2022-23

This dataset was collected near the city of Manaus, Brazil, at the Experimental Station of Tropical Forestry (EEST, aka “ZF2”), inside the BIONTE (BIOmass and NuTrient Experiment). The experiment included three levels of increasing selective logging intensity, along with control, with 1-hectare permanent plots (12 total) located at the center of 4-hectare treatment plots. The vegetation has a high floristic diversity, the soils of the region are poor in nutrients, and the topography is characterized by plateaus (where BIONTE is located), and also valley bottoms and slopes. Three treatments of differing logging intensities were applied in the BIONTE experiment (T1, T2 and T3). The study was conducted in Treatment 3 (Block I – permanent plot), which represents the most intensive logging treatment, with 69% of the basal area (m²∙ha⁻¹) removed in 1988. The present dataset spans the period from May 1, 2022, to May 1, 2023. The data package includes leaf_nutrient_data, litterfall_total_data, leaf_litterfall_species_specific_data, and species_info, all provided in .csv format. These formats allow users to process and analyze the data in various software applications and programming languages, such as Python and R. This dataset was collected to advance knowledge on nutrient cycling in Amazonian forests, specifically distinguishing between species with two distinct functional traits: fast-growing and slow-growing. It also aims to improve Earth System Models, such as the E3SM Functionally Assembled Terrestrial Ecosystem Simulator (FATES). Additionally, it was used in a paper currently in preparation (Carvalho et al., in prep.), which aims to quantify seasonal litter production and foliar nutrient resorption in pioneer (fast-growing) and non-pioneer (slow-growing) tree species in the central Amazon. Specifically, it seeks to answer two key questions: 1) Is there a difference in leaf litter production, leaf nutrient flux and leaf nutrient concentration between pioneers and non-pioneers species? Is there a difference in the efficiency of foliar nutrient resorption between pioneers and non-pioneers species?

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Statistics of Green's functions on a disordered Cayley tree and the validity of forward scattering approximation

The accuracy of the forward scattering approximation for two-point Green's functions of the Anderson localization model on the Cayley tree is studied. A relationship between the moments of the Green's function and the largest eigenvalue of the linearized transfer-matrix equation is proved in the framework of the supersymmetric functional-integral method. The new large-disorder approximation for this eigenvalue is derived and its accuracy is established. Using this approximation the probability distribution of the two-point Green's function is found and compared with that in the forward scattering approximation (FSA). It is shown that FSA overestimates the role of resonances and thus the probability for the Green's function to be significantly larger than its typical value. The error of FSA increases with increasing the distance between points in a two-point Green's function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Radial Steering Displacement Method to Measure Dispersion

The measurement of the dispersion function in a transfer line plays the crucial role for the development of an accurate optics model for the beam tracking simulation. In this technical note, we describe an experimental method to measure the dispersion value at multi-wire location MW063 in the NASA Space Radiation Laboratory (NSRL) transfer line. Finally, we use our MADX model to calculate the dispersion function at the same multi-wire location. The ultimate goal of this measurement and calculation is to understand the right value of dispersion function and the beam size (emittance) at the entrance of NSRL beam line.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

AUTOMATIC GENERATION OF EVENT TREES AND FAULT TREES: A MODEL-BASED APPROACH

In the past few decades, increasing complexity in modern engineering systems has been driven by the integration of a large number of components and by the fact that the system operations involve many disciplines (e.g., thermal-hydraulics, plant operations, cyber-security). Current safety/reliability modeling approaches to such systems are labor intensive, difficult to learn, and rely heavily on simplistic Boolean logic to depict failure propagation and accident progression. While these methods serve well for simple systems (i.e., linear causal systems with limited small inter- and intra-system interactions), their results are difficult to verify when modeling complex systems (typically performed through the extensive use of modeling assumptions). The development of new methods is addressed to meet these challenges through a model-based system engineering (MBSE) lens. Under MBSE philosophy, every aspect of the system (form or function) is represented by a model that completely characterizes its architecture or behavior. MBSE approach greatly improves the management of design, analysis and verification of complex systems. An integration of Dynamic Probabilistic Risk Assessment (DPRA) methods with MBSE models is proposed to perform safety/reliability analyses of engineering systems. In particular, MBSE representation of the system (performed using Systems Modeling Language [SysML]) is coupled with DPRA methods to automatically generate event trees and fault trees.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Meanders as a scaling motif for understanding of floodplain soil microbiome and biogeochemical potential at the watershed scale

Abstract Background Biogeochemical exports from watersheds are modulated by the activity of microorganisms that function over micron scales. Here, we tested the hypothesis that meander-bound regions share a core microbiome and exhibit patterns of metabolic potential that broadly predict biogeochemical processes in floodplain soils along a river corridor. Results We intensively sampled the microbiomes of floodplain soils located in the upper, middle, and lower reaches of the East River, Colorado. Despite the very high microbial diversity and complexity of the soils, we reconstructed 248 quality draft genomes representative of subspecies. Approximately one third of these bacterial subspecies was detected across all three locations at similar abundance levels, and ~ 15% of species were detected in two consecutive years. Within the meander-bound floodplains, we did not detect systematic patterns of gene abundance based on sampling position relative to the river. However, across meanders, we identified a core floodplain microbiome that is enriched in capacities for aerobic respiration, aerobic CO oxidation, and thiosulfate oxidation with the formation of elemental sulfur. Given this, we conducted a transcriptomic analysis of the middle floodplain. In contrast to predictions made based on the prominence of gene inventories, the most highly transcribed genes were relatively rare amoCAB and nxrAB (for nitrification) genes, followed by genes involved in methanol and formate oxidation, and nitrogen and CO 2 fixation. Within all three meanders, low soil organic carbon correlated with high activity of genes involved in methanol, formate, sulfide, hydrogen, and ammonia oxidation, nitrite oxidoreduction, and nitrate and nitrite reduction. Overall, the results emphasize the importance of sulfur, one-carbon and nitrogen compound metabolism in soils of the riparian corridor. Conclusions The disparity between the scale of a microbial cell and the scale of a watershed currently limits the development of genomically informed predictive models describing watershed biogeochemical function. Meander-bound floodplains appear to serve as scaling motifs that predict aggregate capacities for biogeochemical transformations, providing a foundation for incorporating riparian soil microbiomes in watershed models. Widely represented genetic capacities did not predict in situ activity at one time point, but rather they define a reservoir of biogeochemical potential available as conditions change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Defect graph neural networks for materials discovery in high-temperature clean-energy applications

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. In conclusion, it will help to tackle future materials discovery problems in clean energy and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Boosting the Oxygen Reduction Performance via Tuning the Synergy between Metal Core and Oxide Shell of Metal−Organic Frameworks‐Derived Co@CoO x

Abstract Co core@Co oxide shell (Co@CoO x ) catalysts represent a large family with promising oxygen reduction reaction (ORR) catalytic activity. However, inadequate understanding of Co@CoO x synergy prohibits further pursuit of catalytic performance enhancement. Herein, a Co zeolitic−imidazolate framework was converted into metallic Co, followed by controlled air treatment to form Co@CoO x . The composition and structure evolution as a function of air treatment temperature were studied thoroughly through conventional and synchrotron (both ex‐situ and in‐situ) characterizations, confirming the coexistence of CoO and Co 3 O 4 in the shell. The optimal catalyst showed an ORR half‐wave potential of 0.87 V (vs. RHE) in an alkaline half‐cell and delivered high discharge capacity in an aprotic Li−O 2 battery (7,124 mAh g Cat+C −1 ) and an aqueous Zn−air battery (694 mAh g Zn −1 ) with good performance retention after durability test. Modeling simulation and density functional theory calculation confirmed the charge donation from metal core to oxide shell and shed light on new insights of how metal@metal oxide synergy impacted the ORR via tuning the charge conductivity, oxygen affinity and intermediate transfer pathway. This work opens up a venue to boost ORR catalytic activity from an interfacial synergy perspective.

Wang, Hao↗
DOE OSTI2024

Structural insights reveal interplay between LAG-3 homodimerization, ligand binding, and function

Lymphocyte activation gene-3 (LAG-3) is an inhibitory receptor expressed on activated T cells and an emerging immunotherapy target. Domain 1 (D1) of LAG-3, which has been purported to directly interact with major histocompatibility complex class II (MHCII) and fibrinogen-like protein 1 (FGL1), has been the major focus for the development of therapeutic antibodies that inhibit LAG-3 receptor-ligand interactions and restore T cell function. Here, we present a high-resolution structure of glycosylated mouse LAG-3 ectodomain, identifying that cis-homodimerization, mediated through a network of hydrophobic residues within domain 2 (D2), is critically required for LAG-3 function. Additionally, we found a previously unidentified key protein-glycan interaction in the dimer interface that affects the spatial orientation of the neighboring D1 domain. Mutation of LAG-3 D2 residues reduced dimer formation, dramatically abolished LAG-3 binding to both MHCII and FGL1 ligands, and consequentially inhibited the role of LAG-3 in suppressing T cell responses. Intriguingly, we showed that antibodies directed against D1, D2, and D3 domains are all capable of blocking LAG-3 dimer formation and MHCII and FGL-1 ligand binding, suggesting a potential allosteric model of LAG-3 function tightly regulated by dimerization. Furthermore, our work reveals unique epitopes, in addition to D1, that can be targeted for immunotherapy of cancer and other human diseases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Interlayer magnetophononic coupling in MnBi2Te4

Abstract The emergence of magnetism in quantum materials creates a platform to realize spin-based applications in spintronics, magnetic memory, and quantum information science. A key to unlocking new functionalities in these materials is the discovery of tunable coupling between spins and other microscopic degrees of freedom. We present evidence for interlayer magnetophononic coupling in the layered magnetic topological insulator MnBi 2 Te 4 . Employing magneto-Raman spectroscopy, we observe anomalies in phonon scattering intensities across magnetic field-driven phase transitions, despite the absence of discernible static structural changes. This behavior is a consequence of a magnetophononic wave-mixing process that allows for the excitation of zone-boundary phonons that are otherwise ‘forbidden’ by momentum conservation. Our microscopic model based on density functional theory calculations reveals that this phenomenon can be attributed to phonons modulating the interlayer exchange coupling. Moreover, signatures of magnetophononic coupling are also observed in the time domain through the ultrafast excitation and detection of coherent phonons across magnetic transitions. In light of the intimate connection between magnetism and topology in MnBi 2 Te 4 , the magnetophononic coupling represents an important step towards coherent on-demand manipulation of magnetic topological phases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detecting Reactive Products in Carbon Capture Polymers with Chemical Shift Anisotropy and Machine Learning

Aminopolymers are attractive sorbents for CO 2 direct air capture applications due to their high density of amine groups, which can readily react with atmospheric levels of CO 2 to form chemisorbed species. The identity of these chemisorbed species and the functional groups that form upon oxidative degradation depends on both material properties and processing conditions, forming a variety of carbonyl-type sites such as ammonium carbamates, bicarbonates, carbonates, carbamic acids, ureas, and amides. 13 C solid-state nuclear magnetic resonance (NMR) is often used to help elucidate the identity of these reacted species, but it is challenging due to the narrow chemical shift range of carbonyl sites. Herein, we demonstrate the application of a two-dimensional (2D) chemical shift anisotropy (CSA) recoupling pulse sequence (ROCSA) to obtain CSA tensor values at each isotropic chemical shift, overcoming limitations of isotropic peak resolution. CSA tensor values describe the local chemical environment and can readily differentiate between the chemisorbed and degradation products. To aid identification, we also developed a k-nearest neighbor (kNN) classification model to distinguish the functional groups via their CSA tensor parameters. This methodology was demonstrated on poly(ethylenimine) in γ-Al 2 O 3 exposed to CO 2 and showed that the chemisorbed products are ammonium carbamate and a mixed carbamate–carbamic acid species. The sample was analyzed again after desorption at 100 °C inducing mild degradation, and the remaining products were strongly bound carbamate and urea species. In conclusion, the combination of 2D CSA measurements coupled with a kNN classification model enhances the ability to accurately identify chemisorbed or degradation products in complex carbon capture materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Effects of coupling a stochastic convective parameterization with the Zhang–McFarlane scheme on precipitation simulation in the DOE E3SMv1.0 atmosphere model

Abstract. A stochastic deep convection parameterization is implemented into the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1.0 (EAMv1). This study evaluates its performance in simulating precipitation. Compared to the default model, the probability distribution function (PDF) of rainfall intensity in the new simulation is greatly improved. The well-known problem of “too much light rain and too little heavy rain” is alleviated, especially over the tropics. As a result, the contribution from different rain rates to the total precipitation amount is shifted toward heavier rain. The less frequent occurrence of convection contributes to suppressed light rain, while more intense large-scale and convective precipitation contributes to enhanced heavy total rain. The synoptic and intraseasonal variabilities of precipitation are enhanced as well to be closer to observations. The sensitivity of the rainfall intensity PDF to the model vertical resolution is examined. The relationship between precipitation and dilute convective available potential energy in the stochastic simulation agrees better with that in the Atmospheric Radiation Measurement (ARM) observations compared with the standard model simulation. The annual mean precipitation is largely unchanged with the use of the stochastic scheme except over the tropical western Pacific, where a moderate increase in precipitation represents a slight improvement. The responses of precipitation and its extremes to climate warming are similar with or without the stochastic deep convection scheme.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Microbial metabolic dependency and its impacts on the soil carbon cycle (Final Technical Report, DE-SC0016364)

The terrestrial biosphere contains a large fraction of global carbon (C) and nearly 70% of the organic C in terrestrial systems is found in soils. Global changes in atmospheric CO 2 , temperature, precipitation, and ecosystem N inputs, will impact primary production and carbon inputs to soils, but it remains difficult to predict soil responses to these environmental conditions. Our ability to predict C-cycle responses to global change remains limited by unexplained variability within the terrestrial C-cycle. Our difficulty in predicting responses of the soil C-cycle to environmental change suggests a need for a greater understanding of the biotic mechanisms that govern soil C-transformations. Changes in microbial community structure and function have been proposed to impact soil processes, but mechanistic connections between microbial activity and C-cycle transformations remain poorly described. We assembled an interdisciplinary team to conduct cross-scale examination of microbial contributions to C-cycle dynamics and C-fate in soils. We investigated microbial contributions to soil C-transformations using a suite of complementary isotopic techniques including DNA-Stable Isotope Probing (DNA-SIP), metagenomic-SIP, metabolomic-SIP, proteomic-SIP, and microspatial analysis of soil organic C. These techniques were deployed in experiments that tested fundamental assumptions that underlie terrestrial C-cycle models. We characterized the functional characteristics of soil microorganisms responsible for major C-transformations in soil, evaluated the metabolic and ecological interactions that underlie soil C-cycle dynamics, and evaluated the degree to which these interactions impact rates of C mineralization and stabilization in soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Flexible and Generic Functional Mock-up Unit Based Threat Injection Framework for Grid-interactive Efficient Buildings: A Case Study in Modelica

Grid-interactive efficient buildings (GEBs) have been considered as an important asset to support the power grid reliability by utilizing the demand flexibility offered by GEBs. GEBs are enabled by advances in sensors and controls, and the communication between building equipment, whole buildings, and the grid. The integration of different building technologies and network-based communication system makes GEBs vulnerable to passive threats such as equipment failure and active threats such as cyber-attacks. Modeling and simulation is an effective way to evaluate the impact of threats on the system performance. This paper proposes a generic and flexible threat injection framework for commonly-used building energy simulators such as EnergyPlus and Modelica to support threat modeling and evaluation. This framework leverages functional mock-up unit (FMU) to develop a general modeling interface for threat injection and simulation. A numerical case study using Modelica as a building energy simulator is conducted to demonstrate the capability of the framework for supporting single/multiple-order threat modeling and simulation of a GEB. Four threats and their combinations are injected on a Modelica-based threat-free building energy and control system, including operating supply fan at its full speed, remotely cycling the chiller on and off, blocking the chiller from receiving the chilled water supply temperature setpoints, and hijacking the global zone air temperature setpoint. Simulation results show that the cyber-attack that leads to short-term signal blocking has small effects on the system operation due to the "self-healing" feature of the heating, ventilation, and air-conditioning (HVAC) interactive control system. The threat that takes control of resetting the global zone air temperature setpoints has the most adverse impact on the system energy use, peak power demand, thermal comfort and the provision of demand flexibility. The combination of four threats have aggregative effects on the system but the effects are less than the additive effects of the individual threat.

Fu, Yanyang↗
DOE OSTI2021

Impact of ion fluxes across thylakoid membranes on photosynthetic electron transport and photoprotection

In photosynthetic thylakoid membranes the proton motive force (pmf) not only drives ATP synthesis, in addition it is central to controlling and regulating energy conversion. As a consequence, dynamic fine-tuning of the two pmf components, electrical (Δψ) and chemical (ΔpH), is an essential element for adjusting photosynthetic light reactions to changing environmental conditions. Good evidence exists that the Δψ/ΔpH partitioning is controlled by thylakoid potassium and chloride ion transporters and channels. However, a detailed mechanistic understanding of how these thylakoid ion transporter/channels control pmf partitioning is lacking. In this work, we combined functional measurements on potassium and chloride ion transporter and channel loss-of-function mutants with extended mathematical simulations of photosynthetic light reactions in thylakoid membranes to obtain detailed kinetic insights into the complex interrelationship between membrane energization and ion fluxes across thylakoid membranes. The data reveal that potassium and chloride fluxes in the thylakoid lumen determined by the K + /H + antiporter KEA3 and the voltage-gated Cl – channel VCCN1/Best1 have distinct kinetic responses that lead to characteristic and light-intensity-dependent Δψ/ΔpH oscillations. These oscillations fine-tune photoprotective mechanisms and electron transport which are particularly important during the first minutes of illumination and under fluctuating light conditions. By employing the predictive power of the model, we unravelled the functional consequences of changes in KEA3 and VCCN1 abundance and regulatory/enzymatic parameters on membrane energization and photoprotection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Regression of exchangeable relational arrays

Relational arrays represent measures of association between pairs of actors, often in varied contexts or over time. Trade flows between countries, financial transactions between individuals, contact frequencies between school children in classrooms and dynamic protein-protein interactions are all examples of relational arrays. Elements of a relational array are often modelled as a linear function of observable covariates. Uncertainty estimates for regression coefficient estimators, and ideally the coefficient estimators themselves, must account for dependence between elements of the array, e.g., relations involving the same actor. Existing estimators of standard errors that recognize such relational dependence rely on estimating extremely complex, heterogeneous structure across actors. This paper develops a new class of parsimonious coefficient and standard error estimators for regressions of relational arrays. Here we leverage an exchangeability assumption to derive standard error estimators that pool information across actors, and are substantially more accurate than existing estimators in a variety of settings. This exchangeability assumption is pervasive in network and array models in the statistics literature, but not previously considered when adjusting for dependence in a regression setting with relational data. We demonstrate improvements in inference theoretically, via a simulation study, and by analysis of a dataset involving international trade.

97 MATHEMATICS AND COMPUTING↗