Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple solutions”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Final Technical Report: Transport of Complex Mixtures in Ion-Containing Polymer Membranes

Permselective ion-containing membranes are an integral component for many applications from water treatment, fuel cells, and solar fuels devices where the selective transport of molecules and ions is desired. In solar fuels devices, ion-containing polymer membranes are responsible for permitting selective transport of ions between electrodes to maintain overall charge neutrality yet limit transport of reaction products produced at the electrodes. While the transport of single solutes through such membranes has been fairly well described, binary and multicomponent transport is poorly understood due to the myriad of interactions that occur in these systems (i.e. between co-permeants and between permeants and the membrane). Solar fuels devices are just one example of an application where understanding the transport of multiple simultaneous species is critically important to improving device performance as product crossover leads to reductions in overall device performance. The objectives of this research was to improve our understanding of the complex array of factors that influence transport behavior of multiple solutes within ion-containing polymer membranes. This experimental project addressed the lack of fundamental understanding of multicomponent transport behavior by synthesizing ion exchange membranes with varied incorporation of comonomers (ionic and neutral moieties) to investigate fundamental relationships between membrane structure, membrane physiochemical properties, and transport behavior of solutes and complex solute mixtures through dense, hydrated membranes.

25 ENERGY STORAGE↗
DOE OSTI2025

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Slicing Solutions for Wire Arc Additive Manufacturing

Both commercial and research applications of wire arc additive manufacturing (WAAM) have seen considerable growth in the additive manufacturing of metallic components. However, there remains a clear lack of a unified paradigm for toolpath generation when slicing parts for WAAM deposition. Existing toolpath generation options typically lack the appropriate features to account for all complexities of the WAAM process. This manuscript explores the key slicing challenges specific to toolpaths for WAAM geometry and pairs each consideration with multiple solutions to mitigate most negative effects on completed components. These challenges must be addressed to minimize voids, prevent bead collapse, and ensure deposited components accurately approximate the desired geometry. Slicing considerations are grouped into four general categories: geometric, process, thermal, and productivity. Geometric considerations are addressed with overhang compensation, corner-sharpening, and toolpath-smoothing features. Process considerations are addressed with start point configuration and controls for the bead lengths and end points. Thermal and productivity considerations are addressed with island optimization, multi-material printing, and connected insets. Finally, tools for the post-processing of generated G-code are explored. Overall, these solutions represent a critical set of slicing features used to improve generated toolpaths and the quality of the components deposited with those toolpaths.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Genomic characterization of rare earth binding by Shewanella oneidensis

Abstract Rare earth elements (REE) are essential ingredients of sustainable energy technologies, but separation of individual REE is one of the hardest problems in chemistry today. Biosorption, where molecules adsorb to the surface of biological materials, offers a sustainable alternative to environmentally harmful solvent extractions currently used for separation of rare earth elements (REE). The REE-biosorption capability of some microorganisms allows for REE separations that, under specialized conditions, are already competitive with solvent extractions, suggesting that genetic engineering could allow it to leapfrog existing technologies. To identify targets for genomic improvement we screened 3,373 mutants from the whole genome knockout collection of the known REE-biosorbing microorganism Shewanella oneidensis MR-1. We found 130 genes that increased biosorption of the middle REE europium, and 112 that reduced it. We verified biosorption changes from the screen for a mixed solution of three REE (La, Eu, Yb) using Inductively Coupled Plasma Mass Spectrometry (ICP-MS) in solution conditions with a range of ionic strengths and REE concentrations. We identified 18 gene ontologies and 13 gene operons that make up key systems that affect biosorption. We found, among other things, that disruptions of a key regulatory component of the arc system ( hptA ), which regulates cellular response to anoxic environments and polysaccharide biosynthesis related genes ( wbpQ , wbnJ , SO_3183 ) consistently increase biosorption across all our solution conditions. Our largest total biosorption change comes from our SO_4685 , a capsular polysaccharide (CPS) synthesis gene, disruption of which results in an up to 79% increase in biosorption; and nusA, a transcriptional termination/anti-termination protein, disruption of which results in an up to 35% decrease in biosorption. Knockouts of glnA , pyrD , and SO_3183 produce small but significant increases (≈ 1%) in relative biosorption affinity for ytterbium over lanthanum in multiple solution conditions tested, while many other genes we explored have more complex binding affinity changes. Modeling suggests that while these changes to lanthanide biosorption selectivity are small, they could already reduce the length of repeated enrichment process by up to 27%. This broad exploratory study begins to elucidate how genetics affect REE-biosorption by S. oneidensis , suggests new areas of investigation for better mechanistic understanding of the membrane chemistry involved in REE binding, and offer potential targets for improving biosorption and separation of REE by genetic engineering.

Medin, Sean↗
DOE OSTI2023

Quantum protocol for decision making and verifying truthfulness among N ‐quantum parties: Solution and extension of the quantum coin flipping game

Abstract The authors devised a protocol that allows two parties, who may malfunction or intentionally convey incorrect information in communication through a quantum channel, to verify each other's measurements and agree on each other's results. This has particular relevance in a modified version of the quantum coin flipping game. The key innovation of the authors’ work includes the new design of a quantum coin that excludes any advantage of cheating, by which the long‐standing problem of the fair design of the game is, affirmatively, solved. Furthermore, the analysis is extended to N ‐parties communicating with each other, where multiple solutions for the verification of each player's measurement is proposed. The results in the N ‐party scenario could have particular relevance for the implementation of future quantum networks, where verification of quantum information is a necessity.

Ikeda, Kazuki↗
DOE OSTI2024

Optimizing and Exploring Untapped Micro-Hydro Hybrid Systems: a Multi-Objective Approach for Crystal Lake as a Large-Scale Energy Storage Solution

Increasing electricity demand and concerns about climate change and fossil fuel consumption have highlighted the importance of renewable energy resources and storage systems. This paper proposes a method for exploring untapped pumped hydro storage potentials to accommodate intermittent renewable energy generation profiles. Hourly measured data from 2022 in Benzie County, Michigan, United States, were gathered for system sizing and a thorough, realistic analysis. By employing the multi-objective grey wolf optimization algorithm, we formulated optimal sizing and energy-management strategies for three different scenarios. Unlike similar studies, the 3rd with triple objective functions (OFs) scenario aims to maximize both reliability and ecological OFs while minimizing the cost OF. It has shown promising results with multiple solutions, considering economic, environmental, and reliability factors. A case study conducted in Crystal Lake, Michigan, revealed that although Crystal Lake would function only as a micro-hydro power facility, it is a promising and huge storage unit with a substantial storage capacity of around 14.9734GWh. The system investigated is significant in the USA due to its rapid deployment capabilities, minimal construction requirements, and ease of integration with the distribution grid. The fuzzy logic method was employed to identify the best non-dominant solution among the other solutions. Furthermore, these outcomes include a notably low levelized cost of energy at 0.046147$/kWh, a robust index of reliability of 99.705%, and a significant reduction in CO₂ emissions amounting to 7.9142×10 3 tons/year, when considering the triple OFs. The paper’s methodology provides valuable insights for regions aiming to utilize renewable energy from untapped storage sources.

13 HYDRO ENERGY↗
DOE OSTI2021

Mass action model of solution activity via speciation by solvation and ion pairing equilibria

Abstract Solutes and their concentrations influence many natural and anthropogenic solution processes. Electrolyte and solution models are used to quantify and predict such behavior. Here we present a mechanistic solution model based on mass action equilibria. Solvation and ion pairing are used to model speciated solute and solvent concentrations such that they correlate to a solution’s vapor pressure (solvent activity) according to Raoult’s law from dilute conditions to saturation. This model introduces a hydration equilibrium constant ( K ha ) that is used with either an ion dissociation constant ( K id ) or a hydration modifier ( m ) with an experimentally determined ion dissociation constant, as adjustable parameters to fit vapor–liquid equilibrium data. The modeled solvation equilibria are accompanied by molecular dynamics (MD) studies that support a decline in the observed degree of solvation with increased concentration. MD calculations indicate this finding is a combination of a solvent that solvates multiple solutes, and changes in a solute’s solvation sphere, with the dominant factor changing with concentration. This speciation-based solution model is lateral to established electrostatics-based electrolyte theories. With its basis in mass action, the model can directly relate experimental data to the modeled solute and solvent speciated concentrations and structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Insights on the continuous representations of piecewise-smooth nonlinear systems: limits of applicability and effectiveness

Dynamical systems subject to intermittent contact are often modeled with piecewise-smooth contact forces. However, the discontinuous nature of the contact can cause inaccuracies in numerical results or failure in numerical solvers. Representing the piecewise contact force with a continuous and smooth function can mitigate these problems, but not all continuous representations may be appropriate for this use. Here, five representations used by previous researchers (polynomial, rational polynomial, hyperbolic tangent, arctangent, and logarithm-arctangent functions) are studied to determine which ones most accurately capture nonlinear behaviors including super- and subharmonic resonances, multiple solutions, and chaos. The test case is a single-DOF forced Duffing oscillator with freeplay nonlinearity, solved using direct time integration. This work intends to expand on past studies by determining the limits of applicability for each representation and what numerical problems may occur.

42 ENGINEERING↗
DOE OSTI2022

One Stomatal Model to Rule Them All? Toward Improved Representation of Carbon and Water Exchange in Global Models

Abstract Stomatal conductance schemes that optimize with respect to photosynthetic and hydraulic functions have been proposed to address biases in land‐surface model (LSM) simulations during drought. However, systematic evaluations of both optimality‐based and alternative empirical formulations for coupling carbon and water fluxes are lacking. Here, we embed 12 empirical and optimization approaches within a LSM framework. We use theoretical model experiments to explore parameter identifiability and understand how model behaviors differ in response to abiotic changes. We also evaluate the models against leaf‐level observations of gas‐exchange and hydraulic variables, from xeric to wet forest/woody species spanning a mean annual precipitation range of 361–3,286 mm yr −1 . We find that models differ in how easily parameterized they are, due to: (a) poorly constrained optimality criteria (i.e., resulting in multiple solutions), (b) low influence parameters, (c) sensitivities to environmental drivers. In both the idealized experiments and compared to observations, sensitivities to variability in environmental drivers do not agree among models. Marked differences arise in sensitivities to soil moisture (soil water potential) and vapor pressure deficit. For example, stomatal closure rates at high vapor pressure deficit range between −45% and +70% of those observed. Although over half the new generation of stomatal schemes perform to a similar standard compared to observations of leaf‐gas exchange, two models do so through large biases in simulated leaf water potential (up to 11 MPa). Our results provide guidance for LSM development, by highlighting key areas in need for additional experimentation and theory, and by constraining currently viable stomatal hypotheses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Analysis of nonlinear features associated with radio-frequency sheaths

An analytic study is conducted of nonlinear features that arise when radio-frequency (RF) sheath boundary conditions are solved self-consistently with the electromagnetic fields in a magnetized plasma. An instability associated with RF sheath–plasma interactions is identified and shown to be closely linked with sheath–plasma resonance and the occurrence of multiple solutions. It is proposed that this instability serves as a physical mechanism that enables jumping from one solution branch to another, providing a physical basis for the prediction of an RF sheath hysteresis in the time domain. It is further shown that this instability can lead to non-intuitive features such as spontaneous symmetry breaking in an otherwise completely symmetric problem. In plane-stratified geometry, two analytically solvable examples are considered for the cases of an unbounded half-space and a bounded domain

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Robust wind farm layout optimization

Wake interactions in wind farms cause losses in annual energy production (AEP) on the order of 10%. Wind farm designers optimize the layout of the farm to mitigate wake losses, especially in the dominant site-specific wind directions. As wind turbines and wind farms grow in scale, optimization becomes more complex. Offshore wind farms regularly comprise more than 100 wind turbines and are characterized by complex boundaries due to shipping lanes, neighboring wind farms, and other constraints. Layout optimization methods are broadly split between gradient-based and gradient-free approaches. Gradient-based approaches can converge quickly and perform well for smaller, academic problems but are often sensitive to initial conditions and tuning parameters and require expert knowledge to use. On the other hand, gradient-free approaches can be more robust to problem complexities. We present a robust layout optimization approach based on a random search algorithm. The algorithm is intended for those who are not optimization experts and has few tuning parameters that need specification to achieve satisfactory results. Unlike off-the-shelf methods, which use generally available, non-domain-specific optimization routines that accept as inputs an optimization function and constraint definitions, this approach takes advantage of the relative computational costs of the different evaluations by evaluating cheaper computations first (boundary and minimum distance constraints) and running expensive AEP evaluations only if all other checks pass. Moreover, an outer genetic algorithm allows multiple solutions to evolve in parallel, enabling rapid solution development on high-performance computers. We discuss the relative ease of selecting necessary tuning parameters and demonstrate the efficacy of the genetic random search on a complex layout problem consisting of placing 70 turbines in a nonconvex and unconnected boundary region.

17 WIND ENERGY↗
DOE OSTI2023

Shedding light on shadow generalized parton distributions

The feasibility of extracting generalized parton distributions (GPDs) from deeply virtual Compton scattering (DVCS) data has recently been questioned because of the existence of an infinite set of so-called “shadow GPDs” (SGPDs). These SGPDs depend on the process and manifest as multiple solutions (at a fixed scale 𝑄 2 ) to the inverse problem that needs to be solved to infer GPDs from DVCS data. SGPDs therefore pose a significant challenge for extracting GPDs from DVCS data. With this motivation we study the extent to which QCD evolution can provide constraints on SGPDs. This is possible because the known classes of SGPDs begin to contribute to observables after evolution, and can then be constrained (at the input scale 𝑄$^{2}_{0}$) by data that has a finite 𝑄 2 range. The impact that SGPDs could have on determining the total angular momentum, pressure and sheer force distributions, and tomography is also discussed. We find that 𝑄 2 evolution, given empirical data over a wide range of skewness (𝜉 ) and 𝑄 2 scales, constrains the SGPDs studied at low 𝑥 and low 𝜉. If this finding generalizes to all SGPDs—and we have not been able to find any SGPDs with qualitatively different behavior—then this potentially makes extraction of GPDs from DVCS possible over a limited range of 𝑥 and 𝜉.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A Cross-Domain Optimization Framework of PMU and Communication Placement for Multidomain Resiliency and Cost Reduction

Phasor measurement units (PMUs) play a crucial role in real-time monitoring and control of power grids. They rely on a communication network to transfer measurement data to the phasor data concentrator (PDC) for further processing and analysis. In this paper, a resilient cross-domain PMU and communication link placement method for minimizing the overall installation cost of the wide-area measurement system (WAMS) is proposed. Here, the main idea is to break down the barrier between the power grid domain and the communication domain, and consider the impact of one when design the other. The PMU placement in the power grid domain takes into account the cost of communication links by generating multiple solutions with equally minimum PMU costs for communication link placement evaluation. On the other hand, the communication link placement problem reduces the cost by customizing the routing policies based on the different roles of PMUs in grid observability. The proposed WAMS design is capable of withstanding any single component failure in the power domain (PMU failure or power branch failure) or in the communication domain (communication link failure or PDC failure). Numerical study on the IEEE 57-bus system reveals that the developed cross-domain optimization framework can significantly reduce the overall installation cost of WAMS while attaining multi-domain resiliency.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

On the Viability of Quantum Annealers to Solve Fluid Flows

This paper explores the suitability of upcoming novel computing technologies, particularly adiabatic annealing based quantum computers, to solve fluid dynamics problems that form a critical component of several science and engineering applications. For our experiments, we start with a well-studied one-dimensional simple flow problem, and provide a framework to convert such problems in continuum to a form amenable for deployment on such quantum annealers. Since the DWave annealer returns multiple states sampling the energy landscape of the problem, we explore multiple solution selection strategies to approximate the solution of the problem. We analyze the continuum solutions obtained both qualitatively and quantitatively as well as their sensitivities to the particular solution selection scheme.

42 ENGINEERING↗
DOE OSTI2021

Amphiphile conformation impacts aggregate morphology and solution structure across multiple lengthscales

Although the self-assembly of amphiphiles is well-studied in aqueous solutions, much less is understood about the fundamental driving forces and structure property relationships in non-polar media. In recent work [Journal of Physical Chemistry B, 2020, 124, 10822.] the authors have studied a series of malonomide-based amphiphiles that are relevant to liquid-liquid extraction. That work demonstrated that aggregation is largely driven by local dipole-dipole interactions between molecules. Here, we build upon this observation to develop a more detailed understanding of how the balance of dipole-dipole interactions (controlled by conformation) and molecular architecture influences the morphology of the aggregates across lengthscales. Using constrained molecular dynamics about key degrees of freedom, we demonstrate that the conformation of N,N-dimethyl,N,N-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N-dimethyl,N,N-dibutyltetradecyl malonamide (DMDBTDMA) has a significant impact upon self-association - where appropriate conformational sampling is essential. To quantify the aggregate morphology, several graph theoretic and persistent homology based properties are determined. The former examines the patterns of intermolecular interactions within clusters, while the latter examines the 3-dimensional spatial distribution. Based upon these analyses, we find that the morphology of aggregates, particularly at higher malonamide concentration, depends on a balance of dipole alignment and alkyl tail sterics. Dipole alignment encourages linear patterns of the intermolecular interactions within aggregates, while the alkyl tail steric interactions between the malonamide result in noticeably less linear aggregates for DMDOHEMA than DMDBTDMA. This is reflected in the spatial distribution, where more holes or voids exist between extractants within the DMDOHEMA that distribute within the solution in more of a "swiss cheese" arrangement as opposed to the more filamentous distribution of DMBDTDMA. This study links conformation and molecular structure to the morphology of amphiphile assemblies, and serves as a basis for ongoing study of multicomponent amphiphile solutions with polar and other solutes, and how these impact aggregation phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Microbial potentiometric sensor technology for real-time detecting and monitoring of toxic metals in aquatic matrices

Considering that toxic metals can affect metabolic processes in microorganisms adversely, it can be hypothesized that these metals in water matrices would induce a decrease in metabolic activity of the biofilm microorganisms populating the surface of a sensing electrode, which could be registered as a change in the open-circuit potential (OCP) generated by the biofilm microorganisms. The goal of this study was to test this hypothesis and demonstrate the underlying principle that microbial potentiometric sensor (MPS) technology could be used for long-term and real-time monitoring and detection of rapid changes in metal concentrations in realistic aquatic environments. To address the goal, four objective were addressed: (1) a batch reactor with three graphite-based MPS electrodes was fabricated; (2) a set of single-ion solutions and one multiple ion solution were prepared reflecting realistic concentrations of metals found in electroplating wastewaters; (3) the responses of the MPS to the simultaneous presence of multiple toxic metal ions in a single solution were measured; and (4) the changes of the MPS signals to the presence of individual metal ion solutions were examined. While the hypothesis was validated, the study also revealed that the MPS was sufficiently sensitive to not only detect, but also quantify, toxic metal ion concentrations in aqueous solutions. The coefficients of determination, which were R 2 > 0.995, and responsiveness of <1 μmol/L for some toxic metal cations, strongly support the performance of MPS technology in the echelons of expensive analytical tools capable detecting and measuring trace elements.The magnitude of the MPS response was toxic metal specific. When the molar concertation normalizes the inhibition portion of the signal area, the assessed sensitivity order was: Se>Cd>Pb>Ag>Ni> Zn. The study provides valuable information for enforcement agents, environmental professionals, and wastewater treatment operators, so toxic metal pollution and its detrimental impacts can be prevented and mitigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Multi-Fidelity Scheme for Accelerating 4D Finite Element Analysis for Mircoreactors

Next-generation microreactors are currently being designed to be operated terrestrial and extraterrestrial for remote surface power production. These systems will provide an alternative source of carbon-free energy that is versatile and can be utilized for various applications. These applications of microreactors have prompted the usage of high-fidelity unstructured finite element (FE) based approaches to provide time-dependent solutions for multiple physics fields. Computing these high-fidelity (4-D) solutions for multiple design iterations, physics fields, and transient events requires an immense computational cost. These high-fidelity solutions are computational expensive and the cost can be reduced through the implementation of less accurate low-fidelity solutions. Unlike the current fleet of commercial nuclear reactors, these next-gen systems present challenges due to the material and physical limitations required. Due to these constraints, the brute force technique of parameterizing important system characteristics determines whether a design meets the project objective. A system such as a nuclear reactor could have thousands of design parameters that affect system performance. In order to analyze the entire parameter space of such a complex system would require millions of CPU hours and countless design iterations. This costly approach is not practical due to regulatory and budget limitations. In this proposal, an approach that utilizes hybrid high-fidelity and low-fidelity FE models to reduce the computational cost of evaluating these applications in 4-D will be presented. The accuracy of the high-fidelity model and the computational efficiency of the low-fidelity model are taken advantage of to produce a solution that closely resembles the full order high-fidelity solution. By utilizing an unconverged coarse FE mesh, operation limits such as temperature, structural loading, etc., can be evaluated in an accelerated fashion and then can be spatially interpolated onto a finer mesh. The resulting error arising from the coarse mesh can be mitigated with a discrepancy function that actively quantifies and corrects the error in the coarse solution. This discrepancy function can be periodically updated with high-fidelity calculations across the temporal domain thus requiring less iterations on the fine mesh. As a result, the computational cost can be reduced for the evaluation and design iteration of microreactors. The proposal for this research contains three sections: Section 2 provides a literature review, Section 3 outlines the methodology for the proposed multi-fidelity scheme, and Section 4 displays preliminary results of the proposed research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Optimization of Water-Alternating-CO2 Injection Field Operations Using a Machine-Learning-Assisted Workflow

Summary This paper will present a robust workflow to address multiobjective optimization (MOO) of carbon dioxide (CO2)-enhanced oil recovery (EOR)-sequestration projects with a large number of operational control parameters. Farnsworth unit (FWU) field, a mature oil reservoir undergoing CO2 alternating water injection (CO2-WAG) EOR, will be used as a field case to validate the proposed optimization protocol. The expected outcome of this work would be a repository of Pareto-optimal solutions of multiple objective functions, including oil recovery, carbon storage volume, and project economics. FWU’s numerical model is used to demonstrate the proposed optimization workflow. Because using MOO requires computationally intensive procedures, machine-learning-based proxies are introduced to substitute for the high-fidelity model, thus reducing the total computation overhead. The vector machine regression combined with the Gaussian kernel (Gaussian-SVR) is used to construct proxies. An iterative self-adjusting process prepares the training knowledge base to develop robust proxies and minimizes computational time. The proxies’ hyperparameters will be optimally designed using Bayesian optimization to achieve better generalization performance. Trained proxies will be coupled with multiobjective particle swarm Optimization (MOPSO) protocol to construct the Pareto-front solution repository. The outcomes of this workflow will be a repository containing Pareto-optimal solutions of multiple objectives considered in the CO2-WAG project. The proposed optimization workflow will be compared with another established methodology using a multilayer neural network (MLNN) to validate its feasibility in handling MOO with a large number of parameters to control. Optimization parameters used include operational variables that might be used to control the CO2-WAG process, such as the duration of the water/gas injection period, producer bottomhole pressure (BHP) control, and water injection rate of each well included in the numerical model. It is proved that the workflow coupling Gaussian-SVR proxies and the iterative self-adjusting protocol is more computationally efficient. The MOO process is made more rapid by squeezing the size of the required training knowledge base while maintaining the high accuracy of the optimized results. The outcomes of the optimization study show promising results in successfully establishing the solution repository considering multiple objective functions. Results are also verified by validating the Pareto fronts with simulation results using obtained optimized control parameters. The outcome from this work could provide field operators an opportunity to design a CO2-WAG project using as many inputs as possible from the reservoir models. The proposed work introduces a novel concept that couples Gaussian-SVR proxies with a self-adjusting protocol to increase the computational efficiency of the proposed workflow and to guarantee the high accuracy of the obtained optimized results. More importantly, the workflow can optimize a large number of control parameters used in a complex CO2-WAG process, which greatly extends its utility in solving large-scale MOO problems in various projects with similar desired outcomes.

Energy & Fuels↗