Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factor analysis”

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,009 records · Page 56

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗

Modeling and Simulation of Tank Pressure Control using Zero-Boiloff Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗

Advanced Nuclear Technology. Using Technology for Small Modular Reactor Staff Optimization, Improved Effectiveness, and Cost Containment, 3002007071

Through this grant, the U.S. Department of Energy (DOE) will review several functional areas within a nuclear power plant, including fire protection, operations and operations support, refueling, training, procurement, maintenance, site engineering, and others. Several functional areas need to be examined since there appears to be no single staffing area or approach that alone has the potential for significant staff optimization at new nuclear power plants. Several of the functional areas will require a review of technology options such as automation, remote monitoring, fleet wide monitoring, new and specialized instrumentation, human factors engineering, risk informed analysis and PRAs, component and system condition monitoring and reporting, just in time training, electronic and automated procedures, electronic tools for configuration management and license and design basis information, etc., that may be applied to support optimization. Additionally, the project will require a review key regulatory issues that affect staffing and could be optimized with additional technology input. Opportunities to further optimize staffing levels and staffing functions by selection of design attributes of physical systems and structures need also be identified. A goal of this project is to develop a prioritized assessment of the functional areas, and R&D actions needed for those functional areas, to provide the best optimization

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Finalizing Construction of a New Data Center at BNL

Computational science, data management and analysis have been key factors in the success of Brookhaven National Laboratory's scientific programs at the Relativistic Heavy Ion Collider (RHIC), the National Synchrotron Light Source II (NSLS-II), the Center for Functional Nanomaterials (CFN), and in biological, atmospheric, and energy systems science, Lattice Quantum Chromodynamics (LQCD) and Materials Science, as well as our participation in international research collaborations, such as the ATLAS experiment at Europe's Large Hadron Collider (LHC) at CERN (Switzerland) and the Belle II experiment at KEK (Japan). The construction of a new data center is an acknowledgement of the increasing demand for computing and storage services at BNL in the near term and enable the Lab to address the needs of the future experiments at the High-Luminosity LHC at CERN and the Electron-Ion Collider (EIC) at BNL in the long term.

99 GENERAL AND MISCELLANEOUS↗

Decarbonization Scenarios in the United States: Comparing Biofuels Growth in Two Models - GCAM and BSM

Scenarios for deep decarbonization rely on biomass for biofuels, biopower, and bioproducts, often including negative emissions via carbon capture and storage or utilization. Despite the prominence of biomass in many deep decarbonization pathways, critical questions remain about biomass allocation, effects of transportation electrification, the pace of growth, and implications for agriculture and land use. We address these questions through a unique comparison of carbon pricing effects on the growth of biomass utilization and its effects on land use in the United States by comparing results from a multisectoral integrated assessment model, the Global Change Analysis Model [GCAM], with results from a biomass-to-biofuels system dynamics model, the Biomass Scenario Model [BSM]. We contribute to model comparison efforts by analyzing the biomass deployment needed for a scenario consistent with a "Middle of the Road" Shared Socioeconomic Pathway [SSP2] and a representative concentration pathway of 2.6 W/m2. The GCAM scenarios solve for global equilibrium conditions that are consistent with this pathway, including demands for biomass across all economic sectors and representing bioenergy with carbon capture and storage as a technology option. The BSM scenarios assess those biomass and biofuel results for the United States and identify challenges associated with that pace and amount of expansion. In the scenario analysis, we harmonize key factors such as carbon price trajectory, domestic ethanol fuel demand, ethanol blending, and arable land availability, and vary them in both models. In GCAM, we vary the carbon price, transportation electrification, ethanol blending constraints, and arable land availability inputs and the value of the carbon in land; in BSM, in addition to directly inputting certain GCAM results, we vary the maximum rate of biorefinery construction, flexibility of feedstock types across conversion processes, and policy incentives such as tax credits and renewable identification number payments. The selected carbon price trajectory results in a rapid increase in biofuel production in the United States, reaching about 9.4 EJ/year in 2060 in the highest scenario analyzed in GCAM. Results differ between the two models in timing and ultimate quantity of biomass and biofuel production. GCAM biofuel quantities generally exceed BSM amounts because CCS is applied to biofuel pathways in GCAM, and because of differences in capacity expansion and related dynamics of land allocation, biomass production, and price dynamics. These dynamics include rapid biorefinery capacity expansion in high demand cases. To satisfy this biomass demand, GCAM rapidly equilibrates land allocation, but the BSM limits the rate at which this re-allocation can occur. A further contrast with the equilibrium approach in GCAM is that the BSM represents a delay between planting and harvesting woody biomass resources. As a result of these model contrasts, feedstock costs in BSM increase more than in GCAM, and the absence of CCS in the BSM also reduces the relative economic attractiveness of biofuels production. The bottlenecks, lags, and price increases also lead to potential for volatility in feedstock price and land allocation to biomass in the BSM. GCAM has more biomass production than BSM in all scenarios, partly because of the broader, economy-wide coverage of GCAM, in contrast to BSM's exclusive focus on biofuels. In both models, trends like those of biofuels production were observed for biomass production: minimal growth without a carbon price and policy incentives, and increases with a carbon price, particularly with carbon capture and storage, because the inputs assume that biopower and biofuels decrease greenhouse gas emissions. In high policy scenarios, biomass demand is high, and the consequent high biomass prices due to the land re-allocation bottleneck in the BSM limit biofuel production even if the biorefinery capacity is expanded. However, because biomass prices do not increase as much in the low policy scenario, growth is slower and the land-reallocation bottleneck no longer dominates, such that the effect of increased capacity can be seen. Across both the models, a change in assumptions from less to more land availability increases biofuel production in both GCAM and BSM, as the upward pressure on feedstock price and volatility are both reduced.

biofuels↗

Considering the Variability of Soiling in Long-Term PV Performance Forecasting: Preprint

This study presents the development of a methodology for evaluating the variability associated with soiling on long-term PV forecasting. Independent engineering firms typically build P50 forecasts for large PV plants through the use of the PVsyst software, where monthly soiling losses are one of many inputs to the P50 model. Subsequently, long-term performance distributions, or Pvalues, are constructed through a Monte Carlo analysis that includes various factors such as: satellite irradiance modeling uncertainty, uncertainty in the PVsyst model, and long-term irradiance variability. Often the PVsyst model uncertainty is increased to account for sites with significant soiling concerns but no systematic method has been presented in the literature to specifically include soiling variability within Pvalues. In this work soiling information from 16 sites in the U.S. Southwest are combined with 20 years of rainfall data to generate 20 years of energy production with soiling losses and then subsequently generate Pvalues. The results show that the spread of Pvalues (P1-P99) can increase from 0-13% when interannual soiling variability is included.

interannual variability↗

Use of Uncertainty in Modeling DFLAW Treatment at Hanford - 17115

Average or mean is understood by most; standard deviation (SD) is understood by few. Certainty is easy to model; uncertainty is not. The flowsheet simulation, described in this paper, placed uncertainties in the G2 model of the Hanford Tank Waste Treatment and Immobilization Plant's (WTP) Direct Feed Low-Activity Waste (DFLAW) process for performance and sensitivity evaluations. G2 is a dynamic material balance model that starts and proceeds with operating logic until all the feed is processed. Dynamic material balance is achieved by tracking the flows, storage, and change of all materials within the plant as time increments. Processing equipment and limitations of the flowsheet are modeled (e.g., equipment types, piping, volumes, flowrates, efficiencies, and physical and chemical environments that impact separations, phase equilibriums, and chemical reactions). Operating logic represents the rules and strategies of a running plant. The G2 model simulated both certainty and uncertainty scenarios for comparison purposes. For the uncertainty scenarios, process decisions were made using presumed values to simulate various operational measurement errors and process variabilities. Presumed means that the known G2 values are adjusted to replicate measurement errors and variabilities. These include: 1) volume measurement, 2) sample analysis, 3) melter decontamination factor (DF), 4) glass former weighing, and 5) glass former composition variations. The errors were randomly generated in accordance with SD values of measurement devices, laboratory methods, and glass former composition ranges. The G2 model runs give insight on the DFLAW process' controllability and the expected range of glass properties. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Droning on to Delivery: Examining the Energy Impacts of Using Drones for Moving Goods

The demand for fast, localized delivery has grown significantly in recent years. Whether for cheesy snacks, prepared food, medical solutions, or business deliveries, fast and efficient delivery is increasingly a demand and differentiator. Drone delivery in the freight sector offers to revolutionize last-mile logistics and improve services. This research analyzes the impacts of drone energy for delivery operations, aiming to compare various types of drones (large and small, rotary and VTOL) and various types of business methods. The project executes novel open-air and laboratory-based testing to look at the impacts of weights, operations, temperatures and weather conditions. The study combines this real-world experimental data with fleet optimization mathematical models to assess energy consumption of different scenarios as well as examine the minimum fleet size and additional battery requirements. These fleet optimization models compare energy between types of deployments with existing delivery methods. The model was also extended to look at mixed fleets of aerial drones and ground vehicles to accommodate different restrictions on drones or when weather prohibits their use. The analysis provided insights into factors such as drone design, payload weight, flight distance, weather conditions and operational parameters and impacts of each. It showed how to combine different types of vehicles to reduce energy and improve services. And it showed that drone speed, routing restrictions, and unfavorable weather significantly influence energy consumption. Our open experiment data and optimization models can assist stakeholders and industry in understanding drone package delivery and offers keys to improving deployment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Chemistry Based Criticality Safety Evaluation of the HM-Process Mixer Settlers

The criticality safety of Savannah River Site's H Canyon two cycle uranium extraction and purification process, called the HM process, is examined. Explicit Monte Carlo models of the multi-stage mixer-settlers are generated for the first time using the KENO-VI module of the SCALE 6.1 code package and spanning normal operating conditions and a wide range of potential upset conditions with stagewise compositions based directly on SEPHIS Mod-4 chemistry models of the same conditions. To achieve industrial scale processing, mixer-settlers are often not geometrically favorable to prevent criticality so other means must be used for criticality safety, namely control of fissile concentration and mass. Historically, the worst case upset is a reflux condition resulting from a stream flow rate, temperature, or composition that is off from the normal operating conditions. A reflux is considered to be the state in which the fissile material is not proceeding in the desired direction, i.e. is stripping rather than extracting or vice versa. The criticality safety analyses which set the current operational limits were conducted many years ago with less sophisticated codes and only loosely tied criticality to chemistry. Those analyses conservatively assumed that a reflux changed the behavior of the mixer-settler such that the fissile content collected continuously in only three stages of the mixer-settler, at a rate based on the SEPHIS modeling, until criticality occurred. This implies that under some condition the distribution coefficient of the fissile material either goes to infinity or zero, which is not physically achievable. In such a scenario, the operator would have some amount of time to react and correct the problem before sufficient fissile material built up in the mixer-settler and caused a criticality. In reality, an upset in a stream would cause the process to transition to a new, though undesired, equilibrium state of operation. The distribution coefficients would change, but would not diverge. In this work, aqueous and organic streams are modeled for both the mixer and the settler region of each stage of each bank using chemical distribution data from SEPHIS and solving the chemical equations for a wide range of operational parameters spanning normal and upset conditions. Conservative assumptions are made regarding enrichment and reflection and the whole bank multiplication factor is evaluated using discrete stage-wise compositions. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Methane Emissions Produced by Pneumatic Devices and Produced Water Tanks on Natural Gas Wellsites

In recent years, more scrutiny has been placed on the release of greenhouse gasses (GHGs) by the oil and gas industry in North America as climate concerns increase. In the coming years, stricter regulations proposed by the Environmental Protection Agency (EPA) concerning the release of GHGs will further direct attention to GHG emissions quantification, modeling, and overall reduction. Active natural gas wellsites employ production equipment that routinely emits methane during normal operation. Production equipment often includes, but is not limited to, gas production units (GPUs), pneumatic controllers (PCs) and actuators powered by produced natural gas, and produced water storage tanks. This work details direct methane emission measurements taken at active natural gas wellsites originating from intermittent pneumatic controllers and produced water storage tanks. Efforts to model methane flows originating from these sources are also detailed, and those results are compared to direct measurement. The direct methane emissions measurement campaign collected flow data from two active dry wellsites in the Marcellus shale region. These measurement campaigns captured, in total, several weeks of continuous data from these sites. Data included multiple pneumatic device vent flowrate channels, a single channel monitoring total emissions vented from the produced water storage tanks, and corresponding site weather data. Periodic composition measurements were taken from flows vented from produced water tanks. Computational modeling efforts attempted to mimic a given site’s methane emissions behavior given limited input parameters. Direct measurement data and simulation data yielded, through reduction and analysis, apparent methane emissions factors, statistics on pneumatic device actuation events, important relationships between GPU parameters and pneumatic device actuation events, produced water storage tank vent flows and flow compositions, and parameters most contributing to the magnitude and frequency of these fugitive methane emissions. Measurements, considering data from a 72hr period, produced PC emissions factors from 0.008 to 0.03 SCFH per device with an average device emitting 0.026 SCFH at the first site. The second site produced PC emissions factors spanning from 0.05 to 25.47 SCFH with an average device emitting 1.93 SCFH not including the 25.47 SCFH device, which was assumed to be malfunctioning. Two unique water tank emission measurements were made, 618 g CH4 per hour and 32.6 g CH4 per hour. On a water basis, the water tank emissions factors were 0.66 and 0.042 kg CH4 per bbl of water, respectively.

03 NATURAL GAS↗

Unsupervised machine learning based on non-negative tensor factorization for analyzing reactive-mixing

Analysis of reactive-diffusion simulations representing complex mixing processes requires a large number of independent model runs. For each high-fidelity model simulation, the model inputs are varied and the predicted mixing behavior is represented by temporal and spatial changes in species concentration. It is then required to discern how the model inputs (such as diffusivity, dispersion, anisotropy, and velocity field properties) impact the mixing process. This task is challenging and typically involves interpretation of large model outputs representing temporal and spatial changes of species concentration within the model domain. Yet, the task can be automated and substantially simplified by applying Machine Learning (ML) methods. Here, we present an application of an unsupervised ML method (called NTF k ) using Non-negative Tensor Factorization (NTF) coupled with a custom clustering procedure based on k-means to reveal the temporal and spatial features in product concentrations. An attractive and unique aspect of the proposed ML method is that it ensures the extracted features are non-negative, which is important to obtain a meaningful deconstruction of the mixing processes. The ML methodology is applied to a large set of high-resolution finite-element model simulations representing anisotropic reaction-diffusion processes in perturbed vortex-based velocity fields. The applied finite-element method ensures that spatial and temporal species concentration are always non-negative, even in the case of high anisotropic contrasts. The simulated reaction is a fast irreversible bimolecular reaction A+B→C , where species A and B react to form species C. The reactive-diffusion model input parameters that control mixing include properties of the velocity field (such as vortex structures), anisotropic dispersion, and molecular diffusion. We demonstrate the applicability of the ML feature extraction method to produce a meaningful deconstruction of model outputs to discriminate between different physical processes impacting the reactants, their mixing, and the spatial distribution of the product C. The introduced ML analysis allowed us to identify additive temporal and spatial features that characterize mixing behavior. The application of the proposed NTF k approach is not limited to reactive-mixing. NTFk can be readily applied to any observed or simulated datasets that can be represented as tensors (multi-dimensional arrays) and have separable latent signatures or features.

97 MATHEMATICS AND COMPUTING↗

Impact of Deferring Radiation Therapy in Patients With Epidermal Growth Factor Receptor–Mutant Non-Small Cell Lung Cancer Who Develop Brain Metastases

Purpose: To perform a retrospective analysis of patients with epidermal growth factor receptor (EGFR)-mutant lung adenocarcinoma who developed brain metastases (BM) to evaluate our hypothesis that the use of upfront EGFR–tyrosine kinase inhibitors (TKIs), and deferral of radiation therapy (RT), would result in inferior intracranial progression-free survival but similar overall survival (OS). Methods and Materials: Of 202 patients diagnosed with EGFR-mutant NSCLC between July 1, 2008, and December 31, 2014, 71 developed BM. Twenty-one patients were excluded owing to prior EGFR-TKI use, EGFR-TKI resistance mutation, failure to receive EGFR-TKI after whole-brain radiation therapy (WBRT)/stereotactic radiosurgery (SRS) or <6 months' follow-up. Of the remaining 50 patients, 17 received upfront EGFR-TKI followed by SRS or WBRT, 17 WBRT then EGFR-TKI, and 16 SRS followed by EGFR-TKI. Disease-specific-graded prognostic assessment was similar among all 3 groups. Results: The median OS was longer in the upfront RT group compared with the upfront EGFR-TKI group (34.1 vs 19.4 months; P=.01). On subgroup analysis, the SRS group had longer OS than the upfront EGFR-TKI group (58.4 vs 19.4 months; P=.01), but the WBRT group did not (29.9 vs 19.4 months; P=.09). Intracranial progression-free survival was improved in patients receiving upfront RT compared with those receiving upfront EGFR-TKI (37.9 vs 10.6 months; P<.001). Conclusions: The present study suggests that the use of upfront EGFR-TKI, and the deferral of SRS or WBRT, may result in inferior OS in patients with EGFR-mutant NSCLC who develop brain metastases. A prospective, multi-institutional, randomized trial of upfront EGFR-TKI with RT at intracranial progression versus upfront RT followed by EGFR-TKI is urgently needed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Assessing heterogeneity of patient and health system delay among TB in a population with internal migrants in China

Backgrounds The diagnostic delay of tuberculosis (TB) contributes to further transmission and impedes the implementation of the End TB Strategy. Therefore, we aimed to describe the characteristics of patient delay, health system delay, and total delay among TB patients in Shanghai, identify areas at high risk for delay, and explore the potential factors of long delay at individual and spatial levels. Method The study included TB patients among migrants and residents in Shanghai between January 2010 and December 2018. Patient and health system delays exceeding 14 days and total delays exceeding 28 days were defined as long delays. Time trends of long delays were evaluated by Joinpoint regression. Multivariable logistic regression analysis was employed to analyze influencing factors of long delays. Spatial analysis of delays was conducted using ArcGIS, and the hierarchical Bayesian spatial model was utilized to explore associated spatial factors. Results Overall, 61,050 TB patients were notified during the study period. Median patient, health system, and total delays were 12 days (IQR: 3–26), 9 days (IQR: 4–18), and 27 days (IQR: 15–43), respectively. Migrants, females, older adults, symptomatic visits to TB-designated facilities, and pathogen-positive were associated with longer patient delays, while pathogen-negative, active case findings and symptomatic visits to non-TB-designated facilities were associated with long health system delays (LHD). Spatial analysis revealed Chongming Island was a hotspot for patient delay, while western areas of Shanghai, with a high proportion of internal migrants and industrial parks, were at high risk for LHD. The application of rapid molecular diagnostic methods was associated with reduced health system delays. Conclusion Despite a relatively shorter diagnostic delay of TB than in the other regions in China, there was vital social-demographic and spatial heterogeneity in the occurrence of long delays in Shanghai. While the active case finding and rapid molecular diagnosis reduced the delay, novel targeted interventions are still required to address the challenges of TB diagnosis among both migrants and residents in this urban setting.

Sun, Ruoyao↗

Quantum assignments and intensity measures for methane between 1100 and 1800 per cm - A comparison between theory and experiment

The paper analyzes line positions and absolute line strengths of Blatherwick et al. (1979), based on moderately high resolution spectra of methane between 1100 and 1800 per cm, obtained using the Fourier transform spectrometer and the multiple-pass cold cell at the NASA Ames Research Center. Hamiltonian models are used to calculate theoretical relative line strengths, which, in combination with measured line strengths, yield integrated band strengths for the fundamentals v2 and v4. Ratios of calculated intensities to experimental intensities are analyzed, and the systematic deviations in the P- and R-branches of the v2 band are found to be represented by a Herman-Wallis type factor for Coriolis interactions. An analysis of ratios of calculated to experimental intensities of the v4 band yields a small correction to the Herman-Wallis factor.

Lutz, B. L.↗

Uncovering transcription factor and microRNA risk regulatory pathways associated with osteoarthritis by network analysis

Osteoarthritis (OA) is the most common form of joint disease. The development of inflammation have been considered to play a key role during the progression of OA. Regulatory pathways are known to play crucial roles in many pathogenic processes. Thus, deciphering these risk regulatory pathways is critical for elucidating the mechanisms underlying OA. We constructed an OA-specific regulatory network by integrating comprehensive curated transcription and post-transcriptional resource involving transcription factor (TF) and microRNA (miRNA). To deepen our understanding of underlying molecular mechanisms of OA, we developed an integrated systems approach to identify OA-specific risk regulatory pathways. In this study, we identified 89 significantly differentially expressed genes between normal and inflamed areas of OA patients. We found the OA-specific regulatory network was a standard scale-free network with small-world properties. It significant enriched many immune response-related functions including leukocyte differentiation, myeloid differentiation and T cell activation. Finally, 141 risk regulatory pathways were identified based on OA-specific regulatory network, which contains some known regulator of OA. The risk regulatory pathways may provide clues for the etiology of OA and be a potential resource for the discovery of novel OA-associated disease genes.

60 APPLIED LIFE SCIENCES↗

Optimizing Dynamic Timing Analysis with Reinforcement Learning

There are multiple factors involved in successfully manufacturing ASICIVLSI chips, and ensuring operational specifications are maintained throughout the design and manufacturing process is often challenging. Dynamic timing analysis (DTA) is the principal method used to validate that a manufactured chip complies to its design specifications. In DTA functionality of both synchronous and asynchronous designs are verified by applying input signals and checking for correct output signals. In complex designs where the number of input signal permutations is extremely large, the computing resources required to properly verify the functionality of a chip is prohibitive. In this paper, a strategy using reinforcement learning (RL) for reducing DTA time and resources in such cases is discussed. RL assisted DTA holds much promise in ensuring that VLSI chip design and functionality are fully and optimally verified.

42 ENGINEERING↗

Sun-view angle effects on reflectance factors of corn canopies

The effects of sun and view angles on reflectance factors of corn (Zea mays L.) canopies ranging from the six leaf stage to harvest maturity were studied on the Purdue University Agronomy Farm by a multiband radiometer. The two methods of acquiring spectral data, the truck system and the tower systrem, are described. The analysis of the spectral data is presented in three parts: solar angle effects on reflectance factors viewed at nadir; solar angle effects on reflectance factors viewed at a fixed sun angle; and both sun and view angles effect on reflectance factors. The analysis revealed that for nadir-viewed reflectance factors there is a strong solar angle dependence in all spectral bands for canopies with low leaf area index. Reflectance factors observed from the sun angle at different view azimuth angles showed that the position of the sensor relative to the sun is important in determining angular reflectance characteristics. For both sun and view angles, reflectance factors are maximized when the sensor view direction is towards the sun.

Ranson, K. J.↗

A Correction to the Stress–Strain Curve During Multistage Hot Deformation of 7150 Aluminum Alloy Using Instantaneous Friction Factors

Multistage stress–strain curve correction based on an instantaneous friction factor was studied for axisymmetric uniaxial hot compression of 7150 aluminum alloy. Experimental friction factors were calculated based on continuous isothermal axisymmetric uniaxial compression tests at various deformation parameters. Then, an instantaneous friction factor equation was fitted by mathematic analysis. After verification by comparing single-pass flow stress correction with traditional average friction factor correction, the instantaneous friction factor equation was applied to correct multistage stress–strain curves. The corrected results were reasonable and validated by multistage relative softening calculations. This research provides a broad potential for implementing axisymmetric uniaxial compression in multistage physical simulations and friction optimization in finite element analysis.

36 MATERIALS SCIENCE↗