Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Process optimization”

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,513 records · Page 84

FORCE-DISPATCHES Integration - Initial Demonstration

Integrated energy systems (IES) combine, in mutually beneficial ways, power from variable renewable energy sources and nuclear power plants (NPP) to improve economic viability under uncertain market and weather conditions. The open-source Framework for Optimization of Resources and Economics (FORCE) tool suite, developed at Idaho National Laboratory (INL), has enabled comprehensive modeling and simulation of IES. The capabilities within FORCE include grid portfolio optimization through the Holistic Energy Resource Optimization Network (HERON) and the transient process model analysis library HYBRID, among others. Continuous efforts and investments from the IES programs have been made to expand and improve the versatility of the FORCE toolset in fiscal year 2022. Code-coupling and cross-tool communication have been important methods for improving this versatility. This report focuses on an additional workflow in the HERON tool for capacity and dispatch stochastic optimization through integration with the external tool Design Integration and Synthesis Platform to Advance Tightly Coupled Hybrid Energy Systems (DISPATCHES). DISPATCHES was primarily developed by the National Energy Technology Laboratory, in collaboration with other national laboratories, which included INL, universities, and industry partners. It is coupled to a library of algebraic models for specific plant components, and to a framework for stochastic optimization different from that provided in the current Risk Analysis Virtual Environment (RAVEN)-running-RAVEN algorithm in HERON. HERON currently conducts stochastic optimization via an outer-inner loop: it optimizes over variable capacity on the outer loop, and at each step within the capacity parameter space, conducts an inner optimization over scenarios (of market signals, demand, and/or weather patterns) and hourly dispatch throughout a user-specified number of years. On the other hand, DISPATCHES conducts stochastic optimization via an “all-at-once” strategy in which capacity variables are optimized at the same level as dispatch variables, as all scenarios are considered at once. The latter method works especially well for projects of limited size and project length, as the necessary computational power and memory increases with the number of variables and scenarios. The new capability to use the DISPATCHES workflow in HERON enhances standalone simulations by leveraging FORCE tools—namely, the economic metrics from the Tool for Economic Analysis (TEAL) and reduced-order model (ROM) sampling from RAVEN. The initial demonstration of the DISPATCHES workflow simulates an existing nuclear-case flowsheet within the DISPATCHES repository—this models a NPP with a secondary revenue stream for hydrogen production. Electrical output from the plant is converted to hydrogen via a proton-exchange membrane (PEM) electrolyzer, hydrogen tanks are used for storage, and an additional turbine is added for hydrogen combustion. Continued work regarding this FORCE-DISPATCHES integration will include automatic generation of DISPATCHES models from HERON inputs, offering analysts the option of using either the RAVEN-runsRAVEN or DISPATCHES workflow to solve technoeconomic optimization problems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Designing a Robust MEA-Based Post-Combustion Carbon Capture Process with Capture Rate Guarantees

This work presents an application of the nonlinear two-stage robust optimization solver PyROS to the model-based design and operation of a monoethanolamine scrubbing process for CO<sub>2</sub> capture under epistemic uncertainty. Through this application, risk-averse process designs are successfully obtained for CO<sub>2</sub> capture targets ranging from 90% to over 99%. In particular, the risk-averse solutions for CO<sub>2</sub> capture targets of up to 98% are shown to be only marginally more expensive than their nominally optimal counterparts. Thus, the results demonstrate the utility of recently developed nonlinear robust optimization approaches for the solution of large-scale chemical process models under uncertainty.

20 FOSSIL-FUELED POWER PLANTS↗

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Optimization of artificial viscosity in production codes based on Gaussian Regression surrogate models

To accurately model flows with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity has to be introduced to convert kinetic energy into internal energy, thereby increasing the entropy across shocks. Determining the appropriate strength of the artificial viscosity is an art and strongly depends on the particular problem and experience of the researcher. The objective of this study is to pose the problem of finding the appropriate strength of artificial viscosity as an optimization problem and solve this problem using machine learning (ML) tools, specifically using surrogate models based on Gaussian Process regression and Bayesian analysis. We describe the optimization method and discuss various practical details of its implementation. The shock-containing problems for which we apply this method all have been implemented in the LANL code FLAG. First, we apply ML to find optimal values to isolated shock problems of different strengths. Second, we apply ML to optimize viscosity for a 1D propagating detonation problem based on Zel’dovich-von Neumann-Doring (ZND) detonation theory using a reactive burn model. We compare results for default (currently used values in FLAG) and optimized values of artificial viscosity for these problems demonstrating the potential for significant improvement in the accuracy of computations.

42 ENGINEERING↗

Additive manufacturing of AISI M2 tool steel by binder jetting (BJ): Investigation of microstructural and mechanical properties

The presented research demonstrates for the first time the successful processing of AISI M2 tool steel by binder jetting, a promising additive manufacturing technique capable of producing complex shapes with minimal residual stresses and isotropic properties. The optimal printing parameters were explored by varying processing parameters such as the binder saturation (45 %–105 %), binder set time (0 to 10 s), targeted bed temperature (50–60 °C), oscillator (2600–2750 rpm), recoater (20–28 mm/s), and roller speeds (200–300 rpm). Microstructural characterization and evaluation of mechanical properties of binder jetted parts were performed using x-ray diffraction (XRD), scanning electron microscopy (SEM), and energy dispersive spectroscopy (EDS) to study their chemical composition, powder morphology, microstructure, carbide morphologies, relative density, hardness, compressive strength, and ductility. Two powder sizes (5 and 10 μm) were used, and sintering was performed at varying temperatures (1270, 1280, and 1300 °C) and durations (60 and 120 min), followed by a furnace, air, and water cooling. An optimum hardness of ~970 HV was obtained when parts were sintered at 1270 °C for 60 min, followed by water quenching. Impressive compressive strength of ~ 3580 MPa was observed in the sample sintered at 1280 °C for 60 min duration, followed by air cooling. Furnace-cooled parts showed the highest density of ~95 %, whereas the relative density of air- and water-cooled parts varied between ~91 to 93.50 %, respectively. The microstructure of sintered samples revealed the formation of M 6 C stable carbide, M 2 C metastable carbide, MC as a secondary carbide, and a-Fe matrix, which contributed to the observed increase in mechanical properties.

36 MATERIALS SCIENCE↗

Development of a Technical, Economic, and Risk Assessment Tool for the Evaluation of Work Reduction Opportunities

Efficient and cost-effective operation of a nuclear power plant (NPP) is essential to ensuring long-term economical and safe operation. Multiple cost saving opportunities exist, referred to here as work reduction opportunities (WRO). These WROs reduce plant operating costs by employing various cost-effective strategies (e.g., implementation of modern technologies). Identifying and objectively screening WROs is an essential task to help reduce overall costs. However, there is no comprehensive framework for assessing WROs in the nuclear industry and evaluating their impact on plant operations. This report presents a novel framework for systematically evaluating WROs from a technical, economic, and risk perspective. As NPPs continue to add new technology and implement modernization strategies into their current processes, potential WROs are commonly identified. Although most WROs have the potential to reduce costs, not all opportunities will result in significant cost savings due to unforeseen risks, large implementation costs, or benefits that fall short of expectations. Examples of this can be the result of a technology that is not fully developed, uncertainty in the amount of cost reduction, or difficulties introducing a new process into an organization. These uncertainties can manifest several ways and can result in a WRO with limited cost savings or even a loss of investment. The framework developed emphasizes the importance of effectively screening the WROs from a holistic perspective to objectively identify inefficiencies and ensure a positive impact to the organization. This report presents the Technical, Economic, and Risk Assessment (TERA) as a key methodology for the screening and evaluation of potential WROs. The TERA framework begins with a screening phase where the process is examined through a hybrid combination of Lean Six Sigma and Integrated Operations for Nuclear (ION) guiding principles. This framework examines the current processes using the Lean Six Sigma SIPOC (Suppliers, Inputs, Process, Outputs, Consumers) methodology but retains the ION key elements of People, Technology, Process, and Governance as important factors to the nuclear decision-making process. By combining the principles of Lean Six Sigma and ION, the developed screening process is specific to the nuclear industry in that it systematically evaluates WROs in order to implement new technology that is comprehensively evaluated. The TERA begins by mapping current processes as they relate to WROs and examining the inefficiencies. Furthermore, the created process map can be used to identify and evaluate potential solutions. Using key performance indicators (KPIs), the TERA evaluates each area—technology, economics, and risk—for uncertainties and to perform cost-benefit analysis. The results of the TERA are important KPIs that allow for an evaluation of different processes and technology implementations. This assessment enables decision-makers to compare various WROs based on metrics and then make informed decisions for which opportunity to implement first. This research includes not only the creation of the TERA framework, but also the evaluation of its performance. A case study for screening potential WROs at Southern Nuclear Company is presented that utilizes the TERA methodology. Through the use of TERA, various WROs were screened, and the solutions evaluated for cost-benefit expectations. The report concludes by summarizing the overall effort and implications for utility modernization. The performance of the screening and TERA are discussed as well as the impact on the nuclear industry. The TERA process enables utilities to evaluate and inform investment decisions for WROs and mitigate any potential risks. Through this research, we provide utilities with a valuable framework to optimize operations, reduce costs, and drive continuous process improvement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-throughput homogenization of a quasi-Gaussian ultrafast laser beam using a combined refractive beam shaper and spatial light modulator

Efficiently shaping femtosecond, transverse Gaussian laser beams to flat-top beams with flat wavefronts is critical for large-scale material processing and manufacturing. Existing beam shaping devices fall short either in final beam homogeneity or efficiency. Here, we present an approach that uses refractive optics to perform the majority of the beam shaping and then uses a fine-tune device (spatial light modulator) to refine the intensity profile. For the beam that we selected, circularly asymmetric with intensity fluctuations, our method achieved a uniformity of 0.055 within 90% of the beam area at 92% efficiency. The optimization involved an iterative beam shaping process that converged to optimum within 10 iterations.

43 PARTICLE ACCELERATORS↗

Ink-based laser powder bed fusion of barium titanate

Barium titanate (BTO) is a lead-free functional ceramic that is widely used in sensors, transducers, and actuators. The increasing demand for complex geometries and quick design iterations motivates the use of additive manufacturing techniques. Laser powder bed fusion (L-PBF) additive manufacturing technique enables the net-shape or near-net-shape fabrication of metals and alloys. Ceramic L-PBF remains challenging due to their high melting temperature. Here, to address this challenge, a nanoparticle ink feedstock was incorporated into the L-PBF process to allow low power melting of BTO. Key print parameters (i.e. laser power, scanning speed, and hatching spacing) were optimized to fabricate millimeter-scale BTO parts. The printed BTO samples possessed an inhomogeneous equiaxed-columnar-cellular microstructure. Semiconducting current-voltage behavior was observed in as-printed samples, which was attributed to the creation of oxygen vacancies. The piezoelectric d 33 coefficient was measured to be 4.6 pC/N.

Dielectrics↗

Descriptor Aided Bayesian Optimization for Many-Level Qualitative Variables With Materials Design Applications

Abstract Engineering design often involves qualitative and quantitative design variables, which requires systematic methods for the exploration of these mixed-variable design spaces. Expensive simulation techniques, such as those required to evaluate optimization objectives in materials design applications, constitute the main portion of the cost of the design process and underline the need for efficient search strategies—Bayesian optimization (BO) being one of the most widely adopted. Although recent developments in mixed-variable Bayesian optimization have shown promise, the effects of dimensionality of qualitative variables have not been well studied. High-dimensional qualitative variables, i.e., with many levels, impose a large design cost as they typically require a larger dataset to quantify the effect of each level on the optimization objective. We address this challenge by leveraging domain knowledge about underlying physical descriptors, which embody the physics of the underlying physical phenomena, to infer the effect of unobserved levels that have not been sampled yet. We show that physical descriptors can be intuitively embedded into the latent variable Gaussian process approach—a mixed-variable GP modeling technique—and used to selectively explore levels of qualitative variables in the Bayesian optimization framework. This physics-informed approach is particularly useful when one or more qualitative variables are high dimensional (many-level) and the modeling dataset is small, containing observations for only a subset of levels. Through a combination of mathematical test functions and materials design applications, our method is shown to be robust to certain types of incomplete domain knowledge and significantly reduces the design cost for problems with high-dimensional qualitative variables.

Engineering↗

Optimizing hierarchical membrane/catalyst systems for oxidative coupling of methane using additive manufacturing

The advantage of a membrane/catalyst system in oxidative coupling of methane (OCM) compared to conventional reactive systems is that by introducing oxygen to the OCM catalytic sites through a membrane, parasitic gas phase reactions of O 2 (g), responsible for lowering product selectivity, can be avoided. The design and fabrication of membrane/catalyst systems has, however, been hampered by low volumetric chemical conversion rates, high capital cost, and difficulties in codesigning membrane and catalyst properties to optimize the performance. We solve these issues by developing a dual-layer additive manufacturing process, based on phase inversion, to design, fabricate and optimize a hollow fiber membrane/catalyst system for OCM. We demonstrate the approach though a case study using BaCe 0.8 Gd 0.2 O 3-δ (BCG) as the basis of both the catalyst and separation layers. We show that by using the manufacturing approach we can codesign the membrane thickness and catalyst surface area so that the flux of oxygen transport through the membrane and methane activation rates in the catalyst layer match each other. Here, we demonstrate that this “rate matching” is critical for maximizing the performance, with the membrane/catalyst system significantly overperforming conventional reactor designs under identical conditions.

36 MATERIALS SCIENCE↗

Synthesis of High-Quality Mg-MOF-74 Thin Films via Vapor-Assisted Crystallization

The unique features of metal-organic frameworks (MOFs), such as their large surface areas and diversity of structures, makes them suitable to a broad range of applications. One of the most critical enabling technologies for implementing MOF films is the preparation of dense and continuous films that would maximize the performance of the configurations. While formation of crystalline Mg-MOF-74 thin films presents significant challenges, demonstrating large-scale coatings of the MOF, which exhibits the highest CO2 uptake among MOFs at 1 bar and 25 °C, would be particularly beneficial for CO2- related sensing applications. In this work, we demonstrate for the first time high-quality Mg-MOF-74 films synthesized via a vapor-assisted crystallization (VAC) process. The VAC process described herein produces dense and highly crystalline layers of Mg-MOF-74 with a low coefficient of variation (COV) of film thickness below 7%. We first optimized a precursor solution for the VAC process, and then investigated the effects of synthesis temperature, time and droplet volume on the growth, crystallinity and thickness of Mg-MOF-74 films. The porosity of the film was assessed by measuring the CO2 uptake at room temperature and 1 bar. The obtained Mg-MOF-74 films possess a well-defined microporosity, as deduced from CO2 adsorption studies via quartz crystal microbalance (QCM) and comparison with bulk Mg MOF-74 reference data. Furthermore, CO2 cyclic adsorption-desorption experiments on the Mg-MOF-74

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of soil heterogeneity on soybean plant development and crop yield evaluated using time-series of UAV and ground-based geophysical imagery

Understanding the interactions among agricultural processes, soil, and plants is necessary for optimizing crop yield and productivity. This study focuses on developing effective monitoring and analysis methodologies that estimate key soil and plant properties. These methodologies include data acquisition and processing approaches that use unmanned aerial vehicles (UAVs) and surface geophysical techniques. In particular, we applied these approaches to a soybean farm in Arkansas to characterize the soil–plant coupled spatial and temporal heterogeneity, as well as to identify key environmental factors that influence plant growth and yield. UAV-based multitemporal acquisition of high-resolution RGB (red–green–blue) imagery and direct measurements were used to monitor plant height and photosynthetic activity. We present an algorithm that efficiently exploits the high-resolution UAV images to estimate plant spatial abundance and plant vigor throughout the growing season. Such plant characterization is extremely important for the identification of anomalous areas, providing easily interpretable information that can be used to guide near-real-time farming decisions. Additionally, high-resolution multitemporal surface geophysical measurements of apparent soil electrical conductivity were used to estimate the spatial heterogeneity of soil texture. By integrating the multiscale multitype soil and plant datasets, we identified the spatiotemporal co-variance between soil properties and plant development and yield. Our novel approach for early season monitoring of plant spatial abundance identified areas of low productivity controlled by soil clay content, while temporal analysis of geophysical data showed the impact of soil moisture and irrigation practice (controlled by topography) on plant dynamics. Our study demonstrates the effective coupling of UAV data products with geophysical data to extract critical information for farm management.

54 ENVIRONMENTAL SCIENCES↗

A Unique Case of the “Goldilocks Rule” in Solid-State Electrolytes: Two Are Good, Four Are Too Many

We report the syntheses of two new series of methacrylate monomers with different backbones: ureidopyrimidinone (PU) and boron-substituted urea pyrimidine (U), which enhance both the mechanical and electrochemical properties of the solid-state electrolyte (SSE) while improving the cycle life of lithium iron phosphate (LiFePO 4 , LFP) cathodes. The PU backbone is characterized by four hydrogen bonds (H-bonds), while the U backbone bears only two. Importantly, our research reveals that two H-bonds in these monomers are optimal; in contrast, four are excessive. The exceptional mechanical properties and processability of the SSE with the U series additives, resulting from the optimal H-bonds, were unexpectedly achieved. This leads to the establishment of a “Goldilocks rule” for additive design. The key strategies include: 1) reducing hydrogen-bonding (H-bonding) sites by changing pyrimidinone to pyrimidine and 2) shifting from intermolecular to intramolecular H-bonding and π−π bonding. Furthermore, this reduction in H-bonding also offers significant advantages in processability. The advancement can be extended to electrode fabrication, making the manufacturing of all-solid-state batteries more practical and efficient.

hydrogen bonding↗

QuCLEAR

SF-25-008 This software optimizes quantum circuits using a two-step process. The first step, Clifford Extraction, moves the Clifford gates to the end of the circuit and includes circuit optimizations. The second step, Clifford Absorption, then addresses these extracted gates on a classical computer. By reducing the number of gates that need to run on the quantum device, the software achieves a significantly lower gate count.

LIU, JI [Argonne National Laboratory (ANL), Argonn↗

Integrated RMP-based ELM-crash-control process for plasma performance enhancement during ELM crash suppression in KSTAR

The integrated Resonant Magnetic Perturbation (RMP)-based Edge-Localized Mode (ELM)-crash-control process aims to enhance the plasma performance during the RMP-driven ELM crash suppression, where the RMP induces an unwanted confinement degradation. In this study, the normalized beta ($\beta_\textrm{N}$) is introduced as a metric for plasma performance. The integrated process incorporates the latest achievements in the RMP technique to enhance $\beta_\textrm{N}$ efficiently. The integrated process triggers the n = 1 Edge-localized RMP (ERMP) at the L–H transition timing using the real-time Machine Learning (ML) classifier. The pre-emptive RMP onset can reduce the required external heating power for achieving the same $\beta_\textrm{N}$ by over 10% compared to the conventional onset. During the RMP phase, the adaptive feedback RMP ELM controller, demonstrating its performance in previous experiments, plays a crucial role in maximizing $\beta_\textrm{N}$ during the suppression phase and sustaining the $\beta_\textrm{N}$-enhanced suppression state by optimizing the RMP strength. The integrated process achieves $\beta_\textrm{N}$ up to ~2.65 during the suppression phase, which is ~10% higher than the previous KSTAR record but ~6% lower than the target of the K-DEMO first phase ($\beta_\textrm{N}$ = 2.8), and maintains the suppression phase above the lower limit of target $\beta_\textrm{N}$ (= 2.4) for ~4 s (~60$\tau_\textrm{E}$). In addition to $\beta_\textrm{N}$ enhancement, the integrated process demonstrates quicker restoration of the suppression phase and recovery of $\beta_\textrm{N}$ compared to the adaptive control with the n = 1 Conventional RMP (CRMP). The post-analysis of the experiment shows the localized effect of the ERMP spectrum in radial and the close relationship between the evolution of $\beta_\textrm{N}$ and the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cost optimization of high recovery single stage gap membrane distillation

Estimating the cost of membrane distillation (MD) processes at scale is confounded by the vast array of possible process configurations, module configurations, and operating conditions and the extreme sensitivity of MD costs to these parameters. This work develops a comprehensive MD cost optimization framework for minimizing the levelized cost of water (LCOW) as a function of feed salinity and temperature. Our single stage gap MD model fully captures solution properties, heat transfer, and mass transfer relationships in the modules to determine the cost optimal design as a function of multiple decision variables, including process configuration (feed location and connectivity between the hot and cold-side), the gap type (air, permeate, conductive), equipment sizes (membrane, heater, chiller, and heat exchanger units), and operating conditions (flow rates, hot and cold side temperatures, and heating and cooling duties). We estimate that the minimum LCOW of conductive gap MD ranges between $\$$10 and $\$$16 per m 3 for feed salinities of 25 g/L to 200 g/L and water recoveries of 30%–75%. Further, we update several heuristic assumptions for cost optimal system operation, including: increased thermal conductivity of the gap decreases system costs; cost optimal MD operates at the highest allowable temperature; and the LCOW is a strong function of brine salinity. Finally, we find that improving membrane properties and lowering heating and cooling costs will provide only modest reductions in costs (<25%).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Scientist-in-the-Loop Data Analytics Framework for Intelligent Simulation Model Tuning and Validation

This project developed a scientist-in-the-loop data analytics framework for intelligent simulation model tuning and validation, targeting the Weather Research and Forecasting (WRF) model and its solar energy variant, WRF-Solar-BNL. Domain experts, such as climate scientists, depend on large-scale numerical simulations for knowledge discovery and decision-making, yet the complexity of parameter tuning and the disconnect between automated optimization and domain expertise pose significant challenges. We extended an interactive visual analytics framework that enables domain experts to observe and intervene in the computational steering process by identifying disagreements between the simulation model, surrogate model, and the expert’s domain knowledge. Using Bayesian Optimization with Gaussian Process Regression as the surrogate model, our system allows users to probe parameter relationships, analyze correlation patterns, and adjust tuning parameters in real time. We developed use cases for solar irradiance forecasting through sustained collaboration with Brookhaven National Laboratory, resolving critical model configuration challenges and achieving meaningful reductions in prediction error. The project supported one PhD student, one MS student, and eight undergraduate students across three Data Science Capstone projects, resulting in one master’s thesis.

Dasgupta, Aritra [New Jersey Institute of Technolo↗