SEARCH · Engineering Papers

Results for “Combustion simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

784 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Vitrification of Inactive Incineration Ash in a Full Industrial Scale In Can Melter - 20502

In the framework of the H2020 THERAMIN project funded by the European Commission, CEA performed a pilot-scale trial consisting of the vitrification of inactive ash (Al, Ca, Si, Zn, Bi, K, Mg, etc.) using the In- Can Melter technology. These ashes were coming from incineration of simulated technological wastes such as cotton, latex, neoprene, polyethylene. The borosilicate glass feasibility was first evaluated at a laboratory scale (10-100 g of glass). The volatile behavior of fly ash was controlled by pelletizing to avoid pipes clogging when performing the pilot scale test. A 10 wt.% inorganic binder was added to obtain cohesive pellets. The glass frit composition was selected and its ratio was fixed at 50 wt.%. After this step, a bench test (1 kg of glass) was done to verify the absence of container corrosion. Finally, it was followed by a pilot scale test (50 kg of glass) with 75,000 ash pellets used. The same glass microstructure and chemical composition were observed at each scale. Typical crystals of apatites (Ca, P), zincochromites (Zn, Cr, Al) and bismuth alloys (Bi, Sb) dispersed into a borosilicate glass matrix were observed. The chemical durability of the wasteform was measured during 28 days in a 90 deg. C water. Although there were durable crystalline phases, the hydrolysis rate of the vitreous part is relatively high because of the B and Na contents. However, leach rate profiles suggest that an alteration layer will form, leading to an alteration rate decrease. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Materials for Advanced Ultra-Supercritical (A-USC) Steam Turbines --- A-USC Component Demonstration

The U.S. Advanced Ultra-Supercritical (A-USC) Consortium was formed in 2001 as a government/industry program, sponsored by the U.S. Department of Energy (DOE) and the Ohio Coal Development Office (OCDO) and cost shared by industrial and not-for-profit partners. The purpose of the consortium was to advance the state of the art for power generation by evaluating and developing materials that allow the use of advanced steam cycles in coal-based power plants. These advanced cycles, with steam temperatures up to 1400°F (760°C), can increase the efficiency of coal-fired boilers from an average of 35% (current U.S. fleet) to more than 45% higher heating value (HHV) (>49% lower heating value [LHV]). The increase in a plant’s efficiency is limited unless new materials able to withstand these higher operating temperatures and pressures are identified and approved for use. The A-USC Consortium identified these needed materials during earlier phases of the program. It developed the welding and joining techniques along with manufacturing processes for casting and wrought products made from these new high-nickel alloys. It subjected these materials to extensive laboratory and steam loop testing. It then obtained ASME code approval for their use in U.S. boiler systems. The program’s successes leave this last remaining activity (ComTest Phase 2) that the U.S. utility industry has recommended to be accomplished prior to commercialization. The focus of the activity is the evaluation and demonstration of commercial readiness for “full scale” components to be made from these nickel-based alloy materials and provided by a U.S. domestic supply chain that is new to working with these alloys. According to studies completed by the Electric Power Research Institute (EPRI), the cost of an A-USC plant is approximately 20% higher than a non-A-USC plant because of its use of nickel-based alloys needed for the high temperature operating conditions. However, CO 2 reductions of approximately 30% from the current fleet average provide a strong incentive for its consideration. The actual costs and perceived value for CO 2 abatement will determine whether new or retrofitted plants are undertaken, although decisions to build A-USC plants in India would indicate its economic feasibility while also being part of a global carbon emissions strategy. The work by the A-USC Consortium, prior to the start of the ComTest project, has included lab scale and pilot scale materials testing, both in air and oxy-combustion. This testing has included air-cooled and steam-cooled “loops” that were installed into existing operating utility boilers to gain exposure of these materials to realistic conditions of high temperature and corrosion caused by the constituents in the coal ash. The A-USC Consortium also gained ASME Code approval of the Inconel 740 material, has cast and extruded the largest high nickel precipitation hardened alloys, and developed unique welding techniques to avoid problems identified by the competing European program. However, as valuable as these material test loops and accomplishments have been for obtaining information, their scale is below that required to minimize the risk associated for a U.S. utility to build a multibillion-dollar A-USC power plant. To reduce the final identified risk barrier to full-scale commercialization of these advanced materials and systems, the A-USC Consortium (guided by a utility industry advisory committee) has identified the key areas of the technology they desire to see as being capable of full-scale manufacturing and/or fabrication from an identified, capable U.S. domestic supplier base. A significant amount of work was accomplished during Phase 1 to identity the components, as well as the component size, that would be manufactured from advanced alloys such as Inconel 740H or Haynes 282 alloys. Pathways to supply these components for ComTest have been identified, as well as any further development that would be required. The Phase 2 effort used Phase 1 findings for designing these key full-scale components for A-USC boilers and turbines to include large castings; extrusions, forgings, fabrication of water walls and steam loops with headers from advanced materials, raw material (such as pipe extrusion billets) are at the commercial readiness level to permit advancement to a demonstration project. The Phase 2 work scope was addressed by a diverse team, including government, industry, and not-for-profit partners. The work scope under Phase 2 addressed fabrication of components identified as being outside of the proven capabilities of the existing supply chain, including the following: Steam turbine rotor forging and Haynes 282 nozzle carrier casting Superheater and reheater header and tube assemblies Large-diameter pipe extrusions and forgings Test valve articles to support ASME Code approval. In addition, key fabrication steps were completed, including boiler weld overlays and simulated field repairs. Throughout, extensive inspection and quality assurance testing of the components were performed. The team worked to advance ASME Code approval for key components and processes. Although much of the focus of ComTest Phase 2 was the high-temperature nickel-based alloy materials, a broader range of materials were incorporated, which would be representative of the materials used in full-scale A-USC power plant applications and have cross-cutting applicability on other high-temperature power generation options, such as advanced nuclear, supercritical CO 2 cycles, and central solar receivers. This report that has been submitted is organized in the following manner: Section 1 contains an Executive Summary. Section 2 discusses the ComTest project background and organization. Section 3 discusses project management and reporting. Section 4 discusses the procurement of nickel-based alloy and other A-USC materials and components. Section 5 discusses the fabrication of procurement of nickel-based alloy and other A-USC materials and components. Section 6 discusses the fabrication of cast nickel-based A-USC steam turbine components. Section 7 discusses the fabrication of forged nickel-based A-USC steam turbine piping and steam pipe components. Section 8 discusses the qualification of pressure relieve valves (PRVs) for A-USC power plants. Section 9 discusses proposed plans for future evaluation of A-USC components. Section 10 contains the summary and conclusion.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗

DOE OSTI2024

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Distributed fiber sensing systems for 3D combustion temperature field monitoring in coal-fired boilers using optically generated acoustic waves (Final Report)

In this project, we have developed and tested three kinds of fiber optic sensing systems for real time monitoring of temperature variations within an industrial scale boiler furnace. The fiber optic sensing systems target spatial and temporal distributions of high temperature profiles in a boiler furnace in fossil power plants. The reconstructed temperature profile will provide critical input for the control mechanisms to optimize the combustion process. This temperature profile will address the essential problem for fossil power plants in achieving higher efficiency and fewer pollutant emissions. Acoustic pyrometer systems have been used to reconstruct temperature field of power plant boilers based on measuring TOF (times-of-flight) of sound waves along some straight paths in a 2D cross-section of the boiler. In this project, optically generated acoustic signals from a fiber optic sensing system have replaced the acoustic signals generated from an electrical transducer. A 3D reconstruction algorithm replaced the previous 2D model. In this project, three kinds of fiber optic sensing systems have been developed and tested. They are fiber optic sensing system I, fiber optic sensing system II (Distributed Sensing System I) and fiber optic sensing system III (Distributed Sensing System II). For fiber optic sensing system I, the fiber optic ultrasound generator acts as a signal generator. A microphone, hydrophone or other electronic devices serve as a signal receiver. In this system, there are one generator and one receiver. Distance test, water temperature test, air temperature test, air temperature reconstruction, and GE ISBF pilot test were performed by Fiber optic sensing system I. The fiber optic sensing system I successfully detected temperature in all these tests. We got 2D temperature reconstruction results by using the fiber optic sensing system I and it matched the reference data. The fiber optic sensing system I successfully survived in GE ISBF boiler environment (480 °F). For fiber optic sensing system II (Distributed Sensing System I), it is an all optical ultrasound system. The fiber optic ultrasound generator acts as a signal generator. Fiber Bragg Grating (FBG) and Fabry-Perot (FP) sensor act as a signal receiver. In this system, there is one generator and one receiver. Aluminum plate temperature test, furnace high temperature test, and GE ISBF pilot test were performed by the fiber optic sensing system II. Fiber optic sensing system II successfully detected the temperature in all these tests. The fiber optic sensing system II successfully survived at up to 700 °C furnace environment and 320 °C GE ISBF boiler environment. For fiber optic sensing system III (Distributed Sensing System II), it is also an all optical ultrasound system. The fiber optic ultrasound generator acts as a signal generator. Multiple FP fiber sensors act as signal receivers. In this system, there is one generator and three receivers. Three GE ISBF pilot tests were performed by fiber optic sensing system III. The fiber optic sensing system III survived in the cold flow tests in GE’s ISBF pilot test facility. However, we didn’t get high temperature data by using this system since the nanosecond laser issues. During the period of the project, test trials, data simulation and algorithm optimization was performed successfully. For real time temperature field construction, the sampling rate must be fast enough to capture the field variations. The technology of Code-division multiple access (CDMA) is well studied which could allow parallel multiplexing, even if signals overlap in time or frequencies. Moreover, it has been known that extending the length of signal significantly improves SNR. For acoustic signals, these multiplexing techniques have also been widely used, mainly for sonar and acoustic communications. The CDMA modulation technique has been proposed and studied to guarantee high network throughput, low channel access delay and low energy consumption. We have studied the temperature field reconstruction using Gaussian Radial Basis Functions (GRBF)-based approximation approach. Reconstruction of 3D temperature field using Neural Networks with measured TOF and known propagation paths is feasible. 2D and 3D temperature field reconstruction simulation results are achieved. The milestone status is shown in Table 1. We finished milestone 1-8 and milestone 10. For milestone 9, we did three pilot tests by using the fiber optic sensing system III (Distributed Sensing System II) at GE Power. However, due to the failure of the ns laser, we did not get the temperature results. We conducted some additional tasks that were not originally proposed: 1) We fabricated a fiber optic sensing system I and did a pilot test based on this system. 2) In the proposal, we proposed two pilot tests at GE Power. In reality, we finished at least seven pilot tests at GE Power. GE Power has made a lot of efforts for supporting the pilot tests. 3) We got a simulation results based on CDMA. In summary, most of the tasks have been accomplished. The outcome of this project removed a few barriers that hinder the achievement of the final product of the distributed sensing systems. With the successful accomplishment of this project, a prototype of the fiber optic sensing system can be fabricated to attract more interests from companies and other funding agencies.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Numerical Study of PAHs and Soot Emissions from Gasoline–Methanol, Gasoline–Ethanol, and Gasoline– n -Butanol Blend Surrogates

Soot formation is an intricate phenomenon, and soot propensity of a fuel is interwoven with the fuel composition, physical and chemical properties, and combustion environment. Here, the present study examines the hypothesis that in addition to the chemical composition of the fuel, the sooting nature of the fuel is closely coupled with its chemical property known as octane sensitivity (S). With this motivation, the present study numerically investigates the effects of gasoline surrogate composition and its property, octane sensitivity (S), on polycyclic aromatic hydrocarbons (PAHs) and soot emissions. Four-component toluene primary reference fuel (TPRF)–alcohol blends, comprising iso-octane, n-heptane, toluene, and one of the three different alcohols- methanol, ethanol, and n-butanol, are used as gasoline surrogates. A total of 320 TPRF–alcohol mixtures, with S in the range of 1–10, are examined under laminar counterflow diffusion flame conditions. A detailed chemical mechanism coupled with a comprehensive soot model, which includes reactions for soot inception, surface growth, PAH condensation, and oxidation, is adopted. The analysis indicates that the toluene content in the fuel mixture has a prominent effect, while the alcohol content and octane sensitivity of the fuel have a weak correlation with the PAHs and soot. Thus, it is not clear if any of these three variables, namely, toluene content in the fuel, alcohol content in the fuel, and S, are individually sufficient to characterize the PAHs and soot across various blends. For this reason, a new variable (X CHO ) based on the elemental composition of the fuel mixture is identified and it is shown that X CHO along with S of the fuel characterize soot emissions satisfactorily. Further, a reaction path analysis indicates that the efficacy of alcohols in reducing soot emissions follows the order: methanol > ethanol > n-butanol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Quantifying the thermal effect and methyl radical production in nanosecond repetitively pulsed glow discharges applied to a methane-air flame

In this work, we investigated non-equilibrium plasma produced by nanosecond repetitively pulsed glow discharges applied across a lean premixed methane-air flame. The flame is stationary, axisymmetric, and laminar. The discharges are applied on the symmetry axis crossing the reactant gases, flame front, and product gases, allowing phase-locked averaged measurements and comparisons with axisymmetric numerical simulations. The thermal effect and methyl radical production are quantified in the discharge in the reactant gas region. One-dimensional, two-beam, hybrid, femtosecond-picosecond, coherent anti-Stokes Raman scattering is used to acquire spatial and temporal profiles of temperature and oxygen-to-nitrogen concentration ratio. Photo-fragmentation laser-induced fluorescence is used to acquire quantitative two-dimensional profiles of methyl radicals in the discharge providing the first quantitative imaging of methyl produced ahead of a flame by plasma-induced methane dissociation. The spatial profiles of temperature and oxygen-to-nitrogen concentration ratio are in steady state, indicating that individual discharges have an insignificant heating effect. Upper and lower bounds of the produced mole fraction of methyl radicals in the plasma are obtained due to uncertainties in the collisional quenching rates of excited state methylidyne radicals in the plasma. The discharges produce a maximum of 600–1100 ppm of methyl radicals upstream of the flame front within 25 ns. This amount is similar to the predicted methyl mole fraction for the flame without plasma and thus represents a significant chemical perturbation to the reactants upstream of the flame front. The produced methyl follows an exponential decay in the first microsecond after the discharge with a decay constant of 8 µs close to the flame, and 0.8 µs further from the flame. The decay then deviates from the exponential curve and the methyl persists for tens of microseconds. The results suggest that for the tested configuration, the thermal effect of individual discharges through fast gas heating is negligible, while active chemical species are produced in large quantities in the reactant gases, upstream of the flame front.

hybrid fs-ps CARS↗

DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗

DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗

DOE OSTI2021

High-dimensional Data-driven Energy optimization for Multi-Modal Transit Agencies (HD-EMMA) (Final Technical Report)

Public bus transit services in the U.S. are responsible for at least 19.7 million metric tons of CO 2 emission annually. Electric vehicles (EVs) can have a much lower environmental impact than comparable internal combustion engine vehicles (ICEVs), especially in urban areas. Unfortunately, EVs are also much more expensive than ICEVs. As a result, many public transit agencies can afford only mixed fleets of transit vehicles, consisting of EVs, hybrids (HEVs), and ICEVs. Transit agencies that operate such mixed fleets of vehicles face a challenging optimization problem: these agencies need to decide which vehicles are assigned to serving which transit trips. Since the advantage of EVs over ICEVs varies depending on the route and time of day (e.g., the benefit of EVs is higher in slower traffic with frequent stops and lower on highways), the assignment can have a significant effect on energy use and, hence, environmental impact. Through this project, we have developed reference data about energy collections and constructed a set of machine learning models that can accurately predict the energy consumption for the whole fleet at the level of each trip. We have used these models to develop a scheduling and assignment strategy that can rotate the different vehicle types across the transit agencies’ routes. The optimization algorithm ensures that the vehicles are matched to trips considering weather patterns, expected congestion, and road gradients to minimize the overall energy usage. We list the key observations from our project for other practitioners below. Details are available in the report, and the list of source code and our publications are included in the appendix. 1. We have demonstrated the feasibility of collecting, merging and analyzing large volumes of high-resolution real-world telemetry data from a mixed vehicle fleet. To mitigate the inherent noise of the recorded GPS points, the team developed an algorithm that filters data and maps the points onto a street. The algorithm considers previous and subsequent location measurements and different characteristics of nearby streets to determine how likely the vehicle travels on them. Then, the team segmented the time series into disjoint contiguous samples based on adjacent road segments and repeated the outlier detection and removal. For each data point, the team added features corresponding to elevation changes within the samples, weather features, such as temperature, and traffic data, such as speed ratio between actual speed and free-flow speed. 2. We have developed two forms of machine learning models that be used to understand and analyze the energy operations of a mixed vehicle transit fleet. The micro prediction model provides estimates of instantaneous energy prediction for all types of buses (diesel, hybrid, and electric). Such a model is important in evaluating the energy impacts of real-time bus operation strategies, but it is challenging due to diversified driving cycles of transit buses. The model can help the drivers understand the impact of their driving behaviors and short-term congestions. The macro prediction models estimate average energy consumption across the whole trip considering the features: distance traveled, various road-type features, elevation change, day of the week, time of day, various weather features (temperature, humidity, etc.), and traffic features (speed ratio and jam factor). 3. We have demonstrated that it is possible to transfer the machine learning models we have developed in this project to other teams and cities by using inductive transfer learning. We also showed that the performance of the macro energy prediction models can be improved using a multi-task learning approach where the learning parameters are shared between the models being developed for different vehicle types. The advantage of this approach is improved learning performance as the models can exploit common spatio-temporal and environmental characteristics. 4. Finally, we have developed trip and vehicle assignment and scheduling algorithms that use the energy prediction models and develop a trip to vehicle type (diesel, electric, hybrid) assignment for the whole operation to reduce overall emissions and cost. We have shown through simulations that the proposed algorithms can save $\$$ 48,910 in energy costs and 175 metric tons of CO 2 emission annually for CARTA.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗