Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “configuration mixing”

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,351 records · Page 75

Annual Status Report (FY 2020): Performance Assessment for the Environmental Restoration Disposal Facility

DOE O 435.1 and DOE M 435.1-1 require that a determination of continued adequacy of the performance assessment (PA) (CP-60089), composite analysis, and disposal authorization statement (DAS) be made annually, and these guidelines must be used to consider the results of data collection and analysis from research, field studies, and monitoring as well as provide the need to update any radioactive waste management basis documents. Beginning in 1996, the Environmental Restoration Disposal Facility (ERDF) started accepting low-level radioactive, hazardous, and mixed wastes generated during cleanup activities at the Hanford Site. ERDF is composed of a series of cells or disposal areas and can accommodate future design expansions as needed. Currently, there are 10 cells. During this reporting period (fiscal year 2020, which extended from October 1, 2019, through September 30, 2020), approximately 3.39E+04 U.S. tons (3.07E+04 metric tons) of waste was disposed at ERDF. From ERDF inception through September 30, 2020, approximately 18.5 million U.S. tons of waste has been disposed at ERDF, which equates to the consumption of approximately 88% of the disposal volume. As a condition of the DAS, disposal operations within ERDF must be in accordance with the waste acceptance criteria (ERDF-00011) that provide specific radionuclide disposal limits, waste form restrictions, and descriptions of acceptable waste packages in compliance with DOE M 435.1-1 requirements. The ERDF waste acceptance criteria stipulate that waste destined for disposal at ERDF be controlled based on source, physical form, and contaminant concentration and activity levels. There have been no changes to the physical configuration of ERDF or to the waste forms (source, physical form, etc.). No new Unreviewed Disposal Question Screenings or Evaluations have been generated during this reporting period. Therefore, there are no noted impacts to the PA, composite analysis, DAS, or radioactive waste management basis documents resulting from the evaluations and screenings. Sum of fraction analysis shows that the disposed inventory meets both the concentration and inventory threshold requirements. A sum of fractions value is computed for ERDF sensitive radionuclides contributing to the all pathways and air pathway inventory limits. Computed values were 8.85E-02 and 1.78E-01, respectively. The disposed waste inventory remained well under the PA imposed limits. Required monitoring was satisfactorily completed during the fiscal year reporting period (fiscal year 2020). Compliance with performance objectives were met as each of the reported values were well below the established limit. Overall, there are no substantive changes to primary PA assumptions nor changes to the PA analysis conclusion; therefore, compliance with DOE O 435.1 and the DAS is maintained.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ion Exchange Processing of AN-107 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System

The Hanford Site stores an estimated 56 million gallons of mixed radioactive and chemically hazardous waste in large underground tanks. In support of the Direct Feed Low-Activity Waste (DFLAW) Program for expediting Hanford tank waste supernate treatment, laboratory-scale ion exchange processing using prototypic unit operations was conducted on AN-107 tank waste at the Pacific Northwest National Laboratory Radiochemical Processing Laboratory. This report describes the small-scale ion exchange testing with 13.7 L of diluted and filtered supernate from Tank 241-AN-107 (hereafter referred to as AN-107) at 16 °C (62 °F). One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AN-107 tank waste to meet this criterion, only 0.147% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 678. Testing with AN-107 matched current Tank Side Cesium Removal (TSCR) facility prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitanate (CST) bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 40% Cs breakthrough after processing ~1700 BVs of feed; the 50% Cs breakthrough was extrapolated from the breakthrough data to occur at 1873 BVs. Testing compared to previous AP-101 and AP-107 testing at 16 °C showed ~300 BV increases in volume processed to reach the WAC limit for both lead and lag columns. The increase in capacity was determined to be due to the significantly lower K concentration in the AN-107 compared to the other tank waste matrices. A comparison in breakthrough curves for the three tests indicated slightly slower kinetic behavior in the AN-107, with variations in feed matrices (high organic complexants) likely responsible for the deviation. The Cs effluent from the lag column reached the WAC limit after processing 1097 BVs. Anticipating this breakthrough point, the polish column was preemptively installed around 900 BVs. Cs breakthrough from the lag column began at 500 BVs, reaching 3.06×10 0 µCi/mL, or 2.6 % Cs breakthrough, after processing all 1700 BVs of feed. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ion Exchange Processing of AW-105 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System

The Hanford Site stores an estimated 56 million gallons of mixed radioactive and chemically hazardous waste in large underground tanks. In support of the Direct Feed Low-Activity Waste (DFLAW) Program for expediting Hanford tank waste supernate treatment, laboratory-scale ion exchange processing using prototypic unit operations was conducted on AW-105 tank waste at the Pacific Northwest National Laboratory Radiochemical Processing Laboratory. This report describes the small-scale ion exchange testing with 9.2 L of diluted and filtered supernate from Tank 241-AW-105 (hereafter referred to as AW-105) at 16 °C (62 °F). One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AW-105 tank waste to meet this criterion, only 0.225% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 445. Testing with AW-105 matched current Tank Side Cesium Removal (TSCR) facility prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitanate (CST) bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 83% Cs breakthrough after processing ~1500 BVs of feed; the 50% Cs breakthrough was interpolated from the breakthrough data and occurred at 1041 BVs. Despite the AW-105 having a significantly higher K concentration (0.55 M compared to 0.10 M), testing compared to previous AP-107 ion exchange column testing at 16 °C showed no difference in BVs processed to reach the WAC on the lead column and only an approximate ~20 BV decrease in volume processed to reach the WAC limit on the lag column. The negligible differences in capacity despite the 5x concentration differences in K was determined to be due to the significantly lower NO3 concentration in the AW-105 supernate compared to the AP-107 tank waste matrix. A comparison in breakthrough curves for the two tests also indicated slightly faster kinetic behavior in the AW-105, with the variations in feed matrices (lower NO3 concentration) likely responsible for the deviation. The Cs effluent from the lag column reached the WAC limit after processing 772 BVs. Anticipating this breakthrough point, the polish column was preemptively installed around 675 BVs. Cs breakthrough from the lag column began at 300 BVs, reaching 1.10×10 1 µCi/mL, or 14.13 % Cs breakthrough, after processing all 1500 BVs of feed. The polish column processed nominally 830 BVs and reached 2.10×10 -1 µCi/mL, or 0.27 % Cs breakthrough at the conclusion of the test. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Using Ensemble Short-Term Initialized Coupled NASA GEOS5 Climate Model Integrations to Study Convective Bias Growth

The representation of convective processes, particularly deep convection in the tropics, remains a persistent problem in climate models. In fact structural biases in the distribution of tropical rainfall in the CMIP5 models is hardly different than that of the CMIP3 versions. Given that regional climate change at higher latitudes is sensitive to the configuration of tropical forcing, this persistent bias is a major issue for the credibility of climate change projections. In this study we use model output from integrations of the NASA Global Earth Observing System Five (GEOS5) climate modeling system to study the evolution of biases in the location and intensity of convective processes. We take advantage of a series of hindcast experiments done in support of the US North American Multi-Model Ensemble (NMME) initiative. For these experiments a nine-month forecast using a coupled model configuration is made approximately every five days over the past 30 years. Each forecast is started with an updated analysis of the ocean, atmosphere and land states. For a given calendar month we have approximately 180 forecasts with daily means of various quantities. These forecasts can be averaged to essentially remove "weather scales" and highlight systematic errors as they evolve. Our primary question is to ask how the spatial structure of daily mean precipitation over the tropics evolves from the initial state and what physical processes are involved. Errors in parameterized convection, various water and energy fluxes and the divergent circulation are found to set up on fast time scales (order five days) compared to errors in the ocean, although SST changes can be non-negligible over that time. For the month of June the difference between forecast day five versus day zero precipitation looks quite similar to the difference between the June precipitation climatology and that from the Global Precipitation Climatology Project (GPCP). We focus much of our analysis on the influence of SST gradients, associated PBL baroclinicity enabled by turbulent mixing, the ensuing PBL moisture convergence, and how changes in these processes relate to convective precipitation bias growth over this short period.

Cohen, Charlie↗

A Test for Characterizing Delamination Migration in Carbon/Epoxy Tape Laminates

A new test method is presented for the purpose of investigating migration of a delamination between neighboring ply interfaces in fiber-reinforced, polymer matrix tape laminates. The test is a single cantilever beam configuration consisting of a cross-ply laminate with a polytetrafluoroethylene insert implanted at the mid-plane and spanning part way along the length of the specimen. The insert is located between a 0- degree ply (specimen length direction) and a stack of four 90-degree plies (specimen width direction). The specimen is clamped at both ends onto a rigid baseplate and is loaded on its upper surface via a piano hinge. Tests were conducted with the load-application point located on the intact portion of the specimen in order to initiate delamination growth onset followed by migration of the delamination to a neighboring 90/0 ply interface by kinking through the 90-degree ply stack. Varying this position was found to affect the distance relative to the load-application point at which migration initiated. In each specimen, migration initiated by a gradual transition of the delamination at the 0/90 interface into the 90-degree ply stack. In contrast, transition of the kinked crack into the 90/0 interface was sudden. Fractography of the specimens indicated that delamination prior to migration was generally mixed mode-I/II. Inspection of the kink surface revealed mode-I fracture. In general, use of this test allows for the observation of the growth of a delamination followed by migration of the delamination to another ply interface, and should thus provide a means for validating analyses aimed at simulating migration.

Ratcliffe, James G.↗

Characterizing Delamination Migration in Carbon/Epoxy Tape Laminates

A new test method is presented for the purpose of investigating migration of a delamination between neighboring ply interfaces in fiber-reinforced, polymer matrix tape laminates. The test is a single cantilever beam configuration consisting of a cross-ply laminate with a polytetrafluoroethylene (PTFE) insert implanted at the mid-plane and spanning part way along the length of the specimen. The insert is located between a 0-degree ply (specimen length direction) and a stack of four 90-degree plies (specimen width direction). The specimen is clamped at both ends onto a rigid baseplate and is loaded on its upper surface via a piano hinge. Tests were conducted with the load-application point located on the intact portion of the specimen in order to initiate delamination growth onset followed by migration of the delamination to a neighboring 90/0 ply interface by kinking through the 90- degree ply stack. Varying this position was found to affect the distance relative to the load-application point at which migration initiated. In each specimen, migration initiated by a gradual transition of the delamination at the 0/90 interface into the 90- degree ply stack. In contrast, transition of the kinked crack into the 90/0 interface was sudden. Fractography of the specimens indicated that delamination prior to migration was generally mixed mode-I/II. Inspection of the kink surface revealed mode-I fracture. In general, use of this test allows for the observation of the growth of a delamination followed by migration of the delamination to another ply interface, and should thus provide a means for validating analyses aimed at simulating migration.

Ratcliffe, James G.↗

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor

Numerical simulations of reacting and non-reacting flows within a scramjet combustor configuration experimentally mapped at the University of Virginia s Scramjet Combustion Facility (operating with Configuration A ) are described in this paper. Reynolds-Averaged Navier-Stokes (RANS) and hybrid Large Eddy Simulation / Reynolds-Averaged Navier-Stokes (LES / RANS) methods are utilized, with the intent of comparing essentially blind predictions with results from non-intrusive flow-field measurement methods including coherent anti-Stokes Raman spectroscopy (CARS), hydroxyl radical planar laser-induced fluorescence (OH-PLIF), stereoscopic particle image velocimetry (SPIV), wavelength modulation spectroscopy (WMS), and focusing Schlieren. NC State's REACTMB solver was used both for RANS and LES / RANS, along with a 9-species, 19- reaction H2-air kinetics mechanism by Jachimowski. Inviscid fluxes were evaluated using Edwards LDFSS flux-splitting scheme, and the Menter BSL turbulence model was utilized in both full-domain RANS simulations and as the unsteady RANS portion of the LES / RANS closure. Simulations were executed and compared with experiment at two equivalence ratios, PHI = 0.17 and PHI = 0.34. Results show that the PHI = 0.17 flame is hotter near the injector while the PHI = 0.34 flame is displaced further downstream in the combustor, though it is still anchored to the injector. Reactant mixing was predicted to be much better at the lower equivalence ratio. The LES / RANS model appears to predict lower overall heat release compared to RANS (at least for PHI = 0.17), and its capability to capture the direct effects of larger turbulent eddies leads to much better predictions of reactant mixing and combustion in the flame stabilization region downstream of the fuel injector. Numerical results from the LES/RANS model also show very good agreement with OH-PLIF and SPIV measurements. An un-damped long-wave oscillation of the pre-combustion shock train, which caused convergence problems in some RANS simulations, was also captured in LES / RANS simulations, which were able to accommodate its effects accurately.

Fulton, Jesse A.↗

The Multiple Use Plug Hybrid for NanoSats (MUPHyN) Miniature Thruster

The Multiple Use Plug Hybrid (for) Nanosats is a prototype thruster is being developed to fill a niche application for NanoSat-scale spacecraft propulsion. When fully developed, the MUPHyN thruster will provide an effective and low-risk propulsive capability that could enable multiple NanoSats to be independently re-positioned after deployment from a parent launch vehicle. Because the environmentally benign, chemically-stable propellants are mixed only within the combustion chamber after ignition and the flow rate of the fuel is determined by a pyrolysis mechanism that is nearly independent of pressure or fuel grain defects, the system is inherently safe and can be piggy-backed near a secondary payload with little or no overall mission risk increase to the primary payload. The MUPHyN thruster uses safe-handling and inexpensive nitrous oxide (N2O) and acrylonitrile-butadiene-styrene (ABS) as propellants. Fused Deposition Modeling (FDM), a direct digital manufacturing process, is used to fabricate short-form-factor solid fuel grains with multiple helical combustion ports from ABS thermoplastic. This manufacturing process allows for the rapid development and manufacture of complex fuel grain geometries that are not possible to extrude or cast using conventional methods. This technology enables the construction of fuel grains with length-to-diameter ratios appropriate for incorporation into CubeSats while maintaining high surface areas and regression rates that allow the system to maintain a near optimal oxidizer to fuel ratio. The MUPHyN system provides attitude control torques by using secondary-injection thrust vectoring on a truncated aerospike nozzle. This configuration allows large impulse delta V burns and small impulse attitude control firings to be performed with the same system. To ensure survivability during extend duration burns, the MUPHyN incorporates a novel regenerative cooling design where the N2O oxidizer flows through a cooling path embedded in the aerospike nozzle before being injected into the combustion chamber near the nozzle base.

Eilers, Shannon D.↗

Integrating Commercial Off-The-Shelf (COTS) graphics and extended memory packages with CLIPS

This paper addresses the question of how to mix CLIPS with graphics and how to overcome PC's memory limitations by using the extended memory available in the computer. By adding graphics and extended memory capabilities, CLIPS can be converted into a complete and powerful system development tool, on the other most economical and popular computer platform. New models of PCs have amazing processing capabilities and graphic resolutions that cannot be ignored and should be used to the fullest of their resources. CLIPS is a powerful expert system development tool, but it cannot be complete without the support of a graphics package needed to create user interfaces and general purpose graphics, or without enough memory to handle large knowledge bases. Now, a well known limitation on the PC's is the usage of real memory which limits CLIPS to use only 640 Kb of real memory, but now that problem can be solved by developing a version of CLIPS that uses extended memory. The user has access of up to 16 MB of memory on 80286 based computers and, practically, all the available memory (4 GB) on computers that use the 80386 processor. So if we give CLIPS a self-configuring graphics package that will automatically detect the graphics hardware and pointing device present in the computer, and we add the availability of the extended memory that exists in the computer (with no special hardware needed), the user will be able to create more powerful systems at a fraction of the cost and on the most popular, portable, and economic platform available such as the PC platform.

Callegari, Andres C.↗

The striking influence of oxophilicity differences in heterometallic Mo–Mn oxide cluster reactions with water

Mixed-metal oxides have proven to be effective catalysts for the hydrogen evolution reaction, often outperforming either of the binary metal oxides. The reactivity of Mn x MoO y - (x = 1, 2; y = 3, 4) clusters toward H 2 O was investigated via time-of-flight mass spectrometry with clear evidence of cluster oxidation and corresponding H 2 production, specifically for Mn x MoO 3 - (x = 1, 2) clusters. Unlike previously studied MoxOy- clusters, which assumed a broad distribution of stoichiometries (typically x ≤ y ≤ 3x), both MnMoOy- and Mn 2 MoO y - preferentially formed y = 3 and 4 compositions in significant quantities under our source conditions. The electronic and molecular structures of the Mn x MoO y (x = 1, 2; y = 3, 4) anion and neutral clusters were probed with anion photoelectron spectroscopy and analyzed with supporting density functional theory calculations. Our studies suggest that both metal centers are involved in initial cluster–water complex formation, while Mo is the center that undergoes oxidation; hence, reactivity terminates when Mo is saturated in its highest oxidation state of +6. Across these four clusters, Mn remains relatively reduced and is stable in a high-spin electronic configuration. The preferential reactivity of water molecules toward the Mo center rather than Mn is rationalized by the much lower relative oxophilicity of Mn.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atmospheric Dust Modeling from Meso to Global Scales with the Online NMMB/BSC-Dust Model Part 2: Experimental Campaigns in Northern Africa

The new NMMB/BSC-Dust model is intended to provide short to medium-range weather and dust forecasts from regional to global scales. It is an online model in which the dust aerosol dynamics and physics are solved at each model time step. The companion paper (Perez et al., 2011) develops the dust model parameterizations and provides daily to annual evaluations of the model for its global and regional configurations. Modeled aerosol optical depth (AOD) was evaluated against AERONET Sun photometers over Northern Africa, Middle East and Europe with correlations around 0.6-0.7 on average without dust data assimilation. In this paper we analyze in detail the behavior of the model using data from the Saharan Mineral dUst experiment (SAMUM-1) in 2006 and the Bodele Dust Experiment (BoDEx) in 2005. AOD from satellites and Sun photometers, vertically resolved extinction coefficients from lidars and particle size distributions at the ground and in the troposphere are used, complemented by wind profile data and surface meteorological measurements. All simulations were performed at the regional scale for the Northern African domain at the expected operational horizontal resolution of 25 km. Model results for SAMUM-1 generally show good agreement with satellite data over the most active Saharan dust sources. The model reproduces the AOD from Sun photometers close to sources and after long-range transport, and the dust size spectra at different height levels. At this resolution, the model is not able to reproduce a large haboob that occurred during the campaign. Some deficiencies are found concerning the vertical dust distribution related to the representation of the mixing height in the atmospheric part of the model. For the BoDEx episode, we found the diurnal temperature cycle to be strongly dependant on the soil moisture, which is underestimated in the NCEP analysis used for model initialization. The low level jet (LLJ) and the dust AOD over the Bodélé are well reproduced. The remaining negative AOD bias (due to underestimated surface wind speeds) can be substantially reduced by decreasing the threshold friction velocity in the model.

photometers↗

Jet Noise Prediction for Chevron Nozzle Flows with Wall-Modeled Large-Eddy Simulation

This paper presents results from ongoing research on jet noise prediction with wall-modeled large eddy simulations (WMLES) performed with the LAVA computational framework. In particular the focus of this study is on mixing enhancements from a single stream chevron nozzle at Reynolds number of 1×10(exp 6). Although the concept of chevron nozzles to reduce jet noise is not new, our understanding of its impact on the overall noise is still not well understood. As a first step towards predicting noise reduction due to mixing enhancement concepts from first principles with WMLES, we simulate the noise generated by a single stream chevron nozzle SMC001 as well as its equivalent axisymmetric round jet nozzle SMC000. Detailed comparisons are made with a dedicated experiment conducted at NASA’s Glenn Research Center and good agreement was achieved. Two different approaches to introduce a turbulent boundary layer were compared but show no major impact on the results. This is especially important given future work were multi-stream nozzles are considered and extended costs of resolving the internal BL would have a bigger cost impact. A permeable Ffowcs Williams Hawkings (FWH) surface enclosing the jet is used to predict far-field noise from the simulated flow-field and excellent comparison to microphone array measurements is achieved within the resolved frequency band. Sensitivity of far-field noise predictions to grid resolution is systematically documented. Near-field comparisons to PIV data shows great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the shear-layer. The paper shows a cost reduction of an order of magnitude compared to an earlier study of this configuration due to algorithmic and software improvements and demonstrates that WMLES can be used as a cost-competitive approach for jet noise predictions.

CST↗

Critical and Californium Source-Driven Noise Analysis Subcritical Measurements with an Unreflected Cylindrical Tank of Mixed Uranium-Plutonium Nitrate Solution

The reported experiments were performed by Oak Ridge National Laboratory (ORNL) at the Battelle Northwest Laboratory’s (now Pacific Northwest National Laboratory) critical experiments facility at Hanford, Washington in 1981 and used 16 days of critical facility time, not including 10 days for setup and removal of ORNL equipment. These measurements were to assess the capability of the Cf source-driven noise analysis (CSDNA) method to measure the subcriticality (keff) of mixed U-Pu nitrate solutions. In addition to the CSNDA measurements, measurements were also performed near delayed criticality where CSDNA measurements cannot be performed. This report documents the experiments that were not reported at that time by presenting the ORNL experimental results and any online analysis performed during and shortly after the measurements. The mixed nitrate solution had a U concentration at 188 grams per liter (g/L), a Pu concentration of 280 g/L, free acid normality of 2.80, H ion molarity of 5.9 and a specific gravity of 1.754 g/cm 3 , a 240 Pu isotopic content of 7.981 wt. %, and a 235 U isotopic content of 0.724 wt. %. The stainless-steel tank for the solution had an inside diameter of 35.38 cm, an outside diameter of 35.53 cm, a height of 56.72 cm, and bottom thickness of 0.9525 cm. A Zircaloy pipe with a 3.1496 cm outside diameter, a 2.7788 cm inside diameter, and bottom thickness of 0.635 cm was available for insertion of the Cf source in the center of the fissile solution. The Cf source was also located at the outside surface of the tank (solution height varied from 10 to 53 cm) and in the center of the solution (solution height varied from 10 to 60.7 cm). The CSDNA measurements were not analyzed online to determine the subcritical neutron multiplication factors. At all subcritical states, the break frequency noise analysis data was fitted to obtain the prompt neutron decay constant. The neutron multiplication factors were determined for the two configurations of the measurements near delayed criticality. The subcritical neutron multiplication factors from the CNSDA measurements can be obtained with further analysis. However, the near delayed critical configuration, the prompt neutron decay constants, the count rates, and the measured cross and auto power spectral densities can be calculated directly for benchmarking. Much of data presented in this report are from ORNL notes—not in the ORNL logbooks. For the final benchmark analysis, the data from the Battelle Northwest Laboratory (which operated the critical facility in 1981) critical facility logbook should be consulted and be incorporated where appropriate. The purpose of this report is to document the experimental information for the measurements performed so that at a later date researchers could perform the required uncertainty and calculational analyses and documentation to use these data for an International Criticality Safety Benchmark Program (ICSBEP) or Nuclear Energy Agency benchmark. The data from these measurements are available from the ORNL Records Management Services Department, and the logbook is available from ICSBEP at Idaho National Laboratory. Preparation of the present report is part of a larger cooperative effort between Idaho National Laboratory (INL) and Oak Ridge National Laboratory (ORNL) to document more than 15 undocumented critical and subcritical experiments enumerated in ORNL/TM-2019/18 and performed by ORNL at ORCEF and other USDOE critical experiments facilities using more than 500 operational days of critical facility time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Primordial Inflation Polarization Explorer (PIPER)

The Primordial Inflation Polarization Explorer (PIPER) is it balloon-borne instrument designed to search for the faint signature of inflation in the polarized component of the cosmic microwave background (C-N-113). Each flight will be configured for a single frequency, but in order to aid in the removal of the polarized foreground signal due to Galactic dust, the filters will be changed between flights. In this way, the CMB polarization at a total of four different frequencies (200, 270, 350, and 600 GHz) will be, measured on large angular scales. PIPER consists of a pair of cryogenic telescopes, one for measuring each of Stokes Q and U in the instrument frame. Each telescope receives both linear orthogonal polarizations in two 32 x 40 element planar arrays that utilize Transition-Edge Sensors (TES). The first element in each telescope is a variable-delay polarization modulator (VPM) that fully modulates the linear Stokes parameter to which the telescope is sensitive. There are several advantages to this architecture. First, by modulating at the front of the optics, instrumental polarization is unmodulated and is therefore cleanly separated from source polarization. Second, by implementing this system with the appropriate symmetry, systematic effects can be further mitigated. In the PIPER design, many of the. systematics are manifest in the unmeasured linear Stokes parameter for each telescope and this can be separated from the desired signal. Finally, the modulation cycle never mixes the Q and U linear Stokes parameters, and thus residuals in the modulation do not twist the observed polarization vector. This is advantageous because measuring the angle of linear polarization is critical for separating the inflationary signal from other polarized components.

Chuss, David T.↗

Shivalik Plasma Device-I, a glow discharge device to study the collective dynamics of dusty plasma

Using a uniquely configured glow discharge-based Shivalik Plasma Device-I, we demonstrate a variety of collective phenomena in dusty plasma away from the glow discharge region. The cylindrical glass device produces plasma using parallel disc-shaped electrodes with a smaller anode size than the cathode. The dust microparticles are initially sprinkled over the grounded cathode. These particles acquire a significant negative charge upon plasma formation, resulting in their levitation due to the balance between the Coulomb force and gravity. The new device supports the levitation of a big-sized (10 × 8 × 5 cm 3 ) three-dimensional dust cloud over the glass surface. It contrasts the dusty plasma formations in-between electrodes reported earlier. As the discharge voltage varies from high to low, the dust cloud travels from over the glass surface to between the electrodes. A complex interplay of dust void over the cathode, a sharp density gradient, and gravity lead to self-excitation of collective dust phenomena. It includes dust density waves (phase velocity, v ph ~4 cm/s), dust cloud oscillation (frequency, f = 5 Hz), sheared flow (flow velocity, v f ~ 1 cm/s), and multiple-sized dust vortices. These dust vortices provided an excellent platform for studying turbulent mixing phenomena. The power spectrum analysis agreed with two-dimensional Kolmogorov power-law scaling. This is an ideal dusty plasma apparatus where we can create or move the dust cloud to a location of choice from the glass surface to in-between the electrodes and excite one among many collective dust dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/$\textit{t}$ algorithm, in particular the short iteration performance. Further, a proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li 0.5 La 0.5 TiO 3 . This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li 0.5 La 0.5 TiO 3 are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li 0.5 La 0.5 TiO 3 system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A plasma generator utilizing the high intensity ASTROMAG magnets

The magnet configuration for the proposed particle astrophysics magnet facility (ASTROMAG) on the space station includes a cusp magnetic field with an intensity of a few tesla. With these large magnets (or others) located in the outer ionosphere, many quite interesting and unique plasma physics experiments become possible. First there are studies utilizing the magnet alone to examine the supersonic, sub-Alfvenic interaction with the ambient medium; the scale length for the magnet perturbation is approx. 20 m. The magnetic field geometry when combined with the Earth's and their relative motion will give rise to a host of plasma phenomena: ring nulls, x-points, ion-acoustic and lower-hybrid shocks, electron heating (possible shuttle glow without a surface) launching of Alfvenwaves, etc. Second, active experiments are possible for a controlled study of fundamental plasma phenomena. A controlled variable species plasma can be made by using an RF ion source; use of two soft iron rings placed about the line cusp would give an adequate resonance zone (ECH or ICH) and a confining volume suitable for gas efficiency. The emanating plasma can be used to study free expansion of plasma along and across field lines (polar wind), plasma flows around the space platform, turbulent mixing in the wake region, long wavelength spectrum of convecting modes, plasma-dust interactions, etc.

Sullivan, James D.↗