Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Flow Reversal”

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.

282 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Enhancements to Linear Stability-Based, CFD-integrated Transition Prediction for High-Speed Flows

Combining linear stability calculations with computational fluid dynamics (CFD) simulations has great potential for the automated modeling of high-speed flows, especially when adequate information about the configuration and the disturbance environment is available. However, a significant impediment to the applicability of this technique is the lack of an efficient method to calculate the crucial amplification ratio corresponding to the onset of transition in hypersonic flows. This ratio, also known as the "transition N-factor," is dependent upon the freestream disturbance environment as well as the surface properties of the test article. In response to the need for an engineering solution to predict the transition N-factor within conventional hypersonic wind tunnels, this paper presents a data-driven correlation that expands the existing correlations from straight circular cones with a narrow range of half angles to a broader array of axisymmetric configurations. Furthermore, when tested against a chosen dataset that was not used in its calibration, the suggested correlation shows good predictive accuracy with an RMS error of only 6.9%. Although similar accuracy may also be achieved via existing correlations based on similar datasets, predictions based on the proposed correlation have the advantage of not requiring an extensive amount of configuration-specific data. Practical applications often have access to the input parameters needed for this correlation, such as the freestream disturbance intensity, Mach number, and body-based slenderness Reynolds number. Additionally, this correlation outperforms the traditional assumption of a constant N-factor, particularly for configurations with blunted nose geometries. The development of this correlation is grounded in an extensive dataset encompassing conical models with body half-angles varying between 5 degrees and 16 degrees, Mach numbers ranging from 5 to 14, and nosetip-based Reynolds numbers approaching the transition reversal limit for blunt-nosed cones.

CFD↗
NASA NTRS1994

Sensitivity of lag-damping correlations to structural and aerodynamic approximations of isolated experimental rotors in forward flight

The predictions of regressive lag-mode damping levels are correlated with the database of an isolated, soft-inplane, three-blade rotor operated untrimmed. The database was generated at the Army Aeroflightdynamics Directorate at Ames. The correlation covers a broad range of data, from near-zero thrust conditions in hover to high-thrust and highly stalled conditions in forward flight with advance ratio as high as 0.55 and shaft angle as high as 20 degrees. In the experimental rotor, the airfoil or blade portion has essentially uniform mass and stiffness distributions, but the root flexure has highly nonuniform mass and stiffness distributions. Accordingly, the structural approximations refer to four models of root-flexure-blade assembly. They range from a rigid flap-lag model to three elastic flap-lag-torsion models, which differ in modeling the root flexure. The three models of root-flexure are: three root springs in which the bending-torsion couplings are fully accounted for; a finite-length beam element with some average mass and stiffness distributions such that the fundamental frequencies match those of the experimental model; and accurate modal representation in which the actual mass and stiffness distributions of the experimental root-flexure-blade assembly are used in calculating the nonrotating mode shapes. The four models of root-flexure-blade assembly are referred to as the rigid flap-lag model, spring model, modified model and modal model. For each of these four models of the root-flexure-blade assembly, the predictions are based on the following five aerodynamic theories: ear theory, which accounts for large angle-of-attack and reverse-flow effects on lift, and has constant drag and pitching moment; quasisteady stall theory, which includes quasisteady stall lift, drag and pitching moment characteristics of the airfoil section, dynamics stall theory, which uses the ONERA dynamic stall models of lift, drag and pitching moment; dynamic wake theory, which is based on a finite-state three-dimensional wake model and includes all wake effects including both shed and trailing vorticity; and dynamics and wake theory, which combines both dynamic stall theory and dynamic wake theory and is a relatively complete aerodynamic representation.

Gaonkar, G. H.↗
NASA NTRS1997

Cabauw Experimental Results from the Project for Intercomparison of Land-Surface Parameterization Schemes

In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W/m in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (+/- 10 W/sq m). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models' neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of 30 W/sq m and 25 W/sq m, respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W/m for sensible beat flux and 10 W/m for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to ground water, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.

Chen, Tian Hong↗
NASA NTRS1994

Quiet-time intensifications along the poleward auroral boundary near midnight

Radar and optical measurements from Sondrestrom are combined with satellite and Goose Bay data in a study of the poleward edge of the nightside auroral oval during a quiet period. The B(sub y) and B(sub z) components of the interplanetary magnetic field were close to zero, and the B(sub x) component was approximately 8 nT for more than 24 hours. On a large scale, the convection and precipitation patterns remained almost constant during this period; on a small scale, however, the conditions were quite dynamic. At 10- to 20-min intervals the arc that marked the poleward auroral boundary intensified, and a new arc appeared poleward of it. About once per hour, stronger intensifications were observed. One such event is examined in detail. The auroral arcs first appeared to dim, and then they brightened, with a factor of 10 increase in E region electron density. At the time of the brightening a new arc formed poleward of all the arcs. The arcs then drifted southward at velocities of approximately 270 m/s. A plasma drift disturbance, characterized by a doubling of the southward velocity and a reversal in the east-west component, propagated westward at 900 m/s through the fields of view of the Sondrestrom and Goose Bay radars. A simultaneous satellite overpass close to the radars revealed the presence of an energetic ion event similar to the 'velocity dispersed ion structures' observed on the Aureol satellite and presumed to be the signature of fast ion beams within the plasma sheet boundary layer. The stronger arc intensification events observed by the Sondrestrom radar are associated with an increase in plasma flow across the boundary between open and closed magnetic field lines. We interpret this increased flow as the ionospheric signature of abrupt, localized increases in the reconnection rate in the midnight sector.

De La Beaujardiere, O.↗
NASA NTRS2008

How Does Mediterranean Basin's Atmosphere Become Weak Moisture Source During Negative Phase of NAO: Use of AIRS, AMSR, TOVS, & TRMM Satellite Datasets Over Last Two NAO Cycles to Examine Governing Controls on E-P

The Mediterranean Sea is a noted 'concentration" basin in that it almost continuously exhibits positive evaporation minus precipitation (E - P ) properties -- throughout the four seasons and from one year to the next. Nonetheless, according to the ECMWF Era-40 48-year (1958-2005) climate reanalysis dataset, for various phases of the North Atlantic Oscillation (NAO) when the pressure gradient between Portugal and Iceland becomes either very relaxed (large negative NAO-Index) or in transition (small positive or negative NAO-Index), the atmospheric moisture source properties of the basin become weak, at times even reversed for several months (i.e., negative E - P). This behavior poses numerous questions concerning how and why these events occur. Moreover, it begs the question of what it would take for the basin to reach its tipping point in which P would exceed E throughout the rainy season (some six months) on an annually persistent basis -- and the sea would possibly transform to a recurring "dilution" basin. This talk investigates these questions by: (1) establishing over a period from 1979 to present, based on detailed analyses of satellite retrieval products from a combination of NASA-AQUA, NOAA-LEO, NASA/JAXA Scatterometer, and NASA-TRMM platforms, plus additional specialized satellite data products and ancillary meteorological datasets, the actual observation-based behavior of E - P, (2) diagnosing the salient physical and meteorological mechanisms that lead to the weaker E - P events during the analysis period, partly based on analyzing surface and upper air data at discrete stations in the western and eastern Mediterranean -- while at the same time evaluating the quality of the ERA-40 data over this same time period, (3) conducting GCM and high-resolution regional modeling experiments to determine if perturbed but realistic meteorological background conditions could maintain Mediterranean as a "dilution" basin through the October to March rainy season on annually recurring basis, and (4) investigating how such conditions might modify important internal and external climatic processes known to be closely related to the dynamical, thermal, and hydrological properties of the basin (e.g., drought frequency over Iberian peninsula, rainfall accumulation within Sahel, alteration of Levantine branch of east-west aligned open thermohaline cell, and modification of warm-salty intermediate flow through Gibraltar straight into North Atlantic).

Smith, Eric A.↗
NASA NTRS2012

Photogrammetric and Global Positioning System Measurements of Active Pahoehoe Lava Lobe Emplacement on Kilauea, Hawaii

Basalt is the most common rock type on the surface of terrestrial bodies throughout the solar system and -- by total volume and areal coverage -- pahoehoe flows are the most abundant form of basaltic lava in subaerial and submarine environments on Earth. A detailed understanding of pahoehoe emplacement processes is necessary for developing accurate models of flow field development, assessing hazards associated with active lava flows, and interpreting the significance of lava flow morphology on Earth and other planetary bodies. Here, we examine the active emplacement of pahoehoe lobes along the margins of the Hook Flow from Pu'u 'O'o on Kilauea, Hawaii. Topographic data were acquired between 21 and 23 February 2006 using stereo-imaging and differential global positing system (DGPS) measurements. During this time, the average discharge rate for the Hook Flow was 0.01-0.05 cubic m/s. Using stereogrammetric point clouds and interpolated digital terrain models (DTMs), active flow fronts were digitized at 1 minute intervals. These areal spreading maps show that the lava lobe grew by a series of breakouts tha t broadly fit into two categories: narrow (0.2-0.6 m-wide) toes that grew preferentially down-slope, and broad (1.4-3.5 m-wide) breakouts that formed along the sides of the lobe, nearly perpendicular to the down-flow axis. These lobes inflated to half of their final thickness within approx 5 minutes, with a rate of inflation that generally deceased with time. Through a combination of down-slope and cross-slope breakouts, lobes developed a parabolic cross-sectional shape within tens of minutes. We also observed that while the average local discharge rate for the lobe was generally constant at 0.0064 +/- 0.0019 cubic m/s, there was a 2 to 6 fold increase in the areal coverage rate every 4.1 +/- 0.6 minutes. We attribute this periodicity to the time required for the dynamic pressurization of the liquid core of the lava lobe to exceed the cooling-induced strength of the lobe margins. Using DGPS-derived DTMs of the topography before and after pahoehoe lobe emplacement, we observed that the lava typically concentrated within existing topographic lows, with the lobe reaching a maximum thickness of approx 1.2 m above the lowest points of the initial topography and above reverse-facing slopes. Lobe margins were typically controlled by high-standing topography, with the zone directly adjacent to the final flow margin having average relief that is approx 4 cm higher than the lava-inundated region. This suggests that irregularities approx 25% of the height of the smallest breakout elements (i.e., toes) can exert a strong control on the paths of low-discharge pahoehoe lobes, with stagnated toes forming confining margins that allow interior portions of flow to topographically invert the landscape by inflation.

Hamilton, Christopher W.↗
NASA NTRS2004

Turbulent Jet Flames Into a Vitiated Coflow

Examined is the vitiated coflow flame, an experimental condition that decouples the combustion processes of flows found in practical combustors from the associated recirculating fluid mechanics. The configuration consists of a 4.57 mm diameter fuel jet into a coaxial flow of hot combustion products from a lean premixed flame. The 210 mm diameter coflow isolates the jet flame from the cool ambient, providing a hot environment similar to the operating conditions of advanced combustors; this important high temperature element is lacking in the traditional laboratory experiments of jet flames into cool (room) air. A family of flows of increasing complexity is presented: 1) nonreacting flow, 2) all hydrogen flame (fuel jet and premixed coflow), and 3) set of methane flames. This sequence of experiments provides a convenient ordering of validation data for combustion models. Laser Raman-Rayleigh-LIF diagnostics at the Turbulent Diffusion Flame laboratory of Sandia National Laboratories produced instantaneous multiscalar point measurements. These results attest to the attractive features of the vitiated coflow burner and the well-defined boundary conditions provided by the coflow. The coflow is uniform and steady, isolating the jet flame from the laboratory air for a downstream distance ranging from z/d = 50-70. The statistical results show that differential diffusion effects in this highly turbulent flow are negligible. Complementing the comprehensive set of multiscalar measurements is a parametric study of lifted methane flames that was conducted to analyze flame sensitivity to jet and coflow velocity, as well as coflow temperature. The linear relationship found between the lift-off height and the jet velocity is consistent with previous experiments. New linear sensitivities were found correlating the lift-off height to coflow velocity and temperature. A blow-off study revealed that the methane flame blows off at a common coflow temperature (1260 K), regardless of coflow or jet velocity. An explanation for this phenomenon is that entrainment of ambient air at the high lift-off heights prevents autoignition. Analysis of the results suggests that flame stabilization occurs through a combination of flame propagation, autoignition, and localized extinction processes. Proposed is an expanded view of distributed reaction combustion based on analysis of the distributions of probe volume conditions at the stabilization region of the lifted hydrogen and methane flames. Turbulent eddies the size of the flame thickness mix fuel and hot coflow across the flame front, thereby enhancing the reaction zone with autoignition of reactants at elevated temperatures; this is the reverse effect of turbulent flames in ambient air, where intense turbulence in cool mixtures result in localized extinction. Each of the three processes (i.e., flame propagation, autoignition and localized extinction) contributes to flame stabilization in varying degrees, depending on flow conditions.

Holdeman, James D.↗
NASA NTRS1993

Numerical Study of Base Pressure Characteristic Curve for a Four-Engine Clustered Nozzle Configuration

Excessive base heating has been a problem for many launch vehicles. For certain designs such as the direct dump of turbine exhaust in the nozzle section and at the nozzle lip of the Space Transportation Systems Engine (STME), the potential burning of the turbine exhaust in the base region has caused tremendous concern. Two conventional approaches have been considered for predicting the base environment: (1) empirical approach, and (2) experimental approach. The empirical approach uses a combination of data correlations and semi-theoretical calculations. It works best for linear problems, simple physics and geometry. However, it is highly suspicious when complex geometry and flow physics are involved, especially when the subject is out of historical database. The experimental approach is often used to establish database for engineering analysis. However, it is qualitative at best for base flow problems. Other criticisms include the inability to simulate forebody boundary layer correctly, the interference effect from tunnel walls, and the inability to scale all pertinent parameters. Furthermore, there is a contention that the information extrapolated from subscale tests with combustion is not conservative. One potential alternative to the conventional methods is computational fluid dynamics (CFD), which has none of the above restrictions and is becoming more feasible due to maturing algorithms and advancing computer technology. It provides more details of the flowfield and is only limited by computer resources. However, it has its share of criticisms as a predictive tool for base environment. One major concern is that CFD has not been extensively tested for base flow problems. It is therefore imperative that CFD be assessed and benchmarked satisfactorily for base flows. In this study, the turbulent base flowfield of a experimental investigation for a four-engine clustered nozzle is numerically benchmarked using a pressure based CFD method. Since the cold air was the medium, accurate prediction of the base pressure distributions at high altitudes is the primary goal. Other factors which may influence the numerical results such as the effects of grid density, turbulence model, differencing scheme, and boundary conditions are also being addressed. Preliminary results of the computed base pressure agreed reasonably well with that of the measurement. Basic base flow features such as the reverse jet, wall jet, recompression shock, and static pressure field in plane of impingement have been captured.

Wang, Ten-See↗
NASA NTRS1996

Identification of Ser-543 as the major regulatory phosphorylation site in spinach leaf nitrate reductase

Spinach leaf NADH:nitrate reductase (NR) responds to light/dark signals and photosynthetic activity in part as a result of rapid regulation by reversible protein phosphorylation. We have identified the major regulatory phosphorylation site as Ser-543, which is located in the hinge 1 region connecting the cytochrome b domain with the molybdenum-pterin cofactor binding domain of NR, using recombinant NR fragments containing or lacking the phosphorylation site sequence. Studies with NR partial reactions indicated that the block in electron flow caused by phosphorylation also could be localized to the hinge 1 region. A synthetic peptide (NR6) based on the phosphorylation site sequence was phosphorylated readily by NR kinase (NRk) in vitro. NR6 kinase activity tracked the ATP-dependent inactivation of NR during several chromatographic steps and completely inhibited inactivation/phosphorylation of native NR in vitro. Two forms of NRk were resolved by using anion exchange chromatography. Studies with synthetic peptide analogs indicated that both forms of NRk had similar specificity determinants, requiring a basic residue at P-3 (i.e., three amino acids N-terminal to the phosphorylated serine) and a hydrophobic residue at P-5. Both forms are strictly calcium dependent but belong to distinct families of protein kinases because they are distinct immunochemically.

NASA Discipline Plant Biology↗
NASA NTRS2003

Molecular and functional analyses of the contractile apparatus in lymphatic muscle

Lymphatics are necessary for the generation and regulation of lymph flow. Lymphatics use phasic contractions and extrinsic compressions to generate flow; tonic contractions alter resistance. Lymphatic muscle exhibits important differences from typical vascular smooth muscle. In this study, the thoracic duct exhibited significant functional differences from mesenteric lymphatics. To understand the molecular basis for these differences, we examined the profiles of contractile proteins and their messages in mesenteric lymphatics, thoracic duct, and arterioles. Results demonstrated that mesenteric lymphatics express only SMB smooth muscle myosin heavy chain (SM-MHC), whereas thoracic duct and arterioles expressed both SMA and SMB isoforms. Both SM1 and SM2 isoforms of SM-MHC were detected in arterioles and mesenteric and thoracic lymphatics. In addition, the fetal cardiac/skeletal slow-twitch muscle-specific beta-MHC message was detected only in mesenteric lymphatics. All four actin messages, cardiac alpha-actin, vascular alpha-actin, enteric gamma-actin, and skeletal alpha-actin, were present in both mesenteric lymphatics and arterioles. However, in thoracic duct, predominantly cardiac alpha-actin and vascular alpha-actin were found. Western blot and immunohistochemical analyses corroborated the mRNA studies. However, in arterioles only vascular alpha-actin protein was detected. These data indicate that lymphatics display genotypic and phenotypic characteristics of vascular, cardiac, and visceral myocytes, which are needed to fulfill the unique roles of the lymphatic system.

Non-NASA Center↗
NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
NASA NTRSDate not supplied

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↗