Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “operational analysis”

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

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

At least 793 records · Page 44

Evaluation of Amorphous Oxide Coatings for High-Voltage Li-Ion Battery Applications Using a First-Principles Framework

Cathode surface coatings are widely used industrially as a means to suppress degradation and improve electrochemical performance of lithium-ion batteries. However, developing an optimal coating is challenging, as different coating materials may enhance one aspect of performance while hindering another. Furthermore, to elucidate the fundamental thermodynamic and transport properties of amorphous cathode coating materials, here, we present a framework for calculating and analyzing the Li + and O 2- transport and the stability against delithiation in such materials. Our framework includes systematic workflows of ab-initio molecular dynamics calculations to obtain amorphous structures and diffusion trajectories coupled with an analysis of critical changes of the active-ion local environment during diffusion. Based on these data, we provide an estimate of room-temperature diffusivities, including statistical error bars, and the evaluation of the coating suitability in terms of its ability to facilitate Li + transport while blocking O 2- transport. Finally, we add the thermodynamic stability analysis of the coating chemistry within the operating voltage of common Li-ion cathodes. We apply this framework to two commonly used amorphous coating materials, Al 2 O 3 and ZnO. We find that (1) in general, a higher Li + content increases both Li + and O 2- diffusivities in both Al 2 O 3 and ZnO. Also, Li + and O 2- diffuse much faster in ZnO than in Al 2 O 3 . (2) However, neither Al 2 O 3 nor ZnO is expected to retain a significant concentration of Li + at high charge. (3) ZnO performs much more poorly in terms of O 2- blocking, and hence, Al 2 O 3 is preferred for high-voltage cathode applications. These results will help to quantitatively evaluate amorphous materials, such as metal oxides and fluorides, for different performance metrics and facilitate the development of optimal cathode coatings.

25 ENERGY STORAGE↗

Data transfers for full core heterogeneous reactor high- fidelity multiphysics studies

Multiphysics simulations for nuclear reactor analysis are usually performed by resorting to operator splitting and fixed point iterations between single-physics solvers. This enables the separate solution of each physics, such as neutronics, fuel performance, and thermal hydraulics, on meshes tailored to the requirements of the respective numerical discretizations of the equations. As the equations are coupled, several fields must be transferred between single-physics solves. Projecting fields between meshes while preserving order of accuracy, conservation properties, and mapping non-overlapping geometries is a complex endeavor. This conference paper will present the transfers as implemented in MOOSE, which can handle arbitrary meshes, arbitrary mappings, conservation of integral quantities, and are made to scale with distributed simulations on both ends of the transfers. Their adequacy for advanced nuclear reactor multiphysics coupling is shown through examples and numerical studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Results from a synthetic model of the ITER XRCS-Core diagnostic based on high-fidelity x-ray ray tracing

A high-fidelity synthetic diagnostic has been developed for the ITER core x-ray crystal spectrometer diagnostic based on x-ray ray tracing. This synthetic diagnostic has been used to model expected performance of the diagnostic, to aid in diagnostic design, and to develop engineering tolerances. The synthetic model is based on x-ray ray tracing using the recently developed xicsrt ray tracing code and includes a fully three-dimensional representation of the diagnostic based on the computer aided design. The modeled components are: plasma geometry and emission profiles, highly oriented pyrolytic graphite pre-reflectors, spherically bent crystals, and pixelated x-ray detectors. Plasma emission profiles have been calculated for Xe 44+ , Xe 47+ , and Xe 51+ , based on an ITER operational scenario available through the Integrated Modelling & Analysis Suite database, and modeled within the ray tracing code as a volumetric x-ray source; the shape of the plasma source is determined by equilibrium geometry and an appropriate wavelength distribution to match the expected ion temperature profile. All individual components of the x-ray optical system have been modeled with high-fidelity producing a synthetic detector image that is expected to closely match what will be seen in the final as-built system. Particular care is taken to maintain preservation of photon statistics throughout the ray tracing allowing for quantitative estimates of diagnostic performance.

47 OTHER INSTRUMENTATION↗

Transfer learning for analysis of collective and non-collective Thomson scattering spectra

Thomson scattering (TS) diagnostics provide reliable, minimally perturbative measurements of fundamental plasma parameters, such as electron density (⁠n e ) and electron temperature (⁠T e ⁠). Deep neural networks can provide accurate estimates of ⁠n e and T e when conventional fitting algorithms may fail, such as when TS spectra are dominated by noise, or when fast analysis is required for real-time operation. Although deep neural networks typically require large training sets, transfer learning can improve model performance on a target task with limited data by leveraging pre-trained models from related source tasks, where select hidden layers are further trained using target data. We present five architecturally diverse deep neural networks, pre-trained on synthetic TS data and adapted for experimentally measured TS data, to evaluate the efficacy of transfer learning in estimating n e and T e in both the collective and non-collective scattering regimes. We evaluate errors in n e and T e estimates as a function of training set size for models trained with and without transfer learning, and we observe decreases in model error from transfer learning when the training set contains ≲ 200 experimentally measured spectra.

Artificial neural networks↗

Genetic switching by the Lac repressor is based on two-state Monod–Wyman–Changeux allostery

High-resolution NMR spectroscopy enabled us to characterize allosteric transitions between various functional states of the dimeric Escherichia coli Lac repressor. In the absence of ligands, the dimer exists in a dynamic equilibrium between DNA-bound and inducer-bound conformations. Binding of either effector shifts this equilibrium toward either bound state. Analysis of the ternary complex between repressor, operator DNA, and inducer shows how adding the inducer results in allosteric changes that disrupt the interdomain contacts between the inducer binding and DNA binding domains and how this in turn leads to destabilization of the hinge helices and release of the Lac repressor from the operator. Based on our data, the allosteric mechanism of the induction process is in full agreement with the well-known Monod–Wyman–Changeux model.

59 BASIC BIOLOGICAL SCIENCES↗

Results from the Alfvén Eigenmode Active Diagnostic during the 2019-2020 JET deuterium campaign

This paper presents results of extensive analysis of mode excitation observed during the operation of the Alfvén Eigenmode Active Diagnostic (AEAD) in the JET tokamak during the 2019-2020 deuterium campaign. Here, six of eight toroidally spaced antennas, each with independent power and phasing, were successful in actively exciting stable MHD modes in 479 plasmas. In total, magnetic resonances were detected with up to fourteen fast magnetic probes. In this work, we present the calculations of resonant frequencies $f_0$, damping rates γ < 0, and toroidal mode numbers $\textit{n}$, spanning the parameter ranges $f_0$ ≈ 30–250 kHz, –γ ≈ 0–13 kHz, and $\vert n \vert \leq {30}$. In general, good agreement is seen between the resonant and the calculated toroidal Alfvén Eigenmode frequencies, and between the toroidal mode numbers applied by the AEAD and estimated of the excited resonances. We note several trends in the database: the probability of resonance detection decreases with plasma current and external heating power; the normalized damping rate increases with edge safety factor but decreases with external heating. These results provide key information to prepare future experimental campaigns and to better understand the physics of excitation and damping of Alfvén Eigenmodes in the presence of alpha particles during the upcoming DT campaign, thereby extrapolating with confidence to future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exploring the energy landscape of RBMs: reciprocal space insights into bosons, hierarchical learning and symmetry breaking

Deep generative models have become ubiquitous due to their ability to learn and sample from complex distributions. Despite the proliferation of various frameworks, the relationships among these models remain largely unexplored, a gap that hinders the development of a unified theory of AI learning. In this work, we address two central challenges: clarifying the connections between different deep generative models and deepening our understanding of their learning mechanisms. We focus on Restricted Boltzmann Machines (RBMs), a class of generative models known for their universal approximation capabilities for discrete distributions. By introducing a reciprocal space formulation for RBMs, we reveal a connection between these models, diffusion processes, and systems of coupled bosons. Our analysis shows that at initialization, the RBM operates at a saddle point, where the local curvature is determined by the singular values of the weight matrix, whose distribution follows the Marc̆enko-Pastur law and exhibits rotational symmetry. During training, this rotational symmetry is broken due to hierarchical learning, where different degrees of freedom progressively capture features at multiple levels of abstraction. This leads to a symmetry breaking in the energy landscape, reminiscent of Landau’s theory. This symmetry breaking in the energy landscape is characterized by the singular values and the weight matrix eigenvector matrix. We derive the corresponding free energy in a mean-field approximation. We show that in the limit of infinite size RBM, the reciprocal variables are Gaussian distributed. Our findings indicate that in this regime, there will be some modes for which the diffusion process will not converge to the Boltzmann distribution. To illustrate our results, we trained replicas of RBMs with different hidden layer sizes using the MNIST dataset. Our findings not only bridge the gap between disparate generative frameworks but also shed light on the fundamental processes underpinning learning in deep generative models.

97 MATHEMATICS AND COMPUTING↗

Experimental Limit on Nonlinear State-Dependent Terms in Quantum Theory

Linear time evolution is one of the fundamental postulates of quantum theory. Past theoretical attempts to introduce nonlinearity into quantum evolution have violated causality. However, a recent theory has introduced nonlinear state-dependent terms in quantum field theory, preserving causality [D. E. Kaplan and S. Rajendran, Phys. Rev. D 105, 055002 (2022)]. We report the results of an experiment that searches for such terms. Our approach, inspired by the Everett many-worlds interpretation of quantum theory, correlates a binary macroscopic classical voltage with the outcome of a projective measurement of a quantum bit, prepared in a coherent superposition state. Measurement results are recorded in a bit string, which is used to control a voltage switch. Presence of a nonzero voltage reading in cases of no applied voltage is the experimental signature of a nonlinear state-dependent shift of the electromagnetic field operator. We implement blinded measurement and data analysis with three control bit strings. Control of systematic effects is realized by producing one of the control bit strings with a classical random-bit generator. The other two bit strings are generated by measurements performed on a superconducting qubit in an IBM Quantum processor and on a N 15 nuclear spin in a nitrogen-vacancy center in diamond. Our measurements find no evidence for electromagnetic quantum state-dependent nonlinearity. We set a bound on the parameter that quantifies this nonlinearity | ε γ | < 4.7 × 10 - 11 , at 90% confidence level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

POWTEX visits POWGEN

The high-intensity time-of-flight (TOF) neutron diffractometer POWTEX for powder and texture analysis is currently being built prior to operation in the eastern guide hall of the research reactor FRM II at Garching close to Munich, Germany. Because of the world-wide 3 He crisis in 2009, the authors promptly initiated the development of 3 He-free detector alternatives that are tailor-made for the requirements of large-area diffractometers. Herein is reported the 2017 enterprise to operate one mounting unit of the final POWTEX detector on the neutron powder diffractometer POWGEN at the Spallation Neutron Source located at Oak Ridge National Laboratory, USA. As a result, presented here are the first angular- and wavelength-dependent data from the POWTEX detector, unfortunately damaged by a 50 g shock but still operating, as well as the efforts made both to characterize the transport damage and to successfully recalibrate the voxel positions in order to yield nonetheless reliable measurements. Also described is the current data reduction process using the PowderReduceP2D algorithm implemented in Mantid [Arnold et al. (2014). Nucl. Instrum. Methods Phys. Res. A , 764 , 156–166]. The final part of the data treatment chain, namely a novel multi-dimensional refinement using a modified version of the GSAS-II software suite [Toby & Von Dreele (2013). J. Appl. Cryst. 46 , 544–549], is compared with a standard data treatment of the same event data conventionally reduced as TOF diffraction patterns and refined with the unmodified version of GSAS-II . This involves both determining the instrumental resolution parameters using POWGEN's powdered diamond standard sample and the refinement of a friendly-user sample, BaZn(NCN) 2 . Although each structural parameter on its own looks similar upon comparing the conventional (1D) and multi-dimensional (2D) treatments, also in terms of precision, a closer view shows small but possibly significant differences. For example, the somewhat suspicious proximity of the a and b lattice parameters of BaZn(NCN) 2 crystallizing in Pbca as resulting from the 1D refinement (0.008 Å) is five times less pronounced in the 2D refinement (0.038 Å). Similar features are found when comparing bond lengths and bond angles, e.g. the two N—C—N units are less differently bent in the 1D results (173 and 175°) than in the 2D results (167 and 173°). The results are of importance not only for POWTEX but also for other neutron TOF diffractometers with large-area detectors, like POWGEN at the SNS or the future DREAM beamline at the European Spallation Source.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Online Output-Based Inertia Estimation of Modern Power Systems

The overall inertia of modern power systems is drastically impacted by the penetration level of converter-based resources (CBRs). The intermittent output of CBRs and their dependence on the weather conditions affirm the need for a realtime inertia estimation approach. This can help the transmission system operators(TSOs) to be continuously aware of the changes in the system inertia and to take suitable control actions. This paper is a step forward to a totally data-driven framework for power system inertia estimation using the available measurements from the phasor measurement units (PMUs). This inertia estimation is done based on the spectral analysis of the so-called Koopman linear operator. Koopman operator considers the relationship between the system's inertia and its dynamical response to slight disturbances in the load. Simulation results show that Koopman operator can approximate the system behavior in a reconstructed linear space hence inertia constant is estimated accurately.

data-driven↗

Neural Lyapunov Control for Power System Transient Stability: A Deep Learning-Based Approach

We report that power system control and transient stability analysis play essential roles in secure system operation. Control of power systems typically involves highly nonlinear and complex dynamics. Most of the existing works address such problems with additional assumptions in system dynamics, leading to a requirement for a complete and general solution. This paper, therefore, proposes a novel control framework for various power system control and stability problems leveraging a learning-based approach. The proposed framework includes a two-module structure that iteratively and jointly learns the candidate Lyapunov function and control law via deep neural networks in a learning module. Meanwhile, it guides the learning procedure towards valid results satisfying Lyapunov conditions in a falsification module. The introduced termination criteria ensure provable system stability. This control framework is verified through several studies handling different types of power system control problems. The results show that the proposed framework is generalizable and can simplify the control design for complex power systems with the stability guarantee and enlarged region of attraction.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Toward Understanding the Differences between Mesoscale and Large-Eddy Simulations of Tropical Cyclones

In this work, we investigate the ability of mesoscale and large-eddy simulation (LES) model configurations to predict the mean wind speed profile within the boundary layer of tropical cyclones (TCs). To this end, we perform idealized simulations of five hypothetical intense storms ranging from categories 1 to 5 on the Saffir–Simpson scale and extract time-averaged quantities near the eyewall region. We compare the model-generated data against mean wind speed profiles compiled from dropsondes launched from reconnaissance aircraft operating in the North Atlantic basin. Our analysis shows that mesoscale- and LES-generated mean wind fields display important differences in the boundary layer, including the magnitude of shear as well as the height where their low-level wind speed maxima are located. In addition, a comparison between the two model configurations with the dropsonde data shows that both modeling approaches are unable to capture the typical structure of mean winds in the lower part of the TC boundary layer (10–500 m), calling into question the use of simulations of near-axisymmetric storms for investigating the wind structure of past events. To better understand these differences, we conduct a momentum-budget analysis and show that modeled turbulent fluxes are underestimated in the mesoscale boundary layer parameterization compared to the LES model. Based on the analysis of the horizontal turbulent fluxes and their potential impact on mean flow quantities, a TC-specific boundary layer parameterization may be needed.

17 WIND ENERGY↗

High Throughput In-Line Coating Metrology Development for Solid Oxide Fuel Cell Manufacturing

Coatings play key roles in solid oxide fuel cell (SOFC) stack durability. For example, diffusion barrier coatings on Cr-containing interconnect and balance of plant (BOP) components protect electrodes from Cr poisoning over the long operational lifetimes (>10,000 hours) of the fuel cell stack. Common defects in coatings, such as cracks, pinholes, and porosity, result in a failure to protect the electrodes, resulting in shorter operational lifetime and thus higher cost. It is very unlikely, even in the best coating process, that all these defects can be mitigated, hence identifying critical defects in parts, and removing defective parts from production before they can damage the stack, becomes paramount. Furthermore, these quality control techniques must be operational in the production/assembly line (in-line), i.e., high throughput and non-destructive, and cost effective. Redox Power Systems, LLC (Redox) together with the National Renewable Energy Laboratory (NREL) developed much needed high throughput, in-line metrology techniques for protective coatings. The overall goal of the project is to lower cost while increasing robustness, reliability, and endurance of SOFC stacks. To accomplish this, we had several objectives, including: to identify key coating and substrate defects that lead to coating failure through the use of detailed characterization methods (e.g., microscopy, XRD, EDS, electrochemistry); to assess capabilities of in-line metrology techniques, e.g., optical profilometry (Redox) and thermography (NREL), to probe these defects, or evidence thereof; demonstrate long-term performance of “defect-free” protective coatings, as identified by in-line metrology, in solid oxide fuel cell (SOFC) stack operation. In the first part of this project, the ability to identify key defects expected to lead to coating and SOFC degradation using in-line metrology tools were evaluated. Coated interconnect samples with controlled defect types and populations were tested under conditions similar to SOFC operation, followed by detailed post-test analysis to reveal the defects responsible for observed degradation. In the second part of the project, the optimal in-line metrology techniques and methodologies were used to map the defect distribution in full-size interconnects with critical defects intentionally allowed to exist in some cases. These interconnects underwent SOFC testing for extended periods (up to ~3,000 hours) followed by post-test analysis to evaluate the effectiveness of in-line metrology techniques in mitigating MCO coating related degradation. Key accomplishments in this project included the following: Demonstrated ASR of < 0.05 ohm-cm 2 at 650 °C for 1,000 hours with low defect (determined by in-line metrology) interconnect samples (average ASR=37 milliohms-cm 2 after over 1,000 hours). Demonstrated that low defect coatings on interconnects (as screened using in-line metrology) have low volatilization of chromium at ~650 °C for 1,000 hours as detected using Cr-getter material (< 5 at% increase above baseline); 1022 hour duration tests under humidified, elevated temperature (750 °C rather than 650 °C) compared a base case against different coating thicknesses. Demonstrated capability to identify initial key defects of interest with in-line metrology techniques using up to 8 cm by 10 cm having coatings with and without intentional defects of interest using thermal imaging and optical profilometry. Correlated key defects identified using metrology techniques with observed coating performance (e.g., ASR and Cr volatility). Conducted several 4 cm by 4 cm cell tests using MCO-interconnects that were pre-screened using some of the metrology techniques developed in the project (e.g., optical profilometry). An analysis of ASR measurements were able to show that defect-free coatings resulted in the anticipated performance in the cell tests.

01 COAL, LIGNITE, AND PEAT↗

The Energizer Bunny: Dual-Use Photovoltaic and Pasture-Raised Rabbit Farms (Final Report)

This case study project lays the groundwork for the broader realization of agrivoltaics by considering the techno-economic, environmental, and social aspects of system development and assessing a novel rabbit-based concept. Research on the social and political barriers to agrivoltaic development is sparse, which emphasizes the importance of this project as it provides key insight into the consequential socio-political challenges that may hinder the diffusion of agrivoltaic technology. Additionally, the rabbit agrivoltaic concept is novel and has not been applied in experimental or commercial applications to date. This project advances our understanding of the viability of this co-location scheme through empirical assessment of integrating grazing rabbits with solar. A total of five manuscripts (3 published in peer-reviewed journals, 2 under review) have resulted from this project. The techno-economic component of the project included a) a pilot test study investigating the feasibility of the rabbit agrivoltaic concept and b) an economic analysis of the advantages available to solar operators who leverage grazing animals for lawn maintenance and a dual-revenue stream. Our calculations indicate that the co-location of solar and rabbit farms is a viable form of agrivoltaics, increasing overall site revenue by 2.5%-24.0% above projected electricity revenue. The environmental component of the project entailed an ISO 14040 and 14044 compliant life cycle assessment to quantify the impacts of the rabbit agrivoltaic concept. The results indicate pasture-based agrivoltaic systems provide a dual synergy that produces 69.3% less emissions and demands 82.9% less fossil energy per functional unit of cumulative MWh output of electricity and cumulative kg of meat compared to non-integrated production. The social component of the project included conducting two separate interview studies with agriculture and solar industry professionals and the distribution of a public survey. The interview findings identify the key socio-political opportunities for agrivoltaics include: the retention of agricultural land and rural interests, and increased local level acceptance of solar development. The key barriers to agrivoltaics include: ensuring long term agricultural productivity is not compromised, and subnational zoning strategies. The survey results indicate that support for local solar projects increase by 10% when an agricultural function is incorporated. Survey respondents prefer agrivoltaic projects that are a) designed to provide economic opportunities for farmers and the local community b) located on private property or existing agricultural land, and c) do not threaten local interests. Cumulatively, this project can inform agrivoltaic decision making, solar development practices, rabbit-based applications, land use management, and policy making in a way that supports the furtherance of the renewable energy transition, conserves arable land, and utilizes innovative solar photovoltaic technologies.

14 SOLAR ENERGY↗

Comprehensive Assessment of On-And Off-Board Vehicle-To-Grid Technology Performance and Impacts on Battery and the Grid

This project defined, developed, and tested the open standards-based requirements, engineering design and system integration of on and off-vehicle Vehicle-to-Grid technologies both on the grid and on the PEV, to assess their effectiveness in implementing the use cases that create value for the grid, which in turn, can be passed on to the ratepayers and the PEV owners. The project comprised of three key categories of work and four key sub-areas of work. The three key categories were: On-Vehicle V2G, Off-Vehicle V2G, and impact of V2G operation on the battery capacity degradation. The four key sub-areas of work under each of these categories were: defining requirements, implementing the technology, performing system integration testing, and estimating the value of each of the grid services under specific assumptions, on a per-vehicle and on an aggregated basis. The technical implementation demonstrated the validity of standards-based approach to grid integration that led to interoperability, in addition to assessing the operational details with grid constraints imposed at the transformer level. The technology implementation also identified the subtle gaps in the standards definition to close the feedback loop on the accuracy of the standards as written, which are being implemented as revisions to the appropriate protocols. The battery impacts assessment developed a test cycle that the batteries from a real PHEV were subjected to for both the mobility-only and mobility with V2G specific energy cycling. While the battery impact results show a clear incremental degradation in response to additional throughput for V2G application purposes, the amount of additional kWh that can be made available without exceeding the end-of-life capacity definition of the battery at its 10-year warranty period was found to be significant. Given that the PEV battery being exercised during testing was a PHEV battery (smaller energy capacity), and 44% of its usable energy (i.e., SOC) was utilized daily for V2G applications, the EVs with much larger batteries will fare much better, either in terms of extracting more value in terms of incremental kWh for V2G, or in terms of relative incremental degradation of the capacity through the battery lifespan. Further, this battery was the second-generation Lithium-Ion chemistry circa 2015, which is at least two generations ago. Continual improvements in electrochemistry and manufacturing techniques as well as on-board capacity will enable the future batteries to provide even more energy for non-mobility-related services. Finally, the valuation of V2G energy services at the premise, distribution, and the ISO/market level was carried out for both the on-vehicle and off-vehicle cases, to understand both the GHG mitigation and operational efficiency improvements in quantitative terms. The analysis indicates that between $\$400$ and $\$1400$ per year of total value to the grid can be realized. When netted of costs to implement V2G either on or off-vehicle, the net value is available to be shared between the ratepayers and the EV owners. Over the life of the EV, therefore, this value becomes significant and can easily provide both ratepayer and EV owner benefits more than the cost of the equipment, either on or off-board, in addition to providing GHG mitigation benefits both on the grid and on the mobility sides. For these benefits to accrue, the key learnings from both the on- and off-vehicle V2G parts of the project indicate the need for at-scale demonstrations through involvement of real EV owners to validate both the technical feasibility, interoperability, as well as real grid benefits by performing extensive data collection and analysis. This will help inform the grid planners, the program designs, and tariff designers, as well as automotive and equipment manufacturers how best to create their products to maximize grid, ratepayer, and EV owner value.

30 DIRECT ENERGY CONVERSION↗

Experimental overview about LANL contributions to the Electron-Ion Collider

This report is a request from the LANL team to design and build the Forward Silicon Tracker (FST) for the Electron-Ion Collider (EIC). The report highlights past accomplishments of LANL and seeks DOE support for continued work for the EIC, including silicon detector R&D; silicon tracker construction and operation; heavy flavor physics development and data analysis.

43 PARTICLE ACCELERATORS↗