Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed analytics”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Clustering of the Diffuse Infrared Light from the COBE DIRBE Maps: Power Spectrum Analysis and Excess Isotropic Component of Fluctuations - 3

The cosmic infrared background (CIB) radiation is the cosmic repository for energy release throughout the history of the universe. The spatial fluctuations of the CIB resulting from galaxy clustering are expected to be at least a few percent on scales of a degree, depending on the luminosity and clustering history of the early universe. Using the all-sky data from the COBE DIRBE instrument at wavelengths 1.25 - 100 microns we attempt to measure the CIB fluctuations. In the near-IR, foreground emission is dominated by small scale structure due to stars in the Galaxy. There we find a strong correlation between the amplitude of the fluctuations and Galactic latitude after removing bright foreground stars. Using data outside the Galactic plane (absolute value of b > 20 deg) and away from the center (90 deg < l < 270 deg) we extrapolate the amplitude of the fluctuations to cosec absolute value of b = 0. We find a positive intercept of delta.F(sub rms) = 15.5(sup +3.7, sub -7.0), 5.9(sup +1.6, sub -3.7), 2.4(sup +0.5, sub -0.9), 2.0(sup +0.25, sub -0.5) nW/sq m.sr at 1.25, 2.2, 3.5 and 4.9 microns respectively, where the errors are the range of 92% confidence limits. For color subtracted maps between band 1 and 2 we find the isotropic part of the fluctuations at 7.6(sup +1.2, sub -2.4) nW/sq m.sr. Based on detailed numerical and analytic models, this residual is not likely to originate from the Galaxy, our clipping algorithm, or instrumental noise. We demonstrate that the residuals from the fit used in the extrapolation are distributed isotropically and suggest that this extra variance may result from structure in the CIB. We also obtain a positive intercept from a linear combination of maps at 1.25 and 2.2 microns. For 2 deg < theta < 15 deg, a power-spectrum analysis yields limits of (theta/5 deg) x delta.F(sub rms)(theta) < 6, 2.5, 0.8, 0.5 nW/sq m.sr at 1.25, 2.2, 3.5 and 4.9 microns respectively. From 10 - 100 microns, the dominant foregrounds are emission by dust in the Solar system and the Galaxy. There the upper limits on the CIB fluctuations are below 1 nW/sq m.sr and are lowest (< equal 0.5 nW/sq m.sr) at 25 microns.

Kashlinsky, A.↗
NASA NTRSDate not supplied

Center for Isotope Cosmochemistry and Geochronology at NASA Johnson Space Center

The Center for Isotope Cosmochemistry and Geochronology (CICG) at the NASA Johnson Space Center, led by Justin Simon, applies novel elemental and isotopic analytical techniques to understand solar system formation, from processes occurring during the earliest Solar Nebula, to the evolution of planetary bodies during planet formation and differentiation. CICG consists of an integrated suite of labs including dedicated spaces for specialized handling and preparation of meteorites and planetary return samples, metal-free clean chemistry, and mass spectrometry. In addition to a Thermo-Scientific Triton thermal ionization mass spectrometer (TIMS) for high precision isotopic analyses in processed samples, CICG houses two additional mass spectrometry instruments with added laser ablation capabilities for high spatial resolution in-situ measurements: (1) a Thermo-Scientific Element-XR inductively coupled mass spectrometer (ICP-MS) with a Teledyne Photon Machines Analyte Excite+ laser ablation system, and (2) a Nu Sapphire 1700 multi-collector ICP-MS, a large geometry instrument with a multi-pole collision cell for interference removal, with an Applied Spectra J200 tandem laser induced breakdown spectroscopy (LIBS) femto-second laser ablation system. We detail the analytical capabilities of CICG in the context of current and future research projects related to the studies of protoplanetary disk evolution via early-formed Solar System solids, volatile abundances and distribution within the inner solar system, and geochemical and temporal investigations related to planet formation and differentiation. We also highlight opportunities through the relationship with the Center of the Advancement of Space Safety and Mission Assurance Research (CASSMAR) at the University of Texas at El Paso. CICG is open to collaborations and participates in the NASA Facility for Astromaterials Research (NFAR) program, which accepts research proposals for access to the unique facility capabilities.

Jacob B. Setera↗
DOE OSTI2020

Distributed Small-Signal Stability Conditions for Inverter-Based Unbalanced Microgrids

The proliferation of inverter-based generation and advanced sensing, controls, and communication infrastructure have facilitated accelerated deployment of microgrids. A coordinated network of microgrids can maintain reliable power delivery to critical facilities during extreme events. Low-inertia offered by the power-electronics–interfaced energy resources, however, can present significant challenges to ensuring stable operation of the microgrids. In this work, distributed small-signal stability conditions for inverter-based microgrids are developed that involve the droop-controller parameters and the network parameters (e.g. line impedances, loads). The distributed closed-form parametric stability conditions derived in this paper can be verified in a computationally efficient manner, facilitating reliable design and operations of networks of microgrids. Dynamic phasor models have been used to capture the effects of electromagnetic transients. Furthermore, numerical results are presented, along with PSCAD simulations, to validate the analytical stability conditions. Effects of design choices, such as the conductor types, and inverter sizes, on the small-signal stability of inverter-based microgrids are investigated to derive useful engineering insights.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2005

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗
DOE OSTI2023

An Online Joint Optimization–Estimation Architecture for Distribution Networks

Here in this article, we propose an optimal joint optimization-estimation architecture for distribution networks, which jointly solves the optimal power flow (OPF) problem and static state estimation (SE) problem through an online gradient-based feedback algorithm. The main objective is to enable a fast and timely interaction between the OPF decisions and state estimators with limited sensor measurements. First, convergence and optimality of the proposed algorithm are analytically established. Then, the proposed gradient-based algorithm is modified by introducing statistical information of the inherent estimation and linearization errors for an improved and robust performance of the online OPF decisions. Overall, the proposed method eliminates the traditional separation of operation and monitoring, where optimization and estimation usually operate at distinct layers and different time scales. Hence, it enables a computationally affordable, efficient, and robust online operational framework for distribution networks under time-varying settings.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Prediction for Pressure Differential Across HEPA Filter Media Based on Media Characteristics and Particle Size Distribution

A new method for predicting the pressure drop across High Efficiency Particulate Air (HEPA) filter media is proposed based upon mass deposited onto the filter and known physical characteristics of the filter media. Detailed are the methods used in conjunction with current filter loading models to predict the pressure drop, as well as tests conducted to validate the methods. The benefit of a prediction model for practical use lies in the manufacturing and service life of nuclear grade HEPA filters. HEPA filters for use in nuclear facilities have a prescribed expiration date and a maximum allowable operating pressure drop. Therefore, the ability to predict the pressure drop across a filter and relate it to an expected length of service time can enable a reduction in wasted filters. This will allow for more informed decisions to be made based upon the dictated life cycle of the filters. Also, understanding and predicting how the pressure drop of HEPA filter media behaves as a function of physical characteristics and loaded mass can assist in future design and manufacturing of filter media. Currently, most existing pressure drop models are either computationally based or analytical methods relying on data gathered during tests. Neither are practical for prediction; the computational methods are difficult to implement, and the current analytical methods are more useful as tools for analysis. The current analytical model developed by Bergman et al. is based upon the pressure drop from each media fiber and modeling the deposited particles as newly formed fibers. Using Bergman's model and the media properties, air properties, and known aerosol particle size distribution, a pressure drop curve as a function of loaded mass can be created. This curve implies initial loading in the depth of the filter media with a transition to pure surface loading on the filter face, and a sensitivity to an evolving particle size distribution as mass continues to load onto the filter. To implement the predictive model, knowledge of the mean filter fiber diameter and porosity of the filter media is required. Traditionally, the mean fiber diameter is calculated as an effective diameter from prior media testing, however, in this study a Scanning Electron Microscope (SEM) technique is used to acquire this variable. Validation of the predictive model is provided by flat sheet media tests under a controlled environment with a measured particle size distribution of the challenge aerosol and shows reasonable preliminary agreement.

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

Cesium Batch Contact Equilibrium Testing of Crystalline Silicotitanate (CST) Sieve Cuts in SRS Average Simulant and Examination of CST Samples Aged in Caustic and Simulant in Support of Tank Closure Cesium Removal 1A (TCCR-1A)

Batch contact testing to determine cesium equilibrium loading on Crystalline Silicotitanate ion exchange media in Savannah River Site Average Waste Simulant at 25 °C indicated that smaller diameter particles isolated by sieving pretreated CST media may load slightly higher (<10%) amounts of cesium, though the differences are within analytical uncertainty. In addition, ion exchange media sub-samples stored in 2-4 M NaOH and caustic simulant solutions for ~2.5 years were examined by optical microscopy and the 3 M NaOH sample was also analyzed to determine whether changes in the particle size distribution occurred during storage. No visual indications of particle attrition or agglomeration were observed for any sample. Particle size analysis indicated that a slight decrease occurred in the average particle diameter following contact with 3 M NaOH (541 µm average diameter versus 566 µm for the pretreated CST prior to contact). A small increase (from 0 to <0.5 wt. %) in the number of particles ranging from 271 and 322 µm was also observed for the CST sample contacted with 3 M NaOH relative to a sample of the original pretreated material. However, this small change could be due to sub-sampling differences or analytical uncertainty. It does not appear that small particles are formed to a significant degree during CST caustic contact or that small particles which do form (presumably from attrition of larger particles during pretreatment) load significantly more cesium than the bulk material. Minimal other negative consequences were observed associated with CST extended caustic or simulant contact, except for the tendency for more concentrated salt solutions to form some salt crystals which deposit on the media over time.

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

Structural Trends and Vibrational Analysis of N,N,N′,N′−Tetramethylmalonamide Complexes Across the Lanthanide Series

Fundamental understanding of coordination chemistry across the lanthanide series is essential for explaining chemical behavior of rare-earth metals in complex liquid-liquid extraction processes, which in turn affects the distribution ratios and efficacy of separations as a whole. In this work, we explore the structural trends between the lanthanides and a neutral N,N,N',N'-tetramethylmalonamide (TMMA) ligand within four isolated families of solid-state compounds: Ln(trans-TMMA) 2 (NO 3 ) 3 Ln=La-Nd, Sm; Ln(cis-TMMA) 2 (NO 3 ) 3 Ln=Eu-Tb, Er; [Ln(TMMA) 3 (NO 3 ) 2 ][Ln(TMMA)(NO 3 ) 4 ] Ln=Dy-Tm; Ln(Κ 2 -TMMA)(iPrOH)(NO 3 ) 3 and Ln(Κ 1 -TMMA)(Κ 2 -TMMA)(NO 3 ) 3 Ln=Yb, Lu. Moving across the lanthanide series, we note the formation of both discrete charge-neutral complexes, as well as charged molecular anion-cation pairs, with variations in spatial ligand arrangement, coordination numbers, and ligand denticities. IR and Raman spectroscopy paired with DFT frequency calculations were used for an in-depth investigation of vibrational modes unique to each structural family. The collection of isolated model compounds was also discussed in the context of liquid-liquid separations based on reported distribution ratios from malonamide extraction.

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

Analytical Expression for DC Link Capacitor Current in a Cascaded H-Bridge Multi-Level Active Front-End Converter

Medium-voltage grid-tied systems often use a cascaded H-bridge multi-level active front-end. In this converter, dc link bus capacitors play an important role in stabilizing the converter and enabling both active and reactive power injections. The present work provides analytical expressions for the capacitor current, which are essential for optimizing system design (especially capacitor size vs. lifetime). Then, the expression is incorporated into the grid connected bidirectional power system model. Consequently, this work contributes to the guiding principles to choose accurate dc link capacitor ratings against grid-side power delivery requirements. The analytical results have been validated with detailed simulations and hardware results.

cascaded H-bridge, carrier-based PWM, dc-link capa↗
NASA NTRS1994

Testing the gravitational instability hypothesis?

We challenge a widely accepted assumption of observational cosmology: that successful reconstruction of observed galaxy density fields from measured galaxy velocity fields (or vice versa), using the methods of gravitational instability theory, implies that the observed large-scale structures and large-scale flows were produced by the action of gravity. This assumption is false, in that there exist nongravitational theories that pass the reconstruction tests and gravitational theories with certain forms of biased galaxy formation that fail them. Gravitational instability theory predicts specific correlations between large-scale velocity and mass density fields, but the same correlations arise in any model where (a) structures in the galaxy distribution grow from homogeneous initial conditions in a way that satisfies the continuity equation, and (b) the present-day velocity field is irrotational and proportional to the time-averaged velocity field. We demonstrate these assertions using analytical arguments and N-body simulations. If large-scale structure is formed by gravitational instability, then the ratio of the galaxy density contrast to the divergence of the velocity field yields an estimate of the density parameter Omega (or, more generally, an estimate of beta identically equal to Omega(exp 0.6)/b, where b is an assumed constant of proportionality between galaxy and mass density fluctuations. In nongravitational scenarios, the values of Omega or beta estimated in this way may fail to represent the true cosmological values. However, even if nongravitational forces initiate and shape the growth of structure, gravitationally induced accelerations can dominate the velocity field at late times, long after the action of any nongravitational impulses. The estimated beta approaches the true value in such cases, and in our numerical simulations the estimated beta values are reasonably accurate for both gravitational and nongravitational models. Reconstruction tests that show correlations between galaxy density and velocity fields can rule out some physically interesting models of large-scale structure. In particular, successful reconstructions constrain the nature of any bias between the galaxy and mass distributions, since processes that modulate the efficiency of galaxy formation on large scales in a way that violates the continuity equation also produce a mismatch between the observed galaxy density and the density inferred from the peculiar velocity field. We obtain successful reconstructions for a gravitational model with peaks biasing, but we also show examples of gravitational and nongravitational models that fail reconstruction tests because of more complicated modulations of galaxy formation.

Babul, Arif↗
DOE OSTI2024

Digital Multiphase Composites via Additive Manufacturing

Abstract Mechanical properties of traditional engineering materials are typically coupled to each other, presenting a challenge to practitioners with multi‐dimensional material property requirements. In this work, continuous, independent control over multiple mechanical properties is demonstrated in composite materials realized using additive manufacturing. For the first time, composites additively manufactured from rigid plastic, soft elastomer, and liquid constituents are experimentally characterized, demonstrating materials which span four orders of magnitude in modulus and two orders of magnitude in toughness. By forming analytical mappings between relative concentrations of constituents at the microscale and resulting macroscale material properties, inverse material design is enabled; the method is showcased by printing artifacts with prescribed toughness and elasticity distributions. The properties of these composites are placed in the context of biological tissues, showing they have promise as mechanically plausible tissue mimics.

Smith, Lawrence T.↗
DOE OSTI2021

An exploration of benchtop X–ray emission spectroscopy for precise characterization of the sulfur redox state in cementitious materials

The evolution of sulfur chemistry in cements is best known in the bailiwick of failure mechanisms via sulfate attack but is equally important for its contributions to the reduction capacity of cementitious materials often used for immobilizing nuclear waste streams destined for long-term storage, e.g., cementitious waste forms (CWF). The total reduction capacity of CWFs, encompassing contributions from both S and Fe reductants, and its implications toward radionuclide immobilization is most often studied by destructive wet chemistry methods requiring acid digestion in the presence of Ce(IV) and subsequent titration and colorimetric interpretation. Here, we investigate a similarly analytical but nondestructive alternative, benchtop high resolution wavelength-dispersive x-ray fluorescence spectroscopy, most commonly known as x-ray emission spectroscopy (XES), for probing the bulk sulfur oxidation state distribution. We present here an initial investigation of S XES, including an improved experimental protocol for lab XES of inhomogeneous samples, both as a complement to the Ce(IV) test and for new scientific opportunities that it enables for observing changes in sulfur chemistry. We discuss future improvements and opportunities, including: (1) the practical challenges associated with coordinating XES and Ce(IV) liquid extraction for a more comprehensive perspective on reduction capacity and for a high-precision evaluation of uncertainties in the Ce(IV) test; and (2) new opportunities, due to the nondestructive nature of XES, for controlled evolution studies aimed at elucidating specific chemical responses of CWF’s exposed to invasive gas or liquid species or to accelerated aging by radiative dose or thermal treatment.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The impact of gas purity on observed reactivity with NO using inductively coupled plasma tandem mass spectrometry

Interference removal in inductively coupled plasma tandem mass spectrometry (ICP-MS/MS) is strongly dependent on the gas selected for use within the collision/reaction cell. There has been little investigation on the effects that reaction gas impurities may have on the resulting spectra. The reactivity of 60 elements was evaluated using nitric oxide (NO 99.5%) with and without a gas purifier to reduce H 2 O impurities to <100 pptV. Experiments were performed using V, Ce, Tl and Th to investigate the effects of purified NO at various flowrates (0.22–1.49 mL min -1 ). Purified NO was shown to significantly mitigate oxy-hydride interferences, improve total ion sensitivity (notable at high gas flows), and shift product distributions advantageously. The reduction in oxy-hydride species results in a product distribution favoring the major expected products, where signals were shown to increase by an order of magnitude. Reduced background and increased signal for the major expected products provides avenues for improving various analytical applications of ICP-MS/MS.

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

Coupling of electroelastic dynamics and direct contact ultrasonic drying formulation for annular piezoelectric bimorph transducers

A newly developed technique for drying clothes without thermal energy has been developed through the utilization of ultrasonic vibrations from piezoelectric transducers. The novel technique incorporates the actuation of a thin stainless steel disk in contact with wet fabric via annular piezoelectric rings, where water in the liquid form is atomized, transported through microchannels in the disk, and ejected as a mist. In such a system, resonance matching between the actuation portion of the transducer and the portion contacting fabric must be realized, with theoretical results from the developed electromechanical model showing a reduction in energy consumption by 50% when resonance matching is achieved. The electrically coupled distributed parameter model for an annular bimorph piezoelectric transducer is developed for optimization of ultrasonic drying technology. The thickness mode vibrations are shown to dominate the behavior of the system, where the analytically developed model can be optimized to increase the output acceleration of the transducer, thus increasing drying performance. Here, the electromechanical equation developed will be connected to the drying rates of fabrics in contact with said vibrations, where the novelty of the coupled equations and its description of the physics of ultrasonic drying will be discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Emergent quantum mechanics at the boundary of a local classical lattice model

We formulate a model in which quantum mechanics emerges from classical mechanics. Given a local Hamiltonian H acting on n qubits, we define a local classical model with an additional spatial dimension whose boundary dynamics is approximately—but to arbitrary precision—described by Schrödinger's equation and H. The bulk consists of a lattice of classical bits that propagate towards the boundary through a circuit of stochastic matrices. The bits reaching the boundary are governed by a probability distribution whose deviation from the uniform distribution can be interpreted as the quantum-mechanical wave function. Bell nonlocality is achieved because information can move through the bulk much faster than the boundary speed of light. Finally, we analytically estimate how much the model deviates from quantum mechanics, and we validate these estimates using computer simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Hydrodynamic fluctuations near a Hopf bifurcation: Stochastic onset of vortex shedding behind a circular cylinder

Here, we investigate hydrodynamic fluctuations in the flow past a circular cylinder near the critical Reynolds number Re c for the onset of vortex shedding. Starting from the fluctuating Navier-Stokes equations, we perform a perturbation expansion around Re c to derive analytical expressions for the statistics of the fluctuating lift force. Molecular-level simulations using the direct simulation Monte Carlo method support the theoretical predictions of the lift power spectrum and amplitude distribution. Notably, we have been able to collect sufficient statistics at distances Re ⁡/ Re c – 1 = O ⁡(10 –3 ) from the instability that confirm the appearance of non-Gaussian fluctuations, and we observe that they are associated with intermittent vortex shedding. These results emphasize how unavoidable thermal-noise-induced fluctuations become dramatically amplified in the vicinity of oscillatory flow instabilities and that their onset is fundamentally stochastic.

42 ENGINEERING↗
DOE OSTI2025

Grey-Box System Identification of Grid-Forming Inverters

This paper demonstrates the use of grey-box system identification methods for simplifying and understanding the nonlinear power dynamics of grid-forming inverters (GFMs). The power and frequency outputs of complex high-order GFM models are fed into system identification software in order to fit them to a predetermined LTI system and learn system parameters such as (synthetic) inertia and droop constants. The same process is then run for a high-order synchronous generator model, and the outputs are fit to the same set of LTI equations. Simulation of a network of GFM inverters with diverse control architecture is also performed for the same process. The intent is threefold: first, to demonstrate the appropriateness of unified LTI models for describing the power and frequency dynamics of individual resources and connected networks, in order to facilitate analysis of larger heterogeneous networked systems; second, to discover the relationship between internal control parameters of GFMs and their externally observed values; and third, to validate that grey-box data-driven system identification techniques can be a valuable tool to discover the values of important parameters in the absence of explicit vendor models.

analytical models↗
DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗