Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “well data”

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 379 records · Page 21

A phosphorylation-deficient ribosomal protein eS6 is largely functional in Arabidopsis thaliana , rescuing mutant defects from global translation and gene expression to photosynthesis and growth

The eukaryote-specific ribosomal protein of the small subunit eS6 is phosphorylated through the target of rapamycin (TOR) kinase pathway. Although this phosphorylation event responds dynamically to environmental conditions and has been studied for over 50 years, its biochemical and physiological significance remains controversial and poorly understood. Here, we report data from Arabidopsis thaliana, which indicate that plants expressing only a phospho-deficient isoform of eS6 grow essentially normally under laboratory conditions. The eS6z (RPS6A) paralog of eS6 functionally rescued a double mutant in both rps6a and rps6b genes when expressed at approximately twice the wild-type dosage. A mutant isoform of eS6z lacking the major six phosphorylatable serine and threonine residues in its carboxyl-terminal tail also rescued the lethality, rosette growth, and polyribosome loading of the double mutant. This isoform also complemented many mutant phenotypes of rps6 that were newly characterized here, including photosynthetic efficiency, and most of the gene expression defects that were measured by transcriptomics and proteomics. However, compared with plants rescued with a phospho-enabled version of eS6z, the phospho-deficient seedlings retained a mild pointed-leaf phenotype, root growth was reduced, and certain cell cycle-related mRNAs and ribosome biogenesis proteins were misexpressed. The residual defects of the phospho-deficient seedlings could be understood as an incomplete rescue of the rps6 mutant defects. There was little or no evidence for gain-of-function defects. As previously published, the phospho-deficient eS6z also rescued the rps6a and rps6b single mutants; however, phosphorylation of the eS6y (RPS6B) paralog remained lower than predicted, further underscoring that plants can tolerate phospho-deficiency of eS6 well. Our data also yield new insights into how plants cope with mutations in essential, duplicated ribosomal protein isoforms.

60 APPLIED LIFE SCIENCES↗

Census of the Young Stellar Population in the Galactic H ii Region Sh2-242

We present here identification and characterization of the young stellar population associated with an active star-forming site Sh2-242. We used our own new optical imaging and spectroscopic observational data, as well as several archival catalogs, e.g., Pan-STARRS 1, Gaia DR2, IPHAS, Wide-field InfraRed Camera, 2MASS, and Spitzer. Slit spectroscopic results confirm the classification of the main ionizing source BD+26 980 as an early-type star of spectral type B0.5 V. The spectrophotometric distance of the star is estimated as 2.08 ± 0.24 kpc, which confirms the source as a member of the cluster. An extinction map covering a large area (diameter ∼50′) is generated with H and K photometry toward the region. From the map, three distinct locations of peak extinction complexes (A {sub V} ≃ 7–17 mag) are identified for the very first time. Using the infrared color excess, a total of 33 Class I and 137 Class II young objects are classified within the region. The IPHAS photometry reveals classification of 36 Hα emitting sources, which might be class II objects. Among 36 Hα emitting sources, 5 are already identified using infrared excess emission. In total, 201 young objects are classified toward S242 from this study. The membership status of the young sources is further windowed with the inclusion of parallax from the Gaia DR2 catalog. Using the optical and infrared color–magnitude diagrams, the young stellar objects are characterized with an average age of ∼1 Myr and masses in the range 0.1–3.0 M {sub ☉}. The census of the stellar content within the region is discussed using combined photometric and spectroscopic data.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating high-energy proton-induced reactions on spherical nuclei: Implications for the preequilibrium exciton model

A number of accelerator-based isotope production facilities utilize $100-200$ MeV proton beams due to the high production rates enabled by high-intensity beam capabilities and the greater diversity of isotope production brought on by the long range of high-energy protons. However, nuclear reaction modeling at these energies can be challenging because of the interplay between different reaction modes and a lack of existing guiding cross section data. A Tri-lab collaboration has been formed between the Lawrence Berkeley, Los Alamos, and Brookhaven National Laboratories to address these complexities by characterizing charged-particle nuclear reactions relevant to the production of established and novel radioisotopes. In the inaugural collaboration experiments, stacked-targets of niobium foils were irradiated at the Brookhaven Linac Isotope Producer ($E_p=200$ MeV) and the Los Alamos Isotope Production Facility ($E_p=100$ MeV) to measure $^{93}$Nb(p,x) cross sections between $50-200$ MeV. The results were compared with literature data as well as the default calculations of the nuclear model codes TALYS, CoH, EMPIRE, and ALICE. The default code predictions largely failed to reproduce the measurements. Therefore, we developed a standardized procedure, which determines the reaction model parameters that best reproduce the most prominent reaction channels in a physically justifiable manner. Overall, the primary focus of the procedure was to determine the best parameterization for the pre-equilibrium two-component exciton model. This modeling study revealed a trend towards a relative decrease for internal transition rates at intermediate proton energies ($E_p=20-60$ MeV) in the current exciton model as compared to the default values. The results of this work are instrumental for the planning, execution, and analysis essential to isotope production.

43 PARTICLE ACCELERATORS↗

Evaluation of “Tight Oil” Well Performance and Completion Practices in the Powder River Basin - “Time Slice” Analysis

The objective of this paper is to assess how well completion practices and well performance have evolved over time (“time slices”) in the tight sands and shales of the Powder River Basin (PRB). This information can provide a foundation for helping operators define more effective well completion practices and thus optimize well performance in an emerging “tight oil” basin. To start, the authors assembled a “database” containing well completion practices, production data, and geologic information for more than 800 horizontal (Hz) wells targeting three “tight oil” formations -- Turner, Frontier, and Mowry--placed on production in the past dozen years. To control for the impact of geologic and reservoir properties on well performance, the authors defined 12 geologically distinct areas (“partitions”) for the Turner, Frontier, and Mowry Shale in the Powder River Basin (4 partitions in each formation). This paper discusses the methodology for establishing “time slices” for each partition within each of the three “tight oil” formations, including (1) taking out, to the extent practical, the effects of geology by partitioning the three “tight oil” formations based on their geologic parameters; (2) using type curves and estimated ultimate recoveries (EURs) to establish a reliable measure of well performance; and (3) rigorously inspecting the production and well completion data to assure a quality dataset. Partition #2, a high thermal maturity area of the Frontier Sandstone in the western Powder River Basin, helps illustrate the results of the study. Using a dataset of 55 Hz wells, the study found that well performance has steadily improved, with oil EURs increasing from 180 MBbl in 2012-13 to 290 MBbl in 2018-19. Much of this improvement was due to the use of longer Hz laterals, increasing notably from 3,765 ft in 2012-13 to 10,260 ft in 2018-19. More intensive completion practices contributed, as well. For example, the number of frac stages more than doubled, from 17 to 35, and proppant concentrations increased from 1,060 lbs/ft to 1,320 lbs/ft. Finally, the data show diminishing returns to longer laterals and more intensive completion practices. For example, the key performance measure of oil EUR per 1,000 ft of lateral is decreasing. For 2012-13, Hz wells recovered 50 MBbl of oil per 1,000 ft of lateral while the more recent 2018-19 Hz wells only provide 30 MBbl of oil per 1,000 ft of lateral. Considerable insight can be gained by using “time slices” and geologic partitioning to better understand the relationship between changes in well drilling and completion practices and changes in well performance in emerging “tight oil” plays. The results from this study can also serve as a foundation for subsequent, more intense efforts involving data analytics for defining more optimum well completion practices targeting specified geologic settings and formations.

02 PETROLEUM↗

Bemnifosbuvir: An HCV NS5B Inhibitor With Multiple Modes of Action

Bemnifosbuvir (BEM) is a potent, pan-genotypic inhibitor targeting the hepatitis C virus (HCV) NS5B polymerase. Its antiviral activity was evaluated in an ascending dose phase I clinical trial involving 30 patients treated once a day for 7 days. After treatment initiation, plasma HCV RNA declined in a biphasic manner with a mean reduction of 2.3 log 10 IU/mL after 24 hours and 4.4 log10 IU/mL by day 7 for the highest dose. Alanine aminotransferase (ALT) also normalized in most patients. We employed a multiscale mathematical model fitted to the HCV RNA and ALT dynamics to quantify the antiviral activity and evaluate the modes of action of BEM. We found that models in which BEM only acted as a typical HCV RNA polymerase inhibitor and reduced the intracellular production of HCV RNA did not fit the data as well as models in which BEM had multiple modes of action, including suppressing viral assembly and secretion and enhancing intracellular HCV RNA degradation. BEM's effectiveness in inhibiting intracellular HCV RNA production increased with dose (150 mg/day: 88.2%, 300 mg/day: 98.8%, 600 mg/day: 99.5%), while inhibition of viral assembly and release was ~95% effective regardless of dose. We observed a dose-dependent enhancement in the degradation of intracellular HCV RNA, with degradation rates 1.5-fold higher in patients receiving 300 mg/day and 2.7-fold higher in those receiving 600 mg/day than in patients receiving 150 mg/day. No significant differences in antiviral activity were detected between HCV genotypes 1b and 3 or between patients with and without compensated cirrhosis.

59 BASIC BIOLOGICAL SCIENCES↗

Valence-to-core X-ray emission spectroscopy of titanium compounds using energy dispersive detectors

X-ray emission spectroscopy (XES) of transition metal compounds is a powerful tool for investigating the spin and oxidation state of the metal centers. Valence-to-core (vtc) XES is of special interest, as it contains information on the ligand nature, hybridization, and protonation. To date, most vtc-XES studies have been performed with high-brightness sources, such as synchrotrons, due to the weak fluorescence lines from vtc transitions. Here, we present a systematic study of the vtc-XES for different titanium compounds in a laboratory setting using an X-ray tube source and energy dispersive microcalorimeter sensors. With a full-width at half-maximum energy resolution of approximately 4 eV at the Ti Kβ lines, we measure the XES features of different titanium compounds and compare our results for the vtc line shapes and energies to previously published and newly acquired synchrotron data as well as to new theoretical calculations. Finally, we report simulations of the feasibility of performing time-resolved vtc-XES studies with a laser-based plasma source in a laboratory setting. Our results show that microcalorimeter sensors can already perform high-quality measurements of vtc-XES features in a laboratory setting under static conditions and that dynamic measurements will be possible in the future after reasonable technological developments.

36 MATERIALS SCIENCE↗

A practical approach to using the Genomic Standards Consortium MIxS reporting standard for comparative genomics and metagenomics

Comparative analysis of (meta)genomes necessitates aggregation, integration, and synthesis of well-annotated data using standards. The Genomic Standards Consortium (GSC) collaborates with the research community to develop and maintain the Minimal Information about any (x) Sequence (MIxS) reporting standard for genomic data. To facilitate use of the GSC’s MIxS reporting standard, we provide a description of the structure and terminology, how to navigate ontologies for required terms in MIxS, and demonstrate practical usage through a soil metagenome example.

standards, metadata, genome, metagenome, schema, v↗

Asymptotic behavior of meson transition form factors

One of the open issues in evaluations of the contribution from hadronic light-by-light scattering to the anomalous magnetic moment of the muon (g - 2)µ concerns the role of heavier scalar, axial-vector, and tensormeson intermediate states. The coupling of axial vectors to virtual photons is suppressed for small virtualities by the Landau–Yang theorem, but otherwise there are few rigorous constraints on the corresponding form factors. In this paper, we first derive the Lorentz decomposition of the two-photon matrix elements into scalar functions following the general recipe by Bardeen, Tung, and Tarrach. Based on this decomposition, we then calculate the asymptotic behavior of the meson transition form factors from a light-cone expansion in analogy to the asymptotic limits for the pseudoscalar transition form factor derived by Brodsky and Lepage. Finally, we compare our results to existing data as well as previous models employed in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement of the ${\Lambda}_c^{+}$ to D 0 production ratio in peripheral PbPb collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV

We report on a measurement of the ${\Lambda}_c^{+}$ to D 0 production ratio in peripheral PbPb collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV with the LHCb detector in the forward rapidity region 2 < y < 4.5. The ${\Lambda}_c^{+}$ (D 0 ) hadrons are reconstructed via the decay channel ${\Lambda}_c^{+}$ → pK – π + (D 0 → K – π + ) for 2 < p T < 8 GeV/c and in the centrality range of about 65–90%. The results show no significant dependence on p T , y or the mean number of participating nucleons. They are also consistent with similar measurements obtained by the LHCb collaboration in pPb and Pbp collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV. The data agree well with predictions from PYTHIA in pp collisions at $\sqrt{s}$ = 5 TeV but are in tension with predictions of the Statistical Hadronization model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Integrated positron emission particle tracking (PEPT) and X-ray computed tomography (CT) imaging of flow phenomena in twisted tape swirl flow

Abstract A combined positron emission particle tracking (PEPT) and X-ray computed tomography (CT) technique is presented, and its utility is demonstrated through investigation of flow in a pipe with twisted tape swirl insert with varying flow conditions (diameter-based Reynolds numbers 16,300–63,300). A description of this technique is given, as well as data handling practices used to relate geometric information captured by CT to fluid flow data gathered via PEPT. It is found that the CT component is readily capable of capturing the stainless steel insert geometry in this present system, but the use of combined plastic and metal materials leads to artifacts in imaging of the plastic surface. Nonetheless, CT data are related to PEPT flow measurements, and average velocity fields are calculated via a pseudo-framing and interpolation scheme and used to visualize and interrogate key flow phenomena within the system. Radial velocity profiles of the mean flow characteristics are seen to collapse to a nearly common form across all flow conditions considered. Helical vortices are seen propagating through the flow field, generated by bypass flow around the gap between the insert and pipe wall, with additional coherent secondary flow structures seen in the higher Reynolds number cases. These findings enhance the understanding of the mixing mechanisms in these swirl flows and encourage the continued development of PEPT-CT methodologies for 3D flow measurements in optically inaccessible systems.

42 ENGINEERING↗

Uncertainty Quantification of Machine Learning Predicted Creep Property of Alumina-Forming Austenitic Alloys

The development of machine learning (ML) approaches in materials science offers the opportunity to exploit existing engineering and developmental alloy datasets, such as Oak Ridge National Laboratory (ORNL)’s consistently measured creep-rupture dataset for alumina-forming austenitic (AFA) alloys, to accelerate their further development. As a first step toward achieving ML insights for improved alloy design, the potential sources of uncertainty and their impacts on ML output are examined. It is observed that the selection of algorithms and features as well as data sampling significantly affects the performance of ML models, either positively or negatively. Further, the performance of various ML models in predicting the creep properties of AFA alloys is compared, with further evaluation by assessment of a small set of new developmental AFA alloys that were not part of the training dataset. The present study demonstrates that uncertainty quantification (UQ) is essential in materials science for evaluating the performance of ML algorithms with specifically selected feature sets and obtaining a comprehensive understanding of their limitations and the resultant capability of effective prediction in complex materials systems.

36 MATERIALS SCIENCE↗

Tracking Improvement in Simulated Marine Biogeochemistry Between CMIP5 and CMIP6

The changes or updates in ocean biogeochemistry component have been mapped between CMIP5 and CMIP6 model versions, and an assessment made of how far these have led to improvements in the simulated mean state of marine biogeochemical models within the current generation of Earth system models (ESMs). The representation of marine biogeochemistry has progressed within the current generation of Earth system models. However, it remains difficult to identify which model updates are responsible for a given improvement. In addition, the full potential of marine biogeochemistry in terms of Earth system interactions and climate feedback remains poorly examined in the current generation of Earth system models. Increasing availability of ocean biogeochemical data, as well as an improved understanding of the underlying processes, allows advances in the marine biogeochemical components of the current generation of ESMs. The present study scrutinizes the extent to which marine biogeochemistry components of ESMs have progressed between the 5th and the 6th phases of the Coupled Model Intercomparison Project (CMIP).

54 ENVIRONMENTAL SCIENCES↗

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗

Online evolutionary neural architecture search for multivariate non-stationary time series forecasting

Time series forecasting (TSF) is one of the most important tasks in data science. TSF models are usually pre-trained with historical data and then applied on future unseen datapoints. However, real-world time series data is usually non-stationary and models trained offline usually face problems from data drift. Models trained and designed in an offline fashion can not quickly adapt to changes quickly or be deployed in real-time. To address these issues, this work presents the Online NeuroEvolution-based Neural Architecture Search (ONE-NAS) algorithm, which is a novel neural architecture search method capable of automatically designing and dynamically training recurrent neural networks (RNNs) for online forecasting tasks. Without any pre-training, ONE-NAS utilizes populations of RNNs that are continuously updated with new network structures and weights in response to new multivariate input data. ONE-NAS is tested on real-world, large-scale multivariate wind turbine data as well as the univariate Dow Jones Industrial Average (DJIA) dataset. These results demonstrate that ONE-NAS outperforms traditional statistical time series forecasting methods, including online linear regression, fixed long short-term memory (LSTM) and gated recurrent unit (GRU) models trained online, as well as state-of-the-art, online ARIMA strategies. Additionally, results show that utilizing multiple populations of RNNs which are periodically repopulated provide significant performance improvements, allowing this online neural network architecture design and training to be successful.

97 MATHEMATICS AND COMPUTING↗

Prototype energy models for data centers

Data centers in the United States consume about two percent of the nation’s electricity. Because heat gains from IT equipment drive cooling demand, data centers offer unique opportunities for energy savings. However, no prototype energy model for data centers is available in the suite of existing U.S. Department of Energy’s Commercial Prototype Building Models. Here we present the development of two new data center prototype models and their implementation in OpenStudio and EnergyPlus. The small-size data center model represents a computer room in a building served by computer room air conditioners (CRACs); while the large-sized model represents stand-alone data centers served by computer room air handlers (CRAHs) with a central chiller plant. For each data center model, two levels of IT equipment (ITE) load density were considered, to cover the wide range of IT power density of data centers: 40 and 100 W/ft 2 (430 and 1076 W/m 2 ) for the computer room, and 100 and 500 W/ft 2 (1076 and 5382 W/m 2 ) for the stand-alone data center. All other assumptions, such as building envelope, lighting, HVAC efficiencies and schedules, were based on the minimal requirements of ASHRAE Standard 90.1 at various vintages. We introduced a novel concept of supply and return air approach temperatures to capture the essential effects of non-uniform airflow and temperature distribution in data centers. The approach temperatures were pre-computed by computational fluid dynamics (CFD) simulations for various configurations of ITE loads and airflow containment management in data centers. A new feature was developed in EnergyPlus to implement the approach temperature method. A case study was conducted to demonstrate the use of the data center models. The two data center models cover all U.S. climate zones and can be used to evaluate energy saving measures for data centers, as well as to support development of data center energy efficiency codes and standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Isotopic signatures as tools to reconstruct the primordial architecture of the Solar System

The significant reorganization of the early Solar System due to giant planet migration has hampered our understanding of where planetary bodies formed. Previously employed proxies for reconstructing the primordial planetary architecture, such as water content or oxidation state, are complicated by post-accretionary processes. Furthermore we investigate basaltic achondrites for their nucleosynthetic isotope signatures in the elements neodymium (Nd) and zirconium (Zr) and show that they are—similar to previously investigated chondritic meteorites—characterized by a relative deficit in isotopes produced by the s-process of nucleosynthesis. Importantly, these data are well correlated with nucleosynthetic signatures observed in other elements, demonstrating that s-process matter was heterogeneously distributed throughout the early Solar System. By comparing these isotopic signatures with potential proxies for Solar System reconstruction and computer modeling, we here argue that this isotopic heterogeneity in bulk meteoritic materials is linked to the original heliocentric distance of formation. Such scaling of nucleosynthetic signatures with heliocentric distance could permit reconstruction of the primordial architecture of the Solar System by ‘cosmolocating’ the accretion orbits of meteoritic parent bodies as a function of incorporated s-process matter.

58 GEOSCIENCES↗

Magnetothermal properties of Tm x Dy 1- x Al 2 ( x = 0.25, 0.50 and 0.75)

Here, we describe magnetic, thermal, and magnetocaloric properties of rare earth intermetallic compounds Tm x Dy 1-x Al 2 with 0.25, 0.5 and 0.75. Using model Hamiltonian we consider contributions of the crystalline electric field anisotropy in both Tm and Dy magnetic sublattices, disorder in exchange interactions among Tm-Tm, Dy-Dy and Tm-Dy magnetic ions, and the Zeeman effect. Employing earlier reported and new experimental measurements, we first determine a single free variable – the intersublattice magnetic exchange parameter – to properly model the temperature and magnetic field dependencies of heat capacity and magnetization, and then use the modeling results to explain the emergence of an anomalous spin reorientation transition and its influence on the magnetocaloric effect in the title compounds. Theoretical results agree with experimental data reasonably well.

36 MATERIALS SCIENCE↗

Production of polyacrylonitrile/ionic covalent organic framework hybrid nanofibers for effective removal of chromium(VI) from water

Hexavalent Cr(VI) found in industrial wastewater is a proven carcinogen which causes serious health issues in humans around the world. This study presents a novel method to enhance the Cr(VI) oxyanion removal from wastewater by polyacrylonitrile (PAN) nanofibers through incorporation of a guanidinium-based ionic covalent organic framework (BT-DG) in the nanofibers structure. In this work, simple electrospinning technique was employed to produce PAN nanofibers and BT-DG was synthesized through condensation between benzene-1,3,5-tricarbaldehyde and N,N'-diaminoguanidine monohydrochloride. In-situ polymerization of BT-DG onto PAN nanofibers resulted in generation of hybrid PAN-BT-DG nanofibers. This modified PAN-BT-DG was characterized by obtaining its point of zero charge (PZC), differential scanning calorimeter (DSC), scanning electron microscopy (SEM) morphology and surface elements and oxidation states by X-ray photoelectron spectroscopy (XPS). PAN-BT-DG exhibited positive surface charge below pH 4, making it an outstanding adsorbent, for Cr(VI) removal. Cr(VI) adsorption onto PAN-BT-DG followed pseudo second order kinetics and adsorption data fitted well to Freundlich isotherm model. Highest Cr(VI) removal was obtained at 55 °C with a maximum Langmuir adsorption capacity of 173 mg/g at pH 3. Kinetic studies revealed that Cr(VI) adsorption onto PAN-BT-DG is endothermic and thermodynamically feasible. Desorption studies were conducted on PAN-BT-DG using 1 M NaOH as the stripping solvent and PAN-BT-DG exhibited excellent regeneration after five consecutive cycles.

36 MATERIALS SCIENCE↗