Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative capture”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermonuclear reaction rate of 29 Si(p,γ) 30 P

The thermonuclear rate of the 29 Si(p,γ) 30 P reaction impacts the 29 Si abundance in classical novae. A reliable reaction rate is essential for testing the nova paternity of presolar stardust grains. At present, the fact that no classical nova grains have been unambiguously identified in primitive meteorites among thousands of grains studied is puzzling, considering that classical novae are expected to be prolific producers of dust grains. We investigated the 29 Si + p reaction at center-of-mass energies of 200–420 keV, and present improved values for resonance energies, level excitation energies, resonance strengths, and branching ratios. One new resonance was found at a center-of-mass energy of 303 keV. For an expected resonance at 215 keV, an experimental upper limit could be determined for the strength. We evaluated the level structure near the proton threshold, and present new reaction rates based on all the available experimental information. Our new reaction rates have much reduced uncertainties compared with previous results at temperatures of T ≥ 140 MK, which are most important for classical nova nucleosynthesis. Furthermore, future experiments to improve the reaction rates at lower temperatures are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Revised decay properties of the key 93-keV resonance in the Mg 25 ( p , γ ) reaction and its influence on the MgAl cycle in astrophysical environments

The γ-decay properties of an excited state in 26 Al at 6398.3(8) keV have been reexamined using the 11 B+ 16 O fusion-evaporation reaction. This level represents a key 93.1(8)-keV resonance in the 25 Mg+p system and its relative branching to the 26 Al ground state, f 0 , has been determined to be 0.76 ± 0.03 (stat.) ± 0.10 (syst.). Furthermore, this is a significantly higher value than the most recent evaluation and implies a considerable increase in the production of cosmic γ rays from 26 Al radioactivity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F

The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.

6 ≤ A ≤ 19↗
DOE OSTI2023

Measuring the cross section of the N 15 ( α , γ ) 19 F reaction using a single-fluid bubble chamber

15 N(α, γ) 19 F is believed to be the primary means of stellar nucleosynthesis of fluorine. Here we present the use of a single-fluid bubble chamber to measure the cross section of the time-inverse photo-dissociation reaction. The method benefits from a luminosity increase of several orders of magnitude due to the use of a thicker liquid target - compared to thin films or gas targets - and from the reciprocity theorem. We discuss the results of experiments at the Thomas Jefferson National Accelerator Facility, where the cross section of the photodisintegration process 19 F(γ, α) 15 N was measured by bombarding a superheated fluid of C 3 F 8 with bremsstrahlung γ rays produced by impinging a 4 - 5.5 MeV electron beam on a Cu radiator. From the photodissociation yield the cross section was extracted by performing a convolution with a Monte Carlo–generated γ-ray beam spectrum. The measurement produced a cross section that was then time inverted using the reciprocity theorem. The cross section for the 15 N(α, γ) 19 F reaction was determined down to a value in the range of hundreds of picobarns. We report that with further improvements of the experimental setup the technique could potentially push cross section measurements down to the single picobarn range.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Measurement of the 25 Al(𝑑,𝑛)⁢ 26 Si reaction and impact on the 25 Al⁢(𝑝,𝛾)⁢ 26 Si reaction rate

The 25 Al(p,γ) 26 Si reaction plays a key role in nucleosynthesis pathways that influence the galactic abundance of 26 Al. A more precise determination of the proton strength of the lowest ℓ = 0 proton resonance in 26 Si is needed to improve reaction-rate calculations. This work measures the 25 Al(d,n) 26 Si proton-transfer reaction in inverse kinematics using a radioactive ion beam at RESOLUT, determining excitation energies and cross sections for the lowest ℓ = 0 resonance associated with the 3 + 3 state at 5.92(2) MeV. Coupled-reaction-channels calculations (FRESCO) are used to extract the ℓ = 0 spectroscopic factor and the corresponding proton width, yielding Γₚ = 2.19(45) eV and a (p,γ) resonance strength of 26(10) meV. This resonance is found to dominate the 25 Al(p,γ) 26 Si reaction rate above 0.2 GK.

20 ≤ A ≤ 38↗
DOE OSTI2024

Nuclear level density and γ -decay strength of Sr 93

This work presents the first experimentally determined nuclear level density and γ-ray strength function of the short-lived fission product 93 Sr, accomplished using the β-Oslo method. Direct measurement of the 92 Sr(n, γ) 93 Sr cross section is not currently possible, as the half-life of 2.66 hours is too short; instead, 93 Sr was formed through β decay of 93 Rb to excitation energies around the neutron separation energy. The γ-ray spectra were measured using a total absorption spectrometer at the National Superconducting Cyclotron Laboratory (NSCL) at Michigan State University (MSU). The statistical properties of the 93 Sr nucleus were experimentally determined, including the γ-ray strength function and nuclear level density. At low energies, the γ-ray strength function exhibits a constant γ-decay strength, rather than a slightly increasing strength with decreasing γ-ray energy as had been previously observed for several nuclei in this mid-mass region. Finally, these statistical properties were then implemented in the reaction code TALYS1.95 to calculate the 92 Sr(n, γ) 93 Sr cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Impact of the Li 6 asymptotic normalization constant onto α -induced reactions of astrophysical interest

Indirect methods have become the predominant approach in experimental nuclear astrophysics for studying several low-energy nuclear reactions occurring in stars, as direct measurements of many of these relevant reactions are rendered infeasible due to their low reaction probability. Such indirect methods, however, require theoretical input that in turn can have significant poorly quantified uncertainties, which can then be propagated to the reaction rates and have a large effect on our quantitative understanding of stellar evolution and nucleosynthesis processes. Here we present two such examples involving α-induced reactions, 13 C (α,n)⁢ 16 O and 12 C (α,γ)⁢ 16 O, for which the low-energy cross sections have been constrained with ( 6 Li,d) transfer data. In this Letter, we discuss how a first-principle calculation of 6 Li leads to a 21% reduction of the 12 C⁡(α,γ) ⁢ 16 O cross sections with respect to a previous estimation. This calculation further resolves the discrepancy between recent measurements of the 13 C (α,n)⁢ 16 O reaction and points to the need for improved theoretical formulations of nuclear reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Strength measurement of the $E^{lab}_α = 830$ keV resonance in the 22 Ne (α,n) 25 Mg reaction using a stilbene detector

The interplay between the 22 Ne (α,γ) 26 Mg reaction and the competing 22 Ne ⁢(α,n) 25 Mg reaction determines the efficiency of the latter as a neutron source at the temperatures of stellar helium burning. In both cases, the rates are dominated by the α-cluster resonance at 830 keV. This resonance plays a particularly important role in determining the strength of the neutron flux for both the weak and main s process as well as the n process. Recent experimental studies based on transfer reactions suggest that the neutron and γ-ray strengths for this resonance are approximately equal. In this study, the 22 Ne (α,n) 25 Mg resonance strength has been remeasured and found to be similar to the previous direct studies. Finally, this reinforces an 830 keV resonance strength that is approximately a factor of 3 larger for the 22 Ne ⁢(α,n) 25 Mg reaction than for the 22 Ne (α,γ) 26 Mg reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Search for Nova Presolar Grains: $γ$-Ray Spectroscopy of 34 Ar and its Relevance for the Astrophysical 33 Cl($p,γ$) Reaction

The discovery of presolar grains in primitive meteorites has launched a new era of research in the study of stellar nucleosynthesis. However, the accurate classification of presolar grains as being of specific stellar origins is particularly challenging. Recently, it has been suggested that sulfur isotopic abundances may hold the key to definitively identifying presolar grains with being of nova origins and, in this regard, the astrophysical Cl 33 ( p , γ ) Ar 34 reaction is expected to play a decisive role. As such, we have performed a detailed γ -ray spectroscopy study of Ar 34 . Excitation energies have been measured with high precision and spin-parity assignments for resonant states, located above the proton threshold in Ar 34 , have been made for the first time. Uncertainties in the Cl 33 ( p , γ ) reaction have been dramatically reduced and the results indicate that a newly identified ℓ = 0 resonance at E r = 396.9 ( 13 ) keV dominates the entire rate for T = 0.25 – 0.40 GK . Furthermore, nova hydrodynamic simulations based on the present work indicate an ejected S 32 / S 33 abundance ratio distinctive from type-II supernovae and potentially compatible with recent measurements of a presolar grain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Development of a pattern recognition algorithm for reconstructing multi-particle events in the Mu2e detector

Mu2e is an upcoming experiment at Fermilab and its main goal is to search for the Charged Lepton Flavor Violation (CLFV) in the coherent transition of a muon into an electron on an Al target. In Mu2e, multi-particle events can occur simultaneously within the same time region and it is crucial to accurately identify each particle track, including signals, to improve the robustness of track finding methods and enhance reconstruction efficiency. ¯p annihilation is one of the background events and produces multiple particles that can mimic signal events. Additionally, photons from radiative pion captures can produce a γ → $e+e−$ pair, which can be used to calibrate the Mu2e momentum scale and the resolution. The Mu2e track reconstruction sequence begins by grouping hits produced in the tracker based on time and z coordinate information, called TimeCluster, and selected hits are processed to reconstruct helices and determine their momentum. The current pattern recognition algorithms identify a single helix per TimeCluster for single track events. A new pattern recognition algorithm is being developed to reconstruct multi-particle events and its features for finding multiple tracks and the current evaluation results are reported.

Kitagawa, H. [Pisa U.]↗
DOE OSTI2020

Search for $\mu^- \rightarrow e^+$ conversion: what can be learned from the SINDRUM-II positron data on a gold target

In their 2006 paper setting the current limit on \mumemconv\ conversion search on a gold target \cite{sindrum_ii:Bertl2006}, the SINDRUM-II collaboration published, along with the electron momentum distribution, the momentum distribution of reconstructed positrons. Near the positron spectrum endpoint, there is a statistically significant excess of observed events over the expected background. We estimate that in the region 88 MeV/c < p < 95 MeV/c there are 13 events with an expected background of about 1-1.5 event, which has not been discussed by the authors. Those 13 events form a bump with a width consistent with the experimental resolution, making one think of a $\mu^- \rightarrow e^+$ conversion signal. However, the reconstructed position of the bump is about 1 MeV/c, or $\sim\ 4\sigma_p$, lower than the expected $\mu^- Au \rightarrow e^+ Ir $ signal, which strongly discourages the exotic interpretation. The excess, however, could be due to an exclusive dipole radiative muon capture (RMC) transition $^{197}Au(GS) ~\rightarrow~ ^{197}Pt(GS)$ with the branching fraction of about $2\cdot 10^{-4}$. Such a transition would not be resolved by the existing RMC measurements. We conclude that the exclusive RMC transitions could significantly modify the positron spectrum near the kinematic endpoint, to fully exploit the physics potential of the upcoming experiments such as Mu2e and COMET, a better theoretical understanding of the RMC spectrum on nuclei is needed. A high-resolution measurement of the RMC photon spectra has to be carried out. Without that, the sensitivity of the searches for $\mu^- \rightarrow e^+$ and $\mu^- \rightarrow e^-$ might be severely limited by unknown probabilities of RMC transitions to the exclusive low-lying states of the daughter nuclei.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Probing the limits of statistical neutron capture for the r process: Experimental constraints on 141 Cs nuclear level densities

The r-process abundance peaks, particularly near mass number A ∼ 130, reflect underlying nuclear structure effects such as closed neutron shells, yet modeling the nucleosynthesis in this region remains hindered by uncertain neutron-capture rates. These rates are especially sensitive to nuclear level densities (NLDs) and γ-ray strength functions of neutron-rich nuclei, where experimental data are scarce. We present the first experimental constraint on the NLD of 141 Cs using the β-Oslo method, extending sensitivity to the neutron-rich regime near the N = 82 closed shell. Our data allow for critical calibration of microscopic NLD models and reveal that 141 Cs lies near the limit of statistical model applicability. Using this experimental input, we evaluate radiative neutron-capture rates across neighboring isotones using both Hauser–Feshbach (HF) and High Fidelity Resonance (HFR) models. Our results show order-of-magnitude rate increases for nuclei along the N = 86 line, signaling a transition to resonance-dominated capture in this region. These findings underscore the importance of constraining NLDs to improve r-process reaction network predictions, particularly in environments where the validity of statistical models breaks down.

Nuclear level density↗
DOE OSTI2022

Sensitivity of Cloud-Radiative Effects to Cloud Fraction Parametrizations in Tropical, Mid-Latitude and Arctic Kilometre-Scale Simulations

The Regional Atmosphere (RA) configuration of the Met Office Unified Model currently requires different cloud fraction parametrizations (CFP) for tropical and midlatitude simulations. To explore the scope for unification of these two RA configurations, this paper presents a detailed evaluation of simulations over tropical, midlatitude and arctic domains, with two different diagnostic CFPs, a prognostic CFP, and no CFP at all. Furthermore, a novel, hybrid approach was used that treats liquid cloud diagnostically and ice cloud prognostically. Using observations from three U.S. Department of Energy Atmospheric Radiation Measurement supersites, it is shown that none of these CFPs stands out as superior over all domains. Over the frequently overcast Arctic, the all-or-nothing approach best captures the cloud radiative properties. Conversely, CFPs are of benefit in regions with frequent partial cloudiness, such as the midlatitudes and the tropics. However, their improved cloud radiative properties often hide an error compensation. All models underestimate overcast, low-base cloud with small water paths in convective environments. In addition, mid-latitude overcast, low-base, optically thick clouds in the morning, possibly associated with overnight convection, are frequently too broken. Diagnostic schemes compensate for these errors by producing spurious, scattered afternoon cloud, which could be due to a correct cloud response to too eager convective initiation. Winter clouds over the midlatitudes are improved when liquid cloud is represented diagnostically with a bimodal saturation-departure PDF, without error compensation. While it is difficult to unify the RA across the globe around a single CFP scheme, the newly proposed Hybrid scheme performs reasonably well for cloud cover across all regions. It also exhibits SW biases that are smaller than most other configurations and is less affected by excessive liquid water paths and compensating errors than fully diagnostic schemes. Surface precipitation is fairly insensitive to the CFP in the simulations shown here.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Review of the Resuspension of Radioactively Contaminated Particles by Vehicle and Pedestrian Traffic—Current Theory, Practice, Gaps, and Needs

Here, the resuspension of radioactively contaminated particles in a built environment, such as from urban surfaces like foliage, building exteriors, and roadways, is described empirically by current plume and dosimetry models used for hazard assessment and long-term risk purposes. When applying these models to radiological contamination emergencies affecting urban areas, the accuracy of the results for recent contamination deposition is impacted in two main ways. First, the data supporting the underlying resuspension equations was acquired for open, quiescent conditions with no vehicle traffic or human activities, so it is not necessarily representative of the urban environment. Second, mechanical disturbance by winds in urban canyons and during emergency operations caused by vehicle traffic and human activities are not directly considered by the equations. Accordingly, plume and dosimetry models allow the user to input certain compensating values, but the models do not necessarily supply users instructions on what values to use. This manuscript reviews the available literature to comprehensively and consistently pool data for resuspension due to mechanically induced resuspension applicable to urban contamination. Because there are few studies that directly measured radioactive resuspension due to vehicles and pedestrians, this review novelly draws on a range of other studies involving non-radioactive particles, ranging from outdoor air pollution emissions to indoor allergen transport. The results lead to tabulated, recommended values for specific conditions in the emergency phase to help users of plume and dosimetry models maintain the conservativeness needed to properly capture the potential radiation dose posed by mechanically induced resuspension. These values are of benefit to model users until better data are available. The results also suggest the types of data that may result in improved plume and dose modeling.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

Trap-assisted Auger-Meitner recombination in GaN p-i-n diodes

Most properties of semiconductor devices are dominated by shallow impurities. However, deep defects often play an important role, for instance, in recombination processes or high field transport. While a variety of techniques are available to assess the density and energy levels of impurities, other properties, such as the recombination mechanisms of the defects, escape observation. We report on the direct measurement of hot electrons generated by trap-assisted Auger-Meitner recombination (TAAR) in GaN p-i-n diodes. By performing electron emission spectroscopy (EES) on diodes with surfaces activated to negative electron affinity by cesium, we observe the expected overflow electrons of p-i-n diodes under low current injection. However, when operating the devices at higher current densities, as low as ∼25 A/c⁢m 2 , we measure the emission of high-energy electrons. At variance with the observed hot electrons in light-emitting diodes (LEDs) using EES, the hot electrons generated in p-i-n diodes at our tested currents cannot be from eeh Auger-Meitner recombination due to the diodes' significantly lower carrier densities compared to those in LEDs. During our measurements, we observe the emission of accumulated electrons with energies ∼0.42 eV, ∼0.99 eV, ∼1.43 eV, and ∼2.32 eV above the conduction-band minimum (CBM) at various bias conditions, suggesting the existence of conduction-band features in GaN at these energies where electrons can be long-lived, such as satellite-valley minima and inflection points. We also measure incompletely relaxed hot electrons approaching energies 1.97 ± 0.13 eV and 2.94 ± 0.13 eV above the CBM, as the diodes are biased to high currents, suggesting at least some of the TAAR partaking defects have an energy level ≳1.97 eV and ≳2.94 eV away from either the conduction or valence band edges. Additionally, at our highest operating currents, we measure hot electrons with energies 3.28 ± 0.13 eV above the CBM, providing direct evidence of TAAR processes involving shallow impurities. This unexpected observation of TAAR in GaN p-i-n diodes spotlights the importance of further studies of defects in GaN and the necessity to incorporate the multi-phonon emission, radiative, and TAAR capture steps of defect-assisted recombination cycles into device modeling. Furthermore, this experiment demonstrates the applicability of the simplest semiconductor structures, p-i-n diodes, as a test bed to study the rich recombination physics of semiconductor materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A Simplified Numerical Approach to Characterize the Thermal Response of a Moving Bed Solar Reactor

Abstract Concentrated solar thermochemical storage in the form of a zero-emission fuel is a promising option to produce long-duration energy storage. Solar fuel is produced using a cavity reactor that captures concentrated solar radiation from a solar field of heliostats. In this paper, heat transfer model of a tubular plug-flow reactor designed and manufactured for a solar fuel production is presented. Experimental data collected from a fixed bed tubular reactor testing are used for model comparison. The system consists of an externally heated tube with counter-current flowing gas and moving solid particles as the heated media. The proposed model simulates the dynamic behavior of temperature profiles of the tube wall, gas, and particles under various gas flowrates and residence times. The heat transfer between gas–wall, solid particle–wall, and gas–solid particle is numerically studied. The model results are compared with the results of experiments done using a 4 kW furnace with a 150 mm heating zone surrounding a horizontal alumina tube (reactor) with 50.8 mm outer diameter and thickness of 3.175 mm. Solid fixed particles of magnesium manganese oxide (MgMn2O4) with the size of 1 mm are packed within the length of 250 mm at the center of the tube length. Simulation results are assessed with respect to fixed bed experimental data for four different gas flowrates, namely, 5, 10, 15, and 20 standard liters per minute of air, and furnace temperatures in the range of 200–1200 °C. The simulation results showed good agreement with maximum steady state error that is less than 6% of those obtained from the experiments for all runs. The proposed model can be implemented as a low-order physical model for the control of temperature inside plug-flow reactors for thermochemical energy storage applications.

Engineering↗
DOE OSTI2022

A Forward Feedback Control Scheme for a Solar Thermochemical Moving Bed Counter-Current Flow Reactor

Abstract Pelletized thermochemical energy storage media has a potential for long-duration energy storage. Production of solid-state energy storage media can be done within a cavity chemical reactor that captures concentrated solar radiation from a solar thermal field. The temperature stability of a solar reactor is directly influenced by the solar flux intercepted. This paper presents a low-order physical model to simulate the dynamic response of temperature inside a tubular plug-flow reactor prototype. Solid granular particles are fed to the reactor from the top whereas a counter-current flowing gas enters the reactor from the bottom. An in-house code was developed to model transient heat transfer of the reactor wall, gas, and moving particles. The model was preliminarily validated with packed beds for different temperature ranges and two gas flowrates. Dynamic response of the reactor temperature is simulated for different input power and gas/particle flowrates. The results show that the system response can be controlled efficiently by utilizing input power (solar flux) as a control parameter. A conventional proportional integral (PI) controller is designed to control the temperature inside the reactor and to maintain it during the solar flux intermittency. The controller parameters are tuned using the Ziegler–Nichols method to ensure optimal system response. The results show that the feedback control model is successful in tracking different reference reactor temperatures within a reasonable settling time of 30 min and eliminated overshoot. This study can be extended to include a hybrid reactor with a multi-input, multi-output variable system.

Energy & Fuels↗