Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solutions”

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 361 records · Page 20

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

ABSTRACT Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell’s equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly rotating magnetized neutron stars, taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. Rs ∼ 0.5 [R*], we show that the new terms lead to a per cent order correction in the magnetic field orientation and strength compared to the case, with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable, when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.

Torres, R. (ORCID:0000000291820228)↗

Exact solutions for the electromagnetic fields of a flying focus

The intensity peak of a “flying” focus travels at a programmable velocity over many Rayleigh ranges while maintaining a near-constant profile. Assessing the extent to which these features can enhance laser-based applications requires an accurate description of the electromagnetic fields. Here we present exact analytical solutions to Maxwell’s equations for the electromagnetic fields of a constant-velocity flying focus, generalized for arbitrary polarization and orbital angular momentum. Furthermore, the approach combines the complex source-point method, which transforms multipole solutions into beam-like solutions, with the Lorentz invariance of Maxwell’s equations. Propagating the fields backward in space reveals the space-time profile that an optical assembly must produce to realize these fields in the laboratory. Comparisons with simpler paraxial solutions provide conditions for their reliable use when modeling a flying focus.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Testing meson portal dark sector solutions to the MiniBooNE anomaly at the Coherent CAPTAIN Mills experiment

A solution to the MiniBooNE excess invoking rare three-body decays of the charged pions and kaons to new states in the MeV mass scale was recently proposed as a dark-sector explanation. This class of solution illuminates the fact that, while the charged pions were focused in the target-mode run, their decay products were isotropically suppressed in the beam-dump-mode run in which no excess was observed. This suggests a new physics solution correlated to the mesonic sector. We investigate an extended set of phenomenological models that can explain the MiniBooNE excess as a dark sector solution, utilizing long-lived particles that might be produced in the three-body decays of the charged mesons and the two-body anomalous decays of the neutral mesons. Over a broad set of interactions with the long-lived particles, we show that these scenarios can be compatible with constraints from LSND, KARMEN, and MicroBooNE, and evaluate the sensitivity of the ongoing and future data taken by the Coherent CAPTAIN Mills experiment to a potential discovery in this parameter space. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Perturbative stability of non-Abelian electric field solutions

We consider SU(2) gauge theory with a scalar field in the fundamental representation. The model is known to contain electric field solutions sourced by the scalar field that are distinct from embedded Maxwell electric fields. We examine the perturbative stability of the solution and identify a region of parameter space where the solution is stable. In the regime where the scalar field has a negative mass squared, the solution has two branches, and we identify an instability in one of the branches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Learning Optimal Power Flow Solutions using Linearized Models in Power Distribution Systems

Solving nonlinear optimal power flow (OPF) problem is computationally expensive, and poses scalability challenges for power distribution networks. An alternative to solving the original nonlinear OPF is the linear approximated OPF models. Although, these linear approximated OPF models are fast, the resulting solutions may result in significant optimality gap. Lately, the application of machine learning (ML) methods in successfully solving the nonlinear OPF has been reported. These methods learn and estimate the nonlinear control policies using a purely data-driven approach. In this paper, we propose an approach to complements the ML based approach to solving OPF using solutions from known linearized OPF model. Specifically, we use supervised learning to map the solutions of linear OPF to nonlinear control variables. Unlike, the traditional ML based methods for OPF that approximate the full distribution feeder model using function approximation, our approach uses a two-node approximation of radial networks. The proposed approach is validated using IEEE 123 bus test system for OPF solutions obtained using the nonlinear OPF models.

optimal power flow, power distribution systems, su↗

A MultiGPU Performance-Portable Solution for Array Programming Based on Kokkos

Today, multiGPU nodes are widely used in high-performance computing and data centers. However, current programming models do not provide simple, transparent, and portable support for automatically targeting multiple GPUs within a node on application areas of array programming. In this paper, we describe a new application programming interface based on the Kokkos programming model to enable array computation on multiple GPUs in a transparent and portable way across both NVIDIA and AMD GPUs. We implement different variations of this technique to accommodate the exchange of stencils (array boundaries) among different GPU memory spaces, and we provide autotuning to select the proper number of GPUs, depending on the computational cost of the operations to be computed on arrays, that is completely transparent to the programmer. We evaluate our multiGPU extension on Summit (#5 TOP500), with six NVIDIA V100 Volta GPUs per node, and Crusher that contains identical hardware/software as Frontier (#1 TOP500), with four AMD MI250X GPUs, each with 2 Graphics Compute Dies (GCDs)for a total of 8 GCDs per node. We also compare the performance of this solution against the use of MPI + Kokkos, which is the cur-rent de facto solution for multiple GPUs in Kokkos. Our evaluation shows that the new Kokkos solution provides good scalability for many GPUs and a faster and simpler solution (from a programming productivity perspective) than MPI + Kokkos.

Valero Lara, Pedro↗

Application of the Solute-Solvent EMF Cell to Measure Activity of NiF 2 in Molten FLiNaK

Electromotive force (emf) measurements made using a combination of solute- and solvent-based electrodes were used to determine the activity of NiF 2 in molten FLiNaK eutectic at 823 K across a concentration range of x NiF2 = 5.2 × 10 –4 –1.0 × 10 –2 . The solute emf values were measured using electrodes consisting of Ni wires immersed in FLiNaK with dissolved NiF 2 contained in graphite crucibles. The measured emf values were then converted to the FLiNaK Eutectic Potassium Electrode (FEKE) potential and used to quantify the activity of dissolved NiF 2 . This quantification was based upon comparative measurements of a reference solvent electrode consisting of a K-Bi alloy immersed in pure FLiNaK contained in a boron nitride crucible and a solute electrode. Short cell lives were characteristic of the measurements due to the corrosive nature of the fluoride salts. Quantifying the activity of NiF 2 will improve the utility of Ni 2+ /Ni reference electrodes in molten fluoride salts, which are notoriously difficult electrolytes to work with because of their reactivity. In conclusion, this work demonstrates the general nature of the solute-solvent approach as a repeatable, easily employed method for measuring the activity values of electroactive species in a variety of molten salts to improve understanding of the electroactive species behavior in these systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗

Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics

Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code

CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Federated Architecture for Secure and Transactive Distributed Energy Resource Management Solutions (FAST-DERMS)

This document provides system-level specifications for a federated architecture for secure and transactive distributed energy resource management solutions (FAST-DERMS), presents a solution, and describes operational concepts for the proposed solution. FAST-DERMS enables the provision of reliable, resilient, and secure transmission and distribution (T&D) grid services through the scalable aggregation and near-real-time management of utility-scale and small-scale distributed energy resources (DERs). We first present the principles and objectives of FAST-DERMS. Then, after discussing important system concepts, we present the specifications for FAST-DERMS and a solution that employs a distributed and federated control methodology in which the DERs connected to a single point of common coupling with the rest of the system, such as individual substations, are optimized coordinately to provide system-level grid services. FAST-DERMS aims to aggregate and coordinate the operations of DERs to support T&D grid operations. The key optimization and control component of this FAST-DERMS reference implementation is a flexible resource scheduler (FRS) that aggregates the DERs within a substation service area. These FRSs operate at the substation level and perform constrained economic dispatch of DERs, either directly or through a transactive market or aggregator, as shown in Figure ES-1. An FRS Coordinator at the distribution system operator (DSO) level aggregates distribution substations operated by FRSs and interfaces with the transmission system operator (TSO) to provide transmission services. FAST-DERMS also allows for the integration of the FRS Coordinator with an existing distribution utility management system that could be employed by the DSO to enhance distribution grid operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Non-Equilibrium Plasma Interactions with Biomaterials, Biological Solutions and Tissues

Cold atmospheric pressure plasma discharges offer an abundant source of reactive oxygen and nitrogen species (RONS) at room temperature enabling unique interactions with biomaterials, biological solutions and tissues. These interactions particularly with living matter are presently an important intellectual frontier in plasma science with promising potential applications ranging from human health care to advanced biomaterial processing. Exciting case studies have been reported that illustrate the huge potential of cold atmospheric plasma technology in wound healing and cancer treatment. The interaction of plasma with conducting and dielectric biomaterials such as tissue strongly influences the plasma properties. In turn this changes the impact of the plasma on the biomaterial. Particularly in the case of living matter, liquid based solutions are ubiquitous which complicates interfacial processes. The lack of insight into the underlying mechanisms of the interaction of plasma with wounds and tumors is currently a bottleneck for the further development of the technology and gives rise to many interesting scientific questions. This project was focused on plasma properties and kinetics during plasma-biomaterial interactions. Both DC pulsed and RF driven atmospheric pressure plasma jets, extensively used by the plasma community were studied. The bio-interfaces included hydrogel as a tissue model, (saline) solutions, bacteria and virus. Plasma diagnostics used include Thomson scattering, Rayleigh scattering, Raman scattering, (two-photon absorption) laser induced fluorescence, optical emission spectroscopy, absorption spectroscopy, molecular beam mass spectrometry and fast imaging allowing to determine electron densities and temperatures, ionic species, reactive species including radicals, gas temperatures, gas composition, electric fields and solution components transferred to the gas phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of Irradiation and Decay Time on Chemical Speciation in Uranyl Sulfate Solution

A pH-1, 200-g U/L, highly enriched uranyl sulfate solution was neutron-irradiated using an accelerator-driven neutron source (the Universal Neutron Irradiator, UNI) to fission U-235 under conditions similar to those used by SHINE to produce Mo-99 in a subcritical assembly. We completed a series of five uranyl sulfate irradiations to assess the impact of irradiation and decay time on uranium precipitation and iodine speciation in solution. The duration of each irradiation was increased in order, with the shortest irradiation taking 2.9 hours and the longest, 64.5 hours. A total of 9.5E16 U-235 nuclei were fissioned in all irradiations, and a maximum fission power density of 0.26 W/mL was reached in the final irradiation. We observed no evidence of uranyl peroxide or other precipitates, demonstrating that the addition of Fe 2+ to a final Fe concentration of 200 ppm in the target solution was successful at preventing precipitation under representative irradiation conditions. Nonvolatile I - was found to be the dominant (>70%) iodine species in solution, both immediately after irradiation and after days to weeks of decay time.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Options for the Removal of Cesium and Strontium from Acidic Uranyl Sulfate Solutions from Molybdenum-99 Production

A small laboratory-scale column process aimed at removing fission products from prototypical acidic sulfate solution wastes from a prototypical Mo-99 recovery process has been operated in a radioactive materials hood at Savannah River National Laboratory. This equipment has been operated with realistic concentrations of uranium (0.4 wt % 235 U), plutonium ( 239 Pu), neptunium ( 237 Np), non-radioactive surrogates of common fission products (including Cs and Sr) and two promising absorbents AMP-PAN and CST R91290-B. The goal of this work was to demonstrate the removal and concentration of high specific activity isotopes which are expected to dominate the classification of the low-level waste from the production of 99 Mo. Both AMP-PAN and CST R9120-B are attractive absorbents for the removal of 137 Cs from an acidic sulfate waste stream. They do not provide a solution for the removal of 90 Sr from the acidic solutions, but CST should be investigated at a higher pH. Limited removal of transuranic species and almost no removal of uranium was observed in this testing. This could be important due to the relatively large amount of uranium expected to be present in such waste streams. The capacity of these absorbents for 137 Cs will likely be dependent on the Cs concentration in the solutions processed but loading in the range of 0.2 to 5 mg Cs per gram of absorbent appears to be achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Interaction of Concentrated (68%) Aqueous Nitric Acid Solutions with Cheesecloth: A TGA-DSC-MS Study

The exposure of cotton cheesecloth to concentrated aqueous solutions of nitric acid reduced the onset of combustion by 30-40°C (from ~350°C to ~310°C) relative to that of both unexposed cheesecloth and cheesecloth that was rinsed immediately following exposure to the same solution. Cheesecloth exposed to concentrated (68%) nitric acid underwent three exothermic pre-combustion reactions which onset at 125-135°C, 160- 170°C and 200-210°C. These reactions give off mostly water vapor, with lesser amounts of NO and carbon dioxide. Exposure to concentrated nitric acid damages the cellulose polymers, oxidizing them, reducing them in length, and increasing their reactivity in response to thermal stress. Cheesecloth which was rinsed after exposure behaved similarly to unexposed cheesecloth, indicating that NC was not formed on exposure to the concentrated aqueous nitric acid solution. The addition of a 239Pu – 241Am spike in an aqueous nitric – hydrofluoric solution to the cheesecloth reduced the onset temperature of the first exothermic reaction to 126°C and lowered the onset temperature of combustion from ~310°C to ~275°C. When dried on to cellulose substrates in the presence of excess nitrate, metal ions may act as oxidizers, or may catalyze the oxidation of the carbonaceous material, lowering ignition temperatures. No evidence for the formation of nitrocellulose following the exposure of cheesecloth to nitric acid was observed during this study

36 MATERIALS SCIENCE↗

Evaluating a Semi-Quantitative, Spectrophotometric Method for the Rapid Determination of Total Rare Earth Element Concentrations in Citrate Leaching Solutions

Quick, reliable methods to measure the rare earth elements (REEs) in process streams are needed to support real-time parameterization and monitoring for a burgeoning number of REE extraction schemes. Such methods would ideally be fast, reproducible, and field-deployable. This study evaluates the feasibility of using a chromogenic indicator, arsenazo(III) (hereafter abbreviated as Arsz(III)) to detect the concentration of REEs in complex citrate leaching solutions using a spectrophotometer. Arsz(III) forms a chromogenic complex with the lanthanide elements that produces a blue/purple color detectable in the visible light range. Matrices of solutions containing varied concentrations of calcium, iron, and lanthanum were scanned in the presence of an Arsz(III) indicator to generate a modeled relationship that could be used to back calculate REE concentrations. While relatively quick (scan times on the order of seconds), the results indicated that the calcium and iron concentrations of typical mining solutions overwhelm the signal of the REEs, to the extent that sensitivity of the method is no less than 50 µM REEs. These results demonstrate the difficulties of accurately measuring trace metals in complex solutions, although the method may still hold utility for processes with significantly lower calcium and iron concentrations or REE concentrations much greater than 50 µM.

36 MATERIALS SCIENCE↗

Performance Demonstration of an Occupancy Sensor-enabled Integrated Solution for Commercial Buildings

Traditionally, a single-loop fixed-gain controller is applied to supply fan (SF) and cooling coil (CC) valve controls while a fixed-damper position control is applied to outdoor air (OA) damper control at air handling units (AHUs) in commercial buildings. With the increasing application of occupancy sensors, the information generated by occupancy sensors is applied to not only reduce the electricity loads from lighting and controllable plug loads, but also reset OA intake and minimum supply airflow setpoints. Meanwhile, these intermittent operation actions greatly elevate the dynamics of AHU systems, which may introduce unstable SF and CC valve operations and inaccurate OA flow control at AHUs and consequently degrade maximum energy efficiency gains. With virtual fan and valve flow meter technologies, two advanced controls, including cascade control and gain scheduling control, can be implemented on both the SF and CC valve, and an advanced control using a virtual OA flow meter can be implemented on OA damper integrated with occupancy sensors. The goal of this project is to demonstrate the savings, cost, and performance of an integrated solution that integrates the three advanced HVAC controls with occupancy sensors to allow accurate and stable AHU operations in real buildings. The project objectives are to: 1) develop and validate an advance SF control algorithm; 2) develop and validate an advanced CC valve control algorithm; 3) validate an algorithm to implement a virtual OA flow meter; and 4) demonstrate the savings, cost, and performance of the integrated solution in real buildings. The technical approaches are to: 1) select a test system at the University of Oklahoma; 2) develop and implement the algorithms of advanced SF and CC valve controls and validate the performance; 3) develop and implement the advanced OA control using a virtual OA flow meter and validate the performance; 4) demonstrate the savings, cost, and performance of the proposed integrated solution with and without three advanced HVAC controls; and 5) disseminate the project results through publications and presentations. For the SF control, both the gain scheduling and cascade controls can improve the fan energy performance by reducing the fan power during the transient period and the fan control performance at lower speeds by reducing fan speed variation. Moreover, the gain scheduling control provides a simple and low-cost solution and is recommended. The fan power savings can reach 30% in a transient period. For the CC valve control, the gain scheduling control can considerably reduce the supply air temperature oscillation range and frequency under both higher and lower load conditions and the control valve response is much more stable. As a result, the gain scheduling control is recommended. The projected pump energy consumption can be reduced by 68.5%. With the developed virtual OA flow meter, the OA can be accurately controlled at its setpoint, which is determined based on the actual number of occupants in the building provided by occupancy sensors. The RMSE of the proposed OA control is 15.9 L/s. The energy data shows that the fan power and CC cooling energy were significantly reduced. On the other hand, the energy savings majorly results from the occupancy sensors and the energy savings by the advanced HVAC controls is minimal because that the controllers in the test AHU were tuned with very slow response. An annual technical savings potential is estimated as 0.5 quads in the commercial sector.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗