Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “maximum heat transfer coefficient”

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.

101 records · Page 6

Separation of volatile organic contaminants from water using a direct-contact dehumidifier: An experimental study and modeling

A packed-bed, direct-contact dehumidifier (PB-DCD) has been developed to selectively condense water from a gaseous mixture of organic contaminants. In this process, hot air is humidified by mixing with wastewater vapor in a static mixer and is dehumidified by subcooled fresh water passing through the PB-DCD in a counter-current configuration. The goal of this study is to understand the relationship between various operating parameters and how they affect separation of isopropyl alcohol (IPA) from the gaseous mixture with the purpose of expanding the investigation to other contaminants in the future. The process has been successfully used to separate IPA from wastewater. Different parameters including the concentration of IPA in water, feed temperature, air-to-vapor ratio, and cooling-water-to-vapor mass ratio, as well as residence time were examined to determine their influence on contaminant separation and clean water recovery rate. The Buckingham Pi theorem was applied to reduce the number of studying variables by generating non-dimensional groups. Experiments were performed to determine the relationship between the non-dimensional parameters and model mass transfer in the system. Mass transfer coefficients for multicomponent gaseous mixture condensation in a packed-bed direct contact condenser column are calculated and a modified correlation of mass transfer is developed. The parameters of the mass transfer model are obtained from operating conditions of the PB-DCD, stream conditions and thermodynamic properties. A modified Sherwood correlation is developed to predict IPA separation from a gaseous mixture consisting of air, water, and IPA. The modified Sherwood correlation is applied to successfully predict the Sherwood number with a mean absolute error equal to 6.6 %. According to the experimental results, the maximum IPA separation achieved was 76.3 %, which corresponds to a water recovery rate of 58.7 %.

42 ENGINEERING↗

Porous mesh manifold for enhanced boiling performance

High-performance electronics are continuously demanding cooling of higher heat fluxes. Phase-change cooling, including pool boiling, is a useful approach to address this challenge; however, competition between liquid and vapor flows generally limit the heat fluxes that can be dissipated. A range of strategies to control these flows have been investigated previously, including capillary guides. Here a manifold structure formed from a metallic mesh is investigated to control the disposition of liquid and vapor phases above a pool fed boiling surface enhanced with porous structures. Copper mesh forms defined liquid flow paths, using capillary action to guide and distribute liquid evenly over the heated surface, along with open channels to facilitate vapor escape. The mesh provides a novel structure for liquid guidance that imposes low resistance to liquid flow while occluding a minimal area of heated surface underneath. The manifold performance is characterized in boiling fed by a pool of water above a laser-textured aluminum nitride heat dissipation surface with pin–fin structures having heights of 110 µm and spacing of 30 µm with a heated area of 5 mm x 5 mm. A maximum heat flux of 490 W/cm 2 is reached with the manifold in the pool fed configuration, representing an increase of more than 65% over the porous pin fin surface alone. The maximum stable superheat observed for the manifold of 36K is 14K higher than that for the porous surface without the manifold. The factors limiting performance of the manifold are analyzed. High superheat is attributed to partial flooding of the boiling surface as suggested by the reduction in superheat using external suction. Similar systems and structures for enhanced two-phase cooling are compared.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tribological properties and surface chemistry of silicon carbide at temperatures to 1500 C

Silicon carbide surfaces were heated to 1500 C in a vacuum and analyzed at room temperature with X-ray photoelectron spectroscopy (XPS) and Auger electron spectroscopy (AES). The basic unit of the surfaces was considered as a plane of a tetrahedron of either SiC4 and CSi4 composition. AES spectra were obtained from 250-1500 C, with an analysis depth of 1 nm revealed the presence of little Si and mostly graphite. XPS analysis depth was 2 nm or less, and Si was found in the second 1 nm. Sliding friction tests with single-crystal silicon carbide in contact with iron in a vacuum were characterized by a stock-slip value. The coefficient of friction increased with increasing temperature up to 400 C, then decreased with increasing temperature from 400-600 C. Reheating surfaces to 800 C after preheating them to that temperature produced no changes in AES readings. It is concluded that the maximum density of silicon and silicon-carbide is at 800 C, and the higher the sliding temperature, the more metal that is transferred.

Miyoshi, K.↗

Numerical simulation of falling film sensible heat transfer over round horizontal tubes

Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.

42 ENGINEERING↗

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗

An Analysis of the Corridor and Guidance Requirements for Supercircular Entry into Planetary Atmospheres

An analysis is presented of supercircular entry into a planet's atmosphere giving particular attention to the corridor through which spacecraft must be guided in order to accomplish various maneuvers. A dimensionless parameter based on conditions at the conic perigee altitude is introduced for characterizing supercircular entries and conveniently pre-scribing corridor widths associated with elliptic, parabolic, or hyperbolic approach trajectories. The analysis applies to vehicles of arbitrary weight, shape, and size. Illustrative calculations are made for Venus, Earth, Mars, Jupiter, and Titan. For nonlifting vehicles having fixed aerodynamic coefficients, curves are presented of dimensionless parameters from which can be calculated the maximum deceleration, maximum rate of laminar convective heating, and total laminar heat absorbed during single-pass entry at velocities up to twice circular velocity. For lifting vehicles, curves are presented of the maximum deceleration and overshoot boundary of an entry corridor; equations are presented for estimating laminar aerodynamic heating from the maximum deceleration. It is shown that the corridor width is independent of vehicle weight, dimensions, and drag coefficient, provided these are the same at the overshoot boundary as at undershoot. The corridors of certain planets can be broadened markedly by the application of aerodynamic lift; for example, the 10-earth-g corridor width for single- pass, nonlifting, parabolic entry is increased from 0 miles for Jupiter, 7 for Earth, and 8 for Venus, to 52, 51, and 52 miles, respectively, by employing a lift-drag ratio of 1. The use of aerodynamic lift does not increase appreciably the corridors of Mars and Titan. All corridor widths decrease rapidly as the entry velocity is increased. Terminal guidance requirements on accuracy of velocity and flight path angle for successfully entering various corridors are compared with analogous requirements for putting a satellite into orbit, for hitting the moon from the earth, and for achieving ICBM accuracy. Consideration is given to the terminal guidance problem involved in using a planet's atmosphere-rather than rocket fuel-to effect orbital transfers from heliocentric to planetocentric motion, thereby converting a hyperbolic approach trajectory to an elliptic orbit about the target planet. This fuel saving maneuver appears technologically feasible for certain planetary voyages, and implies the possibility of achieving a large reduction in required Earth lift-off weight of chemical propulsion systems.

Chapman, Dean R.↗

Investigation of Heat Transfer from a Stationary and Rotating Conical Forebody

The convective heat transfer from the surface of a conical forebody having a hemispherical nose, an included angle of approximately 30 deg, and. a maximum diameter of 18.9 inches was investigated in a wind tunnel for both stationary and. rotating operation. The range of test conditions included free-stream velocities up to 400 feet per second, rotational speeds up to 1200 rpm, and. angles of attack of 0 deg and 6 deg. Both a uniform surface temperature and a uniform heater input power density were used. The Nusselt-Reynolds number relations provided good correlation of the heat-transfer data for the complete operating range at 0 deg angle of attack with and without spinner rotation, and for 6deg angle of attack with rotation. Rotational speeds up to 1200 rpm had no apparent effect on the heat-transfer characteristics of the spinner. The results obtained at 6 deg angle of attack with rotation were essentially the same as those obtained at 0 deg angle of attack without rotation. The experimental heat-transfer characteristics in the turbulent flow region were consistently in closer agreement with the results predicted for a two-dimensional body than with those predicted. for a cone. For stationary operation at 60 angle of attack, the measured heat-transfer coefficients in the turbulent flow region were from 6 to 13 percent greater on the lower surface (windward. side) than on the upper surface (sheltered side) for corresponding surface locations. The spinner-nose geometry appeared to cause early boundary-layer transition. Transition was initiated at a fairly constant Reynolds number (based on surface distance from nose) of 8.0 x 10(exp 4). Transition was completed at Reynolds numbers less than 5.0 x 10(exp 5) for all conditions investigated.

Ruggeri, Robert S.↗

Heat pulse testing at monitoring wells to estimate subsurface fluid velocities in geological CO 2 storage

Monitoring the injected CO 2 during geological CO 2 storage (GCS) is essential to assure containment and identify CO 2 leakage. Here in this work, a new approach is introduced to estimate the evolution of the downhole fluid velocity at a monitoring well and identify CO 2 arrival time using in-well heat pulse/tracer test. The proposed technique involves using a downhole heater to generate a series of heat pulses and measuring their corresponding temperature response. The surface temperature of the downhole heater is controlled by the supplied electrical power and the heat loss by convection to the surroundings. Convective heat transfer is well described using Newton's law of cooling in which the temperature difference between the heater and the surrounding fluids drives the heat transfer, for which the convection heat transfer coefficient (h) controls the magnitude of heat loss. Among various factors that control h, it depends on the type of the flowing fluid and its velocity. Through analyzing the measured temperature at different heat pulses, the changes in h - due to mobilization of the in-situ brine or CO 2 arrival - can be estimated. Consequently, the velocity of the flowing fluid across the heater can be obtained. Since heat transfer by convection is sensitive to the type of the surrounding fluid, intrusion of CO 2 can be detected from the relatively higher surface temperature obtained at CO 2 arrival. Churchill and Bernstein (1977)'s correlation is adopted to estimate the change of fluid velocity in terms of the change in h. To demonstrate the validity of the proposed technique, the results are applied and validated against those of COMSOL Multiphysics simulation tool for single-phase brine (before CO 2 arrival) and single-phase CO 2 (after CO 2 arrival). The observed temperature heating is sensitive to the flowing fluid velocity and fluid type. The temperature signal observed at CO 2 arrival is large and easily detectable using temperature monitoring tool which provides reliable indication for tracking CO 2 arrival at monitoring wells compared with passive temperature monitoring. The results obtained using the proposed technique agree very well with the numerical results obtained from the simulation tool with a maximum estimation error of 7 percent.

02 PETROLEUM↗

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗