SEARCH · Engineering Papers

Results for “Mixture design”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Chemically Recyclable Polyester Thermosets from Activated Adipic Acid and Renewable Polyols

This study outlines a method for producing chemically recyclable crosslinked polyesters using renewable polyols-glycerol and sorbitol-combined with adipic acid (AA), which is transformed/activated into a polyanhydride mixture prior to use. A three-step procedure has been designed: 1) an acid-catalyzed reaction of AA with nontoxic isopropenyl acetate or acetic anhydride to form a crosslinking mixture (CLM) made of adipic-acetic mixed polyanhydrides; 2) a solvent- and additive-free process where glycerol or sorbitol, or a combination thereof, is reacted with the CLM to achieve a prepolymer, and 3) a casting/molding of the liquid viscous prepolymer to yield a thermoset as the end product. Different thermosets (eight examples) are prepared by changing the reagents ratio. These solids are thoroughly characterized by tensile tests, DMA, high-resolution magic angle spinning and solid-state NMR, thermal gravimetric analysis, DSC, and fourier transformed infra red (FT-IR) spectroscopy. The formation of cross-linked polyesters is confirmed in all cases, but mechanical properties varied significantly from one specimen to another. Interestingly, a tensile strength up to 18 MPa-approximately an order of magnitude higher than similar polymers-is achieved when sorbitol and the CLM are used in a 1:1 wt% ratio. The chemical recycle of the resulting polymers is achieved via methanolysis with quantitative recovery of the monomeric units.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Molecular dynamics simulations of uranyl and plutonyl cations in a task-specific ionic liquid

Ionic liquids (ILs) are a unique class of solvents with potential applications in advanced separation technologies relevant to the nuclear industry. ILs are salts with low melting points and a wide range of tunable physical properties, such as viscosity, hydrophobiciy, conductivity, and liquidus range. ILs have negligible vapor pressure, are often non-flammable, and can have high thermal stability and a wide electrochemical window, making them attractive for use in separations processes relevant to the nuclear industry. Metal salts generally have a low solubility in ILs; however, by incorporating new functional groups into the IL cation or anion that promote complexation with the metal, the solubility can be greatly increased. One such task-specific ionic liquid (TSIL) is 1-carboxy-N, N, N-trimethylglycine bis(trifluoromethylsulfonyl)imide ([Hbet][Tf 2 N]). Water, which is detrimental for electrochemical separations, is a common impurity in ILs and can coordinate with actinyl cations, particularly in ILs containing only weakly coordinating components. Understanding the behavior of actinides in TSIL/water mixtures on a molecular level is vital for designing improved separations processes. Classical molecular dynamics simulations of uranyl(VI) and plutonyl(VI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]) with deprotonated Hbet (betaine) and water have been performed to understand the coordination and dynamics of the actinyl cations. We find that betaine is a much stronger ligand than water and prefers to coordinate the metal in a bidentate manner. Potential of mean force simulations yield a relative free energy for betaine coordination of approximately -120 to -90 kJ/mol in mixtures with water. As the amount of betaine coordinated to the actinide increases, the diffusion coefficient of the actinyl cation decreases. Moreover, the betaine ligand is able to bridge between two metal centers, resulting in dimeric complexes with actinide–actinide distances of ~5 Å. Potential of mean force simulations show that these structures are stable, with relative free energies of up to -40 kJ/mol. The crystal structure for [(UO 2 ) 2 (bet) 6 (H 2 O) 2 ][Tf 2 N] 4 shows that the betaine bridges between two uranium atoms to form dimeric complexes similar to those found in our simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Robust Copper-Based Nanosponge Architecture Decorated by Ruthenium with Enhanced Electrocatalytic Performance for Ambient Nitrogen Reduction to Ammonia

Electrochemical conversion of nitrogen to green ammonia is an attractive alternative to the Haber–Bosch process. However, it is currently bottlenecked by the lack of highly efficient electrocatalysts to drive the sluggish nitrogen reduction reaction (N2RR). In this work, we strategically design a cost-effective bimetallic Ru–Cu mixture catalyst in a nanosponge (NS) architecture via a rapid and facile method. The porous NS mixture catalysts exhibit a large electrochemical active surface area and enhanced specific activity arising from the charge redistribution for improved activation and adsorption of the activated nitrogen species. Benefiting from the synergistic effect of the Cu constituent on morphology decoration and thermodynamic suppression of the competing hydrogen evolution reaction, the optimized Ru 0.15 Cu 0.85 NS catalyst presents an impressive N 2 RR performance with an ammonia yield rate of 26.25 μg h –1 mg cat. –1 (corresponding to 10.5 μg h –1 cm –2 ) and Faradic efficiency of 4.39% as well as superior stability in alkaline medium, which was superior to that of monometallic Ru and Cu nanostructures. Additionally, this work develops a new bimetallic combination of Ru and Cu, which promotes the strategy to design efficient electrocatalysts for electrochemical ammonia production under ambient conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗

DOE OSTI2020

Mixed Material Scintillator Systems Particle ID Modelling Report

AMixed-Material Scintillator System (MMSS) is a radiation detector using a scintillator made of a heterogeneous, structured mixture of two or more materials. MMSSs are designed so that the structured mixing of materials encodes properties of the radiation detected in the scintillator. This allows for new radiation detectors with advantages over detectors using traditional homogeneous scintillators. This report contains results from this project’s second task, to model and quantify the promise of the MMSSs invented earlier in the project. Those inventions have been divided into two categories, and this report focuses on the first category: MMSS inventions the use zones of differently-colored scintillators to detect and characterize neutron sources. We call this category particle-ID, or PID, detectors. Inventions using gradients designed to encode the position of radiation interactions using gradients will be discussed in a future report.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Mixed Material Scintillator Systems Position Reconstruction Modelling Report

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Multicomponent gas mixture parametric CFD study of condensation heat transfer in small modular reactor system safety

Safety is always the primary concern for designing and analyzing nuclear reactor systems. The requirements for the safety margin for advanced small modular reactor (SMR) systems are targeted even higher than the conventional commercial large-scale nuclear reactors incorporating the passive and inherent safety systems. The SMR systems are designed with the condensation passive containment cooling system (PCCS), which plays a critical role in removing reactor heat during a steam release accident case. However, the presence of non-condensable gas (NCG), like air, reduces the heat transfer performance. This physics phenomenon becomes multifactorial for nuclear reactor containment during a fuel failure accident case that releases hydrogen gas. Besides, the mixture component of steam-air-hydrogen varies in reactor accident cases, which need simulation and validation keeping parameters of importance. Reviews showed that previous studies for SMR’s PCCCS did not cover the condensation heat transfer (CHT) in the presence of multicomponent NCG mixture parametric computational fluid dynamics (CFD) simulation and validation, making a research gap in the SMR design safety. A comprehensive CHT parametric CFD study was performed for SMR PCCS to fill this research gap. This study used experimental data as simulation 3D physics domain inlet and outlet boundary conditions. However, the wall boundary conditions were constant temperature, curve-fit, and annular coolant for verifying the turbulence models. Parametric simulations were performed, verified, and optimized for steam-NCGs mixtures. The multicomponent gases, multiphase mixtures, and fluid film condensation models were applied with associated turbulence models. The results of the parametric study were evaluated for realistic reactor conditions. Results showed that parametric study provided critical insight about the dependency of multicomponent gas mixture parameters that supports reactor safety design, analysis, and licensing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Measurements of Density and Sound Speed in Mixtures Relevant to Supercritical CO2 Cycles

Abstract The objective of this research is to validate properties of mixtures relevant to supercritical carbon dioxide (sCO2) power cycles. Direct-fired sCO2 cycles are a promising technology for the future power generation systems. The working fluid of sCO2 cycles will be near and above critical point of CO2. One of the challenges is that the simulation of mixtures should consider real gas behavior. Expected operating conditions of Allam cycles reach up to 300 bar and 1000 °C. Characterizing the mixtures at the extreme conditions is an important issue in current researches and industrial applications. Thermophysical properties of mixtures may be beyond the valid range of the widely used database such as NIST REFPROP. Experimental data of mixture properties in the literature are limited which is necessary to develop high-fidelity design tools for sCO2 power cycles. We measured the density and sound speed of several multi-component mixtures. A temperature-controlled high-pressure test cell was used for the density measurements. Sound speed was measured by resonant frequency detection using an external speaker and a piezoelectric pressure sensor. Mixtures studied in this work include carbon dioxide, methane, oxygen, and water vapor. Properties of pure CO2 were measured to show the validity of our technique. Compositions were selected to be close to frozen mixtures at the inlet, mid-progress, and exhaust conditions of a model sCO2 combustor in the previous numerical simulation work. Corresponding reaction progress variables (RPV) were RPV = 0, 0.5, and 1. Temperature and pressure conditions of experiments are 310–450 K and 0–150 bar. In our study, density and sound speed from the NIST REFPROP database agree with experimental measurements within the range of our measurement uncertainties.

Energy & Fuels↗

DOE OSTI2020

Experimental and Simulation Analysis of Binary Mixtures of Biomass and Inert Material

A literature review was conducted to identify experiments on various scales of biomass combustion specific to circulating fluidized bed (CFB) combustors. Typical feedstock and inert bed material characteristics for biomass combustion in a CFB were gathered. This information was used to design laboratory-scale, non-reacting experiments on biomass-inert bed material mixtures in a 1-in. cylindrical test rig in the Multiphase Flow Analysis Laboratory (MFAL) in Morgantown, WV. The biomass chosen was a hardwood. Pellets from the hardwood were milled and sieved to obtain three size ranges—small, medium, and large. Three mass loadings of each of the biomass samples were fluidized with glass beads as the inert bed material to study the mixing characteristics under varying fluidizing gas velocities. Hydrodynamic data including pressure drop measurements at select vertical locations and video recordings of the mixture behavior were generated in this effort.

09 BIOMASS FUELS↗

DOE OSTI2023

Rational Design of A Chemical Bath Deposition Based Tin Oxide Electron Transport Layer for Perovskite Photovoltaics

Chemical bath deposition is widely used to deposit SnO x as an electron transport layer in perovskite solar cells (PSCs). The conventional recipe uses thioglycolic acid (TGA) to facilitate attachments of SnO x particles onto the substrate. However, nonvolatile TGA has been reported to harm the operational stability of PSCs. In this work, we introduced a volatile oxalic acid (OA) as an alternative to TGA. OA, a dicarboxylic acid, functions as a chemical linker for the nucleation and attachment of particles to the substrate in the chemical bath. Moreover, OA can be readily removed through thermal annealing followed by a mild H 2 O 2 treatment, as shown by FTIR measurements. Synergistically, the mild H 2 O 2 treatment selectively oxidizes the surface of the SnOx layer, minimizing nonradiative interface carrier recombination. EELS (electron-energy-loss-spectroscopy) confirms that the SnO x surface is dominated by Sn 4+ , while the bulk is a mixture of Sn 2+ and Sn 4+ . This rational design of a CBD SnO x layer leads to devices with T85~1,500h, a significant improvement over the TGA-based device with T80~250h. Our champion device reached a power conversion efficiency of 24.6%. This work offers a rationale for optimizing the complex parameter space of CBD SnO x to achieve efficient and stable PSCs.

14 SOLAR ENERGY↗

DOE OSTI2021

Multi-objective optimization of low-GWP mixture composition and heat exchanger circuitry configuration for improved system performance and reduced refrigerant flammability

The recently introduced hydrofluoroolefin (HFO) refrigerants, including R1234yf and R1234ze(E), have significantly lower global warming potentials (GWPs) than traditional hydrofluorocarbon (HFC) refrigerants like R410A. However, prior tests show that direct drop-in of pure R1234yf or R1234ze(E) into equipment designed for R410A results in a decrease in heat exchanger capacity and the system coefficient of performance. The primary reason is the lower in-tube heat transfer performance of R1234yf and R1234ze(E) compared with that of R410A. To address this issue, previous studies have mixed the mildly flammable HFC R32 with HFOs to improve system performance, with HFC R125 also added to suppress flammability. Previous studies selected compositions based on simple cycle analyses and did not consider modifications of the heat exchanger circuitry configuration to adapt to the new refrigerants. This study presents a novel multi-objective optimization approach to design a refrigerant composition that maximizes energy efficiency within flammability and GWP limits. The approach in this work simultaneously optimizes mixture composition and heat exchanger circuitry configuration. A case study on a rooftop unit indicates that, compared with mixture-only optimization, simultaneous optimization of mixture and heat exchanger circuitry yields a 5.9% improvement in cycle efficiency and a 48.6% reduction in refrigerant flammability with a GWP of 268. Circuitry optimization using refrigerants with different temperature glides shows that the larger the temperature glide is, the larger EER improvement is obtained. The results show that zeotropic blends with a large temperature glide are more sensitive to the refrigerant circuitry than pure refrigerants and may suffer significant performance degradation with subpar heat exchanger circuitry design. The proposed optimization approach is generally applicable to mixtures with any number of components. Using this approach to design a HVAC system can yield higher system efficiency within flammability and GWP constraints.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Calculating the Effects of Solids Input and Removal as a Temperature Control in the Advanced Scale Up Reactor Experiment (ASURE) Facility at NETL Using Aspen

The Advanced Scale Up Reactor Experiment (ASURE) facility at NETL is being designed to be a fuel-flexible multi-purpose reactor that can be used for pyrolysis/gasification or evaluation of other high pressure “circulating fluidized bed” (CFB) chemical processes. The initial system design calculations for pyrolysis/gasification are presented in this work showing the expected performance of the ASURE facility when used as a biomass conversion reactor. Several other areas of application include gasification of any carbonaceous fuel including biomass, coal, plastics, and other waste materials. The reactor can therefore be used to produce SYNGAS of various compositions and hydrogen as well as other high value chemicals resulting from a typical tuned gasification process. This paper discusses an ASPEN model of the facility, focusing on the riser of the CFB reactor and the solids recirculation loop. The ASPEN model divides the riser into two sections. A bottom section which receives ash, char and sand which have been recirculated from a return loop. In this section an inert fluidization gas, (N2 or CO2), is introduced which acts as the primary mover of the solids through the system. The bottom section is equipped with a restricted air feed so that the recirculated char can be partially oxidized. This oxidation process along with the inventory of recirculating sand are used to effectively control the temperature in the following two chemical conversion sections of the reactor which are the pyrolysis zone followed by a tar cracking zone. Fresh fuel is added to the pyrolysis zone and undergoes drying and devolatilization. The products ash, char, volatile matter, and water vapor exit the pyrolysis zone and enter the reaction block for tar cracking. Steam and CO2 gasification reactions will be incorporated into the tar cracking zone, however at the design operating temperature, conversion from these reactions is expected to be essentially zero. The unit when completed in 2026 will test mixtures of biomass, plastics, and waste coal. This presentation discusses the basic ASPEN engineering design model for this project and provides preliminary sensitivity studies to determine how the various parts of the reactor will perform.

ASSURE↗

DOE OSTI2025

Calculating the Effects of Solids Input and Removal as a Temperature Control in the Advanced Scale Up Reactor Experiment (ASURE) Facility at NETL Using Aspen

The Advanced Scale Up Reactor Experiment (ASURE) facility at NETL is being designed to be a fuel-flexible multi-purpose reactor that can be used for pyrolysis/gasification or evaluation of other high pressure “circulating fluidized bed” (CFB) chemical processes. The initial system design calculations for pyrolysis/gasification are presented in this work showing the expected performance of the ASURE facility when used as a biomass conversion reactor. Several other areas of application include gasification of any carbonaceous fuel including biomass, coal, plastics, and other waste materials. The reactor can therefore be used to produce SYNGAS of various compositions and hydrogen as well as other high value chemicals resulting from a typical tuned gasification process. This paper discusses an ASPEN model of the facility, focusing on the riser of the CFB reactor and the solids recirculation loop. The ASPEN model divides the riser into two sections. A bottom section which receives ash, char and sand which have been recirculated from a return loop. In this section an inert fluidization gas, (N2 or CO2), is introduced which acts as the primary mover of the solids through the system. The bottom section is equipped with a restricted air feed so that the recirculated char can be partially oxidized. This oxidation process along with the inventory of recirculating sand are used to effectively control the temperature in the following two chemical conversion sections of the reactor which are the pyrolysis zone followed by a tar cracking zone. Fresh fuel is added to the pyrolysis zone and undergoes drying and devolatilization. The products ash, char, volatile matter, and water vapor exit the pyrolysis zone and enter the reaction block for tar cracking. Steam and CO2 gasification reactions will be incorporated into the tar cracking zone, however at the design operating temperature, conversion from these reactions is expected to be essentially zero. The unit when completed in 2026 will test mixtures of biomass, plastics, and waste coal. This paper presents the basic ASPEN engineering design model for this project and provides preliminary sensitivity studies to determine how the various parts of the reactor will perform.

ASSURE↗

DOE OSTI2021

Optimization of Refrigerant Compositions for Low-GWP Refrigerant Mixtures Using Segment-by-segment Heat Exchanger and Detailed System Models

Li, Zhenning↗

DOE OSTI2025

Modeling hydraulic fracture opening and closure with proppant transport and settlement

Hydraulic fracturing is a widely used reservoir stimulation technique for improving fluid circulation in rock formations with extremely low permeability, particularly in enhanced geothermal systems (EGS). To better understand the complex processes involved and improve hydraulic stimulation performance, we have developed ELK (ELectrical fracKing), a MOOSE-based 3D finite element application designed to model the behavior of proppant-fluid mixtures in propagating fractures. ELK integrates both the fluid and proppant components, incorporating particle-driven processes such as gravity settling, particle-particle interactions, and strong density and viscosity contracts, in addition to conventional fluid-driven fracture propagation. In this contribution, we extend ELK to model propped fracture closure, which occurs after the injection phase due to a dramatic drop in the effective stress on the fracture plane. During the shut-in, flowback, and production periods, the fracture width decreases, with the closure behavior depending on proppant concentration.. At low concentrations, closure follows a nonlinear joint law linked to the stiffness of asperities in the fracture walls. While at high concentrations, it is controlled by the properties of packed proppant bed. The extended ELK application is validated against several benchmark examples, including the propagation of an inclined frictional crack, fracture opening and sliding in response to fluid injection, and flowback analysis. We believe that ELK’s enhanced capabilities can serve as a valuable tool for the design and optimization of EGS deployment.

15 - GEOTHERMAL ENERGY↗

DOE OSTI2021

Neutronic Performance of a Compact 10 MWe Nuclear Reactor with Low Enrichment (ThxU1−x)N Fuel

Abstract Recent interest in compact nuclear reactors for applications in space or in remote locations drives innovation in nuclear fuel design, especially non-oxide ceramic nuclear fuels. This work details neutronic modeling designed to support the development of a new nuclear fuel concept based on a mixture of thorium and uranium nitride. A Monte Carlo N-Particle Version 6.2 (MCNP-6) model of a compact 10 MWe reactor design which incorporates (Th x U 1− x )N fuel is presented. In this context, a “compact” reactor is a completely assembled reactor which may be emptied of coolant and transported by specialized commercial vehicle, deployed by a C130J aircraft, or launched into space. Core geometry, reflector barrels, and the heat exchange zones are designed to support reduction of overall reactor volume of core components while maintaining criticality with a fixed total fuel mass of 4500 kg. Dense mixed nitrides of thorium nitride (ThN) additions in uranium nitride (UN) in 5 wt.% increments between $$0.05 \le x \le 0.5$$ 0.05 ≤ x ≤ 0.5 have been considered for calculation of $$k_{\infty }$$ k ∞ and $$k_{{{\text{effective}}}}$$ k effective . ThN additions in UN results in a slight increase in the magnitude of the temperature coefficient of reactivity, which is negative by design. The isotopic distribution of the principal actinide inventory as a function of burnup, time, and initial fuel composition is presented and discussed within the context of the proliferation risk of this core design.

Parker, S. S. (ORCID:000000018597440X)↗

DOE OSTI2020

A new approach for measuring the carbon and oxygen content of atmospherically relevant compounds and mixtures

Abstract. Due to its complexity, gas- and particle-phase organic carbon in the atmosphere is often classified by its bulk physicochemical properties. However, there is a dearth of robust, moderate-cost approaches to measure the bulk chemical composition of organic carbon in the atmosphere. This is particularly true for the degree of oxygenation, which critically affects the properties and impacts of organic carbon but for which routine measurement approaches are lacking. This gap has limited the understanding of a wide range of atmospheric components, including particulate matter, the mass of which is monitored worldwide due to its health and environmental effects but the chemical characterization of which requires relatively high capital costs and complex operation by highly trained technical personnel. In this work, we demonstrate a new approach to estimate the mass of carbon and oxygen in analytes and mixtures that relies only on robust, moderate-cost detectors designed for use with gas chromatography. Organic compounds entering a flame ionization detector were found to be converted with approximately complete efficiency to CO2, which was analyzed downstream using an infrared detector to measure the mass of carbon analyzed. The ratio of the flame ionization detector (FID) signal generated to CO2 formed (FID∕CO2) was shown to be strongly correlated (R2=0.89) to the oxygen-to-carbon ratio (O∕C) of the analyte. Furthermore, simple mixtures of analytes behaved as the weighted average of their components, indicating that this correlation extends to mixtures. These properties were also observed to correlate well with the sensitivity of the FID estimated by structure activity relationships (quantified as the relative effective carbon number). The relationships between measured FID∕CO2, analyte O∕C, and FID sensitivity allow the estimation of one property from another with <15 % error for mixtures and <20 % error for most individual analytes. The approach opens the possibility of field-deployable, autonomous measurement of the carbon and oxygen content of particulate matter using time-tested, low-maintenance detectors, though such an application would require some additional testing on complex mixtures. With some instrumental modifications, similar measurements on gas-phase species may be feasible. Moreover, the potential expansion to additional gas chromatography detectors may provide concurrent measurement of other elements (e.g., sulfur, nitrogen).

42 ENGINEERING↗

DOE OSTI2022

Formulation of the density eigenvalue problem in neutron transport for relevant engineering applications

A new formulation of the density eigenvalue problem for the neutron transport equation is presented. This new eigenvalue, named ζ eigenvalue can be introduced freely in the transport model, acting on a selected portion of the phase space. Despite its broader applications and its connection with the nature of the transport operator, the ζ eigenvalue has been presented here mainly as a design-oriented technique for the efficient evaluation of the critical concentration for a specific nuclide (or mixture of nuclides). This new eigenvalue is particularly adequate to study the definition of the material composition in the criticality design process of a multiplying system. The method is then applied for the study of classical problems such as the critical moderation ratio and the poison concentration to control the reactor. Some results are presented in one dimensional configuration using the multi- group spherical harmonics approach. This eigenvalue formulation proves to be a convenient and useful way to attain criticality, also for complex, realistic systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗