SEARCH · Engineering Papers

Results for “yield point”

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

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Bottlebrush Block Copolymer Thin Films

The self-assembly of block copolymers (BCPs) in thin films is governed by interfacial interactions (enthalpy) between the blocks and interfaces (at both the substrate and the surface). Advantage can also be taken of entropy to control the polymer orientation. Here, we synthesized a series of bottlebrush block copolymers (BBCPs), with deuterated polystyrene (DPS) as the core block and poly(solketal acrylate) (PSA) as the corona block, where the backbone length (N BB) and grafting densities (GDs) were varied. The hydrophobic PSA block was converted to a hydrophilic poly(glyceryl acrylate) (PGA) block by solid-state hydrolysis, bringing the BBCP from a disordered state into an ordered state with a lamellar microdomain morphology. The orientation of the morphology as a function of distance from the interfaces was systematically investigated by atomic force microscopy (AFM), interference microscopy, grazing-incidence small-angle X-ray scattering (GISAXS), and grazing-incidence small angle neutron scattering (GISANS). For N BB = 1 (diblock BCP), a mixed lamellar orientation was found. For N BB = 5 (star-like BBCP), a vertical orientation of the lamellar microdomains was found dominantly at the interfaces that propagated into the bulk of the film due to the entropic penalty associated with constraints on the junction points. Increasing the N BB to 50 (rod-like BBCP) and 100 (worm-like BBCP) yielded a random orientation, as the side chains became overcrowded when aligned in the same direction. By reducing the grafting density (GD) of a N BB = 50 BBCP from 100% to 75% and 50% by the incorporation of spacers into the backbone, the steric hindrance of side chains is reduced markedly, facilitating a vertical orientation of the lamellar microdomains. However, the vertical orientation becomes more random within the bulk of the thin film for different BBCP architectures. By modifying the substrate with a layer of PS, the preferential interactions between the core block of the BBCP and the interfaces induce a reconfiguration of the BBCP, drawing the core block out to contact substrate directly driven by enthalpy. Under these conditions, N BB = 5 showed a predominantly parallel orientation of the lamellar microdomains with the corona block minimizing contact with the substrate. Higher N BB values of 50 and 100 showed a more random orientation of the lamellar microdomains at different depths on the PS-modified substrate. Overall, thin-film morphology was regulated by tuning the configuration of the BBCP, varying the entropic contribution to the microdomain orientation.

Chen, Zhan↗

DOE OSTI2021

Extracting correlation effects from momentum-resolved electron energy loss spectroscopy: Synergistic origin of the dispersion kink in Bi 2.1 Sr 1.9 CaCu 2 O 8 + x

We employ momentum-resolved electron energy loss spectroscopy (M-EELS) on Bi 2.1 Sr 1.9 CaCu 2 O 8 + x to resolve the issue of the kink feature in the electron dispersion widely observed in the cuprates. To this end, we utilize the GW approximation to relate the density response function measured in M-EELS to the self-energy, isolating contributions from phonons, electrons, and the momentum dependence of the effective interaction to the decay rates. The phononic contributions, present in the M-EELS spectra due to electron-phonon coupling, lead to kink features in the corresponding single-particle spectra at energies between 40 and 80 meV, independent of the doping level. We find that a repulsive interaction constant in momentum space is able to yield the kink attributed to phonons in ARPES. Hence, our analysis of the M-EELS spectra points to local repulsive interactions as a factor that enhances the spectroscopic signatures of electron-phonon coupling in cuprates. We conclude that the strength of the kink feature in cuprates is determined by the combined action of electron-phonon coupling and electron-electron interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Projective Integral Updates for High-Dimensional Variational Inference

Variational inference is an approximation framework for Bayesian inference that seeks to improve quantified uncertainty in predictions by optimizing a simplified distribution over parameters to stand in for the full posterior. Capturing model variations that remain consistent with training data enables more robust predictions by reducing parameter sensitivity. This work introduces a fixed-point optimization for variational inference that is applicable when every feasible log density can be expressed as a linear combination of functions from a given basis. In such cases, the optimizer becomes a fixed-point of projective integral updates. When the basis spans univariate quadratics in each parameter, the feasible distributions are Gaussian mean-fields and the projective integral updates yield quasi-Newton variational Bayes (QNVB). Other bases and updates are also possible. Since these updates require high-dimensional integration, this work begins by proposing an efficient quasirandom sequence of quadratures for mean-field distributions. Each iterate of the sequence contains two evaluation points that combine to correctly integrate all univariate quadratic functions and, if the mean-field factors are symmetric, all univariate cubics. More importantly, averaging results over short subsequences achieves periodic exactness on a much larger space of multivariate polynomials of quadratic total degree. The corresponding variational updates require four loss evaluations with standard (not second-order) backpropagation to eliminate error terms from over half of all multivariate quadratic basis functions. Furthermore, this integration technique is motivated by first proposing stochastic blocked mean-field quadratures, which may be useful in other contexts. A PyTorch implementation of QNVB allows for better control over model uncertainty during training than competing methods. Experiments demonstrate superior generalizability for multiple learning problems and architectures.

Gaussian mean-field↗

DOE OSTI2022

The oxygen stable isotope composition of CRM 125-A UO 2 standard reference material

While there is a clear need for standardized reference materials for analytical calibrations and for inter-laboratory comparisons, there are not currently any for the oxygen stable isotopic composition of uranium oxides. In this paper we summarize the results from four laboratories by seven different methods of oxygen stable isotope analyses using fluorination techniques of CRM 125-A UO 2 Standard Reference Material. We synthesize these data and methods to arrive at a consensus oxygen stable isotope composition for CRM 125-A $δ$ 18 O = -9.63‰ (±0.29‰) VSMOW. We discuss methodological differences between analytical approaches, including furnace vs laser heating, fluorination using BrF 5 or ClF 3 , as well as calibration strategies. We highlight the potential effects of calibration scale compression from single-point calibrations using reference material with $δ$ 18 O values having a large relative difference from the sample being analyzed. We demonstrate how calibration scale compression can yield differences in calibrated $δ$ 18 O values up to ~2‰ for samples with ~20‰ difference from a single reference material, if the calibration slope of different analytical systems differs by 0.1. In conclusion, we suggest the use of liquid water calibration standards sealed in silver capillary tubes for multi-point calibrations of fluorination analysis systems.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2023

Engineering the formation of spin-defects from first principles

The full realization of spin qubits for quantum technologies relies on the ability to control and design the formation processes of spin defects in semiconductors and insulators. We present a computational protocol to investigate the synthesis of point-defects at the atomistic level, and we apply it to the study of a promising spin-qubit in silicon carbide, the divacancy (VV). Our strategy combines electronic structure calculations based on density functional theory and enhanced sampling techniques coupled with first principles molecular dynamics. We predict the optimal annealing temperatures for the formation of VVs at high temperature and show how to engineer the Fermi level of the material to optimize the defect’s yield for several polytypes of silicon carbide. Our results are in excellent agreement with available experimental data and provide novel atomistic insights into point defect formation and annihilation processes as a function of temperature.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Design optimization of integrated cooling inserts in modular Fischer-Tropsch reactors

Sustainable production of liquid fuels and feedstocks from atmospheric CO 2 through carbon recycling technologies is necessary to broaden decarbonization efforts and further reduce global emissions. Fischer–Tropsch (FT) reactors can address this challenge by providing storable, high-value liquid hydrocarbon fuels and feedstocks from syngas generated through reductive CO 2 utilization. FT technology, however, is most cost effective at large-scale, fixed-site plants and does not effectively address the smaller, globally distributed CO 2 point sources. Alternatively, smaller, modular reactors can be composed together to the specific scale of the emission source, yielding a flexible solution that enables more widespread deployment. In these modular reactors, thermal management using a finned cooling insert embedded within the catalyst matrix is critical to maintaining the performance and viability. Thus, in this work, topology optimization is used to determine the cooling insert geometry that maximizes reactor productivity while preventing auto-thermal runaway. Optimal designs are generated for varying number of constituent fins and over a range of maximum operating temperatures. Constraints including minimum feature length-scales and prescribed cooling insert material are imposed on the designs to enhance manufacturability. The impact of design features such as insert tapering and increased length scale hierarchy is automatically revealed by the systematic design framework employed. Here, this approach thus generates novel optimal geometries while automating, accelerating, and enhancing the design process compared to traditional heuristic approaches for cooling insert design in modular FT reactors.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Pairwise connected tensor network representation of path integrals

It has been recently shown how the tensorial nature of real-time path integrals (PIs) involving the Feynman-Vernon influence functional can be utilized with matrix product states, taking advantage of the finite length of the bath-induced memory. Tensor networks (TNs) promise to provide a unified language to express the structure of a PI. A generalized TN specifically incorporating the pairwise interaction structure of the influence functional and its invariance with respect to the average forward-backward position or the sojourn value in the form of the blip representation is derived and implemented. This pairwise connected TNPI (PC-TNPI) is illustrated through applications to typical spin-boson problems and explorations of the differences caused by the exact form of the spectral density. The storage and performance scalings are reported, showing the compactness of the representation and the efficiency of the contraction process. Finally, taking advantage of the compressed representation, the viability of using PC-TNPI for simulating multistate problems is demonstrated. The PC-TNPI structure can be shown to yield other TN algorithms currently in use. Consequently, it should be possible to use it as a starting point for deriving other optimized procedures.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Brief: Daily Performance Changes in Metal Halide Perovskite PV Modules

Metal halide perovskite photovoltaic modules deployed outdoors and held at their maximum power point show daily, reversible, relative changes of up to 30% in efficiency between morning and afternoon. Predicting energy yield and quantifying reliability will require properly handling such daily performance changes.

14 SOLAR ENERGY↗

DOE OSTI2020

Evidence for anti-synergism between ion-assisted etching and in-plasma photoassisted etching of silicon in a high-density chlorine plasma

Etching of p-Si in 60 mTorr 10%Cl2/90%Ar Faraday-shielded inductively coupled high density plasmas was investigated under both ionassisted etching (IAE) and photoassisted etching (PAE) conditions. Real-time etching rates and after-etching Si surface chemical compositions were obtained by laser interferometry and vacuum-transfer x-ray photoelectron spectroscopy (XPS), respectively. Precisely controlled ion energy distributions (IEDs) were generated by applying pulsed negative DC bias on the conductive sample stage. Above a 36 eV threshold at a total flow rate of 250 SCCM, the IAE rate increased with the square root of the ion energy. In contrast to the DC bias, etching under RF bias did not exhibit a threshold ion energy because of the wide IED. XPS spectra revealed that the surface layer under PAE conditions had a significantly lower chlorine content, composed of only SiCl. Under IAE conditions, however, silicon dangling bonds (Si•), SiCl 2 , and SiCl 3 were found on the surface, in addition to SiCl, with a relative abundance of SiCl > SiCl 2 , and SiCl 3 . The absence of higher chlorides and Si• under PAE conditions suggested that vacuum ultraviolet photons and above threshold-energy ions interact with the surface very differently. By varying the duty cycle of the pulsed DC bias, it was found that the IAE rate scaled with the energetic ion dose, but only for low duty cycles. For higher duty cycles, the apparent IAE yield fell off with an increasing Cl coverage on the surface, as the duty cycle went up, which pointed to a negative synergy (antisynergism) between PAE and IAE as the explanation. This antisynergism was further supported by the observed decrease of the total etching rate with an increasing period of the pulsed DC bias. A plausible mechanism is that increasing the pulsing period causes more near-surface damage, creating more recombination centers that lead to a higher loss rate of electron-hole pairs through recombination, thereby reducing the PAE rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

ATLAS flavour-tagging algorithms for the LHC Run 2 $pp$ collision dataset

The flavour-tagging algorithms developed by the ATLAS Collaboration and used to analyse its dataset of $\sqrt{s}$ = 13 TeV $pp$ collisions from Run 2 of the Large Hadron Collider are presented. These new tagging algorithms are based on recurrent and deep neural networks, and their performance is evaluated in simulated collision events. These developments yield considerable improvements over previous jet-flavour identification strategies. At the 77% $b$-jet identification efficiency operating point, light-jet (charm-jet) rejection factors of 170 (5) are achieved in a sample of simulated Standard Model $t\bar{t}$ events; similarly, at a $c$-jet identification efficiency of 30%, a light-jet ($b$-jet) rejection factor of 70 (9) is obtained.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Integrating CO2 Electrolysis with Gas Fermentation to Produce Valuable Fuels and Chemicals

Many industrial activities squander CO2, decreasing process yield. We envision a future where this waste carbon is instead captured, upgraded, and valorized directly at the point of emission. Within the CO2 Reduction and Upgrading Consortium (a collaboration of seven US national laboratories and industrial partners), we are pursuing this goal by developing and de-risking new technologies for low temperature CO2 electrolysis, coupled with biological upgrading of intermediates into more valuable compounds. One such process involves electrocatalytic reduction of CO2 to generate carbon monoxide (CO). As both a carbon and energy source, CO represents an attractive feedstock for microbial upgrading by certain syngas-fermenting species, such as the autotrophic bacterium Clostridium autoethanogenum. Our team has developed new genetic tools and optimized cultivation techniques to enhance C. autoethanogenum as a platform host for the biological conversion of syngas into value-added products. For example, we have created novel CRISPR-based genetic engineering techniques to build new, genome-reduced, platform strains of C. autoethanogenum with improved growth rates. Further, we ve introduced heterologous biochemical pathways into C. autoethanogenum to enable the production of high-value compounds from syngas, such as the isoprenoid precursor mevalonic acid. With the tools of electrochemistry and synthetic biology, there is virtually no limit to the spectrum of products that could be sustainably manufactured from CO2.

09 BIOMASS FUELS↗

DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗

DOE OSTI2021

Command and Control of Utilities Systems

The utilities systems that encompass the Y-12 National Security Complex are an integral part of the entire campus, and are essential to maintaining critical mission focus for the site. Although the utilities systems each differ in function, the central production and distribution model creates a common thread on which each system is run. Several incidents in FY 2018 and FY 2019 documented a lack of overall systems control and highlighted the need to create a better command system that could maintain system availability and decrease possible impacts to the surrounding environment. An initial system was started in September of 2019 to better control mechanical utility distribution systems. The measurement of the system can be modeled against a standard command and control related incident. The corresponding work hours from a typical control-related incident were used as key component to this system model. System effectiveness is modeled against utility’s incidents that demonstrate a loss of command and control. The increased simulations of the model also document control can only be effectively contained within the stated time (3-5 years) by yielding more engineers to complete alignments, or drastically reducing the amount of time for alignments. The total control points of the model document a clear correlation between alignments, control points, and the amount of man-hours per incident. Current data shows that actual incident man-hours (1.88 manhours/ incident since initial system implementation) are near expected incident man-hours (1.81 manhours/ incident) with command and control for October 2021, which documents a 33.3% reduction in man-hours/incident since initial system implementation. A conclusion can be drawn that focus on increasing the amount of controlled points in the system can reduce incidents and the decrease the amount of non-value added work to employee’s day-to-day lives, while also achieving a culture that enhances and sustains utilities reliability.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2022

3D simulations of inertial confinement fusion implosions part 2: systematic flow anomalies and impact of low modes on performances in OMEGA experiments

We present the first 3D radiation-hydrodynamic simulations of directly driven inertial confinement fusion implosions with an inline package for polarized crossed-beam energy transfer, which were used to assess the impact of the current distributed polarization rotators (DPRs) on OMEGA as well as other known sources of asymmetry. Applied to OMEGA implosions, the simulations predict bang times with no need for ad hoc multipliers, as well as yields—if you separately account for the impacts of imprint and fuel age. The magnitude of the flow is well reproduced when the low mode sources are large, whereas the modeling of the stalk is thought to be required to match the flow magnitude in the remaining cases. For the cases explored in more detail, polarized cross-beam energy transfer (CBET)—the only known systematic drive asymmetry, brought the results closest to the measured flow vectors. The remaining discrepancies are shown to be stemming from the limited knowledge of the laser pointing modes. For typical current levels of beam mispointing, power imbalance, target offset, and asymmetry caused by polarized CBET, low modes degrade the yield by more than 40%. The current strategy of attempting to compensate the mode-1 asymmetry with a preimposed target offset recovers only about one-third of the losses caused by the low modes due to the dynamic nature of the multiple asymmetries and the presence of low modes other than l = 1. Furthermore, addressing the root causes of the drive asymmetries is apt to be more beneficial. To that end, one possible solution to the specific issue of polarized CBET (10 µm DPRs) is shown to work well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Soft pions and transport near the chiral critical point

During the expansion of a heavy ion collision, the system passes close to the $\textit{O}$(4) critical point of QCD, and thus the fluctuations of the order parameter $(σ, \overrightarrow{π}$) are expected to be enhanced. In this work, our goal is to compute how these enhanced fluctuations modify the transport coefficients of QCD near the pseudocritical point. We also make a phenomenological estimate for how chiral fluctuations could effect the momentum spectrum of soft pions. We first formulate the appropriate stochastic hydrodynamic equations close to the $\textit{O}$(4) critical point. Then, working in mean field, we determine the correlation functions of the stress tensor and the currents which result from this stochastic real-time theory, and use these correlation functions to determine the scaling behavior of the transport coefficients. The hydrodynamic theory also describes the propagation of pion waves, fixing the scaling behavior of the dispersion curve of soft pions. We present scaling functions for the shear viscosity and the charge conductivities near the pseudocritical point and estimate the absolute magnitude of the critical fluctuations to these parameters and the bulk viscosity. Using the calculated pion dispersion curve, we estimate the expected critical enhancement of soft pion yields, and this estimate provides a plausible explanation for the excess seen in experiment relative to ordinary hydrodynamic computations. Our results motivate further phenomenological and numerical work on the implications of chiral symmetry on real-time properties of thermal QCD near the pseudocritical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Accuracy of predictions made by machine learned models for biocrude yields obtained from hydrothermal liquefaction of organic wastes

Hydrothermal liquefaction (HTL) has potential for converting abundant wet organic wastes into renewable fuels. Because HTL consists of a complex reaction network, deterministic, physics-based prediction of its biocrude yield is prohibitively difficult. Data-driven methods provide an alternative to the physics-based approach; however, rigorous testing must be performed to ensure the accuracy of predictions made by data-driven methods. To this end, a data set was assembled consisting of 570 data points appearing in the open literature. The data set was divided into training, validation, and test sub-sets and used for evaluating different machine learning regression approaches to predict biocrude yield. Among the tested algorithms, Random Forest and eXtreme Gradient Boosting (XGBoost) predicted biocrude yields in a test set that had not been used for training with the greatest accuracy, with root mean square errors (RMSE) of 8.34 and 8.57, respectively. Further refinement of the Random Forest model reduced its RMSE to 8.07. In comparison, predictions of a series of literature models resulted in RMSE ranging from 9.16 in the most accurate case to 27.6 in the least accurate; most literature models yielded RMSE values > 10. Using biocrude yield predictions from the most accurate Random Forest model and a probabilistic economic analysis found that the model accuracy is sufficient to prioritize allocation of resources based on projected minimum fuel selling price. In our report the models and analysis represent a major advance in the ability to use readily available data to predict biocrude yields on new feedstocks that have not previously been studied.

42 ENGINEERING↗

DOE OSTI2022

Predicting tungsten erosion and leakage properties for the new V-shaped small angle slot divertor in DIII-D

Impurity transport modeling of the new tungsten (W)-coated, V-shaped small angle slot (SAS) divertor in the DIII-D tokamak was conducted using the SOLPS-ITER plasma edge code package and the DIVIMP impurity tracking code. The inboard baffle of the current SAS divertor will be shifted closer to the outboard baffle, creating a V-corner at the slot vertex. In addition, the outboard baffle will be coated with 10–15 μm of W for experiments studying high-Z sourcing and leakage in a closed divertor. Modeling of the 'SAS-VW' divertor predicts that these changes to the inner baffle will reduce W gross erosion by 40× relative to the existing SAS divertor when the outer strike point (OSP) is at the V-corner and the ion B × ∇B drift is towards the divertor, driven primarily by significant cooling near the slot vertex. Most W erosion in SAS-VW is expected to occur near the slot entrance, which may pose a higher risk to core contamination than W eroded deeper in the slot. Furthermore, adding a new sheath-based prompt redeposition model increases the sensitivity of redeposition estimates to near-target plasma conditions and may provide more accurate predictions of net erosion. Moving the OSP outboard from the slot vertex ~4 cm onto the W-coated region yields a 40× increase in the gross erosion rate and a 50% decrease in the core leakage fraction. Thus slight variations in strike point location may counteract the potential benefits of the tightly-baffled V slot on minimizing erosion. This impurity transport modeling provides useful guidance for future experiments on the SAS-VW divertor focused on high-Z erosion/redeposition, scrape-off layer transport, and core leakage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Impact of duration and missing data on the long-term photovoltaic degradation rate estimation

Accurate quantification of photovoltaic (PV) system degradation rate (R D ) is essential for lifetime yield predictions. Although R D is a critical parameter, its estimation lacks a standardized methodology that can be applied on outdoor field data. The purpose of this paper is to investigate the impact of time period duration and missing data on R D by analyzing the performance of different techniques applied to synthetic PV system data at different linear R D patterns and known noise conditions. The analysis includes the application of different techniques to a 10-year synthetic dataset of a crystalline Silicon PV system, with emulated degradation levels and imputed missing data. Here, the analysis demonstrated that the accuracy of ordinary least squares (OLS), year-on-year (YOY), autoregressive integrated moving average (ARIMA) and robust principal component analysis (RPCA) techniques is affected by the evaluation duration with all techniques converging to lower R D deviations over the 10-year evaluation, apart from RPCA at high degradation levels. Moreover, the estimated R D is strongly affected by the amount of missing data. Filtering out the corrupted data yielded more accurate R D results for all techniques. It is proven that the application of a change-point detection stage is necessary and guidelines for accurate R D estimation are provided.

14 SOLAR ENERGY↗