Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Temperature sensitivity”

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 739 records · Page 41

A metal-ceramic coaxial cable with multipoint Fabry-Pérot interferometers for monitoring distributed high temperature

A metal ceramic coaxial cable (MCCC) sensor containing multiple Fabry-Pérot interferometers (FPIs) has been designed, fabricated, and evaluated for measuring distributed high temperatures. The MCCC sensor demonstrated here consists of stainless steel tube and wire conductors, alumina tube insulator, and three pairs of air gap reflectors distributed along the cable. The multipoint sensor operation is realized by the joint time-frequency domain technique, which combines the functions of time-domain reflectometry and frequency division multiplexing to enable the measurement of spatially distributed temperature. The 3-point MCCC-FPI sensor has been demonstrated for distributed temperature measurements between 250 and 550 °C and operability up to 1000 °C. The sensor offers high sensitivity of >0.18 MHz/°C with a measurement resolution of 1 Hz and reasonably small high temperature measurement error range of ±3 °C.

42 ENGINEERING↗

Hotspot formation due to shock-induced pore collapse in 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX): Role of pore shape and shock strength in collapse mechanism and temperature

The shock to detonation transition in heterogeneous high energy density solids starts with the spatial localization of mechanical energy into so-called hotspots that form due to the interaction between the leading wave and microstructural features and defects. We used large-scale molecular dynamics to characterize the hotspots resulting from the shock-induced collapse of cylindrical voids and elongated cracks focusing on the effect of shock strength, defect shape, and size. The temperature fields resulting from the collapse of cracks elongated along the shock direction show significantly higher sensitivity to both shock strength and size than cylindrical voids. Cracks 80 nm in length result in temperatures almost three times higher than voids 80 nm in diameter, reaching values corresponding to the ideal case of isentropic recompression of a gas. The molecular dynamics trajectories reveal the atomic origin of this contrasting behavior. While circular voids undergo a transition from viscoelastic pore collapse to a hydrodynamic regime with increasing shock strength, shock focusing in elongated cracks results in jetting and vaporization which, upon recompression, leads to increased heating.

36 MATERIALS SCIENCE↗

Wireless Passive Ceramic Sensor for Far-Field Temperature Measurement at High Temperatures

A passive wireless high-temperature sensor for far-field applications was developed for stable temperature sensing up to 1000 °C. The goal is to leverage the properties of electroceramic materials, including adequate electrical conductivity, high-temperature resilience, and chemical stability in harsh environments. Initial sensors were fabricated using Ag for operation to 600 °C to achieve a baseline understanding of temperature sensing principles using patch antenna designs. Fabrication then followed with higher temperature sensors made from (In, Sn) O2 (ITO) for evaluation up to 1000 °C. A patch antenna was modeled in ANSYS HFSS to operate in a high-frequency region (2.5–3.5 GHz) within a 50 × 50 mm2 confined geometric area using characteristic material properties. The sensor was fabricated on Al2O3 using screen printing methods and then sintered at 700 °C for Ag and 1200 °C for ITO in an ambient atmosphere. Sensors were evaluated at 600 °C for Ag and 1000 °C for ITO and analyzed at set interrogating distances up to 0.75 m using ultra-wideband slot antennas to collect scattering parameters. The sensitivity (average change in resonant frequency with respect to temperature) from 50 to 1000 °C was between 22 and 62 kHz/°C which decreased as interrogating distances reached 0.75 m.

47 OTHER INSTRUMENTATION↗

Experimental validation of an organic rankine-vapor compression cooling cycle using low GWP refrigerant R1234ze(E)

There is a significant global opportunity to capture and utilize low grade waste heat to reduce fossil fuel consumption, greenhouse gas emissions, and improve energy efficiency across a wide range of industries. In this work, an advanced type of thermally activated cooling system, an organic Rankine-vapor compression cycle (ORVC) with novel heat integration strategies, was designed and tested at a relevant scale for industrial waste heat recovery (300 kW th cooling capacity). The ORVC linked an organic Rankine power cycle and a vapor compression cooling cycle using a turbine and compressor that shared a single shaft. The ORVC test facility absorbed waste heat from a liquid stream at 91 °C to simulate engine coolant in diesel generator sets, rejected heat to a glycol stream at 30 °C, and generated chilled water at 7 °C. A cooling capacity of 264 kW ± 3.5 kW was experimentally validated with a COP of 0.56 ± 0.01 during steady-state operation at the design temperatures. The thermal efficiency, accounting for pump work, of the Rankine cycle was 7.7% ± 0.22% and the COP of the vapor compression cycle was 5.25 ± 0.09. The centrifugal turbo-compressor operated at 31.5 kRPM ± 0.3 kRPM, with a turbine and compressor isentropic efficiencies of 76.7% ± 0.90% and 84.8% ± 0.54%, respectively, with near-perfect power transmission between these components. The pressure drop in the piping and heat exchangers were significantly larger than expected which had a detrimental impact on the performance of the ORVC. In addition, the condenser on the cooling cycle could not deliver the subcooling as specified from the design point modeling. Furthermore, the results from the sensitivity analysis showed that the higher condenser glycol outlet temperature had the largest impact on performance, which is consistent with other analytical models in the literature. When the ORVC simulations were updated with experimental values for isentropic efficiencies of the turbomachinery, the thermal COP was 0.66 which represents an estimate of the predicted performance if test facility limitations are overcome.

30 DIRECT ENERGY CONVERSION↗

Flow boiling transient critical heat flux tests with stainless steel and FeCrAl: Transient correlation implementation, model calibration, and sensitivity analysis

In this study, steady-state and transient internal-flow critical heat flux (CHF) experiments were carried out under two flow conditions, at atmospheric pressure, on 316 stainless steel (316-SS) and iron-chromium-aluminum (FeCrAl) tubes. Slow power transients at a low mass flow (300 kg/m 2 - s) with very low subcooling (X e =-0.0054) generated premature CHFs, which were prevented by faster power transients and high mass flow (1000 kg/m 2 - s). The measured transient CHFs increased linearly with increasing power transients, compared with the steady-state CHF. Yet, the wall superheat at the CHFs decreased with faster power transients. Transient CHF correlations highlighting heterogeneous spontaneous nucleation were calibrated to the measured CHFs and compared with other existing correlations. Transient CHF multipliers were acquired from pool and flow boiling empirical CHF correlations that were generated. The multipliers were applied to the RELAP5-3D nuclear system code to analyze the discrepancies between the measured data and the predicted CHF and post-CHF behavior, which improved peak cladding temperature predictions by 24.6%. A variance-based global sensitivity study perturbing the experimental uncertainties and cladding material thermal properties showed the diminishing influence of flow boiling heat transfer with increasing power transients, highlighting the significance of the volumetric heat capacity for the cladding integrity during transients. Transient CHF correlations were applied to the most limiting design basis accident: a hot-zero-power reactivity-initiated accident on a generic pressurized water reactor RELAP5-3D model.

36 MATERIALS SCIENCE↗

Stress Relaxation Cracking of Alloys at Temperatures Higher Than 540°C

Type 347H stainless steel (347H SS), used in commercial concentrating solar power (CSP) thermal energy storage to store solar-salt at a temperature of 565°C, has been reported in the literature to be susceptible to stress-relaxation cracking (SRC). The welded heat-affected zone (HAZ) and fusion zone (FZ) of 347H SS, particularly in thick sections, are known to be susceptible to failure during post-weld heat treatment (reheat cracking). SRC could also occur after months or years under an elevated-temperature service environment. Two conditions must be present for failure to occur in the HAZ and/or FZ: 1) a modified or sensitized microstructure and 2) sufficiently high tensile residual stresses present at the elevated service temperature. The overarching goal of this project is to recommend SRC mitigation protocols to avoid susceptibility to fail through SRC at temperatures relevant for CSP. Post weld heat treatment conditions and alternative alloys are investigated as potential mitigation solutions to SRC. We used Gleeble thermomechanical simulation tests and finite element (FE) models to understand the susceptibility of 347H SS to SRC as a function of temperature, stress, and microstructure. We started with a literature review of weldability issues with 347H SS and techniques to mitigate SRC in 347H SS weldments. Next, we used Gleeble experiments to determine reheat cracking susceptibility in the simulated HAZ of 347SS weldments and an alternative alloy, 316L SS (NUCL 167 SPH) with boron added. We also performed Gleeble experiments to compare the reheat cracking susceptibility of 347H cross welded with two different fillers: E347, which is used in some commercial CSP TES tanks, and E16.8.2. We validated the experimental results with FE models of the residual stresses. We found that although the 316L (NUCL 167 SPH) is less susceptible to reheat cracking, it is slightly weaker than 347H and ASME BP&V codes limit its use to a service condition of 565°C. We also found that welds using E16.8.2 as the weld filler are less susceptible to failure than those using E347, likely due to the higher creep ductility of E16.8.2, and that it may be used as an alternative filler for repair welding of 347H welds or as the primary choice of filler for newly developed weld joints. We also found that post weld heat treatment could be a viable solution for mitigating stress in E347-347H SS thick, constrained welds, like those found in CSP tanks, and propose several options for mitigating SRC in existing and future TES tanks.

14 SOLAR ENERGY↗

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Using 3D position sensitivity to reveal response non-uniformities in CdZnTe, TlBr, and CsPbBr 3 detectors

Position sensitivity enables the correction of response non-uniformities in room-temperature semiconductor detectors caused by crystal defects and other factors. It can also be used to pinpoint the exact location of crystal defects responsible for the response variations. This work describes a technique for revealing and visualizing the detector regions affecting the charge collection efficiency in CdZnTe (CZT), TlBr, and CsPbBr 3 detectors configured as position-sensitive virtual Frisch-grid (VFG) devices. The technique correlates the photopeak events in energy spectra with their spatial distributions inside the detectors using the position information. By selecting the events from narrow energy intervals within a photopeak, we can visualize the detector volumes with particular charge collection efficiencies, which, in turn, correlate with the locations of electrode and crystal defects. Here, we demonstrate this technique in several examples. Columnar structures in the volume plots (position distribution maps) are consistent with signal losses near or at the anode in selected samples of CZT and TlBr. Structures exhibiting a distinct depth dependence are consistent with grain boundaries or other crystal defects.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Assessing thermal comfort and participation in residential demand flexibility programs

Residential space-conditioning-based demand flexibility (DF) has become an increasingly sought-after method for demand-side load management to enhance grid reliability and facilitate integration of renewable energy generation. However, predicting the effectiveness and flexibility of residential DF resources is challenging due to the variability in household energy use behaviors. Current estimates show that only 50 % of projected savings from DF resources are actualized due to regulatory, technological, and social barriers. From a household perspective, concerns over thermal comfort during space conditioning-based DF events significantly impact participation decisions. Currently, there is a very limited understanding of how thermal comfort during space-conditioning-based DF events in real-world settings impacts household energy use behaviors and, consequently, the success of DF programs in achieving targeted savings. This paper proposes a method to comprehensively assess the thermal comfort implications of DF strategies and presents results of their impacts on DF event participation decisions and demand savings. Here, the proposed method was applied to a heat pump DF field study in Cordova, Alaska. The study’s key findings are: 1) DF event setpoint offsets that maintain indoor operative temperatures between 18 to 22 °C (65 to 71°F) may be preferred in Cordova, Alaska; 2) Household-level thermal comfort is more sensitive to the duration of the DF event than to the degree of temperature offset from baseline conditions; 3) The delayed impact of changes in indoor operative temperature in response to setpoint offsets, both during and after a DF event, influences occupants’ thermal comfort perceptions and willingness to persistently participate in events. The findings from application of the proposed method can help inform future larger-scale occupant-centric DF programs as it can capture information not readily available through utility and device-level energy use data. Thus, it can supplement these sources and help program administrators develop occupant-centric DF strategies, enabling more accurate predictions of participation rates and savings estimates for space-conditioning-based DF programs.

Demand side management↗

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗

Impacts of the Thermal Gradient on Inland Advecting Sea Breezes in the Southeastern United States

Sea breezes are frequently observed in the South Carolina/Georgia region of the Southeastern United States (SEUS) and can reach upwards of 150km inland. This region is unique among the places frequently affected by sea breeze due to it being a continental location with relatively flat topography. The thermal gradient between land and water environments is a factor in introducing the sea breeze, but its role in the inland extent of sea breeze propagation isn’t as well known. We investigate the role of the thermal gradient in previously catalogued sea breeze events observed at the Savannah River Site (SRS) by taking differences of temperature measurements at inland and coastal weather stations for the days that the events occurred. We saw that the temperature differences for those days were much higher than in the non-sea breeze days during the mornings and afternoon. Numerical models were also used to conduct a sensitivity study on a sea breeze case, using simple modifications of the temperature gradient. We found that while the modifications did not stop the generation of a sea breeze circulation, the extent of the inland propagation was dependent on the magnitude of the thermal gradient.

54 ENVIRONMENTAL SCIENCES↗

Wind Turbine Main Bearing Rating Lives as Determined by IEC 61400-1 and ISO 281

This presentation studies the rating lives of wind turbine main bearings, as determined by the IEC 61400-1 and ISO 281 standards. A brief review of relevant bearing life theory and turbine design requirements is provided. This includes a discussion on possible shortcomings, including the existence (or not) of the bearing fatigue load limit and the validity of assuming linear damage accumulation. A detailed study is then undertaken to determine rating lives for two models of main bearing in a 1.5 MW wind turbine. Rating life assessment is carried out under different conditions, including various combinations of main bearing temperature, wind field characteristics, lubricant viscosity and contamination levels. Rating lives are found to be sufficiently above the desired 20 year design life for both bearing models under expected operating conditions. For the larger bearing, operational loads are shown to be below or close to the bearing fatigue load limit a vast majority of the time. Key sensitivities for rating life values are shown to be the temperature and contamination. Overall, the results of this study suggest that rating life assessment does not account for reported rates of main bearing failures in 1 to 3 MW wind turbines. In future work, it is recommended that efforts be undertaken to identify principal root causes of main bearing failures, possibly leading to a new application standard specific to this component. It is also recommended that impacts of partial wake impingement on main bearing rating lives are investigated.

ENGINEERING,WIND ENERGY↗

PMT-004 Validation Testing for Thermal Scattering Laws and Some Interesting Temperature-Dependent Implications [Slides]

This presentation finds that PMT-004 has four new benchmark cases highly sensitive to PE TSL (2 cases) and PMMA TSL (2 Cases). Additionally discovered was that the PE cases were well predicted using MCNP6.2 and ENDF/B-VIII.0. PMMA cases overpredicted by approximately 0.6-0.7% in k eff likely due to the lack of integral benchmarks sensitive to PMMA. The NCSU ENDF/B-VIII.1 PMMA TSL is an improvement over VIII.0. This was Accepted at last ICSBEP meeting pending comment resolution, should be published in 2024 version of the Handbook. It was found that temperature had a large impact on reactivity of the critical configuration. This has Implications for validation work for thermal cases- need to adjust TSL data to correct temperature as it can have hundreds of pcm effects for a few °C. Future thermal experiments should try and measure reactivity at multiple temperatures to add in data testing and benchmark adjustment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SILICON CARBIDE THERMOMETRY USING RAMAN SPECTROSCOPY ON IRRADIATED TRISO PARTICLES

Silicon carbide (SiC) passive thermometry has emerged as a promising post-irradiation examination (PIE) technique for estimating irradiation temperature near the end of irradiation. While dilatometry techniques have been traditionally used to analyze prismatic samples after irradiation, Raman spectroscopy has recently been shown to provide comparable results by analyzing Raman-active phonon modes. In this work, Raman spectroscopy has been applied to the SiC layer of cross-sectioned irradiation tristructural isotropic (TRISO) particles from the AGR-5/6/7 experiment to evaluate the feasibility of particle-scale passive thermometry. Two-dimensional Raman mapping was used to measure the position of the SiC longitudinal optical (LO) phonon, which was then converted to an apparent irradiation temperature using a previously established empirical correlation. Particles from AGR-5/6/7 Compacts 2-2-1 and 5-1-3 were selected as their calculated time-averaged, volume-averaged (TAVA) temperatures (828°C and 706°C, respectively) fall within the range of the sensitivity of the experimental approach. The use of SiC thermometry was anticipated to confirm or highlight potential deviations from calculated end-of-life TAVA temperature across compacts. Four particles from Compact 2-2-1 and two particles from Compact 5-1-3 were selected based on 110mAg inventory (either some measurable activity or below the minimum detection limit) which is commonly used as an indicator of in-pile temperature variation of particles within the same compact. Across every particle selected it was determined that the average LO peak position was located around 968 cm-1 to 969 cm-1. Using the previously determined empirical correlation, this corresponds to an irradiation temperature around 950°C to 966°C, which does not align with the reported TAVA temperature values. This discrepancy in apparent irradiation temperatures likely reflects a combination of uncertainties in the calculated particle temperatures and differences in irradiation history between the present specimens and those used to establish the empirical Raman calibration such as neutron flux (damage rate) and SiC microstructure (as fabricated and irradiated).

Vawdrey, Josh [ORNL]↗

Measurement accuracy and spatial resolution of a distributed temperature sensor based on a two-pulse differential coherent reflectometer

We present a model and numerical simulation of a distributed temperature sensor based on a two-pulse differential coherent optical time-domain reflectometer (COTDR). The differential phase measured using a phase-sensitive Rayleigh reflectometer is shown to have a regular component, which is a linear function of temperature, and a random component, which is related to a random distribution of scattering centres in the fibre and restricts the accuracy of measurements of variations in temperature. Measurement accuracy can be improved by reducing the relative contribution of the random component via a decrease in pulse duration and/or an increase in the time delay between pulses. The spatial resolution of a differential two-pulse phase-sensitive reflectometer is shown to be determined by the time delay between pulses and to vary little with pulse duration. At a typical pulse duration (200 ns) and delay time (300 ns), the accuracy in measurements of variations in temperature in the 0.1-K range is 2 % and the spatial resolution is about 30 m. (fibre-optic sensors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Perovskite-Sensitized Upconversion under Operando Conditions

Perovskite-sensitized upconversion (UC) has resulted in near-infrared-to-visible UC at solar-relevant fluxes. However, the successful implementation of UC devices into operating solar cells will result in exposure to similar environmental stressors as for the commercial photovoltaics (PVs), mainly elevated temperatures, and continuous irradiation. In this article, we investigated the effects of these two stressors, heat and light, on the triplet generation process at the perovskite/rubrene interface. Following exposure to both stressors, local discrepancies across the upconversion device were discovered. Here, the first region showed changes to the morphology, and no detectable upconverted emission was observed. Through the combination of optical microscopy and spectroscopy, crystallization of the organic semiconductor layer, degradation of dibenzotetraphenylperiflanthene, and concurrent degradation of the perovskite sensitizer were found. These effects culminate in a reduction in both triplet generation and triplet–triplet annihilation. In the second region, no changes to the morphology were present and visible UC emission was observed following exposure to both stressors. To probe the triplet sensitization process at elevated temperatures, transient absorption spectroscopy was performed. The presence of the excited spin-triplet state of rubrene at 60 °C highlighted successful triplet generation even at elevated temperatures. This work emphasizes the challenges and continued potential for the integration of perovskite-sensitized UC into commercial photovoltaic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis of nanodiamonds encapsulated by zeolitic imidazole framework-8 for quantum sensing applications

Nitrogen vacancy (NV)-containing nanodiamonds are widely used in quantum sensing applications due to their high sensitivity to magnetic fields, relatively low cost, and ability to be initialized, manipulated, and read out at room temperature. Quantum sensing techniques such as optically detected magnetic resonance (ODMR) and spin relaxometry have exploited the sensitivity of the NV nanodiamonds to magnetic fields to detect a range of analytes, such as pH, metal ions, and biomolecules. However, diversifying the sensing targets accessible by NV diamond quantum sensors typically requires careful engineering of the diamond surface chemistry with stimuli-responsive functional groups. Here, a simple protocol for coating NV nanodiamonds with the zeolitic imidazole framework 8 (ZIF-8), a widely used metal-organic framework, is presented. ZIF-8 is a highly porous material that has been used as a selective sensor for gasses, metal ions, and other analytes. The material is well-characterized by x-ray diffraction, transmission electron microscopy, scanning electron microscopy, x-ray photoelectron spectroscopy, and luminescence spectroscopy. Encapsulation of NV nanodiamonds with a porous scaffold such as ZIF-8 provides a promising method for improving the selectivity for the quantum sensing of various analytes. Importantly, the ZIF-8 coating does not impact the luminescence properties of the NV diamond, which is a key readout in ODMR and spin relaxometry sensing approaches. Indeed, the ODMR spectra with and without the ZIF-8 shell is nearly identical. Moreover, the ZIF-8 coating increases the longitudinal spin relaxation time of the NV nanodiamond by a factor of 4 relative to aggregated diamond, a desirable outcome for spin relaxation-based quantum sensing. Metal-organic framework composites with nanodiamonds thus are an exciting strategy for enhancing NV nanodiamond performance in applications such as quantum sensing and quantum-enhanced nuclear magnetic resonance spectroscopy.

nitrogen vacancy nanodiamond↗

Skeletal Kinetics Reduction for Astrophysical Reaction Networks

A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to Type Ia supernovae (SNe Ia). The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of 56 Ni in SNe Ia is produced for different temperatures, densities, and proton-to-neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of α-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and 44 Ti and 56 Ni abundances by each model. The consequences of ye ≠ 0.5 in the initial composition are also explored where ye is the electron fraction. The simulated results show that 56 Ni production decreases by decreasing ye as expected, and that the 43 Sc is a key isotope in proton and neutron channels toward 56 Ni production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the 56Ni abundance in SNe Ia for ye ≲ 0.5 initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗