Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel simulation”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Physicochemical control of solvation and molecular assembly of charged amphiphilic oligomers at air-aqueous interfaces

Hypothesis : Understanding the rules that control the assembly of nanostructured soft materials at interfaces is central to many applications. We hypothesize that electrolytes can be used to alter the hydration shell of amphiphilic oligomers at the air-aqueous interface of Langmuir films, thereby providing a means to control the formation of emergent nanostructures. Experiments : Three representative salts – (NaF, NaCl, NaSCN) were studied for mediating the self-assembly of oligodimethylsiloxane methylimidazolium (ODMS-MIM + ) amphiphiles in Langmuir films. The effects of the different salts on the nanostructure assembly of these films were probed using vibrational sum frequency generation (SFG) spectroscopy and Langmuir trough techniques. Experimental data were supported by atomistic molecular dynamic simulations. Findings : Langmuir trough surface pressure – area isotherms suggested a surprising effect on oligomer assembly, whereby the presence of anions affects the stability of the interfacial layer irrespective of their surface propensities. In contrast, SFG results implied a strong anion effect that parallels the surface activity of anions. These seemingly contradictory trends are explained by anion driven tail dehydration resulting in increasingly heterogeneous systems with entangled ODMS tails and appreciable anion penetration into the complex interfacial layer comprised of headgroups, tails, and interfacial water molecules. These findings provide physical and chemical insight for tuning a wide range of interfacial assemblies.

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

Large Scale MD to Predict Epitope Regions in HIV Env [Slides]

Highly dense carbohydrates located on the surface of the HIV Env protein play a key role in immune evasion. Such evolutionary adaptation hampers any attempt to obtain a full mechanistic understanding of the role played by the glycans in protecting the virus against an effective immune response. Moreover, and due to their chemical variability, an accurate molecular understanding of the so called “glycan shield” is still limited by the lack of effective resolution of state-of-the-art experimental technics. Here, we have used extensive computational modelling in order to fill this gap, addressing the presence of a large glycan variability as observed experimentally. Based on an automated pipeline, we were able to assemble, set-up and simulate via Molecular dynamics hundreds of different glycosylated Env variants at nearly atomic resolution, recapitulating the glycosylation distributions observed experimentally. Results from these simulations were subjected to machine learning and very accurate prediction of simulation derived glycan shielding areas of each glycan as a function of static sequence features. Such predictive models of per-glycan shielding, incorporating both glycan dynamics and heterogeneity, were used to develop a novel sequence-based glycan shield mapping strategy. Parallel to these studies, we also developed an accurate machine learning approach to predict glycan heterogeneity data using sequence features and found good prediction accuracy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Fully Kinetic Simulations of Proton-beam-driven Instabilities from Parker Solar Probe Observations

The expanding solar wind plasma ubiquitously exhibits anisotropic nonthermal particle velocity distributions. Typically, proton velocity distribution functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by the Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam population that drift with respect to each other. This configuration triggers a proton beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing perpendicular heating at the expense of the parallel temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Machine Learning-Driven Conservative-to-Primitive Conversion in Hybrid Piecewise Polytropic and Tabulated Equations of State

We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.

conservative-to-primitive conversion↗
DOE OSTI2020

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Membrane-mediated dimerization of spherocylindrical nanoparticles

Here we present a numerical investigation of the modes of adhesion and endocytosis of two spherocylindrical nanoparticles (SCNPs) on planar and tensionless lipid membranes, using systematic molecular dynamics simulations of an implicit-solvent model, with varying values of the SCNPs' adhesion strength and dimensions. We found that at weak values of the adhesion energy per unit of area, ξ, the SCNPs are monomeric and adhere to the membrane in the parallel mode. As ξ is slightly increased, the SCNPs dimerize into wedged dimers, with an obtuse angle between their major axes that decreases with increasing ξ. However, as ξ is further increased, we found that the final adhesion state of the two SCNPs is strongly affected by the initial distance, d 0 , between their centers of mass, upon their adhesion. Namely, the SCNPs dimerize into wedged dimers, with an acute angle between their major axes, if d 0 is relatively small. However, for relatively high d 0 , they adhere individually to the membrane in the monomeric normal mode. For even higher values of ξ and small values of d 0 , the SCNPs cluster into tubular dimers. However, they remain monomeric if d 0 is high. Finally, the SCNPs endocytose either as a tubular dimer, if d 0 is low or as monomers for large d 0 , with the onset value of ξ of dimeric endocytosis being lower than that of monomeric endocytosis. Dimeric endocytosis requires that the SCNPs adhere simultaneously at nearby locations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Correlated fission fragment spin dynamics

Here, this study explores the role of nucleon exchange for the generation of the fission fragment angular momenta. For a number of typical fission cases, samples of 10 4 shape evolutions are generated by Langevin simulation and, subsequently, for each such evolution, the nucleon exchange transport theory previously developed for damped nuclear reactions is used to obtain the development of the fragment spin-spin distribution within the Fokker-Planck transport framework. The characteristic evolution of both parallel and perpendicular spin components is discussed. A common feature is that the rotational modes fall out of equilibrium before scission when the temperature rises rapidly while the concurrent shrinking of the neck suppresses further exchange. A number of fission observables are extracted from the event ensembles: the distribution of the magnitude of the fragment spin and its orientation relative to the fission axis, as well as the correlation between the two spins and the distribution of their opening angle. The dependence of these observables on the mass asymmetry is also examined.

fission↗
DOE OSTI2022

Modeling mesoscale fission gas behavior in UO2 by directly coupling the phase field method to spatially resolved cluster dynamics

Abstract Fission gas release within uranium dioxide nuclear fuel occurs as gas atoms diffuse through grains and arrive at grain boundary (GB) bubbles; these GB bubbles grow and interconnect with grain edge bubbles; and grain edge tunnels grow and connect to free surfaces. In this study, a hybrid multi-scale/multi-physics simulation approach is presented to investigate these mechanisms of fission gas release at the mesoscale. In this approach, fission gas production, diffusion, clustering to form intragranular bubbles, and re-solution within grains are included using spatially resolved cluster dynamics in the Xolotl code. GB migration and intergranular bubble growth and coalescence are included using the phase field method in the MARMOT code. This hybrid model couples Xolotl to MARMOT using the MultiApp and Transfer systems in the MOOSE framework, with Xolotl passing the arrival rate of gas atoms at GBs and intergranular bubble surfaces to MARMOT and MARMOT passing evolved GBs and bubble surface positions to Xolotl. The coupled approach performs well on the two-dimensional simulations performed in this work, producing similar results to the standard phase field model when Xolotl does not include fission gas clustering or re-solution. The hybrid model performs well computationally, with a negligible cost of coupling Xolotl and MARMOT and good parallel scalability. The hybrid model predicts that intragranular fission gas clustering and bubble formation results in up to 70% of the fission gas being trapped within grains, causing the increase in the intergranular bubble fraction to slow by a factor of six. Re-solution has a small impact on the fission gas behavior at 1800 K but it has a much larger impact at 1000 K, resulting in a twenty-times increase in the concentration of single gas atoms within grains. Due to the low diffusion rate, this increase in mobile gas atoms only results in a small acceleration in the growth of the intergranular bubble fraction. Finally, the hybrid model accounts for migrating GBs sweeping up gas atoms. This results in faster intergranular bubble growth with smaller initial grain sizes, since the additional GB migration results in more immobile gas clusters reaching GBs.

Kim, Dong-Uk↗
DOE OSTI2025

Wettability of Two-Dimensional Carbon Allotropes from Molecular Simulations

Force-field Monte Carlo and Molecular Dynamics simulations are used to compare wetting behaviors of model carbon sheets mimicking neat graphene, its saturated derivative, graphane, and related planar allotropes penta-graphene, γ-graphyne, and ψ-graphene in contact with aqueous droplets or an aqueous film confined between parallel carbon sheets. Atomistic and area-integrated surface/water potentials are found to be essentially equivalent in capturing moderate differences between the wetting free energies of tested substrates. Despite notable differences in mechanical and electric properties of distinct allotropes, the predicted allotrope/water contact angles span a narrow window of weakly hydrophilic values. Contact angles in the range of 80 ± 10° indicate modest hydration repulsion incapable of competing with van der Waals attraction between carbon particles. Poor dispersibility in neat water is hence a common feature of studied materials.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Fermilab's controls development with virtual accelerator

Control Systems development is often the last thing considered when designing and building new equipment, e.g. a new detector or superconducting RF LINAC; however when the new equipment is installed, it is the first thing desired to be operational for testing. Due to frequent delays in building new equipment and project deadlines, control system development and testing is often curtailed. A way to alleviate this problem is to simulate the control system, though this will be challenging for complex systems.The Fermilab PIP-II (proton improvement plan - II) project is being constructed at Fermilab to deliver $800\,MeV$ protons of $>1\,MW$ beam power to replace the present LINAC for the remainder of the existing accelerator complex. The new LINAC consists of a warm front end (WFE), 23 superconducting RF cryomodules (of 5 types), and a beam transfer line (BTL) to the existing complex.The accelerator physics group has a parallel project to create a digital twin (DT) of the PIP-II accelerator. We have coupled the EPICS controls to this DT and are developing both the DT and EPICS software in parallel. This will allow us to develop the EPICS software framework, the HMIs, sequences, high level physics applications, and other services for use in a fully functional control system.This presentation will detail the work that we have performed to date and show demonstrations of controlling and monitoring the status of the accelerator, as well as future plans for this work.

Hanlet, Pierrick [Fermilab]↗
DOE OSTI2020

Determining Compression Characteristics of Honeycomb Material - 19674

The objective is to determine compression test characteristics of stainless steel honeycomb material to be able to represent honeycomb structures accurately in analytical models used to simulate hypothetical accident scenarios of shipping packages. Honeycomb is a material used primary in the aerospace industry due to its high strength-to-weight ratio. Because it doesn't have a shelf life and can withstand high heats, it is an excellent candidate for a structural material in package designs. Honeycomb has orthotropic material properties. The honeycomb currently being evaluated is made of metal ribbons spot welded together to form hexagon pattern between metal plates. The hexagon ribbons are brazed to the metal plates. The complexity of the honeycomb significantly lengthens the simulation time, which is further complicated when the brazing, welding, and imperfections of the material is considered. This means to affectively represent the honeycomb in simulation programs, such as Abaqus, the material needs to be approximated as a uniform orthotropic material. This requires material properties in each of the three directions, T, W, and L. The T direction is defined as the 'strong' direction, perpendicular to the honeycomb sheet. The L direction is parallel the ribbon and the W direction is perpendicular to the ribbon. Honeycomb has 3 different compression stages. Stage 1 is the initial compression where the honeycomb maintains its structural integrity and does not permanently deform from loads from normal operations. Stage 2 is when the initial buckling causes irreversible damage to the honeycomb. Stage 2 is the one we are most interested in because it absorbs the most energy from a hypothetical accident scenario. During Stage 2, the honeycomb fails layer-by-layer, indicating that the more layers, the longer material crushing is sustained, and thus the more energy that is absorbed. Stage 3 is final compression, similar to compressing solid metal. During stage 3 the stress increases with a diminishing rate of elongation and the honeycomb is completely failed where the plates between the honeycomb core sandwich the crushed honeycomb ribbon. The stress-strain graph bellow shows compression test of 3 different samples in the T direction. The red lines divide the three stages. The far left is Stage 1, the middle stage 2, and the right stage 3. The data collected is displayed in the table below. The expressions in the table represent the regression formula that represents each stage of the compression on a stress-strain graph. Stage 1 is assumed to intersect with the origin; data for stage 3 for the crush in W and L directions where unable to be gathered due to the nature of failure for those directions. Stage 1 and 2 are linear regressions while stage 3 is a degree 2, to best match the curve. Linear regression is used on stage 2 because ultimately that represents the energy absorbed. With more data, a degree n x 2 regression would be a more appropriate regression for stage 2, where n is the number of layers. The table below is the approximation of each scenario and stage. Each stage starts and ends at the intersection of the next stage. There where two main objectives for these tests: firstly to understand how honeycomb performs under extreme compression, and secondly to be able to numerically represent the honeycomb structure. Both of these where completed to varying degrees. It is important understood how the honeycomb fails. If it is crushed in the T direction it retains its integrity even after being crushed, however when crushed in the L and W crush direction, if it fails, it disintegrates, and loses all integrity. Also, the more layers the more time the material spends in stage 2. Failure is started by buckling, thus if there is any imperfection, the stress will not spike but transition straight into stage 2. The data collected and aggregated can be used for initial simulation of honeycomb used in packages. The initial testing has set the ground work for more data to be collected in order to verify results and to allow more confidence in the simulation results. This is only the initial data collected. The next steps is to continue to collect more data to verify results and to test more variants of honeycomb, with different brazing, layers, and shape.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nonvolatile Magnetoelectric Switching of Magnetic Tunnel Junctions with Dipole Interaction

The magnetoelectric effect is technologically appealing because of its ability to manipulate magnetism using an electric field rather than magnetic field or current, thus providing a promising solution for the development of energy-efficient spintronics. Although 180° magnetization switching is vital to spintronic devices, the achievement of 180° magnetization switching via magnetoelectric coupling is still a fundamental challenge. Here, in this work, voltage-driven full resistance switching of a magnetic tunnel junction (MTJ) with dipole interaction on a ferroelectric substrate through switchable parallel/antiparallel magnetization alignment is demonstrated. Parallel magnetization alignment along the y direction is obtained under a bias magnetic field. By rotating the magnetic easy axis via strain-mediated magnetoelectric coupling, the parallel magnetizations in the MTJ reorient to the x axis with opposite paths because of dipole interaction, thus resulting in antiparallel alignment. Moreover, this voltage switching of MTJs is nonvolatile owing to variations in dipole interaction and can be well understood via phase field simulations. The results provide an avenue to realize electrical switching of MTJs and are significant for exploring energy-efficient spintronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows

We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.

compressible flows↗
DOE OSTI2023

Multiscale Ecosystem for solving Maxwell-Schrodinger equations of open quantum systems (OpenMS)

Light-matter interactions play an important role in many branches of physics, chemistry, energy, and materials science. In the strong coupling regime, light-matter interactions are able to tune the materials properties via the formation of new quasiparticles, such as plasmons and polaritons. However, theoretical and numerical modeling of the light-matter interaction-mediated processes remain a big challenge because light-matter interactions are fundamentally multiscale and multiphysics problems involving multiple interactions between electrons, nuclei, and photons at different time/length scales. In the current community, light-matter interactions were treated at different levels of theoretical complexity in quantum chemistry and quantum optics. In quantum optics or quantum photonics, the matter is usually simplified as a few-level system, and the light is treated quantum-mechanically. On the other hand, quantum chemistry explores first-principles methods, including both single-particle and many-body-based techniques, to describe the electronic properties of matter in detail. However, light is usually prescribed as a classical electromagnetic field, and the light-matter interaction is taken into account as an external potential via classical approximations. This software is designed to fill current modeling shortcomings by delivering a first-ever scalable multiscale platform for simulating light-matter interactions in realistic electromagnetic environments. The software solves Maxwell and Schrodinger equations self-consistent on the heterogeneous platforms. It implements HF/DFT, TDDFT, and coupled-cluster counterparts for light-matter interactions and adopts modular programming to offload massively parallel algorithms on a large number of CPU/GPUs.

Zhang, Yu↗
DOE OSTI2020

Separate Effects Validation for Subcooled Boiling in CTF

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution modeling and simulation (M&S) tools for modeling of light water reactors (LWRs). This hub, the Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called the Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including CTF (the CASL version of Coolant Boiling in Rod Arrays— Three Field (COBRA-TF)), Michigan Parallel Characteristics Transport (MPACT), BISON, and Materials Performance and Optimization (MPO) Advanced Model for Boron Analysis (MAMBA). As part of its M&S efforts, CASL has identified various challenge problems, including Crud Induced Power Shift (CIPS), Crud-Induced Localized Corrosion (CILC), Pellet-Cladding Interaction (PCI), and Departure from Nucleate Boiling (DNB). This work addresses CASL milestone L2:VVI.P19.03, which focuses on uncertainty quantification of crud, which is relevant to both CIPS and CILC. This is achieved through an analysis and separate effects validation of the thermal hydraulic phenomenon known as subcooled boiling. As part of this work, various sources of experimental data are examined and compared to different options for empirical modeling of subcooled boiling. Through this analysis, a complete understanding of the underlying models and their implementation details are understood. A subset of these data are incorporated into a separate effects validation study of CTF. The Westinghouse Advanced Loop Tester (WALT) and Rohsenow experiments are modeled, and it is shown that the newly-implemented Gorenflo correlation is more accurate than the existing Chen and Thom correlations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗