Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal network model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

A Novel Approach for Modeling Chemical Reaction in Generalized Fluid System Simulation Program

The Generalized Fluid System Simulation Program (GFSSP) is a computer code developed at NASA Marshall Space Flight Center for analyzing steady state and transient flow rates, pressures, temperatures, and concentrations in a complex flow network. The code, which performs system level simulation, can handle compressible and incompressible flows as well as phase change and mixture thermodynamics. Thermodynamic and thermophysical property programs, GASP, WASP and GASPAK provide the necessary data for fluids such as helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, water, a hydrogen, isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, several refrigerants, nitrogen trifluoride and ammonia. The program which was developed out of need for an easy to use system level simulation tool for complex flow networks, has been used for the following purposes to name a few: Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump Secondary Flow Circuits, Axial Thrust Balance of the Fastrac Engine Turbopump, Pressurized Propellant Feed System for the Propulsion Test Article at Stennis Space Center, X-34 Main Propulsion System, X-33 Reaction Control System and Thermal Protection System, and International Space Station Environmental Control and Life Support System design. There has been an increasing demand for implementing a combustion simulation capability into GFSSP in order to increase its system level simulation capability of a liquid rocket propulsion system starting from the propellant tanks up to the thruster nozzle for spacecraft as well as launch vehicles. The present work was undertaken for addressing this need. The chemical equilibrium equations derived from the second law of thermodynamics and the energy conservation equation derived from the first law of thermodynamics are solved simultaneously by a Newton-Raphson method. The numerical scheme was implemented as a User Subroutine in GFSSP.

Sozen, Mehmet↗

FORCE Development Status Update: Vertical Integration and Benchmarking of System Dynamics

Recent efforts to establish effective models for grid energy analysis, especially given the increase in variable renewable energy (VRE) sources and the economic challenges faced by traditional nuclear energy, have generated new technological considerations. One effort to improve the economic viability of nuclear power involves investigating integrated energy systems (IES) which include secondary energy systems that introduce flexibility and secondary market possibilities to existing and perceived future nuclear energy generation technologies. To analyze the technical and economic potential of IES, the Framework for Optimization of Resources and Economics (FORCE) tool suite was developed through a collaboration among national laboratories. Within the FORCE tool suite, the Holistic Energy Resource Optimization Network (HERON) provides algorithms for analyzing the long-term viability of potential IES technologies, while HYBRID provides algorithms and models to achieve high-resolution analysis of coupling physics over a short time period. Continued maturing of the FORCE tool suite requires further interconnections between the various tools in the suite in order to ensure consistent analysis. Analyses performed by applying HYBRID to transient process modeling should be easily harvestable as inputs to HERON analyses. The first item in this status update is a demonstration of an automated data pipeline for loading data from HYBRID into HERON analyses. Application of HYBRID results to HERON, as part of using the FORCE tool suite, relies on robust modeling assumptions for the various models included in HYBRID. The second result of this status update is the benchmarking and validation of cost and operational data, with a particular focus on natural gas energy generators. These generators are benchmarked with a focus on contrasting them with proposed thermal energy storage (TES) technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Analysis of Defects in Metal Additive Manufacturing with Augmented Data Generation

Laser powder bed fusion (LPBF) is a method of additive manufacturing (AM) that selectively melts and fuses together microscopic metallic powder. LPBF offers the benefit of producing custom structures out of high strength metals that can be difficult to fabricate with conventional methods. The challenge of LPBF is that 3D printed structures often have internal pores due to process flaws. Pulsed thermal tomography (PTT) is a method for reconstructing the depth profile of materials, allowing the visualization internal voids in solids. In prior work, we developed a convolutional neural network (CNN) which, having been trained on simulated 2D PTT images of subsurface elliptical defects, was able to classify the semi-major radii, semi-minor radii, and angular orientation of the best-fit ellipses in previously unseen PTT images. The unseen PTT images contained subsurface irregular defect shapes imported from scanning electron microscopy (SEM) images of metallic LPBF-printed specimens. Training the CNN on irregular defect shapes instead of on elliptical shapes would make the resulting classifications more descriptive of actual defect shapes. However, this requires a much higher volume of SEM images of material defects, which are difficult to obtain because of random occurrence of defects in LPBF. To address this challenge, we developed a generative adversarial network (GAN) to augment the existing dataset of SEM defect images. The GAN model is demonstrated to create novel yet realistic defect shapes that can be used as input for simulated PTT images to train CNN.

36 MATERIALS SCIENCE↗

Embedded Error Bayesian Calibration of Thermal Decomposition of Organic Materials

Organic materials are an attractive choice for structural components due to their light weight and versatility. However, because they decompose at low temperatures relative to tradiational materials they pose a safety risk due to fire and loss of structural integrity. To quantify this risk, analysts use chemical kinetics models to describe the material pyrolysis and oxidation using thermogravimetric analysis. This process requires the calibration of many model parameters to closely match experimental data. Previous efforts in this field have largely been limited to finding a single best-fit set of parameters even though the experimental data may be very noisy. Furthermore the chemical kinetics models are often simplified representations of the true de- composition process. The simplification induces model-form errors that the fitting process cannot capture. In this work we propose a methodology for calibrating decomposition models to thermogravimetric analysis data that accounts for uncertainty in the model-form and experimental data simultaneously. The methodology is applied to the decomposition of a carbon fiber epoxy composite with a three-stage reaction network and Arrhenius kinetics. The results show a good overlap between the model predictions and thermogravimetric analysis data. Uncertainty bounds capture devia- tions of the model from the data. The calibrated parameter distributions are also presented. In conclusion, the distributions may be used in forward propagation of uncertainty in models that leverage this material.

36 MATERIALS SCIENCE↗

Sneutrino tribrid inflation, metastable cosmic strings and gravitational waves

Here, we present a successful realization of sneutrino tribrid inflation model based on a gauged U(1)_B-L extension of Minimal Supersymmetric Standard Model (MSSM). A single interaction term involving the B-L Higgs field and the right-handed neutrinos serves multiple purposes. These include the generation of heavy Majorana masses for the right-handed neutrinos to provide an explanation for the tiny neutrino masses via the seesaw mechanism, a realistic scenario for reheating and non-thermal leptogenesis with a reheat temperature as low as 10 6 GeV, and a successful realization of inflation with right-handed sneutrino as the inflaton. The matter parity which helps avoid rapid proton decay survives as a Z 2 subgroup of a U(1) R-symmetry. Depending on the choice of model parameters yields the following predicted range of the tensor to scalar ratio, 3 x 10 -11 ≲ r≲ 7x 10 -4 (6 x 10 -7 ≲ r ≲ 0.01), and the running of the scalar spectral index, -0.00022 ≲ dn_s/dln k ≲ -0.0026 (-0.00014 ≲ dn_s/dln k ≲ 0.005), along with the B-L breaking scale, 3 x 10 14 ≲ M/ GeV≲ 5 x 10 15 (6 x 10 15 ≲ M/ GeV≲ 2 x 10 16 ), calculated at the central value of the scalar spectral index, n_s =0.966, reported by Planck 2018. The possibility of realizing metastable cosmic strings in a grand unified theory setup is briefly discussed. The metastable cosmic string network admits string tension values in the range 10 -8 ≲ Gμ s ≲ 10 -6 , and predicts a stochastic gravitational wave background lying within the 2-σ bounds of the recent NANOGrav 12.5-yr data.

supersymmetry and cosmology↗

Defect Production and Microstructural Feature Impact for Radiation Damage in Additively Manufactured 316 Stainless Steel

This milestone presents multi-scale modeling research results for additively manufactured 316 stainless steel. A combination of phase field, cluster dynamics, molecular dynamics, and density functional theory with machine learning is used, allowing for predictions of radiation-driven microstructural evolution in additively manufactured 316 stainless steel over a range of temperatures, damage rates, neutron spectra, and microstructures, and supporting the development of combined ion and neutron irradiation for material qualification. Informed by ion irradiation and neutron irradiation results across the Advanced Materials and Manufacturing Technologies program, we investigate the unique aspects of radiation-driven microstructure evolution in additively manufactured 316 stainless steel. In particular, we focus on understanding the impact of carbon concentration (varying, for example, between the 316L and 316H standards) on void formation; radiation-induced segregation at dislocation cells and grain boundaries; and the evolution of dislocation loops and network dislocation populations. We find that the unique characteristics of the additively manufactured microstructures must be accounted for in understanding the evolution of dislocation populations under thermal and irradiation conditions, such as the variation in sink strengths arising due to the variation in dislocation density. We also find that the radiation-induced segregation of Cr and Ni to grain boundaries and cell walls differs due to the differences in their defect sink biases. We also find that increasing the Ni content can slow vacancy diffusion, which may provide a mechanism for the observed reduction in transient swelling rate for austenitic Fe-Cr-Ni alloys with increasing Ni content. In addition, ion irradiations have shown that increasing carbon content in 316 SS results in a larger population of smaller voids, suggesting reduced vacancy diffusion. Our results show that the carbon content of additively manufactured 316 SS has a significant impact on the migration rate of defect clusters. The presence of carbon atoms results in carbon-vacancy trapping, significantly reducing the diffusion rate of vacancies. Carbon atoms may also be trapped near the surface of a void, which may reduce void growth by trapping vacancies that diffuse toward the void.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Three-dimensional compressible convection at low Prandtl numbers

Numerical simulations are used to study fully compressible thermal convection at large Rayleigh numbers. Results are presented from a sequence of three-dimensional simulations that reveal a transition from gradually-evolving laminar convection to nearly turbulent convection as the Prandtl number is reduced from a value of unity to one-tenth. The convective flows form irregular cellular patterns near the upper surface, possesing a network of fast downflow at cell peripheries and gentler upflow at cell centers. At greater depths the curving sheets of downflow collapse into plumes which may twist and possess substantial vertical vorticity. For the lowest Prandtl number, the convection near the bottom of the layer appears to be turbulent, yet the rapidly varying small-scale flow structure there is accompanied by more ordered sites of wavering upflow, with the latter able to penetrate all the way to the upper boundary. Thus a significant component of the flow is able to extend over multiple density scale heights, in contrast to what is argued in formulating mixing-length models for stellar convection. Results are also shown from two-dimensional simulations carried out with very high spatial resolution, which reveal that supersonic convection with fluttering shock systems can be realized.

Toomre, Juri↗

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Thermal-Conductivity Fabrics

Heat management with common textiles such as nylon and spandex is hindered by the poor thermal conductivity from the skin surface to cooling surfaces. This innovation showed marked improvement in thermal conductivity of the individual fibers and tubing, as well as components assembled from them. The problem is centered on improving the heat removal of the liquid-cooled ventilation garments (LCVGs) used by astronauts. The current design uses an extensive network of water-cooling tubes that introduces bulkiness and discomfort, and increases fatigue. Range of motion and ease of movement are affected as well. The current technology is the same as developed during the Apollo program of the 1960s. Tubing material is hand-threaded through a spandex/nylon mesh layer, in a series of loops throughout the torso and limbs such that there is close, form-fitting contact with the user. Usually, there is a nylon liner layer to improve comfort. Circulating water is chilled by an external heat exchanger (sublimator). The purpose of this innovation is to produce new LCVG components with improved thermal conductivity. This was addressed using nanocomposite engineering incorporating high-thermalconductivity nanoscale fillers in the fabric and tubing components. Specifically, carbon nanotubes were added using normal processing methods such as thermoplastic melt mixing (compounding twin screw extruder) and downstream processing (fiber spinning, tubing extrusion). Fibers were produced as yarns and woven into fabric cloths. The application of isotropic nanofillers can be modeled using a modified Nielsen Model for conductive fillers in a matrix based on Einstein s viscosity model. This is a drop-in technology with no additional equipment needed. The loading is limited by the ability to maintain adequate dispersion. Undispersed materials will plug filtering screens in processing equipment. Generally, the viscosity increases were acceptable, and allowed the filled polymers to still be processed.The novel feature is that fabrics do not inherently possess good thermal conductivity. In fact, fabrics are used for thermal insulation, not heat removal. The technology represents the first material that is a wearable fabric, based on company textiles and materials that will significantly conduct heat.

Chibante, L. P. Felipe↗

Thermo-hydraulic steam pipe models for district heating simulations: Simplifications to balance accuracy and simulation speed

Steam piping networks are essential for optimizing performance in industrial processes and district heating systems. However, dynamic models that balance thermo-hydraulic accuracy with computational efficiency remain limited. In response, this paper presents a new discretized steam pipe model based on the plug flow approach, capturing key thermo-hydraulic behaviors while simplifying steam phase change processes. Implemented in Modelica, the model accurately calculates temperature and pressure distributions along steam pipelines. To improve computational efficiency for district-scale simulations, five model simplifications are introduced: lumped thermo-hydraulic functions, empirical correlations, fluid state approximations, steady-state dynamics and inclusion of flow derivatives. These simplified models achieve 85%-98% accuracy in predicting pressure drop and condensation losses, including dynamic condensate behavior during pipe warm-up—a factor often overlooked in existing models. The models support diverse network configurations, scaling effectively to systems with multiple distribution pipes and connected building loads. Discrete models provide detailed insights but exhibit a cubic increase in simulation time as the network scales by N connected building O(N 2.42 ). In contrast, lumped models simulate 10–28 times faster than discrete, offering quadratic scaling of simulation time O(N 1.73 ). However, they still require 6 times more computation time than a lossless network, highlighting the inherent computational challenges of modeling compressible fluid flow. In conclusion, the steady-state lumped variant, with its near-linear scalability in computational time O(N 1.01 ), emerges as an efficient solution for preliminary design evaluations and extensive parametric studies.

15 GEOTHERMAL ENERGY↗

Formation and dissociation of shear-induced high-energy dislocations: insight from molecular dynamics simulations

Solid-phase processing (SPP) allows one to create complex microstructures, not achievable via thermal processing alone. The resulting structures exhibit a rich palette of defects, both thermal and non-thermal, including defect substructures, such as dislocation networks. It is essential to understand the mechanisms of deformation and defect structure formation to guide SPP towards achieving desired microstructures and material properties. In this study, large-scale molecular dynamics simulations are used to investigate the effects of inhomogeneous strain distribution, that mimics deformation conditions of tribological tests, on the evolution of defects under severe shear deformation in polycrystalline Al. Analysis of defect nucleation and reaction pathways reveals that strong geometric constraints suppress the nucleation and slide of low energy dislocation 1/2$\langle110\rangle${111} but promote the nucleation and slide of high energy dislocations, such as [$1\bar1$0] (001) and 1/2 [$1\bar1$$\bar2$]($1\bar1$1). Further, a rough contact surface, characteristic to tribological tests, imposes an inhomogeneous stress field leading to inhomogeneous defect substructures due to location-dependent activation of slip systems. The results suggest that high-energy dislocations can dominate the evolution of grain structures in highly constrained environments, which should be considered in modeling plastic deformation and grain refinement during SPP.

42 ENGINEERING↗

Piezocomposite Actuator Arrays for Correcting and Controlling Wavefront Error in Reflectors

Three reflectors have been developed and tested to assess the performance of a distributed network of piezocomposite actuators for correcting thermal deformations and total wave-front error. The primary testbed article is an active composite reflector, composed of a spherically curved panel with a graphite face sheet and aluminum honeycomb core composite, and then augmented with a network of 90 distributed piezoelectric composite actuators. The piezoelectric actuator system may be used for correcting as-built residual shape errors, and for controlling low-order, thermally-induced quasi-static distortions of the panel. In this study, thermally-induced surface deformations of 1 to 5 microns were deliberately introduced onto the reflector, then measured using a speckle holography interferometer system. The reflector surface figure was subsequently corrected to a tolerance of 50 nm using the actuators embedded in the reflector's back face sheet. Two additional test articles were constructed: a borosilicate at window at 150 mm diameter with 18 actuators bonded to the back surface; and a direct metal laser sintered reflector with spherical curvature, 230 mm diameter, and 12 actuators bonded to the back surface. In the case of the glass reflector, absolute measurements were performed with an interferometer and the absolute surface was corrected. These test articles were evaluated to determine their absolute surface control capabilities, as well as to assess a multiphysics modeling effort developed under this program for the prediction of active reflector response. This paper will describe the design, construction, and testing of active reflector systems under thermal loads, and subsequent correction of surface shape via distributed peizeoelctric actuation.

piezoelectric actuation↗

Evaluation of single-footprint AIRS CH4 Profile Retrieval Uncertainties Using Aircraft Profile Measurements

We evaluate the uncertainties of methane optimal estimation retrievals from single footprint thermal infrared observations from the Atmospheric Infrared Sounder (AIRS). These retrievals are primarily sensitive to atmospheric methane in the mid-troposphere through the lower stratosphere (~2 to ~17 km). We compare to in situ observations made from aircraft during the Hiaper Pole to Pole Observations (HIPPO), the NASA Atmospheric Tomography Mission (ATom) campaigns, and from the NOAA ESRL aircraft network, between the surface and 5-13 km, across a range of years, latitudes between 60 S to 80 N, and over land and ocean. After a global, pressure dependent bias correction, we find that the land and ocean have similar biases and that the reported observation error (combined measurement and interference errors) of ~27 ppb is consistent with the standard deviation between aircraft and individual AIRS observations. A single measurement has measurement (noise related) uncertainty of ~17 ppb, a ~20 ppb uncertainty from radiative interferences (e.g. from water, temperature, etc.), and ~ 30 ppb due to “smoothing error”, which is partially removed when making comparisons to in situ measurements or models in a way that account for this regularization. We estimate a 16 ppb validation error because the aircraft typically did not measure methane at altitudes where the AIRS measurements have some sensitivity, e.g. the stratosphere. Daily averaged AIRS measurements of at least 9 observations over spatio-temporal domains of < 1 degree and 1 hour have a standard deviation of ~17 ppb versus aircraft, likely because the observation errors from temperature and water vapor (for example) are only partly reduced through averaging. Seasonal averages can reduce this ~17 ppb uncertainty further to ~10 ppb, as determined through comparison with NOAA aircraft, likely because uncertainties related to radiative effects of temperature and water vapor can be reduced when averaged over a season.

Susan Kulawik↗

Anomalous softening of 3D printed elastomeric foam irradiated under compressive strain

Elastomeric foam is an essential component in many industrial and technological settings, primarily as thermal insulators and as positional/mechanical support cushions. In particular, silicone foam is utilized in harsh environments due to exceptional thermal and chemical stability. Under service conditions within certain applications such material gets exposed to a high dosage of gamma radiation, which can permanently alter the material’s structural and mechanical response properties. Most studies on gamma-exposure under inert or oxidative atmosphere indicate hardening of silicone foam, which is attributed to radiation-induced enhancement in chemical cross-linking. Here we report two contrasting effects depending on whether (non-oxidative) radiation exposure is carried out with the foam under zero or finite compressive strain. While in the former case we observe radiation-hardening consistent with previous studies, in the latter case (50% porous foam under 30% uniaxial compression) we see a monotonic decrease in Young’s modulus with increasing dosage, although solvent swelling experiments on the constituent rubber indicate a net increase in cross-link density independent of the state of strain. We quantitatively model all dose-dependent data using the Ogden Hyperfoam strain-energy function within the framework of Tobolsky two-network scheme and attribute the above anomaly to a combined effect of radiation-induced thickness change (compression set) and inherent nonlinearity in the foam’s stress-strain response.

Coarse-grained models↗

NASA Tech Briefs, October 2013

Topics include: A Short-Range Distance Sensor with Exceptional Linearity; Miniature Trace Gas Detector Based on Microfabricated Optical Resonators; Commercial Non-Dispersive Infrared Spectroscopy Sensors for Sub-Ambient Carbon Dioxide Detection; Fast, Large-Area, Wide-Bandgap UV Photodetector for Cherenkov Light Detection; Mission Data System Java Edition Version 7; Adaptive Distributed Environment for Procedure Training (ADEPT); LEGEND, a LEO-to-GEO Environment Debris Model; Electronics/Computers; Millimeter-Wave Localizers for Aircraft-to-Aircraft Approach Navigation; Impedance Discontinuity Reduction Between High-Speed Differential Connectors and PCB Interfaces; SpaceCube Version 1.5; High-Pressure Lightweight Thrusters; Non-Magnetic, Tough, Corrosion- and Wear-Resistant Knives From Bulk Metallic Glasses and Composites; Ambient Dried Aerogels; Applications for Gradient Metal Alloys Fabricated Using Additive Manufacturing; Passivation of Flexible YBCO Superconducting Current Lead With Amorphous SiO2 Layer; Propellant-Flow-Actuated Rocket Engine Igniter; Lightweight Liquid Helium Dewar for High-Altitude Balloon Payloads; Method to Increase Performance of Foil Bearings Through Passive Thermal Management; Unibody Composite Pressurized Structure; JWST Integrated Science Instrument Module Alignment Optimization Tool; Radar Range Sidelobe Reduction Using Adaptive Pulse Compression Technique; Digitally Calibrated TR Modules Enabling Real-Time Beamforming SweepSAR Architectures; Electro-Optic Time-to-Space Converter for Optical Detector Jitter Mitigation; Partially Transparent Petaled Mask/Occulter for Visible-Range Spectrum; Educational NASA Computational and Scientific Studies (enCOMPASS); Coarse-Grain Bandwidth Estimation Scheme for Large-Scale Network; Detection of Moving Targets Using Soliton Resonance Effect; High-Efficiency Nested Hall Thrusters for Robotic Solar System Exploration; High-Voltage Clock Driver for Photon-Counting CCD Characterization; Development of the Code RITRACKS; and Enabling Microliquid Chromatography by Microbead Packing of Microchannels.

Source record↗

An Integrated Hydrological and Water Management Study of the Entire Nile River System - Lake Victoria to Nile Delta

The Nile basin River system spans 3 million km(exp 2) distributed over ten nations. The eight upstream riparian nations, Ethiopia, Eretria, Uganda, Rwanda, Burundi, Congo, Tanzania and Kenya are the source of approximately 86% of the water inputs to the Nile, while the two downstream riparian countries Sudan and Egypt, presently rely on the river's flow for most of the their needs. Both climate and agriculture contribute to the complicated nature of Nile River management: precipitation in the headwaters regions of Ethiopia and Lake Victoria is variable on a seasonal and inter-annual basis, while demand for irrigation water in the arid downstream region is consistently high. The Nile is, perhaps, one of the most difficult trans-boundary water issue in the world, and this study would be the first initiative to combine NASA satellite observations with the hydrologic models study the overall water balance in a to comprehensive manner. The cornerstone application of NASA's Earth Science Research Results under this project are the NASA Land Data Assimilation System (LDAS) and the USDA Atmosphere-land Exchange Inverse (ALEXI) model. These two complementary research results are methodologically independent methods for using NASA observations to support water resource analysis in data poor regions. Where an LDAS uses multiple sources of satellite data to inform prognostic simulations of hydrological process, ALEXI diagnoses evapotranspiration and water stress on the basis of thermal infrared satellite imagery. Specifically, this work integrates NASA Land Data Assimilation systems into the water management decision support systems that member countries of the Nile Basin Initiative (NBI) and Regional Center for Mapping of Resources for Development (RCMRD, located in Nairobi, Kenya) use in water resource analysis, agricultural planning, and acute drought response to support sustainable development of Nile Basin water resources. The project is motivated by the recognition that accurate, frequent, and spatially distributed estimates of the water balance are necessary for effective water management. This creates a challenge for watersheds that are large, include data poor regions, and/or span multiple nations. All of these descriptors apply to the Nile River basin, yet successful management of the Nile is critical for development and political stability in the region. For this reason, improved hydrological data to support cooperative water management in the Nile basin is a priority for USAID, the US State Department, the World Bank and other international organizations. In this project, the U.S. based research team is working with partners at RCMRD, Nile Basin Initiative (NBI), and their member national-level agencies to develop satellite-based land cover maps, satellite-derived evapotranspiration estimates (using the ALEXI algorithm), and NASA's Land Data Assimilation System (LDAS) customized to match identified information needs. The cornerstone applied sciences product of the project is the development of a customized "Nile LDAS" that will produce optimal estimates of hydrological states and fluxes, as vetted against the in situ observations of NBI and RCMRD member organizations and independent satellite-derived hydrological estimates. Nile LDAS will be applied to improve the reliability of emerging Decision Support Systems in applications that include drought monitoring, reservoir management, and irrigation planning. The end-users such as RCMRD, NBI, Ethiopian and Kenya Meteorological and Famine Early Warning System Network (FEWSNet) will be the eventual benefactors of this work. There will be a capacity building process involving the above end-user organizations and transfer the models and the results for these organizations to execute for future use. The team has already initiated this study and the early results of first years' work are shown. The plan is to complete this work by late 2013.

Habib, Shahid↗

Biomethanation to Upgrade Biogas to Pipeline Grade Methane

NREL is working closely with DOE, Electrochaea GmbH, and Southern California Gas Company (SoCalGas) to reduce costs of a biomethanation process capable of megawatt-scale deployment that upgrades organic biogas waste streams to produce pipeline quality renewable natural gas (RNG). Biomethanation is a two-step process using a single-celled methanogenic archaea that converts low-carbon low-cost hydrogen (H2) and waste carbon dioxide (CO2) to produce renewable methane (CH4). The process upgrades the biogenic CO2 - while allowing the CH4 to pass through - from biogas sources like dairies, wastewater treatment plants, and landfills. The CH4 produced is a drop-in direct replacement fuel and producers can participate in the growing number of carbon markets; like California's Low Carbon Fuel Standard and the Federal Renewable Fuel Standard. NREL and Argonne National Laboratory have completed a life cycle analysis using the GREET model to show that the biomethanation process produces RNG that is carbon negative even when H2 production via low-temperature water electrolysis is driven by the existing carbon intensity of California's electricity grid. And of course, even further carbon negative (-233 kg CO2e/kWh) when the electricity is produced from low-carbon sources like wind and solar. Leveraging lessons learned from operating SoCalGas' 700L 18-bar bioreactor system, NREL is designing and building a flexible RD&D platform that will enable field trials at biogas and other CO2 sources. A custom 16' long trailer will house a 20L 18-bar bioreactor, 3 - 25 kW proton exchange membrane electrolyzer, and dosing, thermal, and controls systems to support operations with only power, biogas, and water feedstocks required by the field locations. The end-of-project goal is to demonstrate pipeline quality RNG production (> 95% CH4, < 4% H2, <1% CO2, < 0.2% O2 and < 4 parts per million H2 sulfide) using real biogas feedstocks - thereby recycling both greenhouse gases for injection into the natural gas network or to be used onsite.

biogas upgrading↗

Coherency-Constrained Spectral Clustering for Power Network Reduction

This paper presents a methodology for reducing the complexity of large-scale power network models using spectral clustering, aggregation of electrical components, and cost function approximation. Two approaches are explored using unconstrained and constrained spectral clustering to determine areas for effective system reduction. Once the system areas are determined, both loads and generators by type are aggregated, and their new cost function is approximated through polynomial curve-fitting or statistical methods. The performance of reduced networks is evaluated in terms of their ability to follow the true daily cost of the original system over a 24-hour period considering a set of several days. Two test systems are taken as test beds. Application of the methodology to a modified version of the IEEE 39-bus system reduces it from 17 generators to a 4-bus system and 9 generators with about 93% of accuracy. Similarly, the IEEE 118-bus system is reduced from 19 generators to a 3-bus system with three aggregated units achieving over 99% of accuracy. These findings address scalability challenges and enhance accuracy for high and mid-loading level conditions, and by aggregating thermal units with similar cost functions.

42 ENGINEERING↗