Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Current measurement”

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 289 records · Page 16

Measurement of Charged-Current $ν_μ$ and $ν_e$ Interactions Towards a Search for Low-Energy $ν_Ε$ Excess Using Wire-Cell In Microboone

This technote summarizes the existing work in searching for $ν_e$ low-energy excess (eLEE) in MicroBooNE Booster Neutrino Beam (BNB) data stream based on the Wire-Cell event reconstruction paradigm. The charged-current $ν_μ$ and $ν_e$ events are selected from the 5.3e19 POT open data from the BNB beam, 6.37e20 POT far sideband from the BNB beam data, and 2.10e20 POT data from Neutrinos at the Main Injector (NuMI) beam. The charged-current $ν_e$ selection results from the BNB data that are sensitive to the eLEE search are not included. Various comparisons between data and Monte Carlo predictions are performed to validate the overall model and demonstrate the power of the analysis techniques. Physics sensitivities in terms of the exclusion and the discovery potential are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of charged-current $\nu_{\mu}$ and $\nu_e$ interactions with wire-cell in MicroBooNE towards a search for low-energy $\nu_e$ excess

This technote summarizes the existing work in searching for $\nu_e$ low-energy excess (eLEE) in MicroBooNE Booster Neutrino Beam (BNB) data stream based on the Wire-Cell event reconstruction paradigm. The charged-current $\nu_{\mu}$ and $\nu_e$ events are selected from the 5.3e19 POT data from the BNB beam and 2.06e20 POT data from Neutrinos at the Main Injector (NuMI) beam. The charged-current $\nu_e$ selection results from the BNB data that are sensitive to the eLEE search are not included. Various comparisons between data and Monte Carlo predictions are performed to validate the overall model and demonstrate the power of the analysis techniques.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Assessment of Existing OpenStudio Measures: Reviews, Interviews, and Future Developments

OpenStudio Measure development is continuously in progress and greatly propelled by the collaborative efforts within the building energy modeling community. To ensure the widespread adoption and benefit of OpenStudio Measures, developers must understand the current status of Measure development and the needs of OpenStudio Measure users. The first step involves a comprehensive review of existing content to prevent redundancy and gain insights into how OpenStudio Measures are used in the building energy modeling community. This knowledge can then be integrated into the Measure development process, and the expertise of practitioners and OpenStudio Measure users can be leveraged to shape future Measures. Numerous OpenStudio Measures have been created and shared on the Building Component Library (BCL). The BCL is an open-source repository housing various OpenStudio-related resources, including building component blocks, descriptive metadata, and Measures describing modifications to building energy models. The OpenStudio Measures in BCL encompass a wide range of energy conservation Measures from basic lighting power reduction to complex HVAC model transformation. They also enable users to generate customized reports and facilitate the integration of energy simulation with other analytical processes. This report presents review of 272 currently available OpenStudio Measures in BCL. The OpenStudio Measures were reviewed by category and subcategory. These Measures are summarized by their functionalities and keywords. To gain insights into how OpenStudio Measures are used in building energy modeling community, interviews were conducted. A total of 12 interview responses were collected from 6 individuals in the industry and 6 individuals in academia. The knowledge acquired from reviewing the existing Measures and interview results will be integrated into the Measure development process, and the expertise of practitioners and OpenStudio Measure users will be leveraged to shape future Measures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Measurements at peak operational beam current in the SNS beam test facility

Work at the SNS beam test facility has focused on high dimensional and high dynamic range measurements of the medium energy (2.5 Mev) beam distribution. This is motivated by the need to understand and predict beam losses down to one-part-per-million. The initial demonstration of full-and-direct 6D phase space measurement was done at a current of 40 mA transported through the RFQ. Since that demonstration, more detailed studies have been performed at lower transported currents (in the range 30 mA and below). This is due to a hardware change - recent runs utilize the original SNS RFQ, which after a decade of service in the SNS achieves transmission significantly below design (50-60%, vs >80%). A short run in 2023 with a newly-commissioned RFQ enables maximum transmission. Preliminary results from beam distribution measurements during this run are discussed.

Ruisard, Kiersten↗

Direct current response of a thin scCVD diamond detector under increased applied field to 14.1 MeV neutrons

An avalanche effect yielding inherent gain can be exploited in thin, single-crystal chemical vapor deposition (scCVD) diamond. It occurs when a high enough bias is applied across the diamond thickness while avoiding breakdown. This charge multiplication effect was studied previously with alpha particles and heavy ions either by using the transient current technique or by measuring the energy spectrum. The measurements we obtained to evaluate the charge multiplication performance of a 10 μm thick scCVD diamond detector used a novel approach—we employed an electrometer to characterize the response of the detector by performing directly coupled current measurements (time-averaged charge, at 1 Hz sampling) when exposed to 14.1 MeV neutrons from deuterium-tritium fusion. We measured both the dark and irradiated currents from the detector over a range of applied displacement field values from 2 to 75 V/μm. A histogram method with central mean and standard deviation width was used to determine the current over each measurement duration typically from 100 to 300 seconds. The dark-subtracted irradiated current (i.e., contrast) was used to evaluate the gain of the detector at each applied displacement field. The contrast at an applied displacement field between 15 and 20 V/μm was higher than the expected linear increase in contrast proportional to the increased applied bias, indicating the possible presence of avalanche events in the diamond. The detector response also indicated possible polarization and charge depletion effects. These results provide an opportunity to further explore the use of thin scCVD diamond as a fast neutron current mode detector with inherent gain.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurement of Inclusive Charged-Current $\barν_μ$ Scattering on C, CH, Fe, and Pb at $\langle E_{\barν}\rangle \sim$ 6 GeV with MINERvA

We report MINERvA's first measurement of inclusive charged-current $\barν_μ$ cross sections on carbon, hydrocarbon, iron, and lead, and their ratios to the cross section on hydrocarbon, as functions of the antimuon transverse momentum, $p_{\mathrm{T}}$. Using a wide-band $\barν_μ$ beam with mean energy $\sim 6~\text{GeV}$, these measurements probe all interaction modes, including the transition from resonance production to deep-inelastic scattering. The total uncertainties are typically $5-10\%$ for the absolute cross sections and $2-5\%$ for the ratios. Comparisons with multiple neutrino interaction models reveal significant discrepancies in the $p_{\mathrm{T}}$ dependence, particularly for heavier nuclei. The disagreements are most pronounced at low $p_{\mathrm{T}}$ but extend across the full $p_{\mathrm{T}}$ range, indicating missing or mis-modelled nuclear effects.

Klustová, A. [Imperial Coll., London] (ORCID:00000↗

Operando Microscopy Characterization of Perovskite Solar Cells

In this work, we developed operando Kelvin probe force microscopy (KPFM) to study the electrostatic potential distribution across perovskite cells under light and forward bias to gain a deeper understanding of device operation physics. As a case study, we selected perovskite cells with a SnO 2 -based electron-selective layer (ESL), which showed great potential for fabricating high-efficiency, hysteresis-free devices due to the deeper conduction band and higher electron mobility of SnO 2 . The as-made device showed a main junction at the perovskite/spiro interface. After light soaking and applying forward bias, the junction quality improved, possibly explained by filling trap states at the interfaces and by the perovskite absorber perhaps having a self-poling effect; the main junction is observed at the ESL/perovskite interface. The results are consistent with current-voltage measurements, device performance improves mainly with fill factor enhancement. The operando KPFM results should more closely reflect the real case during current density-voltage measurements or solar cell operation. The operando KPFM technique that we have developed can be a powerful tool to provide a deeper understanding of the device operation mechanism and to further optimize the device.

14 SOLAR ENERGY↗

Supervisory Control and Data Acquisition for Electrochemical Separation Experimentation

The Python-based program is a laboratory automation tool designed to control and monitor electrochemical systems. The tool was developed for capacitive deionization (CDI) experiments, but it can be used for any system that requires controlled voltage or current segments and multi-parameter monitoring. The program integrates hardware components to run user-defined experimental parameters, providing operational control of a programmable power supply, peristaltic pump, and data acquisition devices. Currently, the program is structured with a workflow that includes an initialization (or pre-run) phase, a main loop, and a post-experiment stabilization (or post-run) phase. The initialization phase prepares and stabilizes the cell, ensuring that the electrodes and solution reach a baseline state before the experiment begins. The main loop consists of multiple voltage segments that repeat, controlling the experiment while recording key parameters such as time, voltage, current, pH, and conductivity. Finally, the post-experiment stabilization phase allows the system to stabilize after the experiment, returning the cell and solution to equilibrium conditions before ending the sequence. The program is designed with four variations, each tailored to different experimental needs. All variations include both the initialization and post-experiment stabilization stages, which run for a set amount of time, voltage, current, and flow rate before and after the main experiment block. The main loop runs for a set number of cycles, as defined by the user input, and each cycle is composed of 2 or 4 segments. The 4 program variations are described as follows: Program 1: The main program includes 2 segments. Each segment is defined to have a set duration, flow rate, voltage, and current. This program measures conductivity, flow rate, voltage, and current. Program 2: The main program expands Program 1 to include 4 segments. Each segment has a specified duration, flow rate, voltage, and current. Like Program 1, it measures conductivity, flow rate, voltage, and current. Program 3: The main program consists of 2 segments, each defined by time, flow rate, voltage, and current. In addition to conductivity, flow rate, voltage, and current, Program 3 collects pH and temperature data through a 4-channel data acquisition device. Program 4: This program independently controls two channels of a multi-channel power supply simultaneously. While conductivity can only be measured for one cell at a time, the dual-channel control makes it possible to operate two cells simultaneously under different voltage/current conditions. The main program includes 2 segments.For each program, all measurements are automatically logged and integrated into a single Excel output file. Data are displayed in numerical format and plotted, both in real time, to track system performance. A key feature of the program is its ability to synchronize all outputs so that every measurement shares a single timestamp, ensuring accurate alignment of voltage, current, pH, conductivity, and pH data.By combining hardware control, real-time monitoring, and unified data collection, this program significantly reduces manual workload and minimizes errors, making it a reliable platform for researchers, engineers, and laboratory technicians conducting CDI experiments, among other electrochemical tests.

Valentino, Lauren [Argonne National Laboratory (AN↗

Laboratory study of the failed torus mechanism in arched, line-tied, magnetic flux ropes

Coronal mass ejections (CMEs) are some of the most energetic and violent events in our solar system. The prediction and understanding of CMEs are of particular importance due to the impact that they can have on Earth-based satellite systems and, in extreme cases, ground-based electronics. CMEs often occur when long-lived magnetic flux ropes (MFRs) anchored to the solar surface destabilize and erupt away from the Sun. One potential cause for these eruptions is an ideal magnetohydrodynamic (MHD) instability, such as the kink or torus instability. Previous experiments on the magnetic reconnection experiment revealed a class of MFRs that were torus-unstable but kink-stable, which failed to erupt. These “failed-tori” went through a process similar to Taylor relaxation, where the toroidal current was redistributed before the eruption ultimately failed. Herein we have investigated this behavior through additional diagnostics that measure the current distribution at the foot points and the energy distribution before and after an event. These measurements indicate that ideal MHD effects are sufficient to explain the energy distribution changes during failed torus events. This excludes Taylor relaxation as a possible mechanism of current redistribution during an event. A new model that only requires non-ideal effects in a thin layer above the electrodes is presented to explain the observed phenomena. This work broadens our understanding of the stability of MFRs and the mechanism behind the failed torus through the improved prediction of the torus instability and through new diagnostics to measure the energy inventory and current profile at the foot points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phased plan for the implementation of the time-resolving magnetic recoil spectrometer on the National Ignition Facility (NIF)

The time-resolving magnetic recoil spectrometer (MRSt) is a transformative diagnostic that will be used to measure the time-resolved neutron spectrum from an inertial confinement fusion implosion at the National Ignition Facility (NIF). It uses a CD foil on the outside of the hohlraum to convert fusion neutrons to recoil deuterons. An ion-optical system positioned outside the NIF target chamber energy-disperses and focuses forward-scattered deuterons. A pulse-dilation drift tube (PDDT) subsequently dilates, un-skews, and detects the signal. While the foil and ion-optical system have been designed, the PDDT requires more development before it can be implemented. Therefore, a phased plan is presented that first uses the foil and ion-optical systems with detectors that can be implemented immediately—namely CR-39 and hDISC streak cameras. These detectors will allow the MRSt to be commissioned in an intermediate stage and begin collecting data on a reduced timescale, while the PDDT is developed in parallel. A CR-39 detector will be used in phase 1 for the measurement of the time-integrated neutron spectra with excellent energy-resolution, necessary for the energy calibration of the system. Streak cameras will be used in phase 2 for measurement of the time-resolved spectrum with limited spectral coverage, which is sufficient to diagnose the time-resolved ion temperature. Simulations are presented that predict the performance of the streak camera detector, indicating that it will achieve excellent burn history measurements at current yields, and good time-resolved ion-temperature measurements at yields above 3 × 10 17 . The PDDT will be used for optimal efficiency and resolution in phase 3.

47 OTHER INSTRUMENTATION↗

Dark matter trigger for early dark energy coincidence

Current cosmological measurements present a persistent tension in the value of the current cosmic expansion rate, the Hubble constant, as inferred from cosmic microwave background (CMB) and large-scale structure (LSS) data compared to that inferred from the classical distance ladder. Early dark energy (EDE), whose cosmological role is localized in time around the epoch of matter-radiation equality just prior to the release of the CMB photons, is designed to resolve this "Hubble tension". However, the model introduces a new coincidence problem: Why should the EDE dynamics occur near matter-radiation equality if EDE is decoupled from both matter and radiation? The resolution of this problem may lie in an early dark sector (EDS), wherein the dark matter mass is dependent on the EDE scalar field. In this work, we construct such an EDS model and show that it naturally resolves the EDE coincidence problem at the background level without any fine-tuning of the coupling to dark matter or of the initial conditions. When fitting to current cosmological data, including that from the local distance ladder, CMB, and LSS, our EDS maximum-likelihood model performs comparably to EDE for resolving the Hubble tension. However, fitting the Planck CMB data requires a specific range of initial field positions to balance the scalar field fluctuations that drive acoustic oscillations, providing testable differences with other EDE models and a platform for future model-building.

79 ASTRONOMY AND ASTROPHYSICS↗

Heterogeneous Integration Technologies for High-temperature, High-density, Low-profile Power Modules of Wide Bandgap Devices in Electric Drive Applications (Final Technical Report)

The goal of this project is to develop packaging technologies for making high-temperature, high-density, and low-profile wide-bandgap (WBG) power electronics modules for electric drives. These modules are aimed at enabling the DOE’s University Consortium to reach its 2025 inverter targets of ≥ 100 kW/L and ≤ 2.7 $/kW. The specific objectives are to: design and fabricate SiC half-bridge power modules with double-sided cooling and parasitic inductances < 5 nH, heat flux density > 400 W/cm 2 , and working junction temperature of 200 o C; design, fabricate, and deliver a gate driver with double-sided cooled modules for the construction of a 100 kW/L inverter at Oak Ridge National Lab; and design and prototype intelligent gate drivers with integrated current sensor and a low-profile DC-DC power supply with air-core transformer for testing power modules at 200 o C junction temperature. We followed an iterative technical approach of design, simulation, fabrication, and testing of various versions of modules, current sensors, and power supply. The state-of-the-art silicon carbide devices rated at 1.2 kV and 149 A were packaged by sintered-silver bonding on an aluminum nitride direct-bond-copper substrate for high thermal conductivity, high working temperature, and high joint reliability. Porous silver posts were used to interconnect the device’s source pads to the other direct-bond-copper substrate for low mechanical stresses, ease of manufacturing, and double-sided cooling. A current sensor based on package parasitic inductance was developed to measure switching current. A dynamic feedback scheme was developed to compensate the effect of parasitic resistance and temperature variation. A constant-current class-E dc-dc converter with air-core transformer was developed. Air-core transformer was used due to the unavailability of magnetic core at high temperatures. Gate driver and power supply were integrated with the double-side cooled, high temperature SiC power modules for testing the modules at 200 o C junction temperature. Double-pulse and continuous testing of the integrated technologies validated the design and fabrication of the three component technologies. Throughout the project, we overcame the challenge for design verification caused by low prototyping yield, which then helped train the graduate students, the future workforce, to learn the engineering know-how for low-cost manufacturing of reliable products. Below is a summary of the major accomplishments of this project: development of a prototyping process for fabricating double-side cooled (1200 V, 149 A) SiC phase-leg modules capable of working to 200 o C Tj; simulation and experimental verification of the improvement of thermo-mechanical reliability of the double-side cooled SiC phase-leg module by using rigid encapsulant; design and experimental validation of a current sensor based on package parasitic inductance and a compensation solution for eliminating the effect of parasitic resistance; design and experimental validation of a low-profile power supply with six-output air-core transformer for gate driver; functional demonstration of a SiC phase-leg module integrated with its gate driver, current sensor, and an air-core power supply at 200 o C Tj in a double-pulse switching test setup and Buck converter continuous test setup; successful completion of six PhD and two MS students who are or will work at Apple Inc., Tesla Inc., Wolfspeed Inc., Microchip Inc., Monolithic Power Systems Inc., and LG Magna Inc.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Opening the Door to Grid-Interactive Efficient Buildings with Energy Codes

Building codes represent standard design practice in the construction industry and must continually evolve to account for advancing technologies and innovative practices. Energy codes have historically focused on energy efficiency within buildings and across their systems. However, many of tomorrow's technologies go beyond efficiency and target increased flexibility. These demand flexibility (DF) measures can postpone or reduce building electric load based on price or other grid signals, and include smart appliances, connected lighting, and connected mechanical systems. Expanding codes to enable such grid-interactive efficient building (GEB) has the potential to influence their realization at scale, which will support renewable-energy grid-integration and building-sector decarbonization. However, considering DF measures in code development creates a new set of challenges for codes and the practices contained therein. This paper considers the role of building energy codes in enabling GEB. Specifically, it reviews the status of demand flexibility (DF) measures in current commercial building energy code. It examines the national model code development process and identifies components barring the consideration and inclusion of DF measures - including code scope, characterization and analysis of proposed new prescriptive measures, and accounting for time-of-day and geographic differences in their benefits. The paper presents findings from code development analyses that indicate the cost benefit of DF measures and the limitations associated with current code development conventions. To encourage building flexibility and improved energy resilience moving forward, recommendations are made for removing code development barriers and sanctioning DF measure consideration in future model codes.

Building Energy Codes and Standards, Grid-Interact↗

Modeling Ring Current Proton Fluxes Using Artificial Neural Network and Van Allen Probe Measurements

Abstract Terrestrial ring current dynamics are a critical part of the near‐space environment, in that they directly drive geomagnetic field variations that control particle drifts, and define geomagnetic storms. The present study aims to specify a global and time‐varying distribution of ring current proton using geomagnetic indices and solar wind parameters with their history as input. We train an artificial neural network (ANN) model to reproduce proton fluxes measured by the Radiation Belt Storm Probes Ion Composition Experiment instrument onboard Van Allen Probes. By choosing optimal feature parameters and their history length, the model results show a high correlation and a small error between model specifications and satellite measurements. The modeled results well capture energy‐dependent proton dynamics in association with geomagnetic storms, including inward radial diffusion, acceleration and decay. Our ANN model produces proton fluxes with their corresponding 3D spatiotemporal variations, capturing the latitudinal distribution and local time asymmetry that are consistent with observations and that can further inform theory.

Li, Jinxing↗

Status Update on the High Precision Isotopic Measurements on High Burnup LWR Fuel in 2020

The US Department of Energy (DOE) Office of Nuclear Energy (NE) is currently investigating the feasibility of directly disposing dual-purpose (storage and transportation) canisters (DPCs) in a spent nuclear fuel (SNF) repository. Criticality during the repository performance period (10,000 years or more) is one of the major concerns related to direct disposal of DPCs, specifically as the system undergoes degradation in the repository environment and timeframe. Oak Ridge National Laboratory (ORNL) is developing an as-loaded criticality analysis methodology using full (actinides + fission products) burnup credit that exploits the inherent criticality margin associated with actual canister-specific loading configuration. Burnup credit criticality analysis requires validation of the depletion/decay codes used to generate the burned isotopic inventory of an assembly by comparing the code-predicted inventory with experimentally determined isotopic data. Currently, isotopic measurement data for boiling water reactors (BWR) SNF are limited, and additional measurements will be highly beneficial for BWR burnup credit analysis, which is essential to demonstrate disposability of BWR DPCs. Moreover, new pressurized water reactor (PWR) samples of isotopic measurements will expand the PWR sample population and consequently will greatly improve sample statistics. In turn, this will reduce uncertainty in the computational determination of the isotopic composition of commercial SNF and eliminate additional penalties currently used for lack of data. As part of the isotopic analysis task, eight diverse, high burnup (HBU) samples from several PWR rods are being dissolved and isotopically analyzed to provide high-quality measurement data for PWR SNF to reducing the uncertainties associated with PWR isotopic depletion validation. The isotopics of interest focus on nuclides important to burnup credit, shielding, and decay heat in PWR SNF. Additionally, eight high burnup BWR fuel samples from the Limerick nuclear power plant have been identified for measurement to support reducing uncertainties associated with BWR isotopic depletion validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Observing differential spin currents by resonant inelastic X-ray scattering

Controlling spin currents, that is, the flow of spin angular momentum, in small magnetic devices, is the principal objective of spin electronics, a main contender for future energy-efficient information technologies. A pure spin current has never been measured directly because the associated electric stray fields and/or shifts in the non-equilibrium spin-dependent distribution functions are too small for conventional experimental detection methods optimized for charge transport. Here we report that resonant inelastic X-ray scattering (RIXS) can bridge this gap by measuring the spin current carried by magnons—the quanta of the spin wave excitations of the magnetic order—in the presence of temperature gradients across a magnetic insulator. This is possible due to the sensitivity of the momentum- and energy-resolved RIXS intensity to minute changes in the magnon distribution under non-equilibrium conditions. Furthermore, we use the Boltzmann equation in the relaxation time approximation to extract transport parameters, such as the magnon lifetime at finite momentum, essential for the realization of magnon spintronics.

36 MATERIALS SCIENCE↗

EWTN: Quantifying Mass Transport to Enable Water Electrolyzer Architectures with Low Flow-Rate Sensitivity

To develop cost-effective and high-performing polymer electrolyte water electrolyzers (PEWEs) for gigawatt-scale applications, researchers have focused on reducing precious metal catalyst loadings and optimizing porous transport layers. However, the performance of PEWEs is also affected by system architecture. Mass transport losses are dependent on localized architecture and material interactions. In-situ measurements, such as current density distribution maps have demonstrated advantages in understanding the intricate characteristics and influence of two-phase flow within PEWEs. This study proposes the parameter of effective water transport number (EWTN) as a quantitative tool to investigate such current density distribution (CDD) measurements for PEWEs. Results show that higher flow-rates have EWTN values of 0.95 and above, indicating no mass transport limitations; while lower flow-rates with large gradients CDD have EWTN values between 0.6-0.8, indicating mass transport limited conditions. The new analysis also identified a correlation between mass transport losses due to bubble accumulation, membrane hydration, and ohmic overpotentials. To address these limitations, an unitized pin-type LGDL/flow-field design was developed, which effectively prevents local gas phase accumulation, resulting in improved mass transport characteristics. The results of this work show reduced flow-rate sensitivity with the pin-type architecture and ∼13% increased performance at 0.24 ml / min / cm 2 .

Electrochemistry↗

Real time trajectory optimization for hybrid energy management utilizing connected information technologies

The disclosure is directed to solving a full trajectory optimization problem in real-time for a hybrid electric vehicle (HEV) such that future driving conditions and energy usage may be fully considered in determining optimal engine energy usage and battery energy usage in real-time during a trip. An electronic control unit of the HEV may be configured to: receive route information for a route to be driven by the HEV; and after receiving the route information, iterating the operations of: measuring a current state of charge (SOC) of the battery; using at least the measured SOC and an initial co-state value stored in a memory, performing a process to iteratively update the co-state value to obtain an updated co-state value; using at least the updated co-state value, computing an updated control value; and applying the updated control value to control a usage of the battery and the internal combustion engine.

Huang, Mike X.↗