Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat transfer performance”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
NASA NTRS1995

Better Gas-Gap Thermal Switches For Sorption Compressors

Gas-gap thermal switches associated with sorption compressors of some heat pumps and cryogenic systems designed for higher performance, according to proposal, by introducing controlled turbulent flows into gas gaps. Utilizes convection in turbulent flow to transfer heat at greater rate. Design takes advantage of flow of working fluid. Working fluid also serve as heat transfer medium in gas gap.

Bhandari, Pradeep↗
NASA NTRSDate not supplied

Analysis of Cryogenic Propellant Liquefaction Rates in Cooled Constant-Wall-TemperatureTanks

NASA has been focused on developing technology that would allow the production of cryogenic propellants on the Lunar and Martian surfaces. Utilizing Lunar/Martian resources, the produced gaseous propellants must first be liquefied and stored prior to use on the Moon or Mars ascent vehicle. Liquefaction of cryogenic propellants is a necessary technical development to enable NASA’s future spaceflight goals. This paper presents an overview of a proposed model for a propellant liquefaction system, and the effect of tank wall temperature and the ullage pressure control band on cryogenic propellant liquefaction rates. The propellant system was assumed to be a receiver of a gas from an In-Situ Resource Utilization (ISRU) harvester and condensed the gas by cooling the tank walls. First, the heat transfer principles pertinent to condensation are summarized, followed by an overview of the model which performed the analysis, including a summary of the Fortran algorithm implementation. The sensitivity of the condensation rate relative to varying tank wall temperatures is also discussed.

cryogenics↗
NASA NTRS1989

Prediction of the aerodynamic environment and heat transfer for rotor-stator configurations

A numerical study of the aerodynamic and thermal environment associated with axial turbine stages is presented. Computations were performed using a modification of the unsteady viscous code, ROTORI, and an improved version of the steady inviscid cascade system, MERIDL-TSONIC, coupled with boundary layer codes, BLAYER and STAN5. Two different turbine stages were analyzed: the first stage of the United Technologies Research Center Large Scale Rotating Rig (LSRR) and the first stage of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump turbine. The time-averaged airfoil midspan pressure and heat transfer profiles were predicted for numerous thermal boundary conditions including adiabatic wall, prescribed surface temperature, and prescribed heat flux. Computed solutions are compared with each other and with experimental data in the case of the LSRR calculations. Modified ROTORI predictions of unsteady pressure envelopes and instantaneous contour plots are also presented. Relative merits of the two computational approaches are discussed.

Griffin, L. W.↗
NASA NTRS1992

Thermal response of rigid and flexible insulations and reflective coating in an aeroconvective heating environment

Described here is the thermal performance of rigid and flexible thermal protection systems considered for potential use in future Aeroassist Space Transfer Vehicles. The thermal response of these materials subjected to aeroconvective heating from a plasma arc is described. Properties that were measured included the thermal conductivity of both rigid and flexible insulations at various temperatures and pressures and the emissivity of the fabrics used in the flexible insulations. The results from computerized thermal analysis models describing the thermal response of these materials subjected to flight conditions are included.

Kourtides, D. A.↗
NASA NTRS1995

Heat cascading regenerative sorption heat pump

A simple heat cascading regenerative sorption heat pump process with rejected or waste heat from a higher temperature chemisorption circuit (HTCC) powering a lower temperature physisorption circuit (LTPC) which provides a 30% total improvement over simple regenerative physisorption compression heat pumps when ammonia is both the chemisorbate and physisorbate, and a total improvement of 50% or more for LTPC having two pressure stages. The HTCC contains ammonia and a chemisorbent therefor contained in a plurality of canisters, a condenser-evaporator-radiator system, and a heater, operatively connected together. The LTPC contains ammonia and a physisorbent therefor contained in a plurality of compressors, a condenser-evaporator-radiator system, operatively connected together. A closed heat transfer circuit (CHTC) is provided which contains a flowing heat transfer liquid (FHTL) in thermal communication with each canister and each compressor for cascading heat from the HTCC to the LTPC. Heat is regenerated within the LTPC by transferring heat from one compressor to another. In one embodiment the regeneration is performed by another CHTC containing another FHTL in thermal communication with each compressor. In another embodiment the HTCC powers a lower temperature ammonia water absorption circuit (LTAWAC) which contains a generator-absorber system containing the absorbent, and a condenser-evaporator-radiator system, operatively connected together. The absorbent is water or an absorbent aqueous solution. A CHTC is provided which contains a FHTL in thermal communication with the generator for cascading heat from the HTCC to the LTAWAC. Heat is regenerated within the LTAWAC by transferring heat from the generator to the absorber. The chemical composition of the chemisorbent is different than the chemical composition of the physisorbent, and the absorbent. The chemical composition of the FHTL is different than the chemisorbent, the physisorbent, the absorbent, and ammonia.

Jones, Jack A.↗
NASA NTRS2017

Ablation and Heating During Atmospheric Entry and Its Effect on Airburst Risk

Large meteoroids and asteroids entering the atmosphere endure tremendous heating from the shock heated air, and thereby lose a significant fraction of their mass during atmospheric entry a process known as ablation. The predicted evolution of the asteroids mass as it passes through the atmosphere can affect both the predicted energy deposition profile relevant to an airburst event, or the residual mass that strikes the ground in the case of an impact event. This presentation is divided roughly into two parts. In the first part, an overview of traditional models for heat transfer and ablation that are historically used in the meteor physics community is presented, and the validity in the asteroid entry regime discussed. Sensitivity analyses performed using the recently developed Fragment-Cloud Model (FCM) will be presented which show illustrate the range of sizes and entry parameters for which the predicted asteroid threat is most sensitive to the models for ablation and heat transfer. The second part of the presentation shall focus on recent work done under NASAs Asteroid Threat Assessment Project (ATAP) to develop new models for heat transfer and ablation using high-fidelity numerical simulation in concert with state-of-the-art experiments. Coupled computational fluid dynamics (CFD)radiation transport simulations preformed using the state-of-the-art entry modeling tools at NASA show that, for large meteoroids and asteroids, there can significant attenuation of the heat transfer to the surface (95 in some cases) by the products of ablation. In addition to the heat transfer, new models for the material response and ablation of asteroidal material have been developed [cite]. In the current work, we present finding from recent novel experiments performed in the arc jet facility at NASA Ames, which allows us to, in part, simulate the extreme environment experienced by the asteroid during entry. Briefly, the experimental set-up was comprised of a 1.5 conical article of machined H5 chondrite, which was exposed to a high-enthalpy flow resulting in approximately 4 kWcm2 of heating to the surface. A still frame capture from high-speed video taken during this experiment can be seen in Figure 1. In this figure, we can observe some of the major mechanisms for meteoroid ablation, such as melt flow, spallation (mechanical removal of material), and vaporization. Major findings from this, and other experiments will be discussed, as well progress on utilizing the data from the experiments to inform and develop improved models for ablation.

Asteroids↗
NASA NTRS2010

Analysis of Unsteady Tip and Endwall Heat Transfer in a Highly Loaded Transonic Turbine Stage

In a previous study, vane-rotor shock interactions and heat transfer on the rotor blade of a highly loaded transonic turbine stage were simulated. The geometry consists of a high pressure turbine vane and downstream rotor blade. This study focuses on the physics of flow and heat transfer in the rotor tip, casing and hub regions. The simulation was performed using the Unsteady Reynolds-Averaged Navier-Stokes (URANS) code MSU-TURBO. A low Reynolds number k-epsilon model was utilized to model turbulence. The rotor blade in question has a tip gap height of 2.1 percent of the blade height. The Reynolds number of the flow is approximately 3x10(exp 6) per meter. Unsteadiness was observed at the tip surface that results in intermittent "hot spots". It is demonstrated that unsteadiness in the tip gap is governed by inviscid effects due to high speed flow and is not strongly dependent on pressure ratio across the tip gap contrary to published observations that have primarily dealt with subsonic tip flows. The high relative Mach numbers in the tip gap lead to a choking of the leakage flow that translates to a relative attenuation of losses at higher loading. The efficacy of new tip geometry is discussed to minimize heat flux at the tip while maintaining choked conditions. In addition, an explanation is provided that shows the mechanism behind the rise in stagnation temperature on the casing to values above the absolute total temperature at the inlet. It is concluded that even in steady mode, work transfer to the near tip fluid occurs due to relative shearing by the casing. This is believed to be the first such explanation of the work transfer phenomenon in the open literature. The difference in pattern between steady and time-averaged heat flux at the hub is also explained.

Shyam, Vikram↗
NASA NTRS1978

SSV Generic OFT first stage ascent base convective heating environments

Space Shuttle Vehicle (SSV) base convective heating environments during the recirculated flow time segment of first stage ascent were determined for the Generic OFT trajectory and performance conditions. These environments are presented in the attachments as cold wall convective heat transfer coefficients for individual or groups of design points for all affected SSV base surfaces. Two gas recovery temperatures applicable to upper and lower base region surfaces are also included. The environments were based upon magnitudes and trends from analysis of data taken during Base Heating Tests IH-39, IH-75, and IH-83. The plume heating trajectory and performance data were supplied by Rockwell International for both the OFT-1 Generic Mission A and Mission C.

Bender, R. L., Jr.↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1987

Effect of element size on the solution accuracies of finite-element heat transfer and thermal stress analyses of space shuttle orbiter

The effect of element size on the solution accuracies of finite-element heat transfer and thermal stress analyses of space shuttle orbiter was investigated. Several structural performance and resizing (SPAR) thermal models and NASA structural analysis (NASTRAN) structural models were set up for the orbiter wing midspan bay 3. The thermal model was found to be the one that determines the limit of finite-element fineness because of the limitation of computational core space required for the radiation view factor calculations. The thermal stresses were found to be extremely sensitive to a slight variation of structural temperature distributions. The minimum degree of element fineness required for the thermal model to yield reasonably accurate solutions was established. The radiation view factor computation time was found to be insignificant compared with the total computer time required for the SPAR transient heat transfer analysis.

Ko, William L.↗
DOE OSTI2020

A microchannel shell-and-tube absorber for ammonia-water absorption

A shell-and-tube absorber with microscale geometric features that yields significant volumetric heat duties is investigated. A vertically-oriented baffled shell-and-tube heat exchanger with a tightly packed tube bundle is designed to serve as the absorber for a 10.5 kW nominal cooling capacity absorption chiller. The heat and mass transfer in this absorber is investigated experimentally at operating conditions spanning a wide range of working fluid and coupling fluid flow rates and ambient temperatures. An experimental facility is fabricated to closely simulate the operating conditions in an absorption refrigeration system. Absorption heat duties in the range of 5-15 kW are measured. In addition to operating the absorber in the baseline co-current flow configuration with premixed inlet, fluid ports are provided to enable operation in counter-current flow, where the dilute solution enters from the top and vapor enters from the bottom. The performance of these two flow configurations is compared. A preliminary heat and mass transfer model is developed to quantify the local temperatures, concentrations, and mass flow rates of the solution and vapor streams along the length of the absorber. The model predictions for total heat duty agree with the experimental results within 10% for almost all the data points, with an AAD of ~3.4%. Finally, the results of this study, along with the modeling framework, can guide further development of highly compact and efficient absorbers absorption systems.

42 ENGINEERING↗
DOE OSTI2020

Iron-Doped Sodium Vanadium Oxyflurophosphate Cathodes for Sodium-Ion Batteries—Electrochemical Characterization and In Situ Measurements of Heat Generation

Sodium-ion batteries (NaIBs) are increasingly being envisioned for grid-scale energy-storage systems because of cost advantages. However, implementation of this vision has been challenged by the low-energy densities delivered by most NaIB cathodes. Toward addressing this challenge, the authors report the synthesis and characterization of a new iron-doped Na 3 Fe 0.3 V 1.7 O(PO 4 ) 2 F 2 cathode using a novel facile hydrothermal route. The synthesized material was characterized using scanning electron microscopy, X-ray diffraction, and Mössbauer spectroscopy techniques. The obtained discharge capacity in the half-cell configuration lies from 119 to 125 to 130 mA h/g at C/10 while tested using three different electrolyte formulations, dimethyl carbonate–ethylene carbonate (EC)–propylene carbonate (PC), diethyl carbonate–EC, and EC–PC, respectively. The synthesized cathodes were also evaluated in full-cell configurations, which delivered an initial discharge capacity of 80 mA h/g with NaTi 2 (PO 4 ) 3 MWCNT as the anode. Ionic diffusivity and interfacial charge transfer kinetics were also evaluated as a function of temperature and sodium concentration, which revealed that electrochemical rate performances in this material were limited by charge-transfer kinetics. To understand the heat generation mechanism of the Na/Na 3 Fe 0.3 V 1.7 O(PO 4 ) 2 F 2 half-cell during charge and discharge processes, an electrochemical isothermal calorimetry measurement was carried out at different current rates for two different temperatures (25 and 45 °C). The results showed that the amount of heat generated was strongly affected by the operating charge/discharge state, C-rate, and temperature. Altogether, this work provides a new synthesis route for the development of iron-doped Na 3 Fe 0.3 V 1.7 O(PO 4 ) 2 F 2 -based high-performance sodium cathode materials aimed at providing a viable pathway for the development and deployment of large-scale energy-storage based on sodium battery systems.

25 ENERGY STORAGE↗
NASA NTRS1998

Conceptual Design Method Developed for Advanced Propulsion Nozzles

As part of a contract with the NASA Lewis Research Center, a simple, accurate method of predicting the performance characteristics of a nozzle design has been developed for use in conceptual design studies. The Nozzle Performance Analysis Code (NPAC) can predict the on- and off-design performance of axisymmetric or two-dimensional convergent and convergent-divergent nozzle geometries. NPAC accounts for the effects of overexpansion or underexpansion, flow divergence, wall friction, heat transfer, and small mass addition or loss across surfaces when the nozzle gross thrust and gross thrust coefficient are being computed. NPAC can be used to predict the performance of a given nozzle design or to develop a preliminary nozzle system design for subsequent analysis. The input required by NPAC consists of a simple geometry definition of the nozzle surfaces, the location of key nozzle stations (entrance, throat, exit), and the nozzle entrance flow properties. NPAC performs three analysis "passes" on the nozzle geometry. First, an isentropic control volume analysis is performed to determine the gross thrust and gross thrust coefficient of the nozzle. During the second analysis pass, the skin friction and heat transfer losses are computed. The third analysis pass couples the effects of wall shear and heat transfer with the initial internal nozzle flow solutions to produce a system of equations that is solved at steps along the nozzle geometry. Small mass additions or losses, such as those resulting from leakage or bleed flow, can be included in the model at specified geometric sections. A final correction is made to account for divergence losses that are incurred if the nozzle exit flow is not purely axial.

Nadell, Shari-Beth↗
NASA NTRS2015

High Thermal Conductivity NARloy-Z-Diamond Composite Liner for Advanced Rocket Engines

NARloy-Z (Cu-3Ag-0.5Zr) alloy is state-of-the-art combustion chamber liner material used in liquid propulsion engines such as the RS-68 and RS-25. The performance of future liquid propulsion systems can be improved significantly by increasing the heat transfer through the combustion chamber liner. Prior work1 done at NASA Marshall Space Flight Center (MSFC) has shown that the thermal conductivity of NARloy-Z alloy can be improved significantly by embedding high thermal conductivity diamond particles in the alloy matrix to form NARloy-Z-diamond composite (fig. 1). NARloy-Z-diamond composite containing 40vol% diamond showed 69% higher thermal conductivity than NARloy-Z. It is 24% lighter than NARloy-Z and hence the density normalized thermal conductivity is 120% better. These attributes will improve the performance and life of the advanced rocket engines significantly. The research work consists of (a) developing design properties (thermal and mechanical) of NARloy-Z-D composite, (b) fabrication of net shape subscale combustion chamber liner, and (c) hot-fire testing of the liner to test performance. Initially, NARloy-Z-D composite slabs were made using the Field Assisted Sintering Technology (FAST) for the purpose of determining design properties. In the next step, a cylindrical shape was fabricated to demonstrate feasibility (fig. 3). The liner consists of six cylinders which are sintered separately and then stacked and diffusion bonded to make the liner (fig. 4). The liner will be heat treated, finish-machined, and assembled into a combustion chamber and hot-fire tested in the MSFC test facility (TF 115) to determine perform.

Bhat, Biliyar↗
DOE OSTI2021

Evaluation of external tubular configurations for a high-temperature chloride molten salt solar receiver operating above 700°C

Next-generation concentrating solar power (CSP) tower technologies target operating temperatures exceeding 700°C to increase the thermal-to-electric conversion efficiency. Molten chloride salts are one possible alternative to current commercial molten nitrate salts to enable the higher operating temperature. This paper analyzes the predicted optical, thermal-fluids, and structural performance of traditional external tubular solar receiver configurations applied with a chloride salt heat transfer fluid (HTF) and inlet/outlet temperatures of 500°C/735°C, and considers sensitivity analysis and optimization relative to receiver sizing, tube sizing, number of panels, flow circuit configurations, and solar flux concentration under constraints on internal velocity, pressure drop, wall thickness, and required creep-fatigue lifetime. The high temperature conditions increase the significance of inelastic deformation mechanisms such as creep relative to that expected in commercial 565°C nitrate salt designs. High-temperature creep and creep-fatigue damage in the metal alloy tubes are the key factors that limit allowable solar flux concentration and achievable receiver thermal efficiency at the near-800°C wall temperature conditions. For a traditional external cylindrical receiver configuration, the design parameters and conditions capable of satisfying all constraints produced, at best, a design point receiver efficiency of 78.2%, or 80.5% when excluding receiver intercept efficiency. Variation in the optimal receiver performance relative to uncertainty in the binding maximum velocity, minimum wall thickness, and minimum lifetime constraints is presented.

14 SOLAR ENERGY↗
DOE OSTI2024

Material selection and manufacturing for high‐temperature heat exchangers: Review of state‐of‐the‐art development, opportunities, and challenges

Abstract Many energy systems demand heat transfer at high temperatures to keep up with high demand for power, so high‐temperature material that can perform and last under these harsh conditions is needed for heat exchangers. The engineering requirements for these high‐temperature heat exchanger material call for high thermal conductivity, high resistance to fracture, high resistance to creep deformation, environmental stability in environments associated with the application, and high modulus of elasticity while maintaining low cost to make and maintain. Naturally, ceramics are a good solution for this endeavor. In the past, high‐temperature heat exchangers made from ceramics have been used. We provide examples of ceramics in relevant heat exchange applications and provide motivation where additive manufacturing (AM) can improve efficiency. AM for the relevant material is under development, and we provide insight on the AM of ceramic materials and examples of AM heat exchangers keeping cost in mind. The motivation of the review paper is to provide a framework for material and manufacturing selection for high‐temperature heat exchangers for AM to keep up with the demand for better efficiency, better material, better manufacturing, and cost moving forward with AM technology in high‐temperature ceramic heat exchangers.

36 MATERIALS SCIENCE↗
NASA NTRS1976

TAP 1: A Finite Element Program for Steady-State Thermal Analysis of Convectively Cooled Structures

The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature dependent thermal parameters is performed using the Newton-Raphson iteration method. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. A companion plotting program for displaying the finite element model and predicted temperature distributions is presented. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗