Search NASA⌕ Search

SEARCH · Search NASA

Results for “pattern separation”

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.

1,635 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling Efforts to Gain Insight into Historical Leak Events from the Single-Shell Tank A-105 in the Hanford 241-A Tank Farm - 20113

The 241-A Tank Farm is a single-shell tank (SST) farm constructed to store process waste from Hanford nuclear operations. Millions of gallons of nuclear waste were stored in the 241-A Tank Farm SSTs and some of the SSTs leaked in the past. In addition, spills and pipeline leaks during transfers and storage and intentional discharges to cribs and trenches resulted in releasing waste to the ground. Tank liner leaks are referred to as 'leaks' and all other discharges to the soil are referred to as 'releases.' Liquid waste that could be removed by pumping has been removed from all of the SSTs to reduce the potential for future leaks. The most significant historical leak event in 241-A Tank Farm occurred in 1965 when tank A-105 experienced a pressurized steam event. This event resulted in damage to the inner steel liner flooring of A-105, which was separated from the sidewalls over part of its circumference and buckled up to 2.5 m, vertically. During this event up to 7,570 L of contaminated liquid may have leaked to the adjacent soil. This leak occurred under extreme temperature and pressure conditions. Modeling efforts to better understand the nature and extent of historical leaks from subsurface SSTs that were used to store highly radioactive, self-boiling liquid wastes at 241-A Tank Farm have recently been undertaken. These efforts have compiled the available historical information to support a preliminary non-isothermal, multiphase flow and chemical transport modeling effort to re-create the conditions under which the leaks occurred and to formulate a conceptual model as to the extent and distribution of leaked radioactive contaminants in the adjacent soil material. The focus of this work was on developing an understanding of the key features, processes and bounding conditions related to tank A-105 leak events that occurred in the 1960's. The activities include estimation of leak composition, non-isothermal multiphase flow and transport modeling, and geochemical modeling. The STOMP{sup C} Water-Air-Energy modeling code was utilized to implement the three-dimensional representation of the subsurface tank and the surrounding flow field. Once constructed, the model was used to evaluate the conceptual understanding of leaks originating from different parts of the tank under the elevated temperature boundary conditions that were imposed by surrounding tanks at the time of the historical leaks. The chemical evolution of the liquid waste was also evaluated as it leaked from the tank under transient pressure and temperature gradients. The results of the preliminary evaluation suggested the presence of a heat-pipe effect beneath tank A-105 in which water vapor at an elevated temperature is driven away from the base of the tank to a position where water vapor cools, condenses and is then drawn back toward the tank by the strong capillary attraction of the dry soil. A non-sorbing contaminant (Tc-99) was introduced into this flow field to better understand the potential distribution patterns of leaked contaminants that may have occurred during the historical leak events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2003

High-Performance Acousto-Ultrasonic Scan System Being Developed

Acousto-ultrasonic (AU) interrogation is a single-sided nondestructive evaluation (NDE) technique employing separated sending and receiving transducers. It is used for assessing the microstructural condition and distributed damage state of the material between the transducers. AU is complementary to more traditional NDE methods, such as ultrasonic cscan, x-ray radiography, and thermographic inspection, which tend to be used primarily for discrete flaw detection. Throughout its history, AU has been used to inspect polymer matrix composites, metal matrix composites, ceramic matrix composites, and even monolithic metallic materials. The development of a high-performance automated AU scan system for characterizing within-sample microstructural and property homogeneity is currently in a prototype stage at NASA. This year, essential AU technology was reviewed. In addition, the basic hardware and software configuration for the scanner was developed, and preliminary results with the system were described. Mechanical and environmental loads applied to composite materials can cause distributed damage (as well as discrete defects) that plays a significant role in the degradation of physical properties. Such damage includes fiber/matrix debonding (interface failure), matrix microcracking, and fiber fracture and buckling. Investigations at the NASA Glenn Research Center have shown that traditional NDE scan inspection methods such as ultrasonic c-scan, x-ray imaging, and thermographic imaging tend to be more suited to discrete defect detection rather than the characterization of accumulated distributed micro-damage in composites. Since AU is focused on assessing the distributed micro-damage state of the material in between the sending and receiving transducers, it has proven to be quite suitable for assessing the relative composite material state. One major success story at Glenn with AU measurements has been the correlation between the ultrasonic decay rate obtained during AU inspection and the mechanical modulus (stiffness) seen during fatigue experiments with silicon carbide/silicon carbide (SiC/SiC) ceramic matrix composite samples. As shown in the figure, ultrasonic decay increased as the modulus decreased for the ceramic matrix composite tensile fatigue samples. The likely microstructural reason for the decrease in modulus (and increase in ultrasonic decay) is the matrix microcracking that commonly occurs during fatigue testing of these materials. Ultrasonic decay has shown the capability to track the pattern of transverse cracking and fiber breakage in these composites.

Roth, Don J.↗
NASA NTRS2000

Thermocapillary Flow and Aggregation of Bubbles on a Solid Wall

During the electrolytic evolution of oxygen bubbles forming on a vertically oriented transparent tin oxide electrode, bubbles were found to be mutually attractive. The mechanism of the aggregation had never been explained satisfactorily until Guelcher et al. attributed it to thermocapillary flow. The gradient of surface tension of the liquid at the bubble's surface, which was established because of reaction heat and ohmic heat loss at the electrode wall, drives flow of the liquid adjacent to each bubble; the bubble "pumps" fluid along its surface away from the wall. Fluid flows toward the bubble to conserve mass and entrains nearby bubbles in the flow pattern. The same logic would apply when two bubbles of equal size are adjacent to each other on a warm wall. Each bubble drives thermocapillary flow and hence entrains the other in its flow pattern, which drives the aggregation. Our objective here is to perform experiments where the temperature gradient at the wall is well known and controlled. The theory can be quantitatively tested by studying aggregation of bubble pairs of equal size, and by varying system parameters such as temperature gradient, bubble size and fluid viscosity. The results are then compared with the theory in a quantitatively rigorous manner. We demonstrate that the theory without adjustable parameters is capable of quantitatively modeling the rate of aggregation of two bubbles. The equations governing the thermocapillary flow around a single stationary bubble on a heated or cooled wall in a semi-infinite domain were solved. Both Reynolds number and Marangoni number were much less than unity. The critical result is that liquid in the vicinity of a warm wall flows toward a stationary collector bubble. Consequently the thermocapillary flow around the stationary bubble entrains another bubble toward itself. The bubbles undergo hindered translation parallel to the wall with velocity U while the fluid flow field is described with u. Two velocities were equated by using a wall hindrance parameter q: U = qu [1] which shows the velocity of bubble is proportional to the entraining velocity. The hindrance parameter q can experimentally be measured independently. q can also be calculated by solving the equations of motion for a bubble translating parallel to a solid wall. The experimental cell is cylindrical with an ID of 10 cm and consists of a 1 cm deep main cell filled with silicone oil and flanked by two thermal reservoirs. The upper thermal reservoir was heated and the lower thermal reservoir was cooled so that the bubbles aggregate. Two types of silicone oil (eta = 0.02 and 0.50 Pa s) were used. Two equal sized air bubbles were injected into the cell with a syringe. The center-to-center distance of bubbles was observed through a microscope. Bubble radius ranged from 0.40 mm to 0.65 mm and the temperature gradients along with the cell ranged from 1400 to 5000 K/m. The bubbles aggregated when heat flows from the wall to the fluid. The velocities of bubbles were in the range of 1 - 10 microns/s. The separation r decreased more quickly when the temperature gradient was higher, bubble size was larger, and the oil viscosity was lower. r decreased more rapidly as the bubbles approached each other. Dimensionless time was arbitrarily set to be zero when the dimensionless center-to-center distance between the bubbles was 4. All the bubble trajectories fall onto one line, especially in the range of dimensionless distance from 4 to 3. This means the relative movement of the bubble pair is proportional to the temperature gradient and bubble size and it is inversely proportional to the viscosity of the oil. This result strongly suggests that the thermocapillary flow-based aggregation mechanism is correct. A value of q can be estimated by fitting the scaled data to Eq. [1]. A best fit value of q was obtained as q = 0.26 with a standard deviation of 0.03. Independent experimental results for q for a 0.5 mm radius bubble, give values of q in the range 0.11 to 0.23. The value of q obtained from solving the equations of motion reveals q has values in the range 0.23 - 0.30. Since the full scale of possible values of q is zero to one, the maximum deviation of independently determined values of q from the best fit value was 15% of this full scale. Thus reasonable quantitative agreement between theory and experiment has been obtained.

Kasumi, Hiroki↗
NASA NTRS1999

Carbon Fiber Reinforced/Silicon Carbide Turbine Blisk Testing in the SIMPLEX Turbopump

A program designed to implement a ceramic matrix composite integrally bladed disk (blisk) into rocket engine style turbomachinery has successfully completed testing. The Marshall Space Flight Center (MSFC) program, utilizing the MSFC turbomachinery design, analysis, and testing capabilities along with materials development capabilities from both Glenn Research Center (GRC) and MSFC, has tested two carbon fiber reinforced silicon carbide blisks in the Simplex Turbopump at MSFC's Test Stand 500. One blisk contained a polar woven fiber preform, while the second blisk tested utilized a quasi-isotropic preform. Vhile earlier papers have chronicled the program's design, material testing, and torque testing efforts, this paper focuses on the testing of the blisks in the Simplex turbopump. Emphasis will be placed on the actual condition of the blisks before and after the testing test program design methodology, and conclusions that can be drawn from the test data and blisk final conditions. The program performed three separate test series. The first series was needed to validate that the Simplex turbopump was correctly re-built following a major incident to the turbopump. The turbopump had two major differences from the original design. The most obvious difference was the sleeve required throughout the bore of the main housing. The second major difference was modifications to the pump diffuser to improve performance. Several areas were burnt during the incident and were either repaired by weld repair (pump inlet housing) or simply smoothed out (turbine nozzle discharge). The test series was designed to weed out any turbopump design and manufacturing flaws or fatigue issues prior to putting the C/SiC blisks into it. The second and third series were the C/SiC blisk test series. The primary goal of these series was to expose the blisks to as much fatigue causing dynamic stress as possible to examine the material's capability. Initially, the test plan was to put equal time on the two blisks, however, as the test series progressed, the funding allowed additional testing to occur. The additional test time was placed on the polar weave blisk. The total test time accrued on the polar blisk was 2550 seconds with 860 seconds near the turbopump design speed of 25,000 rpm. This testing included 6 tests / 775 seconds pumping liquid nitrogen and 7 tests / 1775 seconds pumping liquid oxygen. The drive gas for all of the tests was gaseous nitrogen due to the lack of hot gas source for the Simplex turbopump. The quasi-isotropic blisk was tested for XX total tests and XXXX seconds with X tests/XXXX seconds pumping liquid nitrogen and X tests/XXXX seconds pumping liquid \oxygen.During the test series, the blisks were inspected following each test. Inspections initially were viewed from the downstream side of the blisks only. Midway through the testing, a method of borescoping the leading edges of the blades was devised, and subsequently, both sides of the blades were inspected following each test. The leading and trailing edges of the polar blisk held up better than the quasi-isotropic blisk. This was a known possibility due to the varying fiber direction in the blades as the rectangular preform weave is cut in a circular pattern. The surprising fact about the testing was that there was no measurable performance loss due to the inaccuracies in the blade manufacturing in the C/SiC blisks, the surface roughness C/SiC of the blades, or the loss of the material in the polar blisk. A performance shift was seen in the quasi-isotropic blisk as portions of the leading and trailing edges were lost. After the testing was completed, detailed inspections of the blisks were performed. The largest surprise was the polar blisk had a obvious crack in a single blade that was located nearly midspan which was not detected in test. The crack ran completely through the blade circumferenciary and through the radial length of the blade. However, the crack does not appear to extend into the blisk hub. Although the cause of the crack is still under investigation, the material appears to be tolerant of this crack, and other hairline cracks discovered under higher magnification. This bodies well for eventual use of this material in actual flight turbopumps where monolithic fracture toughness issues limit its use.

Genge, Gary G.↗
DOE OSTI2023

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗
NASA NTRS2010

Tool for Forecasting Cool-Season Peak Winds Across Kennedy Space Center and Cape Canaveral Air Force Station

The expected peak wind speed for the day is an important element in the daily morning forecast for ground and space launch operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). The 45th Weather Squadron (45 WS) must issue forecast advisories for KSC/CCAFS when they expect peak gusts for >= 25, >= 35, and >= 50 kt thresholds at any level from the surface to 300 ft. In Phase I of this task, the 45 WS tasked the Applied Meteorology Unit (AMU) to develop a cool-season (October - April) tool to help forecast the non-convective peak wind from the surface to 300 ft at KSC/CCAFS. During the warm season, these wind speeds are rarely exceeded except during convective winds or under the influence of tropical cyclones, for which other techniques are already in use. The tool used single and multiple linear regression equations to predict the peak wind from the morning sounding. The forecaster manually entered several observed sounding parameters into a Microsoft Excel graphical user interface (GUI), and then the tool displayed the forecast peak wind speed, average wind speed at the time of the peak wind, the timing of the peak wind and the probability the peak wind will meet or exceed 35, 50 and 60 kt. The 45 WS customers later dropped the requirement for >= 60 kt wind warnings. During Phase II of this task, the AMU expanded the period of record (POR) by six years to increase the number of observations used to create the forecast equations. A large number of possible predictors were evaluated from archived soundings, including inversion depth and strength, low-level wind shear, mixing height, temperature lapse rate and winds from the surface to 3000 ft. Each day in the POR was stratified in a number of ways, such as by low-level wind direction, synoptic weather pattern, precipitation and Bulk Richardson number. The most accurate Phase II equations were then selected for an independent verification. The Phase I and II forecast methods were compared using an independent verification data set. The two methods were compared to climatology, wind warnings and advisories issued by the 45 WS, and North American Mesoscale (NAM) model (MesoNAM) forecast winds. The performance of the Phase I and II methods were similar with respect to mean absolute error. Since the Phase I data were not stratified by precipitation, this method's peak wind forecasts had a large negative bias on days with precipitation and a small positive bias on days with no precipitation. Overall, the climatology methods performed the worst while the MesoNAM performed the best. Since the MesoNAM winds were the most accurate in the comparison, the final version of the tool was based on the MesoNAM winds. The probability the peak wind will meet or exceed the warning thresholds were based on the one standard deviation error bars from the linear regression. For example, the linear regression might forecast the most likely peak speed to be 35 kt and the error bars used to calculate that the probability of >= 25 kt = 76%, the probability of >= 35 kt = 50%, and the probability of >= 50 kt = 19%. The authors have not seen this application of linear regression error bars in any other meteorological applications. Although probability forecast tools should usually be developed with logistic regression, this technique could be easily generalized to any linear regression forecast tool to estimate the probability of exceeding any desired threshold . This could be useful for previously developed linear regression forecast tools or new forecast applications where statistical analysis software to perform logistic regression is not available. The tool was delivered in two formats - a Microsoft Excel GUI and a Tool Command Language/Tool Kit (Tcl/Tk) GUI in the Meteorological Interactive Data Display System (MIDDS). The Microsoft Excel GUI reads a MesoNAM text file containing hourly forecasts from 0 to 84 hours, from one model run (00 or 12 UTC). The GUI then displays e peak wind speed, average wind speed, and the probability the peak wind will meet or exceed the 25-, 35- and 50-kt thresholds. The user can display the Day-1 through Day-3 peak wind forecasts, and separate forecasts are made for precipitation and non-precipitation days. The MIDDS GUI uses data from the NAM and Global Forecast System (GFS), instead of the MesoNAM. It can display Day-1 and Day-2 forecasts using NAM data, and Day-1 through Day-5 forecasts using GFS data. The timing of the peak wind is not displayed, since the independent verification showed that none of the forecast methods performed significantly better than climatology. The forecaster should use the climatological timing of the peak wind (2248 UTC) as a first guess and then adjust it based on the movement of weather features.

Barrett, Joe H., III↗
NASA NTRS2015

Risk of Crew Adverse Health Event Due to Altered Immune Response

Determining the effect of space travel on the human immune system has proven to be extremely challenging. Limited opportunities for in-flight studies, varying mission durations, technical and logistical obstacles, small subject numbers, and a broad range of potential assays have contributed to this problem. Additionally, the inherent complexity of the immune system, with its vast array of cell populations, sub-populations, diverse regulatory molecules, and broad interactions with other physiological systems, makes determining precise variables to measure very difficult. There is also the challenge of determining the clinical significance of any observed immune alterations. Will such a change lead to disease, or is it a transient subclinical observation related to short-term stress? The effect of this problem may be observed by scanning publications associated with immunity and spaceflight, which began to appear during the 1970s. Although individually they are each valid studies, the comprehensive literature to date suffers from widely varying sampling methods and assay techniques, low subject counts, and sometimes a disparate focus on narrow aspects of immunity. The most clinically relevant data are derived from in-flight human studies, which have demonstrated altered cell-mediated immunity and reactivation of latent herpes viruses. Much more data are available from post-flight testing of humans, with clear evidence of altered cytokine production patterns, altered leukocyte distribution, continued latent viral reactivation, and evidence of dramatically altered virus-specific immunity. It is unknown if post-flight assessments relate to the in-flight condition or are a response to landing stress and readaptation. In-flight culture of cells has clearly demonstrated that immune cells are gravity-sensitive and display altered functional characteristics. It is unknown if these data are related to in vivo immune cell function or are an artifact of microgravity culture. Ground analog testing of humans and animals, as well as microgravity-analog cell culture, has demonstrated utility. However, in all cases, it is not known with certainty if these data would reflect similar testing during space travel. Given their ready availability, ground analogs may be extremely useful for assay development and the evaluation of potential countermeasures. In general, the evidence base suffers from widely disparate studies on small numbers of subjects that do not directly correlate well with each other or spaceflight itself. Also lacking are investigations of the effect of gender on adaption to spaceflight. This results in significant knowledge 'gaps' that must be filled by future studies to completely determine any clinical risk related to immunity for human exploration-class space missions. These gaps include a significant lack of in-flight data, particularly during long-duration space missions. The International Space Station represents an excellent science platform with which to address this knowledge gap. Other knowledge gaps include lack of a single validated ground analog for the phenomenon and a lack of flight-compatible laboratory equipment capable of monitoring astronauts (for either clinical or research purposes). However, enough significant data exist, as described in this manuscript, to warrant addressing this phenomenon during the utilization phase of the ISS. A recent Space Shuttle investigation has confirmed the 31 in-flight nature of immune dysregulation, demonstrating that it is not merely a post-flight phenomenon. Several current studies are ongoing onboard the ISS that should thoroughly characterize the phenomenon. NASA recognizes that if spaceflight-associated immune dysregulation persists during exploration flights in conjunction with other dangers, such as high-energy radiation, the result may be a significant clinical risk. This emphasizes the need for a continued integrated comprehensive approach to determining the effect of prolonged spaceflight, separated from transient launch and landing stresses, on human immunity. After such studies, the phenomenon will be understood, and, hopefully, a monitoring strategy will have been developed that could be used to monitor the effectiveness of countermeasure

Crucian, Brian↗
DOE OSTI2019

Flexible Geothermal Power Generation utilizing Geologic Thermal Energy Storage (Final Seedling Project Report)

This report concludes that there is a cost-effective strategy for seasonal storage of heat that will provide firm, but dispatchable, electrical generating capacity in times when other renewable energy is not available to meet demand. Deployment of the technology appears to require no new technology, but instead combines solar, geothermal, and conventional oil and gas drilling technologies in a novel way. The study basis is the use of sedimentary geologic formations as a medium for thermal energy storage (TES), specifically for heat collected in concentrating solar collectors. The study identifies methodologies that could be used to transport this heat into and out of the subsurface in order to produce dispatchable electrical power, and reports on initial optimization results. The GeoTES system (heat input, storage, heat recovery, and heat to electric conversion) described in this analysis has the potential to provide a unique pathway for increasing the grid penetration of renewable energy in large blocks of power and across many states and regions. Further, the system can be used both to meet the nation’s flexible energy needs while also improving grid stability and reliability. The present study evaluated the use of a large number of dedicated wells to store and recover the heat, essentially creating a synthetic geothermal reservoir. The use of sedimentary geology allows the wells to be drilled at low cost. Dedicated hot and cold wells are used, arranged in a five-spot well pattern with each hot or cold well surrounded at an appropriate distance by the opposite type of well. In large numbers this becomes alternating rows of hot and cold wells. Each hot and cold well is operated using a push-pull strategy. This configuration provides the ability to immediately recover stored hot fluid from a GeoTES reservoir, or to store the heat over many months for recovery at low loss when needed. This is a practical approach for managing the system’s fluid inventory and reducing parasitic load. The production and injection power requirements are reduced because the rows of wells operating in “push” mode provide help to the wells operating in “pull” mode, and vice-versa. Initial charging of a GeoTES system increases the heat recovery temperature. Increasing the duration of the charging period decreases the magnitude of the temperature fluctuations that occur following prolonged system operation. Because of direct contact of the heated water with the reservoir formation, the production of both hot water and steam from the TES, and the temperature ranges of the recovered fluid (190 – 230°C or 375 – 445°F), conventional geothermal power cycles were used to convert the stored heat to electricity. A power cycle configuration for the GeoTES system was selected following a screening study of a number of flash, and flash/binary hybrid options. This analysis concluded that, of the configurations evaluated, a dual-stage flash steam cycle provides the lowest capital costs per unit net power generation with an acceptable hot brine inlet fluid flow rate. The evaluation included the power plant cost estimate, the cost and number of wells and the associated parasitic loads. Annual power generation performance was simulated to evaluate capacity factor and LCOE. The LCOE calculated for the inherently high capacity GeoTES system was $\$$0.13/kWhe. This value was calculated for the case where the solar thermal collector was sized in such a way that the solar collectors permitted an annual power plant capacity factor of up to 97%. The power cycle was able to provide power to the grid every night of the year, and flexible base-load power during the winter, if needed. This LCOE value compares favorably with reported values for solar photovoltaic plus battery energy storage (PV+BES) systems in the open literature, i.e. $\$$0.148/kWhe for a PV+BES system with 4 hours of electrochemical battery energy storage capacity. Addition of battery energy storage with more hours of storage would further increase PV+BES system LCOE and increase the separation between GeoTES and PV+BES. A GeoTES system would therefore provide superior economics for high capacity and long duration solar energy storage.

14 SOLAR ENERGY↗
DOE OSTI2021

Phasor-Measurement-Unit-Based Data Analytics Using Digital Twin and PhasorAnalytics Software

A major objective of this project was to apply GE’s commercial machine learning and data analytics toolsets to large-scale, real-world, anonymized Phasor Measurement Unit (PMU) datasets in order to extract signatures, correlated and/or causal factors, and precursor patterns associated with significant power system phenomena. The project had a particular emphasis on extraction of insights relevant to asset health monitoring, real-time load modeling and cybersecurity monitoring. Additionally, the team was directed to undertake a comprehensive data quality analysis for the provided datasets and encouraged to estimate the ‘machine-learning readiness’ of the datasets by documenting any major obstacles to the application of commercial machine learning algorithms. To accomplish the aforementioned objectives, the project team’s work centered around the identification of key event signatures and application of the identified event signatures for event detection and event classification. The industry-validated, semi-supervised machine learning strategy employed for event signature identification involved several major tasks, including data-preprocessing, generation of an overabundance of features, normal data identification, normality modeling, and event signature identification through a methodical, quantitative ranking of features in order of relevance to each studied event type. Throughout the project, data quality issues and mitigation techniques were investigated. In this report, insights are provided regarding the readiness of the provided synchrophasor datasets for application of machine learning and data analytics. The methodologies employed for this technical strategy are summarized in this report. With regards to data preprocessing and feature generation, the provided Training and Test Datasets were ingested into GE’s big data environment. Subsequently, the team applied bad data cleansing and data imputation scripts, event detection scripts, and application programming interfaces (APIs) to the datasets for convenient data access. The project team completed development and validation of dozens of physics-based, statistics-based and transformation-based feature functions used for the extraction of over 60 synchrophasor features. Using a new parallel feature generation technology developed on this project, over 60 features have been rapidly generated for the full two years’ worth of Training and Test Dataset data associated with both the Eastern and Western interconnects. Even accommodating for temporal down-sampling inherent to the feature extraction procedure, this parallel feature generation activity resulted in a massive feature set with a storage requirement approximately equal to that of the raw training dataset itself. With regards to normal data identification and normality modeling, a normality model was built using the feature data extracted from the Training Dataset and iteratively refined subsequent to incremental adjustments and expansions of the Training Dataset feature data. With respect to event characterization and signature identification, an event signature identification pipeline was developed and used in conjunction with the normality model to identify over 15 event signatures for key event categories within the Training Dataset. The identified event signatures were used to characterize hundreds of key events in terms of relative severity, duration, and location of the event. An investigation was undertaken to identify correlated and causal factors involved in transformer events. A separate investigation into temporal trends in ring-down analysis results was undertaken to determine possible associations between system dynamics and various other factors such as loading, season or year. To validate the identified event signatures, additional work was undertaken to develop signature-based anomaly detection and classification tools suitable for convenient application to the synchrophasor datasets. The anomaly detection and classification tools, suitable for online application, were then applied to the entirety of the Eastern Interconnect Training and Test Datasets. Performance of the event detection and classification tools was evaluated upon receipt of the Test Dataset event logs (i.e., the labels for events contained in the Test Dataset), and promising results were obtained despite several challenges (documented herein) associated with application of supervised or semi-supervised machine learning methods to large-scale, anonymized datasets. Finally, the detection and classification tools were used to detect, classify, and characterize thousands of new events not included in the original event logs provided by the DOE within both the Training and Test Datasets.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Continuous Counter‐Current Microfluidic Liquid–Liquid Extraction Achieved Using a Pair of Wettable Screen Meshes

Continuous counter‐current microfluidic liquid–liquid extraction performs separations by flowing immiscible liquids in opposing directions within a single flow channel. In principle, this flow arrangement enables a large number of theoretical separation units in a small footprint, without using interstage valving, pumping, and phase separation. Despite its potential for excellent separation performance, this microfluidic scheme rarely appears in literature due to the requirement for capillary forces to be greater than hydrodynamic forces for stable flow. We present a novel microfluidic device and flow approaches that overcome this force‐balance challenge, enabling stable, long‐duration continuous counter‐current flow. Additionally, we cover a suite of methodologies for quantifying the performance of the microfluidic device, revealing the number of theoretical equilibrium stages achieved. The enabling technologies include a woven mesh screen‐based microfluidic device architecture that is easily fabricated outside of a clean room, surface functionalization strategies to promote conjugate (organic/aqueous) wettability, flow approaches to eliminate bubbles and carryover, and computer‐aided flow automation with optical measurement of extraction performance. The reported experiments lasted for over 36 h, terminated only at experiment conclusion, where the device still exhibited good performance. Automated Raman spectroscopy was used for solute quantitation of the ternary system tert‐butanol in a toluene/water matrix, a ternary system that was specifically chosen to analyze the device's performance with a small solute partition ratio and to enable in‐line Raman measurements of solute concentrations in both phases. The microfluidic device possessed a 55 mm contact length and a 38.5 µL internal volume. During counter‐current flow, we observed approximately 37 equilibrium stages (37 ± 13) based on a best‐fit of the solute fraction remaining in the aqueous phase using a Kremser Group Method analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

AI-powered municipal solid waste management: a comprehensive review from generation to utilization

The accumulation of municipal solid waste (MSW) continues to rise due to burgeoning population, rapid global urbanization and economic growth, intensifying ecological concerns associated with landfills and greenhouse gas (GHG) emissions. Over the past 2 decades, global waste generation has surged by 50%, with one-third remaining uncollected and about 70% sent to landfills. This review examines the critical role of integrating emerging technologies, such as advanced sensors and artificial intelligence (AI), into end-to-end MSW management to alleviate landfill burdens. The suitability of various AI tools for different stages of MSW management is assessed, alongside the deployment of advanced sensors including hyperspectral cameras, computer vision systems, and internet of things (IoT) devices for material identification. Applications of genetic algorithms and reinforcement learning for optimizing collection routes, reducing costs, and lowering emissions are highlighted. Life cycle assessment (LCA) across all stages of MSW management is also reviewed, along with future trends in leveraging generative AI, natural language processing (NLP), and agent-based AI systems to analyze waste generation patterns and public sentiment. Efficient collection and handling can be enhanced through route optimization with geographic information systems and real-time bin-level monitoring. Furthermore, sensor-embedded, real-time object detection systems paired with robotics enable material characterization and automated sorting, thereby lowering costs and diverting waste from landfills into value-added products for diverse industrial sectors including packaging, chemicals, textiles, metals and glass, transportation, and electronics industries. Without intervention, global waste is projected to reach 4.54 billion tons by 2050, contributing direct economic costs of $\$$400 billion and roughly 2.38 billion tons of CO 2 -equivalent emissions annually. This review demonstrates how AI-driven, end-to-end solutions for MSW management can mitigate economic and environmental challenges, while directly supporting the United Nations Sustainable Development (UNDP) goals related to innovation and infrastructure (SDG 9), sustainable cities (SDG 11), responsible consumption and production (SDG 12), and climate action (SDG 13).

09 BIOMASS FUELS↗
NASA NTRS1999

A Study of Direct and Cloud-Mediated Radiative Forcing of Climate Due to Aerosols

The Intergovernmental Panel on Climate Change (IPCC) has reported that in the southeastern US and eastern China, the general greenhouse warming due to anthropogenic gaseous emissions is dominated by the cooling effect of anthropogenic aerosols. To verify this model prediction in eastern China and southeastern US, we analyzed regional patterns of climate changes at 72 stations in eastern China during 1951- 94 (44 years), and at 52 stations in the southeastern US during 1949-94 (46 years) to detect the fingerprint of aerosol radiative forcing. It was found that the mean rates of change of annual mean daily, maximum, minimum temperatures and diurnal temperature range (DTR) in eastern China were 0.8, -0.2, 1.8, and -2.0 C/100 years respectively, while the mean rates of change of annual mean daily, maximum, minimum temperatures and DTR in the southeastern US were -0.2, -0.6, 0.2, and -0.8 C/100 years, respectively. This indicates that the high rate of increase in annual mean minimum temperature in eastern China results in a slightly warming trend of daily temperature, while the high rate of decrease in annual mean maximum temperature and low rate of increase in annual mean minimum temperature lead to the cooling trend of daily temperature in the southeastern US. We found that the warming from the longwave forcing due to both greenhouse gases and aerosols was completely counteracted by the shortwave aerosol forcing in the southeastern US in the past 46 years. A slightly overall warming trend in eastern China is evident; winters have become milder. This finding is explained by hypothesizing that increasing energy usage during the past 44 years has resulted in more coal and biomass burning, thus increasing the emission of absorbing soot and organic aerosols in eastern China. Such emissions, in addition to well-known Asia dust and greenhouse gases, may be responsible for the winter warming trend in eastern China that we have reported here. The sensitivity of aerosol radiative properties to aerosol composition, size distribution, relative humidity (RH) is examined for the following aerosol systems: inorganic and organic ions (Cl-, Br-, NO3 -, SO4 2-, Na+, NH4 +, K+, Ca2+, Mg2+, HCOO-, CH3COO-, CH3CH2COO-, CH3COCOO-, OOCCOO2-, MSA-1); water-insoluble inorganic and organic compounds (elemental carbon, n-alkanes, SiO2, Al2O3, Fe2O3 and other organic compounds). The partial molar refraction method was used to calculate the real part of the refractive index. It was found that the asymmetry factor increased by approximately 48% with the real part varying from 1.40 to 1.65, and the single scattering albedo decreased by 24% with the imaginary part varying from -0.005 to -0.1. The asymmetry factor increased by 5.4 times with the geometric standard deviation varying from 1.2 to 3.0. The radiation transmission is very sensitive to the change in size distribution; other factors are not as significant. To determine the aerosol direct radiative forcing (ADRF), the aerosol optical depth (AOD) values at the three operational wavelengths (415, 500 and 673 nm) were determined at a regionally representative site, namely, Mt. Gibbs (35.78 deg N, 82.29 deg W, elevation 2006 m) in Mt. Mitchell State Park, NC, and a site located in an adjacent valley (Black Mountain, 35.66 deg N, 82.38 deg W, elevation 951 m) in the southeastern US. The two sites are separated horizontally by 10 km and vertically by 1 km. It was found that the representative total AOD values at 500 nm at the valley site for highly polluted (HP), marine (M) and continental (C) air masses were 0.68 +/- 0.33, 0.29 +/- 0.19 and 0.10 +/- 0.04, respectively. A search-graph method was used to retrieve the columnar size distribution (number concentration N, effective radius reff and geometric standard deviation=?g) from the optical depth observations at three operational wavelengths. The ground albedo, single scattering albedo and imaginary part of the refractive index were calculated using a mathematically unique procedure involving a Mie code and a radiative transfer code in conjunction with the retrieved aerosol size distribution, AOD, and diffuse-direct irradiance ratio. It was found that N, r(eff) and sigma(g) were in the ranges of 10 to 1.7 x 10(exp 4)/cubic cm, 0.09 to 0.68 micrometers and 1.12 to 2.95, respectively. The asymmetry factor and single scattering albedo were in the ranges of 0.63 to 0.75 and 0.74 to 0.97 respectively. The ground albedo for the forested terrain and imaginary part of refractive index were found to be in the ranges of 0.06 to 0.29 and 0.005 to 0.051 respectively. On the basis of these aerosol radiative properties obtained at the research sites and computations using the Column Radiation Model (CRM) of National Center of Atmospheric Research (NCAR) Community Climate Model (CCM3), it was found that the average cloud-free 24-hour ADRF values were -13 +/- 8, -8 +/- 3, -33 +/- 16 W/square m for marine, continental, and polluted air masses, respectively. On the assumption that the fractional coverage of clouds is 0.61, it was estimated that the annual mean ADRF was 7 +/- 2 W/square m in the southeastern US. The review with respect to the current knowledge of organic acids shows that aerosol formate and acetate concentrations range from 0.02 to 5.3 nmol/cubic m and from 0.03 to 12.4 nmol/cubic m respectively, and that between 34% to 77% of formate and between 21% to 66% of acetate are present in the fine fraction of aerosols. It was found that although most (98-99%) of these volatile organic acids were present in the gas phase, their concentrations in the aerosol particles were sufficient to make them a good candidate for cloud condensation nuclei (CCN). It is hypothesized that organic acids are at least one of the primary sources of CCN in the atmosphere due to their ubiquitous presence in the troposphere, especially over the continental forested areas. The results of our measurements at Palmer Station, Antarctica show that the daily average CCN concentrations at 0.3% and 1% supersaturations ranged from 0.3 to 160/cubic cm and from 4 to 168/cubic cm, respectively, during the period from 17 January to 26 February, 1994. New evidence for substantial and definitive CCN enhancement near and within cloud has been observed at Mt. Mitchell, North Carolina. The results show that the average monthly CCN concentrations were 460 +/- 217, 386 +/- 286, 429 +/- 228 and 238 +/- 134/cubic cm for in-cloud, overcast, clear and rainy conditions, respectively. The typical CCN spectra show that there were a lot of small CCN produced and the ion concentrations (especially H+ and SO4 2-) were very high during the CCN enhancement period. The significantly positive correlation between black carbon (BC) and CCN at 1% supersaturation indicates that a percentage of the BC measured at the site may be in the form of an internal mixture and participated in the formation of CCN.

Yu, Shao-Cai↗
NASA NTRSDate not supplied

Hydrodynamics of Two-Phase Flows through Porous Media in Microgravity: Packed Bed Reactor Experiment onboard of the International Space Station

The objective of the Packed Bed Reactor Experiment was to develop a fundamental understanding of the hydrodynamics of two-phase flow with no phase transition through porous media in microgravity. This work describes two experimental campaigns flown on the International Space Station (ISS) that were designed to achieve the objective. This work presents the results on flow patterns observed, two-phase flow pressure drop in porous media, and the impact of capillary effects on pressure drop at various gas and liquid flow rates. This work is the first to present predictive correlations of the two-phase friction factor for the different regimes identified based on the superficial liquid and gas velocity in microgravity, the first to hypothesize the different regimes based on the change in slope of the pressure gradient versus flow rate plots, the first to address the effects of the capillarity on the pressure gradient and the first to measure column holdup and assess the effective porosity based on the packing and the trapped gas bubbles. Experimental data on pressure drop in gas (N2) – liquid (water) flow show the presence of four different flow regimes in microgravity. Dispersed bubble flow (DB) and pulse (P) flow regimes are detected at high liquid flow rates, whereas at low liquid flow rates “large bubble” or elongated bubble (LB) and “gas channeling” (GC) regimes are observed. For these different flow regimes, different two-phase flow friction factor and pressure gradient correlations are presented for the first time as a function of the gas and liquid modified Reynolds numbers and Suratman number. Within the viscous-capillary (V-C) regime, it is found that the capillary contribution is the dominant force that contributes to the pressure drop for the wetting case (glass). However, for the non-wetting packing (Teflon), the viscous contribution dominates. It was found that the gas hold-up and pressure drop are functions of bed history at low liquid and gas flow rates with the magnitude of the hysteresis decreasing with increasing flow rates. PBRE-2 results show that the capillary force is a strong function of the superficial liquid velocity but is a much weaker function of the superficial gas velocity and varies inversely with the particle diameter. Within the Viscous-Capillary (V-C) regime, over 90% of the pressure gradient is attributed to the capillary contribution in the gas continuous regime. However, in the large bubble regime, the viscous and capillary contributions were comparable. After the completion of PBRE-2 with glass spherical beads, another bed packed with alumina was installed. Pressure gradient data for the alumina packed bed (PBRE-Water Recovery) were transferred to the PI and were not analyzed by the authors of this work. The Packed Bed Reactor Experiment concluded its on-orbit operation after two successful campaigns of testing in 2017 and in 2021. The flight hardware was brought back from orbit and is being reconditioned for a series of future experiments referred to as PBRE-WRS in support of water recovery (WRS). The WRS series consists of testing different two-phase fluid system components and packed beds to assess their pressure gradient characteristics, which will be used for designing packed bed reactors for various applications relevant to life support. Although there is so much relevance of the experimental results obtained from the two PBRE campaign in microgravity to life support, these findings also apply to other applications that involve two phase flow in porous medium such as fuel cells, transport of nutrient to plants in space and other chemical and materials processing that involve two phase flows. These systems operate differently in microgravity because, due to the lack of buoyancy, the density difference between the phases becomes irrelevant and no longer leads to phase separation.

Packed Bed↗
NASA NTRS2024

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↗
DOE OSTI2017

In-situ X-Ray Analysis of Rapid Thermal Processing for Thin-Film Solar Cells: Closing the Gap between Production and Laboratory Efficiency

For materials synthesis, it is well known that the material final state may not reach equilibrium and depends on the synthetic process. In particular, processes that quickly remove the available energy from the material may leave it in a metastable state and the metastability may actually impart desirable functional properties. By its very nature, Rapid thermal processing (RTP) is ideally suited to produce such metastable materials. However, metastability and the dynamics of reaching a metastable state are poorly understood, since this is best accomplished through in situ monitoring. In this regard, RTP is particularly challenging as the processing time are very short (seconds to minutes). As a result, there is only poor understanding, and hence use, of RTP in industry. This is potentially a cost-increasing limitation, because RTP can decrease cost by decreasing processing time, and as such, increase throughput and decrease the total thermal budget of processing - a significant cost. RTP is already being used for key processing steps in PV technologies. With silicon wafer PV, it is used for establishing electrical contact between the Ag metal grid and the silicon (known as firing). In this process, a silicon wafer with deposited metal/frit in a grid pattern is heated rapidly to temperatures between 750 and 800 ºC. The processing time when the temperature is held above 600ºC is short (<5 seconds). This process has historically been optimized empirically and it is unclear how the thermal processing affects formation of the final contact between the metal and the silicon. In the case of thin-film PV, RTP has been demonstrated in the process of making absorber layers, i.e. CIGS and CZTS. Use of RTP can reduce the processing time from 10s of minutes to seconds, reducing the thermal budget and increasing the throughput significantly. The conversion from precursor material to final PV material is not well understood, and most of the process optimization is done empirically. Again, the proposed in situ monitoring will provide insight into the conversion process. We have developed a unique instrument for in situ monitoring of the reaction process during RTP. The most fundamental property to monitor is the atomic arrangement (for solids, crystalline structure) and thus we have focused on this property. The chamber was designed to collect X-ray diffraction (XRD) for crystal phase determination, although it has lent insight in the formation of liquids. We designed and constructed both inert (air) and reactive (S/Se) versions of this chamber. We have applied the RTP-XRD methodology to: silicon solar cell electrical contact formation; the reaction sequence for forming CIGS PV thin films from copper selenide and indium gallium selenide precursors; the phase sequence for forming HC(NH2)2PbI3 (FAPbI3) and CH3NH3PbI3 (MAPbI3) hybrid-perovskite PV thin films. The RTP instrument enables real-time collection of X-ray diffraction data with intervals as short as 100 ms, and while heating at ramp rates as high 100 ºCs-1 up to 1200 ºC. The system is portable and can be installed on a synchrotron beamline. The RTP reactor consists of a metallic body (front cap assembly, central block and end cap assembly) that houses a quartz reaction chamber. The sample is secured inside the reaction chamber on a sample holder connected to the front-cap assembly. Modifications of the chamber were made for safely handling sulfur and selenium vapor environments. A Se/S source holder was added to generate a Se/S over pressure in the chamber on heating and the entrance window is made wider in order to accommodate the detachable additional holder. Se/S are toxic in nature and therefore a proper Se/S trapping is used before releasing the gases to exhaust. Screen-printing provides an economically attractive means for making Ag electrical contacts to Si solar cells, but the use of Ag substantiates a significant manufacturing cost, and the glass frit used in the paste to enable contact formation contains Pb. To achieve optimal electrical performance and to develop pastes with alternative, abundant, and non-toxic materials requires understanding the contact formation process during firing. We use the RTP-XRD chamber during firing of Ag-Si contacts to reveal the reaction sequence for contact formation. The findings show that between 500 ºC and 650 ºC PbO in the frit etches the SiNx antireflective-coating on the solar cell, exposing the Si surface. Then, above 650 ºC, Ag+ dissolves into the molten glass frit – key for enabling deposition of metallic Ag on the emitter surface and precipitation of Ag nanocrystals within the glass. This clarifies contact formation mechanisms and can be used to suggest approaches for development of inexpensive, nontoxic solar cell contacting pastes. We used the RTP/XRD chamber to study the formation of CuInGaS/Se (CIGS) PV films by in-situ selenization of pre-deposited metal stacks and annealing of pre-deposited Cu-In-Ga-Se precursors. During the metal stack selenization, we observed that Se first reacts with the metals and forms selenides. We observed formation of intermetallics such as Cu2In and binary selenides during heating. On further annealing in a Se atmosphere, CuInSe2 and CuGaSe2 forms which on heating up to 550 °C reacts further to form the CIGS phase. On prolonged exposure to 550 °C for more than 5 mins or higher temperature, we see separation into Ga-rich and Ga-depleted CIGS phases. On introducing a holding step at 400 °C for 10 min during the annealing ramp, the CIGS phase (forms between 500-550 °C), appears to be less prone to decomposition into high-Ga and low-Ga phases during holding at peak temperature. Lead halide perovskites have emerged as successful optoelectronic materials with high photovoltaic power conversion efficiencies and low material cost. However, substantial challenges remain in the scalability, stability and fundamental understanding of the materials. We demonstrated the application of radiative thermal annealing, an easily scalable processing method, for synthesizing formamidinium lead iodide (FAPbI3) perovskite solar absorbers. Devices fabricated from films formed via radiative thermal annealing have equivalent efficiencies to those annealed using a conventional hotplate. By coupling results from in situ x-ray diffraction using a radiative thermal annealing system with device performances, we mapped the processing phase space of FAPbI3 and corresponding device efficiencies. Our map of processing-structure-performance space suggests the commonly used FAPbI3 annealing time, 10 min at 170 °C, can be significantly reduced to 40 s at 170 °C without affecting the photovoltaic performance. Applying the Johnson-Mehl-Avrami model to the decomposition into PbI2, the activation energy was determined. We also presented a well-controlled, manufacturing relevant annealing method, radiative thermal annealing, for the methylammonium lead triiodide films formed by solvent engineering process, with dimethylformamide (DMF) and dimethyl sulfoxide as solvent and diethyl ether as the antisolvent. Radiative thermal annealing can produce high quality perovskite films, evidenced by high efficiency solar cell devices, in a shorter time than the widely-used hotplate annealing. Using in-situ x-ray diffraction during the radiative annealing, we show that the role of the antisolvent is not to form an intermediate compound (a methylammonium iodide, lead iodide, and dimethyl sulfoxide compound) by washing of the main solvent (DMF) but to achieve a pin-hole free, uniform film of methylammonium lead iodide. We find that having the intermediate compound does not guarantee a high quality (pin-hole free) perovskite film. Our study was extended to reveal the effect of annealing environment, annealing temperature profile, and deposition substrate.

36 MATERIALS SCIENCE↗