Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Testing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,603 records · Page 89

Current Source Gate Drive to Reduce Switching Loss for SiC MOSFETs

With conventional voltage source gate drives (VSG), the switching speed of SiC MOSFETs is difficult to increase due to large internal gate resistance, high Miller voltage, and limited gate voltage rating. This paper analyzes the requirement of current source gate drive (CSG) for SiC MOSFETs and proposes a CSG that can improve the switching speed and reduce switching loss. With the introduction of bi-directional switches, the influence of the large internal gate resistance of the SiC MOSFET can be mitigated, and sufficient gate current can be guaranteed throughout the switching transient. Therefore, the switching time and loss is reduced. The CSG can be controlled to be a VSG during steady state so the current of the gate drive is discontinuous and the stored energy of the inductor can be returned to the power supply to reduce gate drive loss. Double pulse tests are conducted for a SiC MOSFET with both conventional VSG and the proposed CSG. Testing results show that the switching loss of the proposed CSG is less than one third of the conventional VSG at full load condition.

Gui, Handong↗

Fuel Cell Stack Model for Real-Time Simulation of Grid-Connected Applications

Fuel cell stacks coupled with electrolyzers and hydrogen storage sites can be a promising category of distributed energy resources for both grid-connected and stand-alone power systems. However, because of high costs, at-scale hardware testing of fuel cell stacks for grid-connected applications is not economically viable at present. A model-based system that can accurately captures the steady-state and dynamic response of fuel cell stacks over long time periods (hours), is needed. This paper demonstrate a real-time electromagnetic transient model of a megawatt-scale, grid-connected proton exchange membrane fuel cell stack, coupled with a mass-based hydrogen storage system. This model can emulate the electrical steady-state and dynamic response of grid-connected fuel cell stacks. We validate the model using the response of commercial hardware fuel cell stacks and analytical models in the literature - using a digital real-time simulator (RSCAD). The proposed real-time model is then used to simulate cases spanning different time horizons and to design controller-hardware-in-the-loop experiments to evaluate controllers for hydrogen stations.

activation potential↗

MagicPlate-512: A 2D silicon detector array for quality assurance of stereotactic motion adaptive radiotherapy

Purpose: Spatial and temporal resolutions are two of the most important features for quality assurance instrumentation of motion adaptive radiotherapy modalities. The goal of this work is to characterize the performance of the 2D high spatial resolution monolithic silicon diode array named “MagicPlate-512” for quality assurance of stereotactic body radiation therapy (SBRT) and stereotactic radiosurgery (SRS) combined with a dynamic multileaf collimator (MLC) tracking technique for motion compensation. Methods: MagicPlate-512 is used in combination with the movable platform HexaMotion and a research version of radiofrequency tracking system Calypso driving MLC tracking software. The authors reconstruct 2D dose distributions of small field square beams in three modalities: in static conditions, mimicking the temporal movement pattern of a lung tumor and tracking the moving target while the MLC compensates almost instantaneously for the tumor displacement. Use of Calypso in combination with MagicPlate-512 requires a proper radiofrequency interference shielding. Impact of the shielding on dosimetry has been simulated by GEANT4 and verified experimentally. Temporal and spatial resolutions of the dosimetry system allow also for accurate verification of segments of complex stereotactic radiotherapy plans with identification of the instant and location where a certain dose is delivered. This feature allows for retrospective temporal reconstruction of the delivery process and easy identification of error in the tracking or the multileaf collimator driving systems. A sliding MLC wedge combined with the lung motion pattern has been measured. The ability of the MagicPlate-512 (MP512) in 2D dose mapping in all three modes of operation was benchmarked by EBT3 film. Results: Full width at half maximum and penumbra of the moving and stationary dose profiles measured by EBT3 film and MagicPlate-512 confirm that motion has a significant impact on the dose distribution. Motion, no motion, and motion with MLC tracking profiles agreed within 1 and 0.4 mm, respectively, for all field sizes tested. Use of electromagnetic tracking system generates a fluctuation of the detector baseline up to 10% of the full scale signal requiring a proper shielding strategy. MagicPlate-512 is also able to reconstruct the dose variation pulse-by-pulse in each pixel of the detector. An analysis of the dose transients with motion and motion with tracking shows that the tracking feedback algorithm used for this experiment can compensate effectively only the effect of the slower transient components. The fast changing components of the organ motion can contribute only to discrepancy of the order of 15% in penumbral region while the slower components can change the dose profile up to 75% of the expected dose. Conclusions: MagicPlate-512 is shown to be, potentially, a valid alternative to film or 2D ionizing chambers for quality assurance dosimetry in SRS or SBRT. Its high spatial and temporal resolutions allow for accurate reconstruction of the profile in any conditions with motion and with tracking of the motion. It shows excellent performance to reconstruct the dose deposition in real time or retrospectively as a function of time for detailed analysis of the effect of motion in a specific pixel or area of interest.

60 APPLIED LIFE SCIENCES↗

Deep Reinforcement Learning for Automatic Generation Control of Wind Farms

This paper provides a model-free framework for real-time control of wind farms to accurately track a power reference signal. This problem requires tractable dynamical models for capturing the aerodynamic interaction between wind turbines and controllers that can make decisions in realtime given varying atmospheric conditions. In this paper, we propose a deep reinforcement learning framework to provide real-time yaw control of a wind farm. Modifications have been made to FLOw Redirection and Induction in Steady State (FLORIS), a modeling tool that incorporates transient wake behavior. The control problem is formulated to track a synthetic power reference signal based on historical atmospheric (wind speed and direction) information, price signals, and regulation deployment data from U.S. regional transmission operators. Results indicate that a wind farm, with this control paradigm, can achieve good tracking performance when tested with real atmospheric data.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Experimental Validation of a Modular All-Electric Power Take-Off Topology for Wave Energy Converter Enabling Marine Renewable Energy Interconnection

Power electronic converters are an enabling technology for the emerging marine energy applications, such as using ocean waves to produce electricity. This paper outlines the power take-off system and its key components used in a wave energy converter offering modularity and scalability to generate power efficiently. The proposed power take-off system was implemented based on a modular multilevel converter and could be deployed to convert any alternating current electrical energy to a different alternating current for interconnection to grid or non-grid applications. Examples of widespread deployment are supplying electricity to coastal communities or producing clean drinking water. The analysis using both the simulation tests and laboratory experiments verified the design objectives and basic functionality of the developed power take-off system. An acceptable response using a field programmable gate array-based controlled laboratory testbench was achieved, complying with guidelines specified in the prevalent industry standards. Seamless operation during steady-state and transients for the studied wave energy converter was achieved as supported by the obtained results. The key findings of this work were experimentally examined under different load conditions, direct current bus voltage fluctuations, and generator speed–torque regulation. The ability of the power take-off system to generate high-power quality of the waveforms, e.g., against adhering to the IEEE 519-2022 standard for total harmonic distortion limits, is also confirmed.

Engineering↗

WUFI (Wärme and Feuchte Instationär)-Oak Ridge National Laboratory (ORNL)/Fraunhofer IBP

WUFI - Oak Ridge National Laboratory (ORNL)/Fraunhofer IBP is a menu-driven PC program which allows realistic calculation of the transient coupled one-dimensional heat and moisture transport in multi-layer building components exposed to natural weather. It is based on the newest findings regarding vapor diffusion and liquid transport in building materials and has been validated by detailed comparison with measurements obtained in the laboratory and on outdoor testing fields. Together with Oak Ridge National Laboratory (ORNL) Fraunhofer IBP has developed a special version of WUFI ® for North America. WUFI® ORNL is a functionally limited free version of WUFI® Pro for non-commercial purposes. It contains climate data for 62 cities in the USA and Canada which are all available in the free version. http://web.ornl.gov/sci/ees/etsd/btric/wufi/ http://www.WUFI.com/ORNL

Manfred Kehrer, ORNL↗

PHIL Validation of Ultracapacitor Storage for Black Start Application

A deck of slides that present the modeling, and calibration to conduct PHIL test of ultracapacitor energy storage system. The PHIL test has been done while preparing for the black start field demonstration with Idaho Falls Power. Presentation is invited to the 6th International Workshop on Grid Simulator Testing of Wind Turbine Power Trains and Other Renewable Technologies.

13 HYDRO ENERGY↗

The rate of extreme coronal line emitting galaxies in the Sloan Digital Sky Survey and their relation to tidal disruption events

High-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created by high-energy emission from active galactic nuclei and certain types of transients. In cases known as extreme coronal line emitting galaxies (ECLEs), these CLs are strong and fade away on a time-scale of years. The most likely progenitors of these variable CLs are tidal disruption events (TDEs), which produce sufficient high-energy emission to create and sustain the CLs over these time-scales. To test the possible connection between ECLEs and TDEs, we present the most complete variable ECLE rate calculation to date and compare the results to TDE rates from the literature. To achieve this, we search for ECLEs in the Sloan Digital Sky Survey (SDSS). We detect sufficiently strong CLs in 16 galaxies, more than doubling the number previously found in SDSS. Using follow-up spectra from the Dark Energy Spectroscopic Instrument and Gemini Multi-Object Spectrograph, Wide-field Infrared Survey Explorer mid-infrared observations, and Liverpool Telescope optical photometry, we find that none of the nine new ECLEs evolve in a manner consistent with that of the five previously discovered variable ECLEs. Using this sample of five variable ECLEs, we calculate the galaxy-normalized rate of variable ECLEs in SDSS to be $R_\mathrm{G}=3.6~^{+2.6}_{-1.8}~(\mathrm{statistical})~^{+5.1}_{-0.0}~(\mathrm{systematic})\times 10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$. The mass-normalized rate is $R_\mathrm{M}=3.1~^{+2.3}_{-1.5}~(\mathrm{statistical})~^{+4.4}_{-0.0}~(\mathrm{systematic})\times 10^{-17}~\mathrm{M_\odot ^{-1}}~\mathrm{yr}^{-1}$ and the volumetric rate is $R_\mathrm{V}=7~^{+20}_{-5}~(\mathrm{statistical})~^{+10}_{-0.0}~(\mathrm{systematic})\times 10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$. Our rates are one to two orders of magnitude lower than TDE rates from the literature, which suggests that only 10–40 per cent of all TDEs produce variable ECLEs. Additional uncertainties in the rates arising from the structure of the interstellar medium have yet to be included.

79 ASTRONOMY AND ASTROPHYSICS↗

Radio Frequency Transients Correlated with Electron Flux Measured On-Board the STP-Sat6

We analyzed transient electromagnetic signatures, referred to as spacecraft environment discharges (SEDs) or on-board discharges (OBDs), and the in-situ particle environment simultaneously measured by the Radio Frequency Sensor (RFS) and the Space and Atmospheric Burst Reporting System Z-plasma spectrometer (SABRS-ZPS), respectively. These sensors are on-board the U.S. Department of Defense Space Test Program Satellite 6 (STP-Sat6) in geostationary orbit. Radio frequency signatures of SEDs on STP-Sat6 have relatively high peak amplitudes; they often caused amplitude-saturation in the so-called high band (120–140 MHz frequency range) of the RFS likely due to the proximity of their source to the RFS’s dipole antennas. We examined the occurrence rate of SEDs reported by the RFS and its temporal relationship to electron flux in the 7.9–12.2 keV range. Peaks in SED count rates were time-correlated with peaks in electron flux. The electron-flux peaks appeared to lead the SED count rate peaks by 24–45 min for the 7.9–12.2 keV energies. This indicates that incidence of electrons at this relatively low-energy range likely resulted in ramping up of the charging of the spacecraft leading to periods of high occurrence rates of discharges in and around the spacecraft.

54 ENVIRONMENTAL SCIENCES↗

Numerical modeling of transient phenomena in the distant solar wind and in the heliosheath

In this paper we summarize our group's efforts to develop a self-consistent, data-driven model of the solar wind interaction with the local interstellar medium. In this model, we describe the motion of plasma with the MHD approach, while the transport of neutral atoms is described by either kinetic or multi-fluid equations. The model and its implementation in the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) is thoroughly tested and validated by comparing our results with other models and spacecraft measurements. The model appears to be rather accurate in reproducing the Interstellar Boundary Explorer ribbon and the general features of the plasma flow available from the Voyager 1 and 2 in-situ data. In particular, the timing of the Voyager crossing of the termination shock is reasonably accurately reproduced without involving exotically strong interstellar magnetic field values. We added a turbulence model to our governing equations and started treating pickup ions as a separate fluid, which made it possible to greatly improve our ability to reproduce the proton temperature. We discuss the results of our modeling of corotating interaction regions, global merged interaction regions, and solar cycle effects. A possible explanation is provided for the low radial velocity measurements from Voyager 1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Recent developments in catalyst-related PEM fuel cell durability

Cost and durability remain the two major barriers to the widespread commercialization of polymer electrolyte membrane fuel cell (PEMFC)-based power systems, especially for the most impactful but challenging fuel cell electric vehicle (FCEV) application. Commercial FCEVs are now on the road; however, their PEMFC systems do not meet the cost targets established by the U.S. Department of Energy, primarily due to the high platinum loading needed on the cathode to achieve the requisite performance and lifetime. Furthermore, while the activities of a number of commercial Pt-based alloy cathode catalysts exceed the beginning-of-life (BOL) targets, these activities, and the overall cathode performance, degrade via a variety of mechanisms described herein. Degradation is mitigated in current FCEVs by utilizing a cathode catalyst with a lower BOL activity (e.g., much lower transition metal alloy content and larger BOL nanoparticle size), necessitating higher catalyst loadings, and through the utilization of system controls that avoid conditions known to exacerbate degradation processes, such as limiting the fuel cell stack voltage range. The design and development of active and robust materials and eliminating the need for vehicle mitigation strategies would greatly simplify the operating system, allowing for greater transient operation, avoiding large hybridization, and curtailing of fuel cell power. Although system mitigation strategies have provided the near-term pathway for FCEV commercialization, material-specific solutions are required to further reduce costs and improve operability and efficiency. Future material developments should focus on stabilization of the electrode structure and minimization of the catalyst particle susceptibility to dissolution caused by oxide formation and reduction over PEMFC cathode-relevant operating potentials plus minimization of support corrosion. Ex situ accelerated stress tests have provided insight into the processes responsible for material and performance degradation and will continue to provide useful information on the relative stability of materials and benchmarks for robust and stable materials-based solutions not requiring system mitigation strategies to achieve adequate lifetime.

25 ENERGY STORAGE↗

MEXIICO: a new equipment to study the pressure impact on the irradiated fuel behavior

The new MEXIICO experimental loop has recently been started in the LECA-STAR hot-cell facility of the CEA Cadarache. It is dedicated to the understanding of the basic mechanisms occurring in nuclear fuels during nuclear reactor power transients and accidental situations. Under such conditions the nuclear fuels are submitted to thermo-mechanical stresses. These stresses play a major role on fission product release and fuel fragmentation. Evaluating such fuel behavior represents a key issue for reactor safety analyses. In the first part of the paper, we describe the complete MEXIICO device, which main component is a hot isostatic press allowing to perform analytical experiments in the [20 -1600 C. degrees / 10 MPa - 160 MPa] region. As shown by simulations (ALCYONE fuel performance code), the upper value of the pressure (160 MPa) is consistent with the one experienced by fuels under nuclear power transients. We also present the experimental loop starting from the argon pressurized feed-line, the furnace and its associated instrumentation and the exhaust line equipped with specific valves and flowmeters controlling the pressure and flowrate. We then focus on the gamma spectrometry measurement which allows recording the fission gas release kinetics. In the second part of the paper, we detail the capability of the MEXIICO device to provide new experimental data regarding the influence of pressure on the fission gas release and fuel fragmentation, and to determine a pressure threshold above which fission gas release and fuel fragmentation may be inhibited. This device also allows investigating the impact of pressure on the microstructure evolution, coupling it with post-test microanalysis examinations: this procedure will allow characterizing the effect of temperature and pressure on the bubble evolution and interconnection. Evolutions of the fuel modeling, needed to simulate correctly these new experimental data, will be also presented. In the last section of the paper, we present the commissioning tests of the MEXIICO device showing the good reproducibility of the tests performed first in conventional laboratory and then in hot-cell laboratory. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Tribological Characterization of SiC and B{sub 4}C Manufactured by Plasma Pressure Compaction

The objective of the present study was to assess the tribological properties and the wear behavior of silicon carbide and boron carbide ceramics, produced by the plasma pressure compaction sintering process. Reciprocating sliding tests were executed to evaluate the tribological behavior of silicon carbide and boron carbide. Samples of each material with distinct values of surface roughness Ra were tested. Wear mass loss, coefficients of friction and wear mechanism of the boron carbide and silicon carbide were evaluated by means of microscopes and an energy-dispersive spectroscopy probe. The coefficient of friction of the evaluated boron carbide was found to be 0.8, and the surface roughness affects only the starting transient, while for silicon carbide the steady-state value ranges from 0.5 to 0.8 depending on the roughness value Ra.

36 MATERIALS SCIENCE↗

A Unified Workflow for Sensitivity-Based Kinetic Analysis in Microkinetic Models

Degrees of rate control (DRC), apparent activation energies, and apparent reaction orders are established local sensitivity diagnostics for interpreting microkinetic models, but applying them routinely to large mechanisms often requires substantial reaction-specific bookkeeping, perturbation design, and postprocessing. Here, in this study, we present a unified derivative-based workflow that evaluates these quantities from a single compiled reaction-network model and target-rate definition. For any user-provided microkinetic model, the workflow compiles the mechanism into stoichiometrically consistent mass-action rate equations, solves the surface dynamics, and uses automatic differentiation to compute sensitivities with respect to rate constants, temperature, and gas partial pressures. By combining their calculations in the same framework, the workflow clearly demonstrates the relationships between different DRCs and the apparent activation energy. Using existing examples of propylene partial oxidation and methane oxidation on Pd(100), we verify expected transient redistribution of rate control, distinguish net Campbell DRCs from one-sided directional sensitivities, and show how apparent activation energy can be reconstructed either from one-sided DRCs or from state-based DRCs while critical mechanistic insights are obtained consistently. In the methane oxidation case, a pathway-subset test further illustrates how a simplified mechanism preserves key kinetic signatures of a full model, showing the potential of our user-friendly tool for model construction beyond kinetic analysis.

36 MATERIALS SCIENCE↗

Imaging Fracture Networks Using Joint Seismic and Electrical Change Detection Techniques (Final Report)

During the initial phase of this SubTER project, we conducted a series of high resolution seismic imaging campaigns designed to characterize induced fractures. Fractures were emplaced using a novel explosive source, designed at Sandia National Laboratories, that limits damage to the borehole. This work provided evidence that fracture locations could be imaged at inch scales using high-frequency seismic tomography but left many fracture properties (i.e. permeability) unresolved. We present here the results of the second phase of the project, where we developed and demonstrated emerging seismic and electrical geophysical imaging technologies that characterize 1) the 3D extent and distribution of fractures stimulated from the explosive source, 2) 3D fluid transport within the stimulated fracture network through use of a contrasting tracer, and 3) fracture attributes through advanced data analysis. Focus was placed upon advancing these technologies toward near real-time acquisition and processing in order to help provide the feedback mechanism necessary to understand and control fracture stimulation and fluid flow. Results from this study include a comprehensive set of 4D cross-hole seismic and electrical data that take advantage of change detection methodologies allowing for perturbations associated with the fracture emplacement and particulate tracer to be isolated. During the testing the team also demonstrated near real-time 4D electrical resistivity tomography imaging and 4D seismic tomography using the CASSM approach with a temporal resolution approaching 1 minute. All of the data collected were used to develop methods of estimating fracture attributes from seismic data, develop methods of assimilating disparate and transient data sets to improve fracture network imaging resolution, and advance capabilities for near real-time inversion of cross-hole tomographic data. These results are illustrated here. Advancements in these areas are relevant to all situations where fracture emplacement is used for reservoir stimulation (e.g. Enhanced Geothermal Systems (EGS) and tight shale gases).

42 ENGINEERING↗

HFIR High Power HEU Neutronics Analyses

Department of Energy National Nuclear Security Administration Office of Material Management and Minimization’s mission includes the conversion of civilian research reactors from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel. Analyses have shown that the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) will need to operate at 95 MW for the LEU silicide dispersion fuel designs to match key performance metrics obtained with HEU fuel at 85 MW. To prove safe operation of HFIR after installation of plant modifications to increase power, a high power HEU test cycle was proposed. Neutronics model updates and reactor physics analyses are performed to support the development of safety design reports for the high power (HP) HEU test cycle. Reactor physics metrics evaluated herein include fuel depletion, actinide production, cycle length, fission rate density distributions, axial power peaking factors, and reactor kinetic parameters. These reactor physics analyses support the development of future LEU safety design reports by providing key input for future HP HEU HFIR thermal hydraulics and reactor transient safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DEVELOPMENT OF A POPULATION BALANCE MODEL TO SIMULATE FRACTIONATION OF GROUND SWITCHGRASS

The population balance model represents a time-dependent formulation of mass conservation for a ground biomass that flows through a set of sieves. The model is suitable for predicting the change in size and distribution of ground biomass while taking into account the flow rate processes of particles through a grinder. This article describes the development and application of this model to a switchgrass grinding operation. The mass conservation formulation of the model contains two parameters: breakage rate and breakage ratio. A laboratory knife mill was modified to act as a batch or flow-through grinder. The ground switchgrass was analyzed over a set of six Tyler sieves with apertures ranging from 5.66 mm (top sieve) to 1 mm (bottom sieve). The breakage rate was estimated from the sieving tests. For estimating the breakage ratio, each of the six fractions was further ground and sieved to 11 fractions on a set of sieves with apertures ranging from 5.66 to 0.25 mm (and pan). These data formed a matrix of values for determining the breakage ratio. Using the two estimated parameters, the transient population balance model was solved numerically. Results indicated that the population balance model generally underpredicted the fractions remaining on sieves with 5.66, 4.00, and 2.83 mm apertures and overpredicted fractions remaining on sieves with 2.00, 1.41, and 1.00 mm apertures. These trends were similar for both the batch and flow-through grinder configurations. The root mean square of residuals (RSE), representing the difference between experimental and simulated mass of fractions, was 0.32 g for batch grinding and 0.1 g for flow-through grinding. The breakage rate exhibited a linear function of the logarithm of particle size, with a regression coefficient of 0.99.

09 BIOMASS FUELS↗

A Real-Time Vehicle-Specific Eco-Routing Model for On-Board Navigation Applications Capturing Transient Vehicle Behavior

The paper develops a real-time vehicle-specific eco-routing model for use in on-board navigation systems for both conventional fuel- and battery-powered electric vehicles. The model uses a novel link cost function that retains all microscopic transient behavior along a link by capturing within-link speed, acceleration, and road grade variations. Specifically, eight vehicle-agnostic link-specific variables are transmitted to the cloud and fused with information received from other connected vehicles. The fused data are then sent back to the connected vehicles and used to compute vehicle-specific eco-routes (edge computing). The solution produces an energy-efficient dynamic user-equilibrium traffic assignment that is solved incrementally. A numerical experiment is first designed to test the model, demonstrating the benefits associated with vehicle-specific eco-routes compared to those produced by traditional models. The model is then implemented in a micro-simulation framework to evaluate the network-wide impacts of the proposed energy-efficient user-optimum incremental traffic assignment for different traffic demand levels. Results demonstrate that the proposed model produces lower network-wide energy consumption levels compared to traditional energy-efficient and travel time user-optimum traffic assignment methods for both conventional fuel and electric vehicles. Conventional fuel savings relative to travel time routing decrease with increasing traffic demand levels, while electric power savings do not vary monotonically with congestion levels. Energy savings relative to traditional eco-routing are also demonstrated to not vary monotonically with demand levels. Results of the study also demonstrate that the transportation network layout significantly affects the eco-routing system performance. Finally, the feasibility of the proposed model in real-world applications is discussed and presented.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗