Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High Temperature Methods”

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 847 records · Page 47

Comprehensive defect evaluation of advanced nuclear fuels using high-resolution acoustic signals and optimized sensor separation

Graphite pebble composite structures based on TRistructural-ISOtropic (TRISO) particles are being developed as core nuclear fuels in advanced power reactors, promising safe operation at increased temperatures. Ensuring the structural integrity of these nuclear fuels requires comprehensive and accurate non-destructive evaluation (NDE) techniques to characterize defects and damage in the pebbles. However, traditional acoustic evaluation methods face limitations in defect characterization due to the highly attenuative, and geometrically and compositionally complex nature of these structures. This study proposes an improved acoustic NDE technique for accurate detection and classification of anticipated relevant defects and damage in graphite pebbles using high-resolution acoustic signals and optimized transmit-receive sensor networks. The proposed approach utilizes a triangular three-sensor network as the base unit, comprising three transmit-receive sensors. The sensor separation distance, as well as acoustic excitation center frequency, pulse-width, and bandwidth are optimized to enhance spatial resolution and improve signal-to-noise ratio, enabling effective characterization of the smallest size and widest range of defects in pebbles. Furthermore, the use of the triangular sensor configuration instead of a more conventional transmit-receive sensor pair expands the inspection region from a one-dimensional linear path to a two-dimensional area, increasing spatial coverage. To mitigate challenges associated with processing of complex acoustic signals arising from high-frequency, high-bandwidth excitation in these structures, a machine-learning-based signal processing algorithm is integrated with the sensor network. In the machine-learning-based algorithm, multi-domain features are extracted from the acoustic signals to capture intricate signal characteristics, significantly improving defect identification and classification compared to traditional approaches. The proposed acoustic NDE technique offers considerable promise for practical and reliable defect/damage diagnostics of advanced nuclear pebble fuels.

42 ENGINEERING↗

A simplified approach for creep damage and remaining life calculation for materials surveillance in advanced reactors

High temperature reactor materials will experience a combination of mechanical degradation caused by creep and fatigue and environmental degradation caused by neutron irradiation and, in some cases, exposure to corrosive coolants. A materials surveillance program is one option for ensuring the safe, reliable operation of key reactor components under these conditions. Such a program would monitor the degradation in key material properties over time, using this data to predict how changes in the material properties affect components performance. This could then be used to make plant operational decisions and support eventual plant license extensions. Past work de-scribes such a material surveillance program based on passively actuated mechanical test articles that impose creep-fatigue type loading on test materials driven only by changes in temperature, for example those experienced by components under standard operating cycles. This report focuses on two aspects of the proposed material surveillance program: how to design the test articles to mimic the mechanical response of the corresponding component and how to use ex-situ test data to deter-mine the amount of damaged experienced by the test article (and hence the corresponding compo-nent) in service. Specifically, this report derives and validated simplified methods for both tasks, replacing earlier, more complicated approaches based on simulating the test article response and finding the best specimen design or current damage via complex numerical optimization. The sim-plified methods for each task developed here can be implemented in spreadsheet software and are simple enough for practical use in future operating plants. This report derives the methods (one for sizing and two options for damage inference), verifies the simplified approaches versus more so-phisticated methods, and compares the results of applying the simplified methods to previous re-sults using the complex numerical optimization approach. The report also provides a worked sam-pled problem applying the simplified techniques to a realistic high temperature reactor component.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Predicting melt pool depth and grain length using multiple signatures from in-situ single camera two-wavelength imaging pyrometry for laser powder bed fusion

In laser powder bed fusion (LPBF), the in-situ process signatures are known to have a direct correlation with the microstructural properties of the solidified melt pool (MP). It is known that the MP cooling and heating rates, and laser processing parameters can critically determine the grain structure and thereby affect the part properties. The objective of this work is to study the feasibility of using in-process, high-speed imaging pyrometry for evaluating the solidified MP properties “below” the surface, such as depth and microstructural properties. To accomplish this, we employ an in-house single camera-based two-wavelength imaging pyrometry (STWIP) system for monitoring the printing of single-scan tracks with Inconel 718 on a commercial LPBF printer (EOS M290). Further, the lab designed STWIP system is a coaxial high-speed (>10,000 fps) imaging system capable of monitoring MP temperature, morphology, and intensity profiles. The temperature measurements from STWIP are emissivity independent. The STWIP measured MP signatures of the printed tracks are correlated with the ex-situ microscopy characterized MP depth and the average grain lengths. From the data analysis, using support vector machine (SVM)-based regression models, we found that the MP temperature signatures are crucial for an accurate prediction of MP depth and the grain length, thus validating the novelty and necessity of the developed in-situ monitoring methods and analysis.

36 MATERIALS SCIENCE↗

Ion-beam assisted sputtering of titanium nitride thin films

Abstract Titanium nitride is a material of interest for many superconducting devices such as nanowire microwave resonators and photon detectors. Thus, controlling the growth of TiN thin films with desirable properties is of high importance. This work aims to explore effects in ion beam-assisted sputtering (IBAS), were an observed increase in nominal critical temperature and upper critical fields are in tandem with previous work on Niobium nitride (NbN). We grow thin films of titanium nitride by both, the conventional method of DC reactive magnetron sputtering and the IBAS method, to compare their superconducting critical temperatures $$T_{c}$$ T c as functions of thickness, sheet resistance, and nitrogen flow rate. We perform electrical and structural characterizations by electric transport and x-ray diffraction measurements. Compared to the conventional method of reactive sputtering, the IBAS technique has demonstrated a 10% increase in nominal critical temperature without noticeable variation in the lattice structure. Additionally, we explore the behavior of superconducting $$T_c$$ T c in ultra-thin films. Trends in films grown at high nitrogen concentrations follow predictions of mean-field theory in disordered films and show suppression of superconducting $$T_c$$ T c due to geometric effects, while nitride films grown at low nitrogen concentrations strongly deviate from the theoretical models.

36 MATERIALS SCIENCE↗

Residual resistance ratio measurement system for Nb 3 Sn wires extracted from Rutherford cables

Residual resistance ratio (RRR) of superconducting strands is an important parameter for magnet electrical stability. RRR serves as a measure of the low-temperature electrical conductivity of the copper within a conductor that has a copper stabilization matrix. For Nb 3 Sn, due to the need of a reaction heat treatment, the technical requirements for high quality measurements of strands extracted from Rutherford cables are particularly demanding. Quality of wire, cabling deformation, heat treatment temperature, heat treatment atmosphere, sample handling, and measurement methods can all affect the RRR. Therefore, as an integral part of the electrical quality control (QC) of Nb 3 Sn Rutherford cables manufactured at the Lawrence Berkeley National Laboratory, it was prudent that we established a RRR measurement system that can isolate the assessment of cable-fabrication-related impacts from sample preparation and measurement factors. Here we describe a bespoke cryocooler-based measurement system, capable of measuring RRR of over 80 samples in a single cooldown. The samples are mounted on custom-designed printed circuit boards that accommodate the shape of strands extracted from a Rutherford cable without added deformation, which we will show is critical in ensuring that the measurements accurately represent the RRR values of the conductor within the cable. Using this sample mounting solution, we routinely measure the overall RRR of the strand as well as individual intra-strand sections corresponding to both cable edges and cable broad faces with high reproducibility. Such measurements provide valuable information on the variation of RRR along the length of the strands as well as across strand productions and cable runs over time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

In Situ Microstructure Characterization of Potassium Di-Phosphate (KDP) Densification during Cold Sintering

In order for ceramic additive manufacturing (AM) to achieve its full potential, it is increasingly important to develop a more rigorous understanding of fundamental phenomena that govern the kinetics and thermodynamics of ceramic AM processes. In the case of additive build processes, such as direct ink write and ceramic extrusion, methods for densifying the resulting green-body product need to be considered to complement the efficiencies of ceramics AM, itself. One densification route, at least for monolithic components, built layer-by-layer, is offered by the recently developed cold sintering process, whereby high-density final product is achieved through addition of a small amount of liquid solvent and application of modest uniaxial compressive stress at relatively low temperature. In situ small-angle X-ray scattering methods and X-ray diffraction have been applied to characterize and quantify the pore morphology evolution during cold sintering for a model system: potassium di-phosphate, KH 2 PO 4 (KDP). It is shown that both temperature and applied stress affect the densification rate, but stress has a stronger effect on the evolving morphology. A regime with an approximate linear densification rate can be identified, yielding an effective densification activation energy of ≈90 kJ/mol.

36 MATERIALS SCIENCE↗

Cost analysis of hydrogen production by high-temperature solid oxide electrolysis

In this study, we estimate construction and operation costs of gigawatt-scale solid oxide electrolysis (SOE) facilities for producing high purity hydrogen gas from water. Manufacturing and assembly costs for two types of SOE cell stacks are estimated using a detailed design for manufacture and assembly (DFMA®) analysis. Modular balance of plant (BOP) process equipment is designed and sized with Aspen®, and cost estimated using equipment vendor quotes. Factory and on-site assembly and installation costs for SOEC stack and BOP equipment integration into modular SOE process units are calculated using a simplified DFMA® method. Total stack costs on a stack input power (SIP) basis reduce to <$100/kW e DC SIP for >500 MW e DC SIP /year production rates with electrode cermet, interconnects, and high-temperature heat treatments dominating the total cost. Integration of stacks with larger BOP equipment operating at higher pressures offers ~36% cost reduction in total facility capital cost due to an economies of physical scale effect since BOP equipment comprises >50% of facility costs. Optimized H 2 prices decrease from ~$\$$4/kgH 2 to ~$\$$2/kgH 2 for 1 GWe DCSIP facilities using $\$$0.025/kWh electricity price. All costs are reported in 2021 US dollars.

08 HYDROGEN↗

Low Cost Glass-Ceramic Matrix Composite Heat Exchanger

As part of ARPA-E’s High Intensity Thermal Exchange through Materials and Manufacturing Processes (HITEMMP) program, this project sought to develop novel heat exchanger (HX) capabilities to enable efficient and power dense power generation cycles. This class of HX comes under the category of ceramic/composite materials with the higher temperature goal in the program of ≥1100 °C inlet temperature operation. The enabling capability of this effort is the use of glass-ceramic matrix composite (GCMC) material which provides the high temperature durability of a ceramic, the flaw tolerance of a composite, a significantly faster and lower cost manufacturing process than conventional matrix CMCs and very low porosity levels < 0.5%. For thin-walled HX structures and the need to minimize leakage, the low porosity differentiator is particularly important. RTRC has prior experience with this material system and in the current project advanced the component design and manufacturing methods into new territory to produce features required for effective heat exchange under high pressures. In this approach, silicon carbide fiber is fabricated into a fiber preform using various textile processes. Graphite tooling is used both during the build-up of the fiber preform (interior tooling) and after the fiber preform has been completed (exterior tooling). This tooling assembly is heated to high temperature in an environment that has been evacuated and backfilled with inert gas. A reservoir of specialty glass is present and once the desired temperature has been reached to achieve the desired glass viscosity, an actuator distributes the glass throughout the fiber preform using passageways which are part of the tooling design in a process known as glass transfer molding. After the tooling has been removed, the composite is heat treated to convert the amorphous glass to a crystalline ceramic, providing improved properties. The project was divided into three phases focusing on the following: 1) 10 kW HX design and coupon-level tube sheet fabrication, 2) 10 kW HX fabrication, 3) 50 kW HX fabrication. During Budget Period 1 (BP1), additional risks were encountered and the need for additional funds was agreed upon by ARPA-E program leadership. Due to a variety of factors, the contract modification required nominally 18 months to execute at which time the HITEMMP program was effectively concluding. Because of this and the time that would be required to perform BP2 tasks, it was decided to conclude the project at the end of BP1. During the design of the 10 kW HX, manufacturing constraints were learned and incorporated, leading to a revised configuration for the fiber preform and HX. Heat exchange and pressure drop predictions also played a role in modifying the original design concept to be a higher aspect ratio shell-and-tube HX, simplifying the manufacturing process and improving the heat exchanger performance. Good gravimetric and volumetric thermal power densities of 11.2 kW/kg and 10,200 kW/m3 for the entire HX were projected that involved thermo-structural Finite Element Analysis to determine the structural mass needed for the high operation pressures of 250 bar cold inlet and 80 bar hot inlet. Fiber preforms using textile processes were produced for multiple headered tube sheets. Additional challenges were encountered during the glass transfer molding step for which solutions were identified, but programmatics did not allow them to be implemented in BP1. While complete HX test articles were not fabricated, the benefits of this GCMC material for a variety of high temperature applications remain.

30 DIRECT ENERGY CONVERSION↗

Reduced Order Models Generation for HTGRs Pebble Shuffling Procedure Optimization Studies

This report provides an initial study for producing reduced-order models (ROMs) of pebble-bed high temperature gas reactor (HTGR) models for the purposes of design optimization. As an initial study, this work is meant to be exploratory---identifying useful workflows and methods for ROM generation---and not meant to be a catch-all analysis of HTGR ROM generation and usage for optimization. This report summarizes three tasks performed in Fiscal Year 2022: 1) the creation of HTGR model, 2) the sensitivity analysis of model design parameters, and 3) an introduction to ROM generation techniques. The representative HTGR model created in this work is a multiphysics equilibrium-core using the BlueCRAB (comprehensive reactor analysis bundle) reactor analysis application, coupling four physical phenomena: neutronics, streamline depletion, porous flow thermal hydraulics, and pebble heat conduction. Part of the model creation was identifying some design parameters and quantities of interest that are relevant in an optimization analysis and adjustable in the model. The sensitivity analysis utilized a polynomial chaos expansion methodology to compute global sensitivity metrics. This analysis showed that thermal hydraulics parameters and quantities of interest had a relatively small impact on simulation results. Finally, the ROM generation work involved exploring three different ROM methodologies: polynomial regression, a Gaussian process, and artificial neural networks. Using a cross-validation technique to characterize ROM performance, the Gaussian process and single-layer artificial neural networks showed the most promising results. Overall, this study was insightful and the lessons learned will be invaluable for the eventual development of an HTGR design optimization workflow.

97 MATHEMATICS AND COMPUTING↗

Hydrogenolysis of Polypropylene and Mixed Polyolefin Plastic Waste over Ru/C to Produce Liquid Alkanes

Catalytic depolymerization of polyolefins is a promising chemical recycling strategy to create value-added products from waste plastics, which are accumulating in landfills and the natural environment at unsustainable rates. The cleavage of strong C–C bonds in polyolefins can be performed using a noble metal and hydrogen via a hydrogenolysis mechanism. Previously, we identified ruthenium nanoparticles supported on carbon (Ru/C) as a highly active heterogeneous catalyst for the conversion of polyethylene into liquid and gaseous n-alkanes under mild conditions. In the present study, we investigated the catalytic depolymerization of polypropylene (PP) under mild conditions (200–250 °C, 20–50 bar H 2 ). We demonstrate that Ru/C produces C 5 –C 32 iso-alkane yields above 68% in the absence of solvent and identify trade-offs between product yield and temperature, hydrogen pressure, and reaction time. We apply a rigorous analytical method to quantify all liquid and gaseous alkane products. The characterized catalyst was found to be recyclable after the complete depolymerization of high molecular weight PP (M w ~ 340,000 Da) to liquid and gaseous hydrocarbons and after depolymerization of a postconsumer PP centrifuge tube. Finally, the catalyst was shown to be effective in depolymerizing a mixture of high-density polyethylene and PP to produce a mixture of linear and branched liquid alkanes, demonstrating feasibility for the depolymerization of streams of mixed polyolefin waste.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase transition kinetics revealed by in situ x-ray diffraction in laser-heated dynamic diamond anvil cells

We report successful coupling of dynamic loading in a diamond anvil cell and stable laser heating, which enables compression rates up to 500 GPa/s along high-temperature isotherms. Dynamic loading in a diamond-anvil cell allows exploration of a wider range of pathways in the pressure-temperature space compared to conventional dynamic compression techniques. By x-ray diffraction, we are able to characterize and monitor the structural transitions with the appropriate time resolution i.e., millisecond timescales. Using this method, we investigate the γ − ε phase transition of iron under dynamic compression, reaching compression rates of hundreds of GPa/s and temperatures of 2000 K. Our results demonstrate a distinct response of the γ − ε and α − ε transitions to the high compression rates achieved, possibly due to the different transition mechanisms. These findings open up new avenues to study tailored dynamic compression pathways in the pressure-temperature space and highlight the potential of this platform to capture kinetic effects (over ms time scales) in a diamond anvil cell. Published by the American Physical Society 2024

Ricks, Matthew (ORCID:000000015315612X)↗

Report on FY 2021 testing in support of integrating Alloy 800H and 2.25Cr-1Mo for EPP analysis

The use of simplified Elastic-Perfectly Plastic (EPP) design methods avoids the complexities and limitations of the elastic analysis approach in evaluating strain limits and creep-fatigue damage for Class A metallic coolant boundary components in ASME Boiler and Pressure Vessel Code, Section III, Division 5, Subsection HB, Subpart B. No stress classification and linearization are required in the EPP methods. Further, the strain limits and creep-fatigue damage evaluation procedures based on elastic analysis results were developed for the reactor operating conditions where uncoupled, rate-independent plasticity and stationary creep can be used to approximate the deformation in the component adequately. But at very high temperatures, the deformation is rate-dependent, and it can only be approximated adequately by unified viscoplastic model. Thus, for example, the elastic analysis approach cannot be used to evaluate strain limits and creep-fatigue damage for Alloy 617 above 650 °C. The EPP methods remove the temperature restriction. They can be used to evaluate strain limits and creep-fatigue damage for the full range of temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Device and material investigations of GaN enhancement-mode transistors for Venus and harsh environments

This Letter reports the device and material investigations of enhancement-mode p-GaN-gate AlGaN/GaN high electron mobility transistors (HEMTs) for Venus exploration and other harsh environment applications. The GaN transistor in this work was subjected to prolonged exposure (11 days) in a simulated Venus environment (460 °C, 94 bar, complete chemical environment including CO2/N2/SO2). The mechanisms affecting the transistor performance and structural integrity in harsh environment were analyzed using a variety of experimental, simulation, and modeling techniques, including in situ electrical measurement (e.g., burn-in) and advanced microscopy (e.g., structural deformation). Through transistor, Transmission Line Method (TLM), and Hall-effect measurements vs temperature, it is revealed that the mobility decrease is the primary cause of reduction of on-state performance of this GaN transistor at high temperature. Material analysis of the device under test (DUT) confirmed the absence of foreign elements from the Venus atmosphere. No inter-diffusion of the elements (including the gate metal) was observed. The insights of this work are broadly applicable to the future design, fabrication, and deployment of robust III-N devices for harsh environment operation.

Physics↗

Controlled Synthesis of SnO 2 Nanocrystals with Tunable Band Gaps

Tin(IV) oxide nanocrystals (SnO 2 NCs) have significant potential in various applications, with their performance closely related to their band gap. The band gap is influenced by the size and shape of the NCs, which can be precisely controlled by adjusting reaction conditions. In this study, we present deliberately designed synthesis protocols to produce high-quality SnO 2 NCs with tunable band gaps using different methods. Key factors affecting the synthesis include control of the oxidizing agent, reaction temperature, solvent selection, and reaction time optimization. The resulting NCs were characterized by using TEM, XRD, XPS, and optical spectroscopy. Notably, SnO 2 NCs synthesized by controlling the oxidizing agent (air injection) in a hot organic solution were smaller in size and exhibited abundant oxygen vacancies. In contrast, extending the reaction time or using ethanol as a solvent in hydrothermal systems facilitated larger spherical or rod-like SnO 2 NCs with fewer oxygen vacancies. Further analysis of the band gap and valence band maximum energy revealed that the abundant vacancies in SnO 2 NCs synthesized with the air-controlled hot organic solution method resulted in a narrower band gap and an upshifted valence band. These synthetic strategies illustrate the potential for deliberately designing SnO 2 NCs with optimized electronic structures for various applications.

band gap↗

Real-space Green's function approach for x-ray spectra at high temperature

There has been considerable interest in properties of condensed matter at high temperature, including nonequilibrium behavior and extreme conditions up to the warm dense matter regime. Such behavior is encountered, e.g., in experimental time-resolved x-ray absorption spectroscopy in the presence of intense laser fields. In an effort to simulate such behavior, we present an approach for calculations of high-temperature x-ray absorption spectra in arbitrary materials, using a generalization of the real-space Green's function formalism. Here, the method is incorporated as an option in the core-level x-ray spectroscopy code FEFF10. To illustrate the approach, we present calculations for several materials together with comparisons to experiment and with other methods.

36 MATERIALS SCIENCE↗

The seventh blind test of crystal structure prediction: structure ranking methods

A seventh blind test of crystal structure prediction has been organized by the Cambridge Crystallographic Data Centre. The results are presented in two parts, with this second part focusing on methods for ranking crystal structures in order of stability. The exercise involved standardized sets of structures seeded from a range of structure generation methods. Participants from 22 groups applied several periodic DFT-D methods, machine learned potentials, force fields derived from empirical data or quantum chemical calculations, and various combinations of the above. In addition, one non-energy-based scoring function was used. Results showed that periodic DFT-D methods overall agreed with experimental data within expected error margins, while one machine learned model, applying system-specific AIMnet potentials, agreed with experiment in many cases demonstrating promise as an efficient alternative to DFT-based methods. For target XXXII, a consensus was reached across periodic DFT methods, with consistently high predicted energies of experimental forms relative to the global minimum (above 4 kJ mol −1 at both low and ambient temperatures) suggesting a more stable polymorph is likely not yet observed. The calculation of free energies at ambient temperatures offered improvement of predictions only in some cases (for targets XXVII and XXXI). Several avenues for future research have been suggested, highlighting the need for greater efficiency considering the vast amounts of resources utilized in many cases.

Chemistry↗

Laser Crystallization and Dopant Activation of a-Si:H Carrier-Selective Layer in TOPCon Si Solar Cells

Herein, we present a pulsed-laser processing method for crystallization and dopant activation of a highly n-doped amorphous silicon (a-Si:H) carrier-selective layer in a tunnel oxide passivated contact (TOPCon) Si solar cell structure. Here, the laser method provides enhanced conductivity and implied open circuit voltage while reducing emitter saturation current density and surface heating, as opposed to conventional high-temperature furnace annealing of the bulk Si wafer with a TOPCon structure. We identify an appropriate laser wavelength, fluence, and layer thickness using modeling and simulations. Raman and Hall effect measurements demonstrate increased crystallinity and dopant activation, whereas photoconductive decay shows enhanced surface and interface passivation quality. Additionally, we examine the role of subsequent SiNx deposition on further improving the passivation of laser-processed TOPCon layers to achieve a 5.9 kΩ/sq film sheet resistance, 2 ms effective carrier lifetime, 718 mV implied open-circuit voltage (iVOC), and 8.6 fA/cm 2 one-side recombination current density (J0) with the potential for further passivation improvement via laser process and a-Si layer thickness optimization.

14 SOLAR ENERGY↗

Investigating the Combustion Performance of Dual Fuel Combustion with Diesel and Port Injected Hydrogen in a Large Bore Locomotive Engine

The heavy-duty transportation sector has primarily relied on conventional diesel combustion engines given their reliability and high thermal efficiency relative to spark ignition engines, but increased focus on reducing greenhouse gas emissions has led to investigation into alternative fuels. Gaseous hydrogen fuel has garnered a great deal of recent interest in the engine community given it has zero carbon, but hydrogen is not available at the scale and cost that petroleum fuels are currently available, and this is a barrier to adoption for industries that are looking to decarbonize their operations. Because of the fuel flexibility provided, dual fuel technology offers a pathway for some industries to adopt hydrogen as a fuel source while maintaining sufficient flexibility in times and locations where the new fuel is not yet available. This computational study investigates dual fuel combustion in a large bore locomotive engine architecture using direct injected diesel and port injected gaseous hydrogen fuel. With an optimal port fuel injection configuration from previous work, simulations of varying substitution ratio, compression ratio, manifold air temperature, diesel injection timing, and diesel injection pressure were performed to understand their effect on combustion performance. Results indicated that both increased substitution ratio and higher intake air temperature accelerates hydrogen flame propagation and can result in high peak cylinder pressures. Additionally, diesel injection timing and injection pressure were demonstrated as effective methods for controlling dual fuel combustion heat release rates.

ODonnell, Patrick Christopher↗