Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arrays”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Recovery of Oxpure 612C-50 Coconut Carbon Particles Using Bump Arrays

Cyanide is used to leach gold from crushed ore (solid matrix) into a gold cyanide solution. The gold is extracted from the cyanide solution by adsorption onto activated carbon. The gold extraction process occurs when the activated coconut carbon is placed into tanks that contain the gold cyanide solution either in a batch process or into a continuous flow circuit. The coconut carbon is then removed from the solution for gold recovery. The use of the mesofluidic separation system (a deterministic lateral displacement system) can significantly reduce costs and waste generated from the process of removing the coconut carbon from the solution. In the coconut carbon recovery process used at many gold mines in Nevada, the gold is attached to carbon particles that are still in the cyanide liquor. The mesofluidic system can rapidly remove the gold bearing carbon particles from the cyanide liquor quickly and with no operating costs. The benefits include: • Reduce operational costs related to filtration and hydro-cyclones • Reduced acid usage for the elution gold stripping phase as only 25% of the fluids will follow the carbon particles to the final express lane when two mesofluidic systems are used in series This alternative particle removal techniques would be advantageous. A promising technique for removing larger particles from slurries is mesofluidic separation, similar to deterministic lateral displacement arrays or “bump” arrays [1] but operates at much larger flow rates. As described by Pease, et al. [2], “Bump arrays in deterministic lateral displacement devices separate large particles from small particles using arrays of staggered posts. Large particles, defined as those with radii larger than the distance between the edge of a post and the stagnation streamline from the next downstream post, must bump toward one side of the device, whereas particles smaller than this distance slalom from entrance to exit without net lateral displacement.” Unlike filters or sieves, the posts in the arrays that cause separation do not block or occlude particle flow but work because particles go around the posts. Unlike hydrocyclones separation, separation is not driven by particle density, because gravitational forces are unimportant to mesofluidic separation. Burns, et al., and Pease, et al., have shown that particles may be separated from complex suspensions, under turbulent flow conditions, and at industrially important flow rates [3-5]. They have also recently shown that mesofluidic devices may be arranged in series to increase separation performance [6]. In this paper, we evaluate the mesofluidic system for the separation of commercial OxPure GR 612 Charred Coconut shell particles as a proof of principle test for the rapid separation using the 1500 micron cut mesofluidic separator. We first describe the experimental system and conditions. We then present the experimental results. [1] Huang, L.R., E.C. Cox, R.H. Austin, J.C. Sturm. 2004. Continuous particle separation through deterministic lateral displacement. Science, 304, 987–990. https://doi.org/10.1126/science.1094567. [2] Pease L.F., J.A. Bamberger, C.A. Burns, and M.J. Minette. 2021. Large Particle Separation from Non-Newtonian Slurries using Bump Arrays. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-65904, V003T08A023. New York, New York: ASME. doi:10.1115/FEDSM2021-65904 [3] Burns C.A., T.G. Veldman, J. Serkowski, R.C. Daniel, X.-Y. Yu, M.J. Minette, L.F. Pease. 2021. Mesofluidic separation versus dead-end filtration. Separation and Purification Technology, 254, 117256. https://doi.org/10.1016/j.seppur.2020.117256. [4] Pease L.F., J.E. Serkowski, T.G. Veldman, J. Williams, X.-Y. Yu, M.J. Minette, J.A. Bamberger, C.A. Burns. 2021. Can Bump Arrays Separate Particles from Turbulent Flows?. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-67696, V003T08A024. New York, New York: ASME. doi:10.11

mesofluidic separation, slurry, bump array↗
DOE OSTI2022

Narrow-Band Least-Squares Infrasound Array Processing

Infrasound data from arrays can be used to detect, locate, and quantify a variety of natural and anthropogenic sources from local to remote distances. However, many array processing methods use a single broad frequency range to process the data, which can lead to signals of interest being missed due to the choice of frequency limits or simultaneous clutter sources. In this work, we introduce a new open-source Python code that processes infrasound array data in multiple sequential narrow frequency bands using the least-squares approach. We test our algorithm on a few examples of natural sources (volcanic eruptions, mass movements, and bolides) for a variety of array configurations. Our method reduces the need to choose frequency limits for processing, which may result in missed signals, and it is parallelized to decrease the computational burden. Improvements of our narrow-band least-squares algorithm over broad-band least-squares processing include the ability to distinguish between multiple simultaneous sources if distinct in their frequency content (e.g., microbarom or surf vs. volcanic eruption), the ability to track changes in frequency content of a signal through time, and a decreased need to fine-tune frequency limits for processing. We incorporate a measure of planarity of the wavefield across the array (sigma tau, στ) as well as the ability to utilize the robust least trimmed squares algorithm to improve signal processing and insight into array performance. Our implementation allows for more detailed characterization of infrasound signals recorded at arrays that can improve monitoring and enhance research capabilities.

58 GEOSCIENCES↗
DOE OSTI2024

Development and Assessment of SNP Genotyping Arrays for Citrus and Its Close Relatives

Rapid advancements in technologies provide various tools to analyze fruit crop genomes to better understand genetic diversity and relationships and aid in breeding. Genome-wide single nucleotide polymorphism (SNP) genotyping arrays offer highly multiplexed assays at a relatively low cost per data point. We report the development and validation of 1.4M SNP Axiom® Citrus HD Genotyping Array (Citrus 15AX 1 and Citrus 15AX 2) and 58K SNP Axiom® Citrus Genotyping Arrays for Citrus and close relatives. SNPs represented were chosen from a citrus variant discovery panel consisting of 41 diverse whole-genome re-sequenced accessions of Citrus and close relatives, including eight progenitor citrus species. SNPs chosen mainly target putative genic regions of the genome and are accurately called in both Citrus and its closely related genera while providing good coverage of the nuclear and chloroplast genomes. Reproducibility of the arrays was nearly 100%, with a large majority of the SNPs classified as the most stringent class of markers, “PolyHighResolution” (PHR) polymorphisms. Concordance between SNP calls in sequence data and array data average 98%. Phylogenies generated with array data were similar to those with comparable sequence data and little affected by 3 to 5% genotyping error. Both arrays are publicly available.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Physics-constrained Gaussian process model for prediction of hydrodynamic interactions between wave energy converters in an array

To improve the efficiency of wave farms and achieve maximum power generation, the layout of wave energy converters (WECs) in an array needs to be carefully designed so that the hydrodynamic interactions can be positively exploited. For this, the hydrodynamic characteristics of the WEC array in different layouts need to be calculated. However, such calculations using numerical models usually entail significant computational cost, especially for large arrays of WECs. To address the computational challenge, a physics-constrained Gaussian process (GP) model is proposed to replace the original expensive numerical model and predict the hydrodynamic characteristics of the WECs for any array layout. By exploring the relationship between the WEC array (i.e., the input) and different hydrodynamic characteristics (i.e., the output), here we summarize a set of physical constraints/features, including invariance, symmetry, and additivity. This prior knowledge about the input-output relationship is then directly embedded in the constructed GP model through the design of physics-constrained kernels. In particular, a double-sum invariant kernel is first developed to incorporate the invariance and symmetry features, and then an additive kernel is developed to incorporate the additive feature of the problem. The invariant kernel and the additive kernel are then integrated to construct the physics-constrained GP model. Compared to the standard GP model, the proposed physics-constrained GP models require less training data to achieve the desired accuracy in predicting the hydrodynamic characteristics and are also less vulnerable to the curse of dimensionality (i.e., good scalability for large arrays) due to the use of an additive kernel. The efficiency, accuracy, and scalability of the proposed approach are demonstrated through an application to predict the hydrodynamic characteristics for WEC arrays of different sizes and layouts.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Metallic and porous Ti nanorod arrays for visible-IR light absorption and dendrite-free stable lithium-metal batteries

It is still a challenging task to prepare highly porous nanorod arrays of metals formed on substrates for optical and energy storage applications. Herein, we demonstrate the design and synthesis of black color, metallic and highly porous Ti nanorod arrays as novel current collectors for dendrite-free and highly stable Li-metal anodes. The high porosity of metallic nanorod arrays provides numerous heterogeneous nucleation sites and huge contact area and large space for the accommodation of Li metal. The conductive metallic Ti nanorod arrays enhance electrode integration. Effectively, it eliminates formation Li dendrites and demonstrates superior cycling stability over 300 cycles. Additionally, the unique porous structures of the nanorod arrays can decrease the amplitude of forced vibration in narrow space leading to light absorption. Interestingly, the metal is black instead of metallic color. The black metallic nanorod arrays can absorb more than 96% of both visible and infra-red lights. Furthermore, this black color metallic porous nanorod arrays may find additional applications in aerospace, energy, biomedical, defence, and chemical industries.

25 ENERGY STORAGE↗
DOE OSTI2025

Neutron detection efficiency measurements of organic glass scintillator arrays

We characterize the neutron detection efficiency of pixelated organic glass scintillator (OGS) arrays manufactured in two ways: an additive method in which the OGS is cast in place in a pixelated mold, and a subtractive method in which the scintillator is first cast as a large block before being cut in bars and bonded together to form the array. The cast-in-place array utilizes a mixture of “bis(9,9-dimethyl-9H-fluoren-2-yl)diphenylsilane” (P2), “tri([1,1'-biphenyl]-4-yl)(phenyl)silane” (PB3), and “4-bis(2-methylstyryl)benzene” (MSB), and 50% by weight of polystyrene is added to the mixture for subtractive manufactured array. We compare the performance of the OGS arrays to that of an EJ-299-33M plastic scintillator array that was manufactured using a subtractive process similar to that of the second OGS array.

Neutron detectors (cold, thermal, fast neutrons)↗
DOE OSTI2024

Ultrafast Annealing Improves SNR and Long-Term Stability of a Highly Multiplexed Line-by-Line FBG Array Inscribed by Femtosecond Laser in a Coreless Fiber for Extreme-Temperature Applications

This study reports the fabrication of a fourth-order line-by-line fiber Bragg grating (FBG) array using femtosecond laser inscription within a highly multimode coreless optical fiber, with a particular focus on achieving substantial multiplexing capabilities. An ultrafast annealing procedure is employed, resulting in an impressive enhancement of the FBG sensor’s fringe visibility by approximately 13 dB, signifying a notable improvement of approximately ~4 dB. This substantial enhancement contributes to the long-term stability and performance of the multiplexed FBG array in extreme-temperature conditions. The systematic fabrication approach employed for the multiplexed FBG array guarantees a high signal-to-noise ratio (SNR) for each individual FBG within the array. This FBG array is intended for extreme-temperature applications, addressing limitations associated with traditional FBGs based on doped optical fibers, including SNR degradation and temperature-induced fringe drift. Testing at temperatures up to 1120 °C demonstrates the FBG array’s stability without fluctuations in readings. Furthermore, it endures seven heat cycles, spanning from 500 °C to 1120 °C, over 60 h, exhibiting outstanding thermal stability. Furthermore, this highly multiplexed FBG array with an ultrafast annealing approach holds promise for extreme-temperature environments, such as steelmaking, where precise and reliable distributed temperature monitoring is imperative.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Sensitivity of Wind Turbine Array Downstream Effects to the Parameterization Used in WRF

The Weather Research and Forecasting (WRF) Model has been extensively used for wind energy applications, and current releases include a scheme that can be applied to examine the effects of wind turbine arrays on the atmospheric flow and electricity generation from wind turbines. Herein we present a high-resolution simulation using two different wind farm parameterizations: 1) the “Fitch” parameterization that is included in WRF releases and 2) the recently developed Explicit Wake Parameterization (EWP) scheme. We compare the schemes using a single yearlong simulation for a domain centered on the highest density of current turbine deployments in the contiguous United States (Iowa). Pairwise analyses are applied to diagnose the downstream wake effects and impact of wind turbine arrays on near-surface climate conditions. On average, use of the EWP scheme results in small-magnitude wake effects within wind farm arrays and faster recovery of full WT array wakes. This in turn leads to smaller impacts on near-surface climate variables and reduced array–array interactions, which at a systemwide scale lead to summertime capacity factors (i.e., the electrical power produced relative to nameplate installed capacity) that are 2%–3% higher than those from the more commonly applied Fitch parameterization. It is currently not possible to make recommendations with regard to which wind farm parameterization exhibits higher fidelity or to draw inferences with regard to whether the relative performance may vary with prevailing climate conditions and/or wind turbine deployment configuration. However, the sensitivities documented herein to the wind farm parameterization are of sufficient magnitude to potentially influence wind turbine array siting decisions. Thus, our research findings imply high value in undertaking combined long-term high-fidelity observational studies in support of model validation and verification.

17 WIND ENERGY↗
DOE OSTI2023

Dipole-lattice nanoparticle resonances in finite arrays

We investigate how the periodic lattices define the collective optical characteristics of the silicon and titanium nanoparticle arrays. We examine the effects of dipole lattice on the resonances of optical nanostructures, including those made of lossy materials, such as titanium. Our approach involves employing coupled-electric-magnetic-dipole calculations for finite-size arrays, as well as lattice sums for effectively infinite arrays. Our model shows that the convergence to the infinite-lattice limit is faster when the resonance is broad, requiring fewer array particles. Our approach differs from previous works by altering the lattice resonance through modifications in the array period. We observed that a higher number of nanoparticles is necessary to achieve convergence to the infinite-array limit. Additionally, we observe that the lattice resonances excited next to higher diffraction orders (such as second order) converge more quickly toward the ideal case of an infinite array than the lattice resonances related to the first diffraction order. This work reports on the significant advantages of using a periodic arrangement of lossy nanoparticles and the role of collective excitation in enhancing response from transition metals, such as titanium, nickel, tungsten, and so on. The periodic arrangement of nanoscatterers allows for the excitation of strong dipoles, boosting the performance of nanophotonic devices and sensors by improving the strength of localized resonances.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nanolaser arrays: toward application-driven dense integration

Abstract The past two decades have seen widespread efforts being directed toward the development of nanoscale lasers. A plethora of studies on single such emitters have helped demonstrate their advantageous characteristics such as ultrasmall footprints, low power consumption, and room-temperature operation. Leveraging knowledge about single nanolasers, the next phase of nanolaser technology will be geared toward scaling up design to form arrays for important applications. In this review, we discuss recent progress on the development of such array architectures of nanolasers. We focus on valuable attributes and phenomena realized due to unique array designs that may help enable real-world, practical applications. Arrays consisting of exactly two nanolasers are first introduced since they can serve as a building block toward comprehending the behavior of larger lattices. These larger-sized lattices can be distinguished depending on whether or not their constituent elements are coupled to one another in some form. While uncoupled arrays are suitable for applications such as imaging, biosensing, and even cryptography, coupling in arrays allows control over many aspects of the emission behavior such as beam directionality, mode switching, and orbital angular momentum. We conclude by discussing some important future directions involving nanolaser arrays.

Deka, Suruj S.↗
DOE OSTI2024

IEA Wind TCP Task 49: Reference Site Conditions for Floating Wind Arrays

This report, prepared within Work package 1 of IEA Wind Task 49, presents reference site conditions for floating wind arrays to serve as a design basis for the techno-economic design of reference floating wind arrays. Data of the reference sites presented here are publicly available in an open database and thus support fast development and comparable design of floating wind arrays for various relevant conditions. The development of these reference sites drew on existing open access datasets and ongoing research projects of task participants. Six classes were identified that describe relevant key conditions for the design and development of floating wind arrays: met-ocean conditions, seabed conditions, coastal infrastructure, environmental impact, socio-economic impact, as well as regulations and permissions. The reference sites for the techno-economic design of floating wind arrays are based on a concept with building blocks to synthesize purpose-built site representations. In each of the identified classes with influencing design factors, building blocks are used to describe the characteristic properties and their spread. However, the latter three classes (i.e., environmental impact, socio-economic impact, regulations and permissions) are not included in the reference site conditions due to limited knowledge and lack of reliable criteria to quantify their impact on the techno-economic design in numeric parameters. Building blocks with key parameters for the techno-economic design of floating wind arrays are provided for met-ocean conditions, seabed conditions, and coastal infrastructure. For met-ocean conditions, multiple sites were selected for detailed analysis that represent a range of conditions across the pipeline of floating wind projects. Wind conditions and sea states are separated, and each location considers both the severity of wind and waves e.g. one site may have a moderate wave condition but severe wind condition. From this pipeline, eleven representative sites were selected where both site-specific analysis was available within the consortium, and where they represent different parts of the global pipeline. The eleven sites are: Hannibal (Italy), Humboldt (US), Ulsan (South Korea), MoneyPoint One (Ireland), Havbredey (UK), Fukushima (Japan), Utsira Nord (Norway), Gulf of Maine (US), Sud de la Bretagne II (France), Sorlige Nordsjo II (Norway). Each of these sites is summarized in the main report while more details about the studies and analyses behind the datasets are provided in the appendix. For seabed conditions, general information about the geotechnical parameters is provided and a baseline is established for the geotechnical parameters and stratigraphy that may be encountered on the sites. A set of six 'synthetic cases' is defined as building blocks providing the different parameters required for design under each case/soil condition. For the coastal infrastructure, general information about the main requirements is provided that a port should comply with to provide a satisfactory service during the construction of floating offshore wind arrays. Minimum port infrastructural requirements are provided for three types of ports.

17 WIND ENERGY↗
DOE OSTI2020

Techno-Economic Modelling of Tidal Energy Converter Arrays in the Tacoma Narrows

Hydrokinetic tidal energy converter (TEC) technology is yet to become cost competitive with other renewable energy sources. Understanding the interaction between energy production and the costs incurred harvesting that energy may unlock the economic potential of this technology. Although hydrodynamic simulation of TEC arrays has matured over time, including demonstration of how small and large arrays affect the resource, integration of cost modelling is often limited. The advanced ocean energy array techno-economic modelling tool ‘DTOcean’ enables designers to calculate and improve the levelised cost of energy (LCOE) of an array through parametric simulation of the energy extraction, design of the electrical network, moorings and foundations, and simulation of the installation and lifetime operations and maintenance of the array. This work presents a verification of DTOcean’s ability to simulate the techno-economic performance of TEC arrays by reproducing the hypothetical RM1 reference model, a semi-analytical model of a TEC array based in the Tacoma Narrows of Washington state, U.S.A. It is demonstrated that DTOcean can produce a reasonable estimate to the LCOE predicted by the reference model, giving (in Euro cents per kiloWatt hour) 36.69 ¢/kWh against the reference model’s 34.612 ¢/kWh for 10 TECs, while for 50 TECs, DTOcean calculated 20.34 ¢/kWh compared to 17.34 ¢/kWh for the reference model.

Topper, Mathew B. R. (ORCID:0000000347324347)↗
DOE OSTI2020

The Tianlai Cylinder Pathfinder Array: System Functions and Basic Performance Analysis

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure of the Universe and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors are oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio quiet site in Xinjiang, and saw its first light in September 2016. In this first paper on the data analysis of the Tianlai cylinder array, we describe the system functions of the array, the commissioning observations during 2016- 2018, and present its basic performance characteristics. We show examples of the interferometric visibility data, then using the observational data, we derive the actual beam profile in the east-west direction, the bandpass response, and calibrate the complex gains for the array elements either with a strong astronomical point source, or with an artificial calibrator source. Based on the preliminary analysis of the data, we derive the system temperature and sensitivity of the array.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Ultra-wideband mutual coupling compensation of active electronically scanned arrays in multi-channel radar systems

Technologies pertaining to the calibration of ultra-wide-band active-electronically-scanned arrays to compensate for transmit and receive mutual coupling between array elements and/or non-ideal isolated channels are described herein. A plurality of near-field measurements are taken by way of probe or probes that may be moved among a plurality of positions aligned with respective elements in the array. For each position of the probe, each of the elements of the array is stimulated to transmit or receive a calibration signal to or from the probe, respectively. Frequency-domain transfer functions are computed from the received signals by the probe or the element for each of the array elements in each of the positions of the probe. The inverse of the matrix of transfer functions comprise frequency-domain transmit and receive correction factors that are used to modify desired array inputs/outputs such that the modified signals correct for mutual coupling between elements in the array.

42 ENGINEERING↗
DOE OSTI2022

Microwave assisted and low-temperature fabrication of nanowire arrays on scalable 2D and 3D substrates

A method of making a titanium dioxide nanowire array includes contacting a substrate with a solvent comprising a titanium (III) precursor, an acid, and an oxidant while microwave heating the solvent, thereby forming a hydrogen titanate H2Ti2O5.H2O nanowire array. The hydrogen titanate nanowire array is annealed to form a titanium dioxide nanowire array. The substrate is seeded with titanium dioxide before starting the hydrothermal synthesis of the hydrogen titanate nanowire array. The titanium dioxide nanowire array is loaded with a platinum group metal to form an exhaust gas catalyst. The titanium dioxide nanowire array can be used to catalyze oxidation of combustion exhaust.

Gao, Pu-Xian↗
DOE OSTI2024

Shading and Energy Impacts for Farmer-Friendly Agrivoltaics Arrays

When designing agrivoltaic solar arrays, farmers’ natural row tilling directions might be reasonable choices for array alignment, especially in cases where the tilling direction helps prevent soil erosion or has other benefits for the farmer’s operations. We examined the light field and energy production effects when installing arrays that have azimuth values that deviate from the typical ideal for three important array types: fixed-tilt, single-axis tracking, and vertical bifacial arrays. On a fractional annual energy production basis, though, all three array types lose very little output when rotated away from their ideal. This offers farmers and solar developers a wider range of system design to possibly align solar array rows with the farmer’s preferred tilling directions.

14 SOLAR ENERGY↗
DOE OSTI2020

Seismic monitoring and high-frequency noise using arrays in the Arabian Peninsula

Seismic arrays provide useful tools for regional seismic monitoring. Two small-aperture, regionally-oriented arrays, QWAR and HQAR, were deployed in Saudi Arabia and Oman in 2012 and 2016. We utilize a time-domain beampacking method, similar to frequency-wavenumber analysis, to examine the performance of the arrays in terms of slowness and azimuthal bias and event detection capabilities. Additionally, we investigate persistent ambient noise sources recorded by the arrays. We find that the arrays provide slowness vectors with biases comparable with similar-sized arrays of the International Monitoring System in other locations around the world. At QWAR, regional events of magnitude 3.0 and above are detected a majority of the time, and as magnitudes increase to 4.0 and above, the detection rate is greater than 82%. Finally, strong noise generation, primarily with slownesses characteristic of Lg waves, is found in the directions of the northern Arabian/Persian Gulf and southern Red Sea and may be a factor in event detection capabilities.

58 GEOSCIENCES↗
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↗