Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Optimization and modeling”

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

Techno-economic life cycle assessment of CO 2 -EOR operations towards net negative emissions at farnsworth field unit

Optimizations of CO 2 Water Alternating Gas(WAG)- systems with multi-objectives of incremental recovery and maximization of CO 2 storage are challenging. Here, the incorporation of a total Greenhouse gas (GHG) life cycle assessment is mostly ignored leading to inaccurate estimation of overall net carbon emissions of their operations. In this study, the effect of a total GHG life cycle assessment on a multi-objective CO 2 -WAG optimization with integrated techno-economic assessment (TEA) which factors carbon tax credit is conducted. A life cycle assessment (LCA) was conducted utilizing a 20 -year optimized post history matched data from a high fidelity reservoir simulation model. Using data generated from the optimum result, a techno-economic life cycle analysis was further conducted. The first scenario classified as the base model had an estimated 81% of purchased CO 2 sequestered. The results through a comprehensive techno-economic LCA model yielded a net estimate of 73% of purchased CO 2 . The optimized forecasted model which considered key operational and reservoir factors such as WAG ratio, injection rates and periods, and well specification resulted in an improved sequestration of 92% of purchased CO 2 . However, this also dropped to 84% after taking it through LCA. These results clearly indicate a significant amount of net CO 2 is not accounted for when operations are not analyzed through LCA. From the LCA, direct flaring volumes of CO 2 , energy consumption and efficiency of unit equipment were noticed to be the major causes of these reductions. Considering ten main sources of energy as source of energy generation, a comparative techno-eco LCA was conducted. The results confirmed a lower net volume and NPV for energy sources with higher carbon footprints and vice versa. Thus, a total LCA of CO 2 -WAG greatly influences net storage factor of purchased CO 2 and hence project NPV where tax credit/incentives per ton of CO 2 sequestered is considered. Although operational conditions are optimized for best results, there are significant factors that leads to minimization of net storage factor. This study therefore provides an insightful information for optimizing CO 2 -WAG multi-objectives to achieve minimum GHG emission.

02 PETROLEUM↗

A forecast-driven decision-making model for long-term operation of a hydro-wind-photovoltaic hybrid system

Hydro-wind-photovoltaic (PV) hybrid system has the potential to increase the integration of renewable energy sources into an existing grid. For the long-term operation of the system, due to the non-storable nature of wind and PV power, it is essentially to decide the optimal long-term carryover storage of cascade reservoirs. However, it remains a challenge due to high uncertainties of long-term forecasts and complicated hydraulic/electrical relationships between cascade reservoirs. Here in this study, a forecast-driven decision-making model is proposed for the hybrid system, which converts the multi-stage long-term operation process into a two-stage operation problem (including current stage and carryover stage) to avoid using longer-horizon forecast information with lower accuracy. First, the carryover stage energy surfaces (CESs) considering the forecast uncertainties of wind, solar and hydro resources are proposed to characterize carryover stage benefit quantitatively. Then a CESs-based forward decision-making optimization model is developed to guide the long-term operation of a hydro-wind-photovoltaic hybrid system. The applications in a hydro-wind-PV hybrid system of Yalong River basin results show that: compared with conventional operation, 1) power generation increases 9.03%; 2) in terms of the carryover storages control, the reservoir impounding and drawdown timing are delayed, and the drawdown depth is increased, which can be used to formulate better reservoir operation rules.

13 HYDRO ENERGY↗

Quantum Chemistry-Informed Active Learning to Accelerate the Design and Discovery of Sustainable Energy Storage Materials

Here we employed Density Functional Theory (DFT) to compute oxidation potentials of 1,400 homobenzylic ether molecules to search for the ideal sustainable redoxmer design. The generated data were used to construct an active learning model based on Bayesian optimization (BO) that targets candidates with desired oxidation potentials utilizing only a minimal number of DFT calculations. The active learning model demonstrated not only significant efficiency improvement over the random selection approach but also robust capability in identifying desired candidates in an untested set of 112,000 homobenzylic ether molecules. Our findings highlight the efficacy of quantum chemistry-informed active learning to accelerate the discovery of materials with desired properties from a vast chemical space.

25 ENERGY STORAGE↗

Wind Turbine Design Optimization for Hydrogen Production

To help meet the need for inexpensive green fuels, we are working on wind turbine design optimization specifically for hydrogen production. We have thus far achieved a 1.53% decrease in LCOH as compared to a turbine optimized for LCOE using the same code, design variables, and models. We accomplished this by optimizing some components of the wind turbine tower, rotor, and drivetrain design with hydrogen production and costs in the design loop.

hydrogen↗

FutureTense

Protective vaccines and reliable diagnostics are essential tools for controlling viral diseases. However, the efficacy of these tools can be diminished by mutations in viral genomes. The delay between the emergence of new viral strains and the redesign of vaccines and diagnostics allows for continued viral transmission. Is it possible to address this challenge by computationally predicting viral genome sequence evolution? Can we “future-proof” vaccines and diagnostics by targeting both current and anticipated future sequence variants? While predicting viral evolution is still an unsolved, “grand challenge” problem in biology, the large, and rapidly growing, number of SARS-CoV-2 genome sequences provide an opportunity to quantify the ability of machine learning to predict viral genome sequence evolution. Towards this end, we have developed a simple computational model for predicting viral evolution at the level of individual nucleotides. The key metric for quantifying the per-base, prediction accuracy for viral evolution is the Mann-Whitney U statistic (or, equivalently, the area under the receiver operator curve). Since the Mann-Whitney U statistic is not a differentiable function, existing deep leaning packages (like Pytorch and Keras/TensorFlow) are not useful, as they require that the accuracy metric/objective function be analytically differentiable with respect to the model parameters. To overcome this challenge, we have implemented custom software, “FutureTense”, that can train a machine learning model by maximizing the non-differentiable Mann-Whitney U statistic. This software trains a machine learning model by exploring along the direction of the discrete gradient of the Mann-Whitney U statistic in the model parameter space. Parallel computing and genome sequence-specific optimizations are used to accelerate model training. The resulting machine learning model learns the observed high C->U mutation rates in the SARS-CoV-2 genome (which are potentially induced by host defenses) and provides prediction accuracies that are significantly better than one would expect from random chance. While predicting viral evolution is still quite far from a solved problem, the surprising performance of this simple model gives hope that the accuracy of predicting viral genome evolution can be further increased by more sophisticated approaches.

Gans, Jason↗

Performing Numerical Analysis of Cybersecurity Options Using Dynamic Risk Analysis Tool EMRALD

Cyberattacks can have many different attack paths, durations, and goals. There are also many different mitigation options involving hardware, software, and/or humans. Considering a cyber threat should involve defense-in-depth methods and a quantitative or numerical evaluation of overall effectiveness against dynamic, time-dependent attacks to make cost and risk-informed decisions. Typical cyberattack modeling methods only provide a qualitative evaluation. The main areas of cybersecurity are confidentiality, integrity, and availability. For companies with cyber-physical systems such as advanced nuclear reactors, cyber-related safety is a requirement set by North American Electric Reliability and the U.S. Nuclear Regulatory Commission. They are also concerned about availability or reliability as a business case. As cyber threats are evolving to a business-for-hire structure, more attacks may focus on disrupting business success and reliability, causing financial and economic stability risk. Companies want to know business reliability and recovery from those threats, and that requires modeling physical behavior of the targets. Dynamic-state-based and Markov-based modeling provides a method for better cyber scenario modeling with different tools having issues such as state-base explosion. Dynamic modeling enables time and conditional features not found in other numerical evaluation methods. EMRALD (Event Modeling Risk Assessment using Lined Diagrams) is a dynamic risk analysis modeling and simulation tool and has features that reduce modeling issues. It has been used to model different time-dependent events including plant behavior and operator procedures. As a general modeling tool, EMRALD can also be used to model cyberattack scenarios with varying mitigation options and quantify effectiveness, producing numerical data for risk-informed decisions. This paper uses EMRALD to demonstrate that dynamic numerical risk analysis can be used for cyber threat modeling to provide insights for design decision-making and optimize defense strategies. Keywords: cyber modeling; cyber-physical systems; numerical cyber modeling

97 - MATHEMATICS AND COMPUTING↗

Wind Turbine Design Optimization for Hydrogen Production

To help meet the need for inexpensive green fuels, we are working on wind turbine design optimization specifically for hydrogen production. We have thus far achieved a 1.53% decrease in LCOH as compared to a turbine optimized for LCOE using the same code, design variables, and models . We accomplished this by optimizing some components of the wind turbine tower, rotor, and drivetrain design with hydrogen production and costs in the design loop.

design↗

Numerical Optimization of Spray-Guided Spark Assistance for Cold Idle Operation in a Heavy-Duty Gasoline Compression Ignition Engine

This article describes the results of a response surface model (RSM)-based numerical optimization campaign for spray-guided spark assistance at cold operations in a heavy-duty gasoline compression ignition (GCI) engine. On the basis of an earlier work on spark-assisted GCI cold combustion, a space-filling design of experiments (DoE) method was first undertaken to investigate a multitude of hardware design variables and engine operating parameters. The main design variables included the number of injector nozzles, fuel split quantities and injection timings, and spark timing. The objective variables were engine combustion efficiency (ŋ c ), maximum pressure rise rate (MPRR), and engine-out nitrogen oxide (NO x ) emissions. A total of 150 design candidates were automatically generated using the Sobol sequence method provided by the commercial software package, CAESES. Then, closed-cycle computational fluid dynamic (CFD) spark-assisted GCI simulations under cold idling operations were performed. The outcomes from the CFD-DoE design campaign were utilized to construct high-fidelity RSMs that allowed for further design optimization of the spark plug- and fuel injector-related design variables, along with fuel injection strategy parameters. A merit function with respect to objective variables was formulated with an appropriate weight assignment on each objective variable. Finally, the best design candidate was identified from the RSM-based optimization process and further validated in the CFD analysis. The best design candidate showed the potential to significantly improve combustion efficiency (ŋ c > 90%) over the baseline at cold idle while satisfying MPRR and NOx emissions constraints (MPRR < 5 bar/CAD and NO x < 4.5 g/kWh).

33 ADVANCED PROPULSION SYSTEMS↗

Design and Optimization of Processes for Recovering Rare Earth Elements from End-of-Life Hard Disk Drives

In this conference paper, we propose a superstructure-based approach to finding the optimal pathways for recovering rare earth elements in their commercialized rare earth oxide form from end-of-life HDDs. The proposed superstructure was modeled as a MILP optimization problem, selecting the net present value as the objective function. Whenever possible, costing data taken from the literature was used to inform this mode. However, due to the novelty of this research area data were often not available thus requiring the generation of flowsheets that were implemented in Aspen Plus. To establish the base case optimal result, projections for the number of EOL HDDs in the U.S. available for recycling and estimates of the projected rare earth oxide prices over the lifetime of the plant were used to inform the model. The model was then expanded to include the recycling of EOL HDDs generated prior to the beginning of plant production (period ranging from 2006 through 2024).

Laliwala, Chris↗

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗

SPEARS: A Database-Invariant Spectral modeling API

The Spectral Physics Environment for Advanced Remote Sensing (SPEARS) application programming interface (API) is a Python-based, line-by-line, local thermal equilibrium (LTE) spectral modeling code which is optimized for simultaneously synthesizing optical spectra from any combination of fundamental spectroscopic databases. In this article, we contribute two novel spectral modeling techniques to the scientific literature. First we describe how SPEARS integrates a physics-based collisional model for calculating pressure broadening in the absence of available broadening coefficients. With this collisional model implementation, a generalized approach to fundamental spectroscopic databases can be achieved across multiple databases. We also detail our adaptive grid mesh algorithm developed to make the code scalable for simulating large spectral bandwidths at high spectral fidelity using intuitive grid parameters. Here, we present comparisons to other modeling tools, experiments, and provide a discussion on the SPEARS user interface.

47 OTHER INSTRUMENTATION↗

Microgrid design and multi-year dispatch optimization under climate-informed load and renewable resource uncertainty

Microgrids are an increasingly popular solution to provide energy resilience in response to increasing grid dependency and the growing impacts of climate change on grid operations. However, existing microgrid models do not currently consider the uncertain and long-term impacts of climate change when determining a set of design and operational decisions to minimize long-term costs or meet a resilience threshold. In this paper, we develop a novel scenario generation method that accounts for the uncertain effects of (i) climate change on variable renewable energy availability, (ii) extreme heat events on site load, and (iii) population and electrification trends on load growth. Additionally, we develop a two-stage stochastic programming extension of an existing microgrid design and dispatch optimization model to obtain uncertainty-informed and climate-resilient energy system decisions that minimizes long-term costs. Use of sample average approximation to validate our two case studies illustrates that the proposed methodology produces high-quality solutions that add resilience to systems with existing backup generation while reducing expected long-term costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING↗

Optimization and Evaluation of Stochastic Unified Convection Using Single-Column Model Simulations at Multiple Observation Sites

We extend the previously developed stochastic unified convection scheme (UNICON) for shallow convection to deep convection by parameterizing the impact of mesoscale organized flow on updraft properties. The extended stochastic UNICON parameterizes thermodynamic properties of updrafts at the near-surface as a multivariate Gaussian distribution, where the variances of the distribution are the summation of variances from non-organized turbulence and mesoscale organized flow. The distribution of updraft radius is parameterized as a power-law distribution with a scale break which is parameterized as a linear function of the strength of mesoscale organized flow. The proposed parameterization is validated using a series of large-eddy simulations of deep convection. The free parameters introduced in the formulation of stochastic UNICON are optimized using 10 cases of single-column model simulations over the ocean. Stochastic UNICON with the optimized parameters significantly reduces the biases of thermodynamic profiles and surface precipitation rates simulated in the original UNICON for tropical convection cases. The simulation of the variation in anomalies of temperature and moisture associated with the Madden-Julian oscillation is also improved. The overall improvements in simulated thermodynamic profiles are found to be due to the increased heating and drying tendencies by convective processes in stochastic UNICON. An additional simulation of an idealized deep convection case shows that stochastic UNICON produces enhanced cloud variabilities with dependency on updraft radius, indicating its ability to represent the coexistence of shallow and deep convection.

54 ENVIRONMENTAL SCIENCES↗

Reducing AI RAG Hallucination by Optimizing Routing Techniques

Large Language Models (LLMs), such as ChatGPT, tend to “hallucinate”, meaning they confidently generate false information. Retrieval Augmented Generation (RAG) attempts to diminish hallucination by providing context to the LLM from data stores (indexes) containing relevant information. The LLM uses this context to formulate its response. RAG systems can still suffer from hallucination because of bad embeddings or ineffective routing. For example, a router will often return context from an irrelevant index, resulting in a hallucinated answer. In this study, we aim to minimize the frequency of routing hallucinations by optimizing Index Summary Routing.

97 MATHEMATICS AND COMPUTING↗

Scale-Bridging Optimization Framework for Desalination Integrated Produced Water Networks

In this work, we develop a Pyomo-based non-linear optimization strategy that includes rigorous MVR models. The detailed desalination unit is integrated into the multiperiod produced water network problem using the trust region filter (TRF) method. TRF decomposes the integrated problem into a master problem consisting of the network variables and a simplified surrogate model for the detailed desalination unit. The surrogate is updated using zero and first-order corrections from the optimal solution of the detailed models at every iteration. This framework allows us to co-optimize the design of the desalination units and operating policy for the multiperiod network. A common design is ensured across all periods using global capacity constraints. We validate the solution obtained using the TRF method by solving the full integrated problem for small network instances and show our results on real case studies on produced water networks from the Permian and Appalachian basins. In this work, we describe our TRF formulation, give details on our implementation in Pyomo, and analyze the results obtained by solving the optimization problem using IPOPT. We also present a discussion on the computational efficiency and scaling using the TRF approach against a full-scale integration of the rigorous models within the water network.

Naik, Sakshi↗

Simultaneous Optimization of Crystal Plasticity Hardening Parameters

Crystal plasticity models relate macroscopic deformation behavior to the evolution of slip systems strength, but their parameterization is often non-unique, with multiple parameter sets being able to describe the same macroscopic behavior. To address this issue, the present work adopts a Bayesian optimization framework for the parameterization of face-centered cubic plasticity models while simultaneously considering multiple experimental datasets from the literature. For single crystal Cu, parameter optimization was guided by the tensile stress–strain curves along several crystallographic orientations, with an adequate fit being found for five orientations at once. While additional parameters allowed for the consideration of more physical mechanisms, like different slip system interaction strengths or misorientations inherent to the experimental data, the extra dimensionality was found to limit the efficiency of the global minimization procedure. For polycrystalline Ni, multiple grain sizes were considered together in a representative polycrystalline model, with the optimization able to reconcile the model with the data for three grain sizes at once. As meaningful interpretation of parameters relies on the uniqueness of their values, incorporating multiple datasets into this discerning parameterization procedure enables more robust prediction and application of crystal plasticity models.

36 MATERIALS SCIENCE↗

Speedup of UEDGE Parameter Scans Using Machine-Learning Optimized OpenMP Parallelization and a Continuation Solver

This article presents the OpenMP parallelization of the preconditioning Jacobian assembly and right‐hand side residual evaluation in UEDGE. A continuation algorithm, utilizing the internal NKSOL implicit Jacobian‐Free Newton‐Krylov solver to efficiently scan physical parameters, is also presented. The implemented parallelization reduces the computational time for a benchmark scan run on 32 threads by compared to the serial version when using trained random forest regression models to identify the optimal decomposition of the system of equations. Random forest regression models applied to the UEDGE time‐dependent and continuation solver algorithms did not yield meaningful improvement in computational performance. A benchmark DIII‐D gas injection rate scan in the 0.35–0.75 kA interval, performed on a test cluster using the parallelized code and continuation solver, produced 1066 steady‐state solutions with a 22 s average wall‐clock computational time per steady‐state solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗