Search NASA⌕ Search

SEARCH · Search NASA

Results for “Internal and external flows”

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.

87 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Laser Diode Analysis and Verification (Professional Report)

The goal of this project at Lawrence Livermore National Laboratory (LLNL) was focused on analyzing and verifying laser diodes to ensure that the diodes meet requirements for LLNL mission applications. The different stages of this testing were documented in a process flow map which included location color coding so that each step was clearly defined in scope and location of appropriate testing facility. Data were collected and analyzed and compared to minimum acceptable values to see if requirements were met. The process map documents the initial receipt, inspection, and testing of the laser diodes. Initial inspections started with Keyence Microscope imaging and then moved on to High Potential, Ramp, and Burst Testing. Data from the diode testing were processed through MATLAB and Python codes to verify various metrics such as slope efficiency, threshold current, back irradiance, and beam divergence met requirements. These metrics were then recorded in Excel summary reports. Approximately 95% of the laser diodes passed all tests. Presentations were given to Lawrence Livermore’s internal leadership team, an external partner, and to a Lab-wide audience. The data released for this report was constrained by information protection considerations of LLNL’s national security missions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Sensitivity and Uncertainty Quantification of Transition Scenario Simulations

This report documents the first collective attempt at developing and applying capabilities to quantify uncertainties, assess parametric sensitivities, and optimize multiple parameters and metrics in fuel cycle simulations generated by the SA&I Campaign. To do this, external codes that were designed to perform sensitivity analysis and uncertainty quantification (SA&UQ) needed to be coupled to the SA&I Campaign’s nuclear fuel cycle simulators (NFCS). In FY20, two approaches were pursued: 1) coupling Cyclus to an ORNL-internal code called MOT (Metaheuristic Optimization Tool) and 2) coupling DYMOND to the opensource SA&UQ tool kit Dakota. The primary objective of having these NFCS/SA&UQ coupled capabilities is to better inform DOE-NE and other stakeholders on the results generated from the NFCS. For a given set of fuel cycle strategies, policies, and technology assumptions that make up a fuel cycle scenario, these NFCS have traditionally been used by the SA&I Campaign to provide quantitative answers in terms of year-by-year mass flows, infrastructure requirements, costs, etc. With these newly developed coupled capabilities, the SA&I Campaign can now efficiently simulate hundreds or thousands of these scenarios, sample large ranges of parameters and assumptions, and use the unique features of the SA&UQ tools to process the data. This enables providing answers with known and propagated uncertainties, determining the sensitivity of important metrics to different parameters and assumptions, quantifying how much fuel cycle and technology parameters impact each other, and producing optimized fuel cycle strategies for single and multiple variables. To demonstrate these new capabilities, the Cyclus/MOT was used to model several scenarios ranging from simple fleet retirements to transitions to advanced reactors. Specifically, for a transition scenario from LWRs to SFRs and advanced LWRs, uncertainty quantification, sensitivity analysis, and optimization studies were applied to cases involving single and multiple parameter (input) and single and multiple metric (output) variations. In addition, a similar transition scenario was modeled to demonstrate how to optimize the reprocessing capacity parameter to minimize two performance metrics while taking into account uncertainties from two other parameters. Lastly, a depletion module based on SCALE/ORIGEN was added in Cyclus to simulate the third scenario that was designed to quantify the impact of the modeling assumption that all LWR used nuclear fuel have the same burnup. The newly developed DYMOND/Dakota capability was also applied to a transition scenario from the existing fleet to small modular reactors and fast reactors. This particular scenario involves not only explicit isotopic depletion via ORIGEN-2, but also includes multirecycling and utilizing the criticality search feature to determine the fresh fuel composition of recycled fuel, a feature unique to the DYMOND NFCS. A large database of simulations were run with 4 main parameters that were sampled: start date of reprocessing, reprocessing capacity, energy demand growth rate, and advanced reactor share of the fleet. The 4 main metrics were uranium consumption, enrichment requirements, waste generation, and levelized cost of electricity using data from the Cost Basis Report. The demonstrated SA&UQ results include those that inform on how to choose parameters to avoid “failed” scenarios, Sobol’ indices that inform on the importance of various parameters individually and synergistically, and Analysis of Variance (ANOVA) studies that decompose parameter ranges into groups and informs on whether variations are statistically significant.

Feng, B.↗
DOE OSTI2022

Critical heat flux on zircaloy and accident tolerant fuel cladding under prototypical conditions of pressurized and boiling water reactors

First of a kind high spatial and temporal resolution temperature measurements of Pressurized Water Reactor and Boiling Water Reactor simulated fuel pins (high heat flux cosine profile heaters) in annular rod type geometry during nucleate boiling, departure from nucleate boiling and re-wetting has been performed. Distributed fiber optic measurements of the cladding axial temperature with spatial resolutions down to 2.5 mm axially at up to three equally spaced azimuthal locations per rod have demonstrated the similarities and differences of the different accident tolerant fuel cladding materials. These experiments show the performance of accident tolerant cladding with a precise measurement of the location and evolution of dryout followed by the re-wetting phenomena (see supplementary time history of a critical heat flux event). Experimental parameters varied coolant mass fluxes from 1695 to 2712 kg/m 2 s, pressures from 10 to 20 MPa, and inlet subcooling from 10 to 55 °C. It was found that the coatings and different material have little effect on the value of the critical heat flux. Any discrepancies between the local critical heat flux verified among the cladding material are attributed to different departure from nucleate boiling phenomena captured during the experiment with the high-resolution optical fiber temperature data. Substantial differences in the performance of the different cladding material to survive a critical heat flux event were however observed by post-test X-ray and optical inspection. External evident damage and internal cracks to the simulated fuel pin having the bare zircaloy cladding were observed in regions associated with the occurrence of the critical heat flux for longer periods. Such damages were associated with zirconium oxidation, and not noticed on the Cr coated zircaloy or FeCrAl claddings. These results indicate that the use of FeCrAl or a Cr coated zircaloy as cladding in Pressurized Water Reactor and Boiling Water Reactor may increase the survivability of the cladding material for short transients.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Prognostic analysis of high-flow nasal cannula therapy and non-invasive ventilation in mild to moderate hypoxemia patients and construction of a machine learning model for 48-h intubation prediction—a retrospective analysis of the MIMIC database

Background This study aims to investigate the clinical outcome between high-flow nasal cannula (HFNC) and non-invasive ventilation (NIV) therapy in mild to moderate hypoxemic patients on the first ICU day and to develop a predictive model of 48-h intubation. Methods The study included adult patients from the MIMIC III and IV databases who first initiated HFNC or NIV therapy due to mild to moderate hypoxemia (100 < PaO2/FiO2 ≤ 300). The 48-h and 30-day intubation rates were compared using cross-sectional and survival analysis. Nine machine learning and six ensemble algorithms were deployed to construct the 48-h intubation predictive models, of which the optimal model was determined by its prediction accuracy. The top 10 risk and protective factors were identified using the Shapley interpretation algorithm. Result A total of 123,042 patients were screened, of which, 673 were from the MIMIC IV database for ventilation therapy comparison (HFNC n = 363, NIV n = 310) and 48-h intubation predictive model construction (training dataset n = 471, internal validation set n = 202) and 408 were from the MIMIC III database for external validation. The NIV group had a lower intubation rate (23.1% vs. 16.1%, p = 0.001), ICU 28-day mortality (18.5% vs. 11.6%, p = 0.014), and in-hospital mortality (19.6% vs. 11.9%, p = 0.007) compared to the HFNC group. Survival analysis showed that the total and 48-h intubation rates were not significantly different. The ensemble AdaBoost decision tree model (internal and external validation set AUROC 0.878, 0.726) had the best predictive accuracy performance. The model Shapley algorithm showed Sequential Organ Failure Assessment (SOFA), acute physiology scores (APSIII), the minimum and maximum lactate value as risk factors for early failure and age, the maximum PaCO 2 and PH value, Glasgow Coma Scale (GCS), the minimum PaO 2 /FiO 2 ratio, and PaO 2 value as protective factors. Conclusion NIV was associated with lower intubation rate and ICU 28-day and in-hospital mortality. Further survival analysis reinforced that the effect of NIV on the intubation rate might partly be attributed to the other impact factors. The ensemble AdaBoost decision tree model may assist clinicians in making clinical decisions, and early organ function support to improve patients’ SOFA, APSIII, GCS, PaCO 2 , PaO 2 , PH, PaO 2 /FiO 2 ratio, and lactate values can reduce the early failure rate and improve patient prognosis.

Fu, Wei↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
DOE OSTI2024

Identification of Dynamic Force Coefficients for an Additively Manufactured Hermetic Squeeze Film Bearing Support Damper Utilizing a Pass-Through Channel

Abstract The following paper presents breakthrough experimental results for a new hermetic squeeze film damper (HSFD) concept that is integrally designed within an externally pressurized tilting-pad radial gas bearing support. The flexibly damped gas bearing module was designed for a 7.2″ (183 mm) diameter shaft and fabricated using direct metal laser melting (DMLM); also known as additive manufacturing. The bearing and HSFD were sized based on ongoing studies for oil-free supercritical carbon dioxide (sCO2) power turbines in the 8.5 MW–10 MW power range. The development of the new damper concept was motivated by past dynamic testing on HSFD, which generated frequency-dependent stiffness and damping force coefficients. In efforts to eliminate the frequency dependency, a new HSFD architecture was conceived that adds accumulator volumes and a pass-through channel to previously conceived HSFD flow network designs. The other motivation for the work is the need to develop a cost-effective and reliable oil-free bearing technology that is scalable to large power turbomachinery applications. There were several objectives for the following work. The first objective was to successfully design and fabricate a single piece bearing-damper using additive manufacturing, while dimensionally controlling critical design features. The paper discusses the manufacturing steps and shows cut-ups that reveal adequate clearance control capability with internal damper clearances. The second objective was to perform experimental testing with the new HSFD design in efforts to extract stiffness and damping coefficients for excitation frequencies within 20–160 Hz and peak vibration amplitudes between 0.25 mils (6.35 microns) to 1 mil (25.4 microns). The test results for a single HSFD bearing module indicated that the design modifications to the HSFD architecture were successful in eliminating nearly all the frequency dependencies for the stiffness force coefficient. The dynamic tests yielded a stiffness coefficient that varied between 112 klb/in. (19.6 MN/m) and 96 klb/in. (16.8 MN/m). The damping force coefficient however, exhibited relatively more variation with frequency with values residing between 175 lb-s/in. (31 kN-s/m) to 214 lb-s/in. (37 kN-s/m). Finally, the paper advances a three-dimensional fluid-structure interaction (FSI) model using transient finite element analysis (FEA) coupled to a computational fluid dynamics (CFD) model. The FSI analysis performed between 20 Hz and 80 Hz was used to predict the stiffness and damping of the HSFD using a quarter-section model of the damper. The FSI analysis was able to support test results by showing only a 6–7.4% change in the magnitude of force coefficients. Stiffness predictions agree reasonably well with experiments whereas damping is underpredicted.

Engineering↗
DOE OSTI2022

Modeling the topographic influence on aboveground biomass using a coupled model of hillslope hydrology and ecosystem dynamics

Abstract. Topographic heterogeneity and lateral subsurface flow at the hillslope scale of ≤1 km may have outsized impacts on tropical forest through their impacts on water available to plants under water-stressed conditions. However, vegetation dynamics and finer-scale hydrologic processes are not concurrently represented in Earth system models. In this study, we integrate the Energy Exascale Earth System Model (E3SM) land model (ELM) that includes the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), with a three-dimensional hydrology model (ParFlow) to explicitly resolve hillslope topography and subsurface flow and perform numerical experiments to understand how hillslope-scale hydrologic processes modulate vegetation along water availability gradients at Barro Colorado Island (BCI), Panama. Our simulations show that groundwater table depth (WTD) can play a large role in governing aboveground biomass (AGB) when drought-induced tree mortality is triggered by hydraulic failure. Analyzing the simulations using random forest (RF) models, we find that the domain-wide simulated AGB and WTD can be well predicted by static topographic attributes, including surface elevation, slope, and convexity, and adding soil moisture or groundwater table depth as predictors further improves the RF models. Different model representations of mortality due to hydraulic failure can change the dominant topographic driver for the simulated AGB. Contrary to the simulations, the observed AGB in the well-drained 50 ha forest census plot within BCI cannot be well predicted by the RF models using topographic attributes and observed soil moisture as predictors, suggesting other factors such as nutrient status may have a larger influence on the observed AGB. The new coupled model may be useful for understanding the diverse impact of local heterogeneity by isolating the water availability and nutrient availability from the other external and internal factors in ecosystem modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

In-situ Photophysical Monitors and Corrective Algorithms for Photovoltaic Film Deposition and Rapid Thermal Processing in Scalable Roll-to-Roll Manufacturing

This project evaluated three photophysics-based approaches for in situ process monitoring of the scalable deposition of perovskite photovoltaic films. Based on stakeholder input and the results of laboratory trials and literature analysis, two techniques, angle resolved elastic light scattering (ARLS) and time-resolved photoluminescence (TRPL) were prototyped and demonstrated with a mixture of simulation, offline and online tests. Also, a hybrid gas flow and humidity control enclosure was constructed and demonstrated to stabilize the ambient conditions around the slot die deposition process and were demonstrated to meet the target humidity control in situ. Due in part to a delayed start and COVID-related lab restrictions and personnel changes, some of the milestones were not met. In spite of these challenges, however, both techniques were prototyped through several iterations and demonstrated significant promise for rapid in line perovskite characterization within the targeted characterization time windows and these techniques are under further development for use in DE-FOA000235:2357-1586. The TRPL system was fully mounted and demonstrated on the scaled slot die roll to roll and sheet coating systems at the Washington Clean Energy Testbeds (WCET). Motivated by significant internal and stakeholder interest, the successful prototyping achieved here, the ongoing need for in-situ characterization for the scale-up of perovskite manufacturing, and the benefits of in situ photophysical characterization for future funded process development research, we are continuing with the principle aims of the project to develop effect in situ photophysical probes for fast analysis of the optoelectronic and structural quality of sheet based and roll to roll process perovskites. TRPLS and ARLS probes proposed in the project are to be completed with external support from the Washington Research Foundation at the WCET. We are currently actively engaged in the acquisition of components and plan to complete and test the probes within three months. It is our intention that these probes will then be used in ongoing perovskite development with our research and commercialization partners including First Solar and other member of the US MAP consortium in future funded development work.

14 SOLAR ENERGY↗
DOE OSTI2022

Improved Particle Swarm Algorithm Using Rubik’s Cube Topology for Bilevel Building Energy Transaction

Following the rapid growth of distributed energy resources (e.g., renewables, battery), localized peer-to-peer energy transactions are receiving more attention for multiple benefits, such as reducing power loss and stabilizing the main power grid. To promote distributed renewables locally, the local trading price is usually set to be within the external energy purchasing and selling price range. Consequently, building prosumers are motivated to trade energy through a local transaction center. This local energy transaction is modeled in bilevel optimization game. A selfish upper level agent is assumed with the privilege to set the internal energy transaction price with an objective of maximizing its arbitrage profit. Meanwhile, the building prosumers at the lower level will response to this transaction price and make decisions on electricity transaction amount. Therefore, this non-cooperative leader-follower trading game is seeking for equilibrium solutions on the energy transaction amount and prices. Additionally, a uniform local transaction price structure (purchase price equals selling price) is considered here. Aiming at reducing the computational burden from classical Karush–Kuhn–Tucker (KKT) transformation and protecting the private information of each stakeholder (e.g., building), swarm intelligence-based solution approach is employed for upper level agent to generate trading price and coordinate the transactive operations. On one hand, to decrease the chance of premature convergence in global-best topology, Rubik’s Cube topology is proposed in this study based on further improvement of a two-dimensional square lattice model (i.e., one local-best topology-Von Neumann topology). Rotating operation of the cube is introduced to dynamically changing the neighborhood and enhancing information flow at the later searching state. Several groups of experiments are designed to evaluate the performance of proposed Rubik’s Cube topology-based particle swarm algorithm. The results have validated the effectiveness of proposed topology and operators comparing with global-best version PSO and Von Neumann topology-based PSO and its scalability on larger scale applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Rubik’s Cube Topology Based Particle Swarm Algorithm for Bilevel Building Energy Transaction

Following the rapid growth of distributed energy resources (e.g. renewables, battery), localized peer-to-peer energy transactions are receiving more attention for multiple benefits, such as, reducing power loss, stabilizing the main power grid, etc. To promote distributed renewables locally, the local trading price is usually set to be within the external energy purchasing and selling price range. Consequently, building prosumers are motivated to trade energy through a local transaction center. This local energy transaction is modeled in bilevel optimization game. A selfish upper level agent is assumed with the privilege to set the internal energy transaction price with an objective of maximizing its arbitrage profit. Meanwhile, the building prosumers at the lower level will response to this transaction price and make decisions on electricity transaction amount. Therefore, this non-cooperative leader-follower trading game is seeking for equilibrium solutions on the energy transaction amount and prices. In addition, a uniform local transaction price structure (purchase price equals selling price) is considered here. Aiming at reducing the computational burden from classical Karush-Kuhn-Tucker (KKT) transformation and protecting the private information of each stakeholder (e.g., building), swarm intelligence based solution approach is employed for upper level agent to generate trading price and coordinate the transactive operations. On one hand, to decrease the chance of premature convergence in global-best topology, Rubiks Cube topology is proposed in this study based on further improvement of a two-dimensional square lattice model (i.e., one local-best topology-Von Neumann topology). Rotating operation of the cube is introduced to dynamically changing the neighborhood and enhancing information flow at the later searching state. Several groups of experiments are designed to evaluate the performance of proposed Rubiks Cube topology based particle swarm algorithm. The results have validated the effectiveness of proposed topology and operators comparing with global-best version PSO and Von Neumann topology based PSO and its scalability on larger scale applications.

Feng, Xiaochun↗
DOE OSTI2023

Experimental study on the startup of the annular wick type heat pipe using fiber optical temperature measurement technique

This study used optical fiber-distributed temperature sensors to measure the internal and external temperature distributions of a water-cooled heat pipe. The sensor technology used in this study is fiber optical distributed temperature sensing, a distributed sensing technique based on the naturally occurring Rayleigh backscatter in optical fibers. This measurement technique provides maximum spatial resolution for static and semi-static applications. Using this sensor, the temperature distribution of the heat pipe's internal, external, vapor core, and the wick was measured with a spatial resolution of 0.65 mm, a sampling frequency of 40 Hz, and a temperature resolution of 0.1 °C. Through the measured temperature distribution database, the starting phenomenon, the effective length trend, and the limitation onset were observed. From the results, it is found that a high-temperature peak appears at the evaporator if a high initial power (75 W) is imposed on the heat pipe, even after the heat pipe approaches the normal operating status. The peak is not observed in a slower startup (30 W initial power then slowly increased to 75 W). It is also found that the temperature distributions and effective condenser length of the heat pipe highly depend on the cooling conditions. There are variations in the temperature according to the radial direction of the horizontal heat pipe due to gravity. In conclusion, lead and lag of the temperature evolution were observed at the onset of the operating limitations.

42 ENGINEERING↗
DOE OSTI2020

Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel

In this LDRD “Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel” a peer-reviewed lab-scale hypersonic wind tunnel (WT) was planned to be assembled in-house, along with time and space resolved and related flow and materials diagnostics to derive high quality measurements and the science that is able to inform the development of our modeling and simulation codes for the hypersonic regime. The baseline capability to achieve Mach 3 flow was scheduled to be completed in Q4 2020, in the first year of the LDRD. However due to the COVID-19 delays and travel restrictions starting in March, 2020 this task could not be completed to carry out the materials response studies in-house, or those planned with external partners that were engaged during this LDRD (NASA Langley, NASA Ames, Texas A&M university, and Virginia Tech). Still, using our hypersonic computational fluid simulation codes, many of the hypersonic components have been designed (for Mach 3, 5, and 7) and procured (or in procurement phase), and a new dedicated hypersonic laboratory was facilitized, 2 WT test sections have been manufactured, and a Mach 3 nozzle was 3D printed for prototyping. In parallel, unique materials diagnostics have been tested under static conditions and fully specified and awaiting validation in our partner facilities. Due to the new hypersonic initiative announced in January, 2020 we proposed a related scope increase in hypersonics research requiring energy interactions and new diagnostics well beyond what could be supported by an LDRD ER. A decision was made to end the current LDRD ER and formally propose a new LDRD SI that leverages the current effort, which will be starting in FY21. As a result of the uniqueness of this LDRD work, our early efforts have been encouraged and well received in the hypersonic research community, and is now part of multiple unsolicited internal and external project proposals expected to be funded in the coming years, which is encouraging to make LLNL a leader in experimentally probing extreme physics of local energetics critical to hypersonics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗