Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Averages”

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 487 records · Page 27

Effects of residue removal and tillage on greenhouse gas emissions in continuous corn systems as simulated with RZWQM2

Agricultural production is a major source of carbon dioxide (CO 2 ) and nitrous oxide (N 2 O) globally. The effects of conservation practices on soil CO 2 and N 2 O emissions remain a high degree of uncertainty. In this study, soil CO 2 and N 2 O emissions under different residue and tillage practices in an irrigated, continuous corn system, were investigated using the Root Zone Water Quality Model (RZWQM2). Combinations of no/high stover removal (NR and HR, respectively) and no-till/conventional tillage (NT and CT, respectively) field experiments were tested over the four crop-years (Apr. 2011–Apr. 2015). The model was calibrated using the NRCT, and validated with other treatments. The simulation results showed that soil volumetric water content (VWC) in the NR treatments (i.e., NRCT and NRNT) was 1.3%–1.9% higher than that in the HR treatments (i.e., HRCT and HRNT) averaged across the four years. A higher amount of CO 2 and N 2 O emissions were simulated in the NRCT across the four years (annual average: 7034 kg C/ha/yr for CO 2 and 3.8 kg N/ha/yr for N 2 O), and lower emissions were in the HRNT (annual average: 6329 kg C/ha/yr and 3.7 kg N/ha/yr for N 2 O). A long-term simulation (2001–2015) suggested that the CO 2 and N 2 O emissions were closely correlated with the stover removal degree (SRD), tillage, VWC, soil temperature (ST), years in management (Y), and fertilizer application. Stover and tillage practices had cumulative effects on CO 2 emissions. The simulated annual CO 2 emissions in 1st year from NRCT, NRNT, and HRCT were 7.8%, 0.0%, and 7.7% higher than that from HRNT, respectively; then the emissions in 15th year were 63.6%, 47.7%, and 29.1% higher, respectively. Meanwhile, there were no cumulative effects on N 2 O emissions. Furthermore, the results also demonstrated that the RZWQM2 is a promising tool for evaluating the long-term effects of CO 2 and N 2 O emissions on different conservation practices.

54 ENVIRONMENTAL SCIENCES↗

Electrochemical sulfate production from sulfide-containing wastewaters and integration with electrochemical nitrogen recovery

Electrochemical methods can help manage sulfide in wastewater, which poses environmental and health concerns due to its toxicity, malodor, and corrosiveness. In addition, sulfur could be recovered as fertilizer and commodity chemicals from sulfide-containing wastewaters. Wastewater characteristics vary widely among wastewaters; however, it remains unclear how these characteristics affect electrochemical sulfate production. Here, in this study, we evaluated how four characteristics of influent wastewaters (electrolyte pH, composition, sulfide concentration, and buffer strength) affect sulfide removal (sulfide removal rate, sulfide removal efficiency) and sulfate production metrics (sulfate production rate, sulfate production selectivity). We identified that electrolyte pH (3 × difference, i.e., 25.1 to 84.9 μM h –1 in average removal rate within the studied pH range) and sulfide concentration (16 × difference, i.e., 82.1 to 1347.2 μM h –1 in average removal rate) were the most influential factors for electrochemical sulfide removal. Sulfate production was most sensitive to buffer strength (6 × difference, i.e., 4.4 to 27.4 μM h –1 in average production rate) and insensitive to electrolyte composition. Together, these results provide recommendations for the design of wastewater treatment trains and the feasibility of applying electrochemical methods to varying sulfide-containing wastewaters. In addition, we investigated a simultaneous multi-nutrient (sulfur and nitrogen) process that leverages electrochemical stripping to further enhance the versatility and compatibility of electrochemical nutrient recovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A classical molecular dynamics simulation method for the formation of “dry” gels from boro-aluminosilicate glass structures

In contact with water, glass transforms into amorphous and porous structures called gels. A simulation method based on classical molecular dynamics is proposed here to mimic “dry” gels forming from initial oxide glass structures. Six glass compositions were investigated. Two behaviours were evidenced depending on the initial glass composition, and in particular on the quantity of elements removed. If a large quantity of soluble elements (B, Na) was removed, it induced an increase in the average pore size within the gels, and the time needed to stabilise the gel structure increased because more local atomic rearrangements occurred. The gel network displayed a higher proportion of Si-Q4 at the expense of Si-Q3 and a lower average ring size compared to the glass network, irrespective of the glass composition. Surface effects were also highlighted in the dry gels, such as the presence of 3-coordinated Al and a decrease in the average angle Si-O-Si and Al-O-Al. Our findings will be compared to both wet gels and experimental data in further studies, to help find the best procedure to simulate such structures.

36 MATERIALS SCIENCE↗

Radiation induced hardening of beryllium during low temperature He implantation

The effect of ion irradiation on evolution of microstructure and hardening of beryllium with different impurity levels was investigated using TEM and nanoindentation. High purity S-65 grade and less-pure S-200-F grade were implanted by helium ions at temperatures of 50°C and 200°C. 11 different energies were used, so as to create a quasi-homogeneous 3 µm irradiated layer with average radiation damage of 0.1 dpa and average He content of 2000 appm.Nanoindentation experiments demonstrated that before irradiation, the S-200-F and S-65 grades have an average hardness of 3.7$\pm$0.8 GPa and 3.4$\pm$0.8 GPa correspondently. After implantation the hardness of both grades increased by about 60% for the 200°C irradiation and 100% for the 50°C irradiation. The crystallographic analysis of indented grains demonstrated that in the as-received materials the hardness is about 2.5 times higher when the indentation direction is close to the [0001] c-axis of beryllium compared to indentation perpendicular to [0001]. Hardness anisotropy significantly decreased after irradiation: the “soft orientation” was most sensitive to irradiation-induced hardening, with hardness increasing by about 140% after irradiation at 50°C and 100% after irradiation at 200°C, compared to about 15 - 20% for the “hard” orientation at both irradiation temperatures. The higher purity grade had smaller increase of the “soft orientation” hardness: 2.5$\pm$0.3 GPa for the S-65 and 2.9$\pm$0.2 GPa for the S-200-F.At both temperatures in both grades, under TEM investigation the radiation damage appears as “black dots” which are likely to be small dislocation loops with the number density of ~ 10 22 m -3 . No bubbles were observed by TEM inside grains and at grain boundaries. Analysis of the possible hardening contribution demonstrated that the observed “black dots” could be responsible for up to half of the measured hardening, while the rest of the hardening should originate from helium bubbles with the size below the TEM resolution (at or below 1.5 nm).

36 MATERIALS SCIENCE↗

Determining the tensile strength of fuel surrogate TRISO-coated particle buffer, IPyC, and buffer-IPyC interlayer regions

A novel micro tensile sample fabrication technique for determining the tensile strength of the buffer, IPyC, and buffer -IPyC interlayer regions of surrogate (ZrO2) TRISO fuel particle layers was refined and implemented. Copper micro tensile samples served as baseline materials to verify the methods used. Tensile tests performed in this study, while limited in number, were analyzed using standard and Weibull statistics. As expected, the buffer layer was weakest, with an average ultimate tensile strength of 138.70 MPa, and the IPyC layer samples, were strongest, with an average ultimate tensile strength of 189.74 MPa. In the buffer -IPyC interface samples, all breaks occurred in the buffer region, though the average ultimate tensile strength of the samples, 159.80 MPa, was between the pure buffer and IPyC samples. These results suggest the interlayer region has unique properties, perhaps associated with pyrocarbon infiltration into the buffer layer during particle coating. All interlayer samples fractured within the buffer side; however, the stress strain behavior of some of these samples resembled the behavior of the IPyC layer samples. Here, the buffer and IPyC layer strengths had a normal distribution under Weibull analysis, while the interlayer region had a Rayleigh distribution. Further testing is needed to clarify both the standard and Weibull statistical results

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Combined selection of the dynamic model and modeling error in nonlinear aeroelastic systems using Bayesian Inference

Here, we report a Bayesian framework for concurrent selection of physics-based models and (modeling) error models. We investigate the use of colored noise to capture the mismatch between the predictions of calibrated models and observational data that cannot be explained by measurement error alone within the context of Bayesian estimation for stochastic ordinary differential equations. Proposed models are characterized by the average data-fit, a measure of how well a model fits the measurements, and the model complexity measured using the Kullback–Leibler divergence. The use of a more complex error models increases the average data-fit but also increases the complexity of the combined model, possibly over-fitting the data. Bayesian model selection is used to find the optimal physical model as well as the optimal error model. The optimal model is defined using the evidence, where the average data-fit is balanced by the complexity of the model. The effect of colored noise process is illustrated using a nonlinear aeroelastic oscillator representing a rigid NACA0012 airfoil undergoing limit cycle oscillations due to complex fluid–structure interactions. Several quasi-steady and unsteady aerodynamic models are proposed with colored noise or white noise for the model error. The use of colored noise improves the predictive capabilities of simpler models.

42 ENGINEERING↗

Evaluation of the Prompt Fission Gamma Properties for Neutron Induced Fission of 235,238 U and 239 Pu

Herein we present details of the prompt fission gamma property evaluations included in the released ENDF/B-VIII.0 files. The average prompt gamma multiplicity, the average prompt gamma spectrum, the average total prompt gamma energy released and the prompt gamma multiplicity distributions, included in the ENDF/B-VIII.0 files, are presented for 235 U(n,f), 238 U(n,f), and 239 Pu(n,f) reactions, as a function of incident neutron energy up to 30 MeV. The evaluation is based on available experimental data and model calculations. Notable data include recent measurements by experimental groups at Los Alamos and Lawrence Livermore National Laboratories, at Geel and Budapest, at CEA labs, as well as seminal historical gamma production measurements by Drake at Los Alamos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimizing pressurized-water reactor equilibrium cycle using a novel loading pattern encoding and rule-based genetic crossover operators

This work presents an extended multi-batch approach applied in shuffling scheme optimization for equilibrium cycle for pressurized water reactors using Genetic Algorithms (GAs). A new ruled based GA crossover operator called Inherited Location and Batch (ILB) was introduced to enhance offsprings reproduction efficiency specialized for equilibrium cycle optimization problem. This approach was implemented within the Plant ReLoad Optimization (PRLO) framework and validated using a generic reactor model based on the AP1000 design, with core parameters calculated via the CASMO/SIMULATE software package. The ILB approach is then applied for both single and multi-objective problems in maximizing cycle length and core average exposure while minimizing the average enrichment of the 57 fresh fuel assemblies (FAs) per cycle. The optimal solutions are selected based on their dominance to the objectives from all feasible solutions. This research identified three optimal solutions satisfied safety constraints: The first solution minimizes feed enrichment costs with a cycle length of 338.8 days and core exposure of 25.39 MWd/MT; the second solution extends cycle length to 361.2 days, with the highest core exposure of 26.84 MWd/MT, using 3.75 wt% average fuel enrichment; the third solution balances both objectives with a cycle length of 349.6 days, core exposure of 25.82 MWd/MT with a slight enrichment increase compared to the first solution. Collectively, these findings underscore the efficiency and effectiveness of the proposed approach in achieving practical multi-objective optimal equilibrium cycle designs using GAs optimizer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Sensitivity analysis of thermal contact conductance modeling to inform MiniFuel irradiation capsule designs

The MiniFuel irradiation platform has been developed by Oak Ridge National Laboratory as a flexible, high-throughput separate effects testing capability within the High Flux Isotope Reactor (HFIR). Finite element thermal models are relied upon to design MiniFuel experiments to achieve a specific time-averaged irradiation temperature for experimental objectives. A previous study identified that uncertainty in the component heat generation rates and thermal contact conductance (TCC) model are the most significant contributors to predicted fuel temperature variance. To address both sources of uncertainty, this work performs sensitivity analysis on the TCC model to identify high-impact, high-uncertainty parameters that contribute to fuel temperature variance. The TCC model is analyzed in increasing detail, first using a standalone Python code, then again after coupling Python to the BISON fuel performance code. Furthermore, the parameters with the largest contributions to fuel temperature variance which can be reduced through design changes are identified as the initial subcapsule gas pressure, contact pressure between the fuel and dish, and the effective surface roughness of the interface. A set of design recommendations for future capsule designs has been established and applied to reduce the previously quantified average fuel temperature uncertainty ranges of ± 40 °C in the HFIR vertical experiment facilities (VXF) and ± 80 °C in the removable beryllium (RB) reflector to approximately ± 32 °C and ± 53 °C, respectively. This equates to a 21 % and 33 % reduction in the uncertainty range of the average fuel temperature for VXF and RB, respectively.

BISON↗

An electrochemical generator for the continual supply of 213 Bi from 225 Ac for use in targeted alpha therapy applications

Bismuth-213 is a radionuclide of interest for targeted alpha therapy and is supplied via a radiochemical generator system through the decay of 225 Ac. Radionuclide generators employ longer lived “parent” radionuclides to routinely supply shorter-lived “daughter” radionuclides. The traditional 225 Ac/ 213 Bi radiochemical generator relies on an organic cation exchange resin where 225 Ac binds to the resin and 213 Bi is routinely eluted. These resins degrade when they absorb large doses of ionizing radiation (>1 × 10 6 Gy/mg), which has been observed when the loading activity of 225 Ac exceeds 2.59*10 9 Bq (70 mCi). Herein we report the development of an electrochemical generator for the supply of 213Bi that has the potential to overcome this limitation. Bismuth-213 spontaneously electrodeposits onto nickel foils in 0.1 M hydrochloric acid at 70 °C. Using this method, we were able to plate an average of 73 ± 4 % of the 213 Bi in solution and obtain a final 213 Bi recovery of 65 ± 8 % in 0.1 M citrate pH 4.5 via reverse electrolysis using titanium as the cathode. The recovered 213Bi had an average radiochemical purity of >99.8 % and was successfully used to radiolabel DOTATATE with an average radiochemical yield of 85.1 % (not optimized).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-speed synchrotron X-ray imaging of directed energy deposition of titanium: effects of processing parameters on the formation of entrapped-gas pores

Laser based directed energy deposition (DED) is a competitive method for repairing and remanufacturing metallic parts used in numerous industries including aerospace and biomedical. However, the numerous dynamic phenomena associated with the DED process often result in defects such as entrapped-gas pores, lack of fusion, and undesirable anisotropic properties. The entrapped-gas pore, being one of the most common issues, not only influences melt-pool dynamics but also reduces the fabrication quality and mechanical properties of parts fabricated by the DED process. To reduce and further understand this issue, the real-time observation of the pore formation process needs to be studied first. To directly observe the phenomena in the melt pool, high-speed techniques are needed because rapid solidification leads to rapid pore formation and movement. In-situ high-speed X-ray has been proven to be an effective method in investigating the melt pool dynamics and pore formation mechanisms in the laser powder bed fusion process, in which the fabrication process is quite different from that in DED. Here, the high-speed X-ray method is extended to study the formation of entrapped-gas pores. The real-time formation and quantitative analysis of pores under each set of processing parameters (particle velocity, laser power, and spot welding dwelling time of stationary laser) in the DED process are investigated. We found that the DED with a higher particle velocity (3.19 m/s) produced a smaller average pore size of 27.8 µm and a lower pore area fraction of 0.52%. The DED under lower laser power (156 W) generated a smaller average pore size of 20.3 µm and a lower pore area fraction of 1.94%. The shorter dwelling time (10 ms) benefited the decrease of both average pore size and pore area fraction.

Defect↗

Dose conversion factors for absorbed dose in a mobile phone to absorbed dose in critical organs in an anthropomorphic phantom for emergency dosimetry applications: OSL and TL experimental results, and Monte Carlo simulations

Here, an anthropomorphic phantom was equipped with three smartphones and irradiated in five separate exposure geometries in air using a 137 Cs source. The phantom was fitted internally with thermoluminescence dosimeters (TLD-700 chips) to measure the absorbed dose to the individual organs, while the absorbed doses in the phones were evaluated using optically stimulated luminescence (OSL) from the phones’ surface mount resistors (SMRs), along with reference dosimeters, namely Luxel OSL films and TLD-700 chips. The results indicate a strong dependence on the phone location with respect to the source due to shielding and scattering from the body. Dose conversion factors were calculated for several plausible exposure scenarios in order to correlate the absorbed doses measured in the phones to the average body dose received by the phantom. The measured conversion factors were compared with calculated factors using MCNP6.2 Monte Carlo simulations and an ADAM voxel phantom. Experiments were conducted to assess how the phone orientation and angle impacts the dose received by the phone components. The overall results indicate that the doses measured can vary significantly depending upon exposure geometry (including phone position and angle). However, in rotational irradiation symmetry, the phone doses lead to an over-estimate of the average whole-body dose by an average of ~10%, independent of the phone locations tested in this study. This result, if confirmed in future studies, may be acceptable for triage.

61 RADIATION PROTECTION AND DOSIMETRY↗

Clear-sky detection for PV degradation analysis using multiple regression

A method is presented to detect clear-sky periods for plane-of-array irradiance time-averaged data that is based on the algorithm originally described by Reno and Hansen. Here we show this new method improves the state-of-the-art by providing accurate detection at longer data averaging intervals. Moreover, our new method detects clear periods in plane-of-array data, which is novel. The new method is developed by applying a Design of Experiment approach to optimize the parameters used in the Reno method, and Monte Carlo simulations are used to understand the robustness of the found parameters. Clear-sky detection accuracy is compared among four methods: the Reno method, the default clear-sky filter in RdTools, the Ellis method, and the method outlined in this work, using a hand-labeled two-year data set of 1-min plane-of-array irradiance for a fixed tilt system. The RdTools clear-sky filter is marred by excessive false positives. The other methods all perform well at 1-min data intervals; the method developed here provides more accurate detection at longer data averaging intervals. We show that the parameters are directly linked to the data frequency in the hope that these input variables may not have to be optimized for every data frequency and location. However, only a single fixed system in one location was carefully examined. Finally, we illustrate how accurate determination of clear-sky conditions helps to eliminate data noise and bias in the assessment of long-term performance of PV plants.

14 SOLAR ENERGY↗

Generating high-resolution total canopy SIF emission from TROPOMI data: Algorithm and application

Solar-induced chlorophyll fluorescence (SIF) is a rapidly advancing front in modeling global terrestrial gross primary production (GPP). Canopy total SIF emissions (SIF total ) are mechanistically linked to the plant photosynthesis, and can be estimated from satellite observed SIF (SIF obs ) through radiative transfer modeling. However, the current satellite SIF obs and thus SIF total are available only at coarse spatial resolutions from several kilometers to tens of kilometers, inhibiting the application at fine spatial scales. Here, in this work, we proposed an algorithm to generate both global high-resolution SIF total (HSIF total ) and high-resolution SIF obs (HSIF obs ) at 1 km from low-resolution SIF obs (LSIF obs ) from the TROPOspheric Monitoring Instrument (TROPOMI), which has a spatial resolution at nadir of 3.5 km by 5.6–7 km. Our statistical method is based on the law of energy conservation and uses satellite derived fraction of absorbed photosynthetically active radiation, fluorescence efficiency, and the escape probability of fluorescence. We evaluated the accuracy of our HSIF total using the Orbiting Carbon Observatory-2 SIF (R 2 = 0.78). We found that the spatial resolution had clear effects on the relationship between HSIF total and GPP. We also compared HSIF total to 8-day averaged tower GPP from 135 flux sites and found that they were better correlated when HSIF total was averaged over a 1-km radius around the tower than when averaged over a larger radius. Our study provided a unique high-resolution HSIF total product, which will advance the estimation of GPP by extrapolating site-level relationships to the global scale.

54 ENVIRONMENTAL SCIENCES↗

Long-term spatial and temporal solar resource variability over America using the NSRDB version 3 (1998–2017)

The study assesses the long-term spatial and temporal solar resource variability in America using the 20-year National Renewable Energy Laboratory's (NREL's) National Solar Radiation Database (NSRDB). Specifically, the coefficient of variation (COV) is used to analyze the spatial and temporal (interannual and seasonal) variability. Further, both spatial and temporal long-term variability are analyzed using the Köppen-Geiger climate classification. The temporal variability is found that, on average, the continental United States (CONUS) COV reaches up to 5% for global horizontal irradiance (GHI) and 10% for direct normal irradiance (DNI), and that the NSRDB domain's COV is roughly twice that of CONUS. For the seasonal variability analysis, the winter months are found to exhibit higher COV than the other seasons. In particular, December exhibits the highest variability, reaching on average 30% for DNI and 20% for GHI over various areas. On the other hand, the summer months demonstrate significantly lower variability, reaching only less than 20% for DNI and 10% for GHI, on average. Similarly, the spatial variability is analyzed by comparing each pixel to its neighbors. The long-term spatial variability is found to increase with the number of neighboring pixels being considered, which is equivalent to an increase in distance (within a 100-km x 100-km square grid). As expected, the DNI spatial variability is higher than that of GHI. Moreover, the annual solar irradiance anomalies are found to reach ±25% for both GHI and DNI (and even exceed those value in some instances) during each year of the 20-year period.

14 SOLAR ENERGY↗

Responses of particulate and mineral-associated organic carbon to temperature changes and their mineral protection mechanisms: A soil translocation experiment

Mineral protection mechanisms are important in determining the response of particulate organic carbon (POC) and mineral-associated organic carbon (MAOC) to temperature changes. However, the underlying mechanisms for how POC and MAOC respond to temperature changes are remain unclear. Here, by translocating soils across 1304 m, 1425 m and 2202 m elevation gradient in a temperate forest, simulate nine months of warming (with soil temperature change of +1.41 °C and +3.91 °C) and cooling (with soil temperature change of −1.86 °C and −4.20 °C), we found that warming translocation significantly decreased POC by an average of 10.84 %, but increased MAOC by an average of 4.25 %. Conversely, cooling translocation led to an average increase of 8.64 % in POC and 13.48 % in MAOC. Exchangeable calcium (Ca exe ) had a significant positive correlation with POC and MAOC during temperature changes, and Fe/Al-(hydr)oxides had no significant correlation or a significant negative correlation with POC and MAOC. Our results showed that POC was more sensitive than MAOC to temperature changes. Ca exe mediated the stability of POC and MAOC under temperature changes, and Fe/Al-(hydr)oxides had no obvious protective effect on POC and MAOC. Our results support the role of mineral protection in the stabilization mechanism of POC and MAOC in response to climate change and are critical for understanding the consequences of global change on soil organic carbon (SOC) dynamics.

Mineral protection↗

Siting and sizing of public–private charging stations impacts on household and electric vehicle fleets

To facilitate the provision of electric vehicle charging stations (EVCS) in urban areas, this study investigates the benefits of co-locating fleet-owned chargers with public charging stations to enable construction incentives and cord-sharing cost savings. Shared EVCS can serve charging demand from both user types: private (household) EV owners and those managing fleet vehicles – like shared and fully automated EV (SAEV) fleets. Using POLARIS to simulate all person-travel across the 6-county Austin, Texas region, new EVCS were sited and sized with DC fast-charging (DCFC) plugs to lower operating and construction costs while providing public + private (PP) service across an 81-square-mile core geofence (where 200 SAEVs were active) over 24-hour days. When co-location is permitted, 115 DCFC cords were added to the 23 existing (publicly available) stations to enable SAEVs and household EVs (HHEVs) charging access, within the geofence. Each 250-mile-range SAEV was simulated to travel an average of 330 miles per day, serve over 92 person-trips, and recharge 2.7 times a day (for 2.4 h per session). The new DCFC plugs were primarily added to public EVCS at shopping centers and schools, and in residential settings along freeways. The average plug served 4.8 EVs per day. Most co-located PP EVCS permitted immediate (no-wait) charging, except for 2 stations along freeways that averaged 8 min of wait time to begin charging. In conclusion, the co-location strategy lowered fleet owners’ initial EVCS construction costs by 12 % (thanks to cord-sharing to avoid cord duplication), while reducing SAEV wait times to just 3.1 min (versus 10.7 min if SAEV managers had to build and operate their own EVCS).

EV charging modeling↗