Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Systems”

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,135 records · Page 63

Hybrid PDES Simulation of HPC Networks Using Zombie Packets

Although high-fidelity network simulations have proven to be reliable and cost-effective tools to peer into architectural questions for high-performance computing (HPC) networks, they incur a high resource cost. The time spent in simulating a single millisecond of network traffic in the highest detail can take hours, even for static, well-behaved traffic patterns such as uniform random. Surrogate models offer a significant reduction in runtime, yet they cannot serve as complete replacements and should only be used when appropriate. Thus, there is a need for hybrid modeling, where high-fidelity simulation and surrogates run side-by-side. Here, we present a surrogate model for HPC networks in which: packets bypass the network, while the network state is left untouched, i.e., suspended. To bypass the network, we use historical data to estimate the arrival time at which every packet should be scheduled at; to suspend the network, all in-flight packets are scheduled to arrive at their destinations, and are kept in the system to awaken as zombies when switching back to high-fidelity. Speedup for a hybrid model is relative to the proportion of surrogate to high-fidelity. This light-weight surrogate obtained up to 76× speedup. Keeping the zombies in the network showed an increase in the accuracy of the high-fidelity simulation on restart when compared to restarting the network from an empty state.

HPC networks↗

CyBERT: Cybersecurity Claim Classification by Fine-Tuning the BERT Language Model

We introduce CyBERT, a cybersecurity feature claims classifier based on bidirectional encoder representations from transformers and a key component in our semi-automated cybersecurity vetting for industrial control systems (ICS). To train CyBERT, we created a corpus of labeled sequences from ICS device documentation collected across a wide range of vendors and devices. This corpus provides the foundation for fine-tuning BERT’s language model, including a prediction-guided relabeling process. We propose an approach to obtain optimal hyperparameters, including the learning rate, the number of dense layers, and their configuration, to increase the accuracy of our classifier. Fine-tuning all hyperparameters of the resulting model led to an increase in classification accuracy from 76% obtained with BertForSequenceClassification’s original architecture to 94.4% obtained with CyBERT. Furthermore, we evaluated CyBERT for the impact of randomness in the initialization, training, and data-sampling phases. CyBERT demonstrated a standard deviation of ±0.6% during validation across 100 random seed values. Finally, we also compared the performance of CyBERT to other well-established language models including GPT2, ULMFiT, and ELMo, as well as neural network models such as CNN, LSTM, and BiLSTM. The results showed that CyBERT outperforms these models on the validation accuracy and the F1 score, validating CyBERT’s robustness and accuracy as a cybersecurity feature claims classifier.

97 MATHEMATICS AND COMPUTING↗

Stochastic quantum Krylov protocol with double-factorized Hamiltonians

Here we propose a class of randomized quantum Krylov diagonalization (rQKD) algorithms capable of solving the eigenstate estimation problem with modest quantum resource requirements. Compared to previous real-time evolution quantum Krylov subspace methods, our approach expresses the time evolution operator e –i$\widehat{H}$$\tau$ as a linear combination of unitaries and subsequently uses a stochastic sampling procedure to reduce circuit depth requirements. While our methodology applies to any Hamiltonian with fast-forwardable subcomponents, we focus on its application to the explicitly double-factorized electronic-structure Hamiltonian. To demonstrate the potential of the proposed rQKD algorithm on near-term quantum devices, we provide numerical benchmarks for a variety of molecular systems with circuit-based state-vector simulators including the effects of sampling noise, achieving ground-state energy errors of less than 1 kcal mol -1 with circuit depths orders of magnitude shallower than those required for low-rank deterministic Trotter-Suzuki decompositions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Flow Redirection and Induction in Steady State (FLORIS) Wind Plant Power Production Data Sets

This dataset contains turbine- and plant-level power outputs for 252,500 cases of diverse wind plant layouts operating under a wide range of yawing and atmospheric conditions. The power outputs were computed using the Gaussian wake model in NREL's FLOw Redirection and Induction in Steady State (FLORIS) model, version 2.3.0. The 252,500 cases include 500 unique wind plants generated randomly by a specialized Plant Layout Generator (PLayGen) that samples randomized realizations of wind plant layouts from one of four canonical configurations: (i) cluster, (ii) single string, (iii) multiple string, (iv) parallel string. Other wind plant layout parameters were also randomly sampled, including the number of turbines (25-200) and the mean turbine spacing (3D-10D, where D denotes the turbine rotor diameter). For each layout, 500 different sets of atmospheric conditions were randomly sampled. These include wind speed in 0-25 m/s, wind direction in 0 deg.-360 deg., and turbulence intensity chosen from low (6%), medium (8%), and high (10%). For each atmospheric inflow scenario, the individual turbine yaw angles were randomly sampled from a one-sided truncated Gaussian on the interval 0 deg.-30 deg. oriented relative to wind inflow direction. This random data is supplemented with a collection of yaw-optimized samples where FLORIS was used to determine turbine yaw angles that maximize power production for the entire plant. To generate this data, a subset of cases were selected (50 atmospheric conditions from 50 layouts each for a total of additional 2,500 cases) for which FLORIS was re-run with wake steering control optimization. The IEA onshore reference turbine, which has a 130 m rotor diameter, a 110 m hub height, and a rated power capacity of 3.4 MW was used as the turbine for all simulations. The simulations were performed using NREL's Eagle high performance computing system in February 2021 as part of the Spatial Analysis for Wind Technology Development project funded by the U.S. Department of Energy Wind Energy Technologies Office. The data was collected, reformatted, and preprocessed for this OEDI submission in May 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub repository under explore_wind_plant_data_h5.ipynb.

AI↗

Emerging Technologies for Privacy Preservation in Energy Systems

This study explores the intersection of digitalization and privacy within the energy sector, focusing on the emerging challenges and opportunities presented by integrating Distributed Energy Resources (DERs) and advanced metering infrastructure. The need for robust digital privacy measures has become crucial as the energy industry evolves towards a more decentralized, digitalized, and decarbonized future. This study delves into four cutting-edge privacy-preserving technologies—Homomorphic Encryption (HE), Secure Multiparty Computation (SMPC), Differential Privacy (DP), and Federated Learning (FL)—each offering unique solutions to safeguard consumer data by increasing digital connectivity and data exchange. Through a detailed examination of these methods, the study explains how each technology operates, its applications within the energy sector, and the specific privacy challenges it addresses. Homomorphic Encryption allows for secure computations on encrypted data, enabling data analysis without compromising privacy. Secure Multiparty Computation enables collaborative data analysis across different entities while protecting the confidentiality of the inputs. Differential Privacy introduces randomness into the assembled data set, preventing the identification of individual records in statistical databases. Lastly, Federated Learning offers a paradigm shift in data analysis, where machine learning models are trained at the edge, minimizing the centralization of sensitive data. The research underscores the significance of implementing these privacy-enhancing technologies to comply with strict data protection regulations, foster consumer trust, and enhance the security of the energy infrastructure. By providing a comprehensive overview of these methodologies and their practical implications for the energy sector, this study aims to contribute to the ongoing discourse on digital privacy, offering insights into how the energy industry can navigate the complexities of data privacy in the digital age.

Cali, Umit↗

The exceedingly strong two-dimensional ferromagnetism in bi-atomic layer SrRuO 3 with a critical conduction transition

In recent years, few-layer or even monolayer ferromagnetic materials have drawn a great deal of attention due to the promising integration of two-dimensional (2D) magnets into next-generation spintronic devices. The SrRuO 3 monolayer is a rare example of stable 2D magnetism under ambient conditions, but only weak ferromagnetism or antiferromagnetism has been found. The bi-atomic layer SrRuO3 as another environmentally inert 2D magnetic system has been paid less attention heretofore. Here we study both the bi-atomic layer and monolayer SrRuO 3 in (SrRuO 3 ) n /(SrTiO 3 ) m (n = 1, 2) superlattices in which the SrTiO 3 serves as a non-magnetic and insulating space layer. Although the monolayer exhibits arguably weak ferromagnetism, we find that the bi-atomic layer exhibits exceedingly strong ferromagnetism with a T c of 125 K and a saturation magnetization of 1.2 μ B /Ru, demonstrated by both superconducting quantum interference device (SQUID) magnetometry and element-specific X-ray circular dichroism. Moreover, in the bi-atomic layer SrRuO 3 , we demonstrate that random fluctuations and orbital reconstructions inevitably occurring in the 2D limit are critical to the electrical transport, but are much less critical to the ferromagnetism. Our study demonstrates that the bi-atomic layer SrRuO 3 is an exceedingly strong 2D ferromagnetic oxide which has great potentials for applications of ultracompact spintronic devices.

36 MATERIALS SCIENCE↗

A Computational Framework to design 3D stiffness gradient acoustic metamaterials for impedance matching

Acoustic waves play a crucial role in various applications, including medical imaging, non-destructive testing, and sonar systems. One of the significant challenges in these applications is impedance matching, which is essential for minimizing reflections and maximizing the transfer of acoustic energy between different media. Acoustic metamaterials offer a promising solution to this challenge. In addition to impedance control, gradient stiffness can enhance structural efficiency and enable spatial control of wave propagation, making it a valuable feature in acoustic metamaterial design. In this pa- per, we present our developed computational method to design 3D stiffness gradient acoustic metamaterials for impedance matching. The key steps in our approach include generating initial designs using a periodic covariance function to provide unit cells that are both periodic on the boundaries and randomly formed inside the unit cell. Furthermore, we integrated manufacturing constraints into the design process, ensuring that the structures are interconnected for fabrication. We propose two computational optimization algorithms: GenUnit, based on a non-dominated sorting genetic algorithm (NSGA-II), and MLMatch, which leverages differentiable machine learning. The two approaches are not separate contributions but complementary com- ponents of a unified framework. GenUnit requires no training data and directly interfaces with physics-based simulations, making it highly accurate but slower for large-scale exploration. In contrast, MLMatch is data-hungry during training but, once trained, enables near-instantaneous inference and broad design-space coverage. Together, they form a hybrid strategy: ML- Match rapidly explores the global design space, and GenUnit provides local refinement with high-fidelity accuracy. This balance between training cost, inference time, and precision is the motivation for including both methods in the same study. We applied this dual-algorithm framework to generate two metallic-based metamaterial designs that match the acoustic impedance of water while exhibiting a controlled gradient in stiffness (from stiff to soft). The stiffness gradient is particularly advantageous in applications where one side of the structure must interface with soft or sensitive surfaces, such as human tissue or delicate components. Here, this work paves the way for improved materials in various acoustic applications, particularly in ultrasound devices, by providing better impedance.

Metamaterial↗

Excitons in fractionally filled moiré superlattices

Long-range Coulomb forces give rise to correlated insulating states when charge particles populate a moiré superlattice at certain fractional filling factors. Such behavior is characterized by a broken translation symmetry wherein particles spontaneously form a Wigner crystal. Here, focusing on the experimental findings of Xu et al. [Nature (London) 587, 214 (2020)], we present a theory that captures the correlated insulating state of a fractionally filled moiré superlattice through the energy shift and change in oscillator strength of the exciton absorption resonance. The theory shows that the experimental findings can only be supported if the electrons reside in a charge-ordered state (i.e., electrons are not randomly distributed among the sites of the moiré superlattice). Furthermore, we explain why the energy shifts of exciton resonances are qualitatively different in cases where the superlattice is nearly empty compared with a superlattice whose sites are doubly occupied.

excitons↗

GraphAlign: Graph-Enabled Machine Learning for Seismic Event Filtering

This report summarizes results from a 2 year effort to improve the current automated seismic event processing system by leveraging machine learning models that can operated over the inherent graph data structure of a seismic sensor network. Specifically, the GraphAlign project seeks to utilize prior information on which stations are more likely to detect signals originating from particular geographic regions to inform event filtering. To date, the GraphAlign team has developed a Graphical Neural Network (GNN) model to filter out false events generated by the Global Associator (GA) algorithm. The algorithm operates directly on waveform data that has been associated to an event by building a variable sized graph of station waveforms nodes with edge relations to an event location node. This builds off of previous work where random forest models were used to do the same task using hand crafted features. The GNN model performance was analyzed using an 8 week IMS/IDC dataset, and it was demonstrated that the GNN outperforms the random forest baseline. We provide additional error analysis of which events the GNN model performs well and poorly against concluded by future directions for improvements.

58 GEOSCIENCES↗

The coupled SYK model at finite temperature

Sachdev-Ye-Kitaev (SYK) model, which describes N randomly interacting Majorana fermions in 0+1 dimension, is found to be an solvable UV-complete toy model for holographic duality in nearly AdS 2 dilaton gravity. Ref. [1] proposed a modified model by coupling two identical SYK models, which at low-energy limit is dual to a global AdS 2 geometry. This geometry is an “eternal wormhole" because the two boundaries are causally connected. Increasing the temperature drives a Hawking-Page like transition from the eternal wormhole geometry to two disconnected black holes with coupled matter field. To gain more understanding of the coupled SYK model, in this work, we study the finite temperature spectral function of this system by numerical solving the Schwinger-Dyson equation in real-time. We find in the low-temperature phase the system is well described by weakly interacting fermions with renormalized single-particle gap, while in the high temperature phase the system is strongly interacting and the single-particle peaks merge. We also study the q dependence of the spectral function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Machine-learning-aided cognitive reconfiguration for flexible-bandwidth HPC and data center networks [Invited]

This paper proposes a machine-learning (ML)-aided cognitive approach for effective bandwidth reconfiguration in optically interconnected datacenter/high-performance computing (HPC) systems. The proposed approach relies on a Hyper-X-like architecture augmented with flexible-bandwidth photonic interconnections at large scales using a hierarchical intra/inter-POD photonic switching layout. We first formulate the problem of the connectivity graph and routing scheme optimization as a mixed-integer linear programming model. A two-phase heuristic algorithm and a joint optimization approach are devised to solve the problem with low time complexity. Then, we propose an ML-based end-to-end performance estimator design to assist the network control plane with intelligent decision making for bandwidth reconfiguration. Numerical simulations using traffic distribution profiles extracted from HPC applications traces as well as random traffic matrices verify the accuracy performance of the ML design estimator ( < <#comment/> 9 % <#comment/> error) and demonstrate up to 5 × <#comment/> throughput gain from the proposed approach compared with the baseline Hyper-X network using fixed all-to-all intra/inter-portable data center interconnects.

Chen, Xiaoliang (ORCID:0000000278056237)↗

Estimating Spatial Distribution Impacts of Rooftops Solar PV on Dynamic Hosting Capacity Evaluation for a Real Distribution Feeder

This paper models the sensitivity of dynamic hosting capacity analysis to the spatial distribution of distributed PV scenarios, random, close and far from substation deployments, at various penetration levels on a real distribution feeder, in a quasi-static time series (QSTS) power flow simulation using high resolution time-aware metrics. These spatial distribution scenarios were chosen to capture a wide range of potential operation impacts of these deployments using a set of thermaland voltage-based metrics such as instantaneous violations and moving averages for both thermal and voltage constraints. This study uses actual load and PV data up-scaled from 1-hour time step to 1-minute resolution to fully characterize the interaction between the daily changes in load and PV output, and their impacts on distribution system operations.

distribution system↗

Quantum annealing-assisted lattice optimization

High Entropy Alloys (HEAs) have drawn great interest due to their exceptional properties compared to conventional materials. The configuration of HEA system is considered a key to their superior properties, but exhausting all possible configurations of atom coordinates and species to find the ground energy state is extremely challenging. In this work, we proposed a quantum annealing-assisted lattice optimization (QALO) algorithm, which is an active learning framework that integrates the Field-aware Factorization Machine (FFM) as the surrogate model for lattice energy prediction, Quantum Annealing (QA) as an optimizer and Machine Learning Potential (MLP) for ground truth energy calculation. By applying our algorithm to the NbMoTaW alloy, we reproduced the Nb depletion and W enrichment observed in bulk HEA. We found our optimized HEAs to have superior mechanical properties compared to the randomly generated alloy configurations. Our algorithm highlights the potential of quantum computing in materials design and discovery, laying a foundation for further exploring and optimizing structure-property relationships.

36 MATERIALS SCIENCE↗

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗

Assessment of BQ-9000 Biodiesel Properties for 2020 (CRADA CRD-15-593)

Biodiesel producers in the United States and Canada can voluntarily participate in the industry's BQ-9000 quality assurance program. This is the fourth in a series of reports documenting biodiesel quality from participating producers. Participants in the BQ-9000 program were requested to voluntarily provide data to a third-party team. This team anonymized and randomized the data prior to providing to the National Renewable Energy Laboratory (NREL) for analysis and reporting. The critical quality parameters analyzed are: sodium and potassium, calcium and magnesium, phosphorus, flash point and alcohol control, water and sediment, cloud point, acid number, free and total glycerin, monoglycerides, sulfur, oxidation stability, and cold soak filterability test (CSFT). The statistical analysis for these parameters is presented in Table ES-1.

09 BIOMASS FUELS↗

High Thermoelectric Performance in the New Cubic Semiconductor AgSnSbSe 3 by High-Entropy Engineering

We investigate the structural and physical properties of the AgSn m SbSe m+2 system with m = 1-20 (i.e., SnSe matrix and ~5-50% AgSbSe 2 ) from atomic, nano, and macro length scales. We find the 50:50 composition, with m = 1 (i.e., AgSnSbSe 3 ), forms a stable cation-disordered cubic rock-salt p-type semiconductor with a special multi-peak electronic valence band structure. AgSnSbSe 3 has an intrinsically low lattice thermal conductivity of ~0.47 W m -1 K -1 at 673 K owing to the synergy of cation disorder, phonon anharmonicity, low phonon velocity, and low-frequency optical modes. Furthermore, Te alloying on Se sites creates a quinary high-entropy NaCl-type solid solution AgSnSbSe 3-x Te x with randomly disordered cations and anions. The extra point defects and lattice dislocations lead to glass-like lattice thermal conductivities of ~0.32 W m -1 K -1 at 723 K and higher hole carrier concentration than AgSnSbSe 3 . Concurrently, the Te alloying promotes greater convergence of the multiple valence band maxima in AgSnSbSe 1.5 Te 1.5 , the composition with the highest configurational entropy. Facilitated by these favorable modifications, we achieve a high average power factor of ~9.54 μW cm -1 K -2 (400-773 K), a peak thermoelectric figure of merit ZT of 1.14 at 723 K, and a high average ZT of ~1.0 over a wide temperature range of 400-773 K in AgSnSbSe 1.5 Te 1.5 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Performance characterization of freeform finished surfaces of potassium dihydrogen phosphate using fluid jet polishing with a nonaqueous slurry

Abstract Potassium dihydrogen phosphate (KDP) and its deuterated analog (DKDP) are unique nonlinear optical materials for high power laser systems. They are used widely for frequency conversion and polarization control by virtue of the ability to grow optical-quality crystals at apertures suitable for fusion-class laser systems. Existing methods for freeform figuring of KDP/DKDP optics do not produce surfaces with sufficient laser-induced–damage thresholds (LIDT’s) for operation in the ultraviolet portion of high-peak-power laser systems. In this work, we investigate fluid jet polishing (FJP) using a nonaqueous slurry as a sub-aperture finishing method for producing freeform KDP surfaces. This method was used to selectively polish surface areas to different depths on the same substrate with removals ranging from 0.16 μm to 5.13 μm. The finished surfaces demonstrated a slight increase in roughness as the removal depth increased along with a small number of fracture pits. Laser damage testing with 351 nm, 1 ns pulses demonstrated excellent surface damage thresholds, with the highest values in areas devoid of fracture pits. This work demonstrates, for the first time, a method that enables fabrication of a waveplate that provides tailored polarization randomization that can be scaled to meter-sized optics. Furthermore, this method is based on FJP technology that incorporates a nonaqueous slurry specially designed for use with KDP. This novel nonaqueous FJP process can be also used for figuring other types of materials that exhibit similar challenging inherent properties such as softness, brittleness, water-solubility, and temperature sensitivity.

36 MATERIALS SCIENCE↗

Cross-reactive sarbecovirus antibodies induced by mosaic RBD nanoparticles

Broad immune responses are needed to mitigate viral evolution and escape. To induce antibodies against conserved receptor-binding domain (RBD) regions of SARS-like betacoronavirus (sarbecovirus) spike proteins that recognize SARS-CoV-2 variants of concern and zoonotic sarbecoviruses, we developed mosaic-8b RBD nanoparticles presenting eight sarbecovirus RBDs arranged randomly on a 60-mer nanoparticle. Mosaic-8b immunizations protected animals from challenges from viruses whose RBDs were matched or mismatched to those on nanoparticles. Here, we describe neutralizing mAbs isolated from mosaic-8b-immunized rabbits, some on par with Pemgarda, the only currently FDA-approved therapeutic mAb. Deep mutational scanning, in vitro selection of spike resistance mutations, and single-particle cryo-electron microscopy structures of spike–antibody complexes demonstrated targeting of conserved RBD epitopes. Rabbit mAbs included critical D-gene segment RBD-recognizing features in common with human anti-RBD mAbs, despite rabbit genomes lacking an equivalent human D-gene segment, thus demonstrating that the immune systems of humans and other mammals can utilize different antibody gene segments to arrive at similar modes of antigen recognition. These results suggest that animal models can be used to elicit anti-RBD mAbs with similar properties to those raised in humans, which can then be humanized for therapeutic use, and that mosaic RBD nanoparticle immunization coupled with multiplexed screening represents an efficient way to generate and select broadly cross-reactive therapeutic pan-sarbecovirus and pan-SARS-CoV-2 variant mAbs.

Science & Technology - Other Topics↗