Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid strength”

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 235 records · Page 13

Multilayered Cu-Carbon Nanotube Composites for Advanced Conductors

Improving the efficiency of electrical components is critical in reducing energy consumption for various industrial and residential applications, ranging from rotating machinery to all electric devices and electric vehicle (EV) components to power grid systems. Substituting Cu wires with reduced resistance conductors that incorporate carbon nanotubes (CNTs) into Cu─ultraconductive Cu (UCC) composites─has recently been considered a promising strategy to improve energy efficiency, power density, and/or performance across various applications. Here, in this study, we created stable material formulations [CNT-containing polyvinylpyrrolidone (PVP) in dimethylformamide (DMF) solution] and utilized commercially viable fabrication approaches (electrospinning and magnetron sputtering) that produced high-performance multilayered tape-based UCC composite architectures. Increasing the CNT volume fraction by sequential layering of the structure with additional Cu-CNT layers showed a nearly stepwise improved performance in electrical and mechanical properties. This study also provides valuable insight into the effectiveness of nitrogen doping in modifying the conductivity of the CNT matrix. Fabricated prototypes demonstrated a >10% increase in current carrying capacity and >10% improvement in mechanical strength compared to those obtained on pure Cu. We believe that the properties demonstrated here, combined with the scalable manufacturing pathway of our approach, pave the way in designing future advanced conductors for diverse energy efficient and high-performance electrical systems and applications.

Carbon nanotubes↗

An Estimate of Low-Cloud Feedbacks from Variations of Cloud Radiative and Physical Properties with Sea Surface Temperature on Interannual Time Scales

Simulations of climate change have yet to reach a consensus on the sign and magnitude of the changes in physical properties of marine boundary layer clouds. In this study, the authors analyze how cloud and radiative properties vary with SST anomaly in low-cloud regions, based on five years (March 2000 - February 2005) of Clouds and the Earth s Radiant Energy System (CERES) -- Terra monthly gridded data and matched European Centre for Medium-Range Weather Forecasts (ECMWF) meteorological reanalaysis data. In particular, this study focuses on the changes in cloud radiative effect, cloud fraction, and cloud optical depth with SST anomaly. The major findings are as follows. First, the low-cloud amount (-1.9% to -3.4% /K) and the logarithm of low-cloud optical depth (-0.085 to -0.100/K) tend to decrease while the net cloud radiative effect (3.86 W/m(exp 2)/ K) becomes less negative as SST anomalies increase. These results are broadly consistent with previous observational studies. Second, after the changes in cloud and radiative properties with SST anomaly are separated into dynamic, thermodynamic, and residual components, changes in the dynamic component (taken as the vertical velocity at 700 hPa) have relatively little effect on cloud and radiative properties. However, the estimated inversion strength decreases with increasing SST, accounting for a large portion of the measured decreases in cloud fraction and cloud optical depth. The residual positive change in net cloud radiative effect (1.48 W/m(exp 2)/ K) and small changes in low-cloud amount (-0.81% to 0.22% /K) and decrease in the logarithm of optical depth (-0.035 to -0.046/ K) with SST are interpreted as a positive cloud feedback, with cloud optical depth feedback being the dominant contributor. Last, the magnitudes of the residual changes differ greatly among the six low-cloud regions examined in this study, with the largest positive feedbacks (approximately 4 W/m(exp 2)/ K) in the southeast and northeast Atlantic regions and a slightly negative feedback (-0.2 W/m(exp 2)/ K) in the south-central Pacific region. Because the retrievals of cloud optical depth and/or cloud fraction are difficult in the presence of aerosols, the transport of heavy African continental aerosols may contribute to the large magnitudes of estimated cloud feedback in the two Atlantic regions.

Eitzen, Zachary A.↗

A comparative study of Navier-Stokes codes for high-speed flows

A comparative study was made with four different codes for solving the compressible Navier-Stokes equations using three different test problems. The first of these cases was hypersonic flow through the P8 inlet, which represents inlet configurations typical of a hypersonic airbreathing vehicle. The free-stream Mach number in this case was 7.4. This 2-D inlet was designed to provide an internal compression ratio of 8. Initial calculations were made using two state-of-the-art finite-volume upwind codes, CFL3D and USA-PG2, as well as NASCRIN, a code which uses the unsplit finite-difference technique of MacCormack. All of these codes used the same algebraic eddy-viscosity turbulence model. In the experiment, the cowl lip was slightly blunted; however, for the computations, a sharp cowl leading edge was used to simplify the construction of the grid. The second test problem was the supersonic (Mach 3.0) flow in a three-dimensional corner formed by the intersection of two wedges with equal wedge angles of 9.48 degrees. The flow in such a corner is representative of the flow in the corners of a scramjet inlet. Calculations were made for both laminar and turbulent flow and compared with experimental data. The three-dimensional versions of the three codes used for the inlet study (CFL3D, USA-PG3, and SCRAMIN, respectively) were used for this case. For the laminar corner flow, a fourth code, LAURA, which also uses recently-developed upwind technology, was also utilized. The final test case is the two-dimensional hypersonic flow over a compression ramp. The flow is laminar with a free-stream Mach number of 14.1. In the experiment, the ramp angle was varied to change the strength of the ramp shock and the extent of the viscous-inviscid interaction. Calculations were made for the 24-degree ramp configuration which produces a large separated-flow region that extends upstream of the corner.

Rudy, David H.↗

Methods of testing parameterizations: Vertical ocean mixing

The ocean's velocity field is characterized by an exceptional variety of scales. While the small-scale oceanic turbulence responsible for the vertical mixing in the ocean is of scales a few centimeters and smaller, the oceanic general circulation is characterized by horizontal scales of thousands of kilometers. In oceanic general circulation models that are typically run today, the vertical structure of the ocean is represented by a few tens of discrete grid points. Such models cannot explicitly model the small-scale mixing processes, and must, therefore, find ways to parameterize them in terms of the larger-scale fields. Finding a parameterization that is both reliable and plausible to use in ocean models is not a simple task. Vertical mixing in the ocean is the combined result of many complex processes, and, in fact, mixing is one of the less known and less understood aspects of the oceanic circulation. In present models of the oceanic circulation, the many complex processes responsible for vertical mixing are often parameterized in an oversimplified manner. Yet, finding an adequate parameterization of vertical ocean mixing is crucial to the successful application of ocean models to climate studies. The results of general circulation models for quantities that are of particular interest to climate studies, such as the meridional heat flux carried by the ocean, are quite sensitive to the strength of the vertical mixing. We try to examine the difficulties in choosing an appropriate vertical mixing parameterization, and the methods that are available for validating different parameterizations by comparing model results to oceanographic data. First, some of the physical processes responsible for vertically mixing the ocean are briefly mentioned, and some possible approaches to the parameterization of these processes in oceanographic general circulation models are described in the following section. We then discuss the role of the vertical mixing in the physics of the large-scale ocean circulation, and examine methods of validating mixing parameterizations using large-scale ocean models.

Tziperman, Eli↗

Comparison of Four Precipitation Forcing Datasets in Land Information System Simulations over the Continental U.S.

The NASA Short ]term Prediction Research and Transition (SPoRT) Center in Huntsville, AL is running a real ]time configuration of the NASA Land Information System (LIS) with the Noah land surface model (LSM). Output from the SPoRT ]LIS run is used to initialize land surface variables for local modeling applications at select National Weather Service (NWS) partner offices, and can be displayed in decision support systems for situational awareness and drought monitoring. The SPoRT ]LIS is run over a domain covering the southern and eastern United States, fully nested within the National Centers for Environmental Prediction Stage IV precipitation analysis grid, which provides precipitation forcing to the offline LIS ]Noah runs. The SPoRT Center seeks to expand the real ]time LIS domain to the entire Continental U.S. (CONUS); however, geographical limitations with the Stage IV analysis product have inhibited this expansion. Therefore, a goal of this study is to test alternative precipitation forcing datasets that can enable the LIS expansion by improving upon the current geographical limitations of the Stage IV product. The four precipitation forcing datasets that are inter ]compared on a 4 ]km resolution CONUS domain include the Stage IV, an experimental GOES quantitative precipitation estimate (QPE) from NESDIS/STAR, the National Mosaic and QPE (NMQ) product from the National Severe Storms Laboratory, and the North American Land Data Assimilation System phase 2 (NLDAS ]2) analyses. The NLDAS ]2 dataset is used as the control run, with each of the other three datasets considered experimental runs compared against the control. The regional strengths, weaknesses, and biases of each precipitation analysis are identified relative to the NLDAS ]2 control in terms of accumulated precipitation pattern and amount, and the impacts on the subsequent LSM spin ]up simulations. The ultimate goal is to identify an alternative precipitation forcing dataset that can best support an expansion of the real ]time SPoRT ]LIS to a domain covering the entire CONUS.

Case, Jonathan L.↗

An Overview of the NCC Spray/Monte-Carlo-PDF Computations

This paper advances the state-of-the-art in spray computations with some of our recent contributions involving scalar Monte Carlo PDF (Probability Density Function), unstructured grids and parallel computing. It provides a complete overview of the scalar Monte Carlo PDF and Lagrangian spray computer codes developed for application with unstructured grids and parallel computing. Detailed comparisons for the case of a reacting non-swirling spray clearly highlight the important role that chemistry/turbulence interactions play in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that some of the combustion occurs in a predominantly premixed-flame environment and the rest in a predominantly diffusion-flame environment. However, the non-PDF solution predicts wrongly for the combustion to occur in a vaporization-controlled regime. Near the premixed flame, the Monte Carlo particle temperature distribution shows two distinct peaks: one centered around the flame temperature and the other around the surrounding-gas temperature. Near the diffusion flame, the Monte Carlo particle temperature distribution shows a single peak. In both cases, the computed PDF's shape and strength are found to vary substantially depending upon the proximity to the flame surface. The results bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations. Finally, we end the paper by demonstrating the computational viability of the present solution procedure for its use in 3D combustor calculations by summarizing the results of a 3D test case with periodic boundary conditions. For the 3D case, the parallel performance of all the three solvers (CFD, PDF, and spray) has been found to be good when the computations were performed on a 24-processor SGI Origin work-station.

Raju, M. S.↗

Statistical equilibrium calculations for silicon in early-type model stellar atmospheres

Line profiles of 36 multiplets of silicon (Si) II, III, and IV were computed for a grid of model atmospheres covering the range from 15,000 to 35,000 K in effective temperature and 2.5 to 4.5 in log (gravity). The computations involved simultaneous solution of the steady-state statistical equilibrium equations for the populations and of the equation of radiative transfer in the lines. The variables were linearized, and successive corrections were computed until a minimal accuracy of 1/1000 in the line intensities was reached. The common assumption of local thermodynamic equilibrium (LTE) was dropped. The model atmospheres used also were computed by non-LTE methods. Some effects that were incorporated into the calculations were the depression of the continuum by free electrons, hydrogen and ionized helium line blocking, and auto-ionization and dielectronic recombination, which later were found to be insignificant. Use of radiation damping and detailed electron (quadratic Stark) damping constants had small but significant effects on the strong resonance lines of Si III and IV. For weak and intermediate-strength lines, large differences with respect to LTE computations, the results of which are also presented, were found in line shapes and strengths. For the strong lines the differences are generally small, except for the models at the hot, low-gravity extreme of our range. These computations should be useful in the interpretation of the spectra of stars in the spectral range B0-B5, luminosity classes III, IV, and V.

Kamp, L. W.↗

Quantum Computing in Next-Generation Transportation Optimization

We explore how quantum computing (QC) can advance transportation optimization, with a focus on two high-impact areas: traffic signal control and vehicle electrification with grid integration. As transportation systems grow in complexity, classical optimization methods increasingly struggle to deliver scalable and efficient solutions, particularly for real-time, data-rich environments. This work identifies key challenges within these two domains where QC may offer advantages, particularly in handling combinatorial decision spaces and dynamic constraints. We begin by outlining the limitations of classical approaches for traffic signal control optimization and electric vehicle charging coordination, highlighting where computational limitations arise. Previous quantum formulations are presented and new formulations are proposed to demonstrate how emerging quantum algorithms, including quantum annealing and the Quantum Approximation Optimization Algorithm, could be leveraged to reformulate and address these problems. We also evaluate the suitability of current quantum hardware and discuss recent trends that indicate when QC may become a viable tool for transportation applications. While acknowledging the present limitations of QC technologies, this poster emphasizes the importance of preparing quantum-compatible models today. By reviewing and establishing formulations that align with the strengths of quantum algorithms, researchers and practitioners can better position themselves to take advantage of QC advancements as they occur. This work aims to provide a practical, forward-looking perspective on the near-term potential of quantum computing in transportation optimization.

33 ADVANCED PROPULSION SYSTEMS↗

Improved Estimates of Pentad Precipitation through the Merging of Independent Precipitation Datasets

Three independent, quasi-global, gridded datasets of precipitation (a rain gauge-based dataset, the satellite-only component of the NASA Integrated Multi-satellitE Retrievals for Global Precipitation Measurement mission [IMERG] Final Run precipitation product, and precipitation estimates derived from NASA Soil Moisture Active Passive [SMAP] soil moisture retrievals), are objectively combined into a single pentad precipitation dataset at 36-km resolution using a unique approach based on extended triple collocation. The quality of each of the four datasets is then evaluated against independent observations. When a global land surface model at 36-km resolution is integrated four times, once utilizing the merged precipitation forcing and once with each of the three contributing datasets, the near-surface soil moisture variations produced with the merged forcing validate best against independent satellite-based soil moisture fields. In addition, the merged dataset is found to be more consistent, relative to each contributor, with estimates of air temperature variations across the globe. The merged dataset thus appears to draw successfully on the complementary strengths of each contributor: the particularly high quality of the rain gauge-based dataset in areas of high gauge density, the more uniform accuracy across the globe of the IMERG data, and the moderate accuracy, particularly in semi-arid regions, of the soil moisture retrieval-based data. Plain Language Summary Obtaining measurements of precipitation across the globe can be challenging. Rain gauges in some ways provide the most accurate measurements, but gauges are absent in many parts of the world, and even where they exist, they only measure precipitation at the gauge itself and therefore may not provide an accurate large-scale average. Satellite-based estimates of precipitation largely overcome these problems, but such data have their own issues, notably a “snapshot” (rather than a time-average) character of the measurements and difficulty associated with interpreting the measured radiances in the presence of complex land surfaces. In the present paper, we use a novel approach to generate a “merged” dataset, one that optimally combines the gauge precipitation information and the satellite-based precipitation information with a third set of estimates derived from soil moisture retrievals. The merged precipitation dataset and each of the three contributors (aggregated here to 5-day averages at a spatial resolution of about 36-km) are then evaluated for consistency with independent geophysical fields. The merged dataset is found to perform best, a clear indication that it takes proper advantage of the complementary strengths of each contributor and, accordingly, that the presented approach for merging the different contributors is indeed viable.

Precipitation↗

An ML-based terrestrial data fusion and augmentation framework to enable advanced understanding of the terrestrial carbon and water interactions

Soil moisture is essential to the terrestrial carbon and water cycles and land–atmosphere interactions. There are various types of soil moisture data, and each type has the distinct spatiotemporal strengths and limitations, depending on the diverse applications and retrieval methodologies of different data types (Li et al., in review; The PNNL-82151 FY23 Report). However, the limitations of different soil moisture data in terms of accuracy and spatiotemporal coverage hinder our ability to further understand the soil moisture dynamics across scales. To have a gap free soil moisture data product with a fine spatiotemporal coverage and vertical profiles, we train extreme gradient boosting (XGBoost) models by using (1) in-situ soil moisture measurements from the International Soil Moisture Network (ISMN), (2) soil moisture from the ECMWF reanalysis (ERA) at the 9 km and sub-daily spatiotemporal resolution, (3) the Daymet meteorological fields, and (4) data products that characterize surface conditions, including soil texture, organic content, topography, vegetation type, and rooting depth. We use the trained XGBoost models that have consistent performance across seven soil layers, i.e., 0–5 cm, 5–10 cm, 10–20 cm, 20–40 cm, 40–60 cm, 60–100 cm, and 100–200 cm, and the gridded model predictors to generate a soil moisture data at the 1 km and daily spatiotemporal resolution for the Continental United States (CONUS) from 2001–2020. This dataset can be broadly used for Earth system model benchmark, monitoring extreme weathers, making informed decisions regarding agriculture, water resource management, climate change mitigation, and ecosystem preservation.

58 GEOSCIENCES↗

Biomaterials Out of Thin Air: in Situ, On-Demand Printing of Advanced Biocomposites

Upmass is the single most significant limitation of our current space mission capability. Although biomaterials and biocomposites have mass, strength, flexibility, and self-healing properties that could significantly reduce upmass, their use is limited by the following drawbacks: Expensive, specific production. Many biomaterials can only be produced as part of significant support ecosystem; Inaccessible functional customization. The grain of wood, the porosity of bone, and so on are an integral part of the materials' desired mechanical properties, but are not deterministic when the material is naturally grown; Limited compositions. Most biomaterials (unlike metal, plastic, etc.) cannot be easily combined or modified to produce new materials. This project builds on recent advances in: Synthetic biology. Libraries of standardized genetic parts which can be used for controlled cellular material production, delivery, and binding; 3D printing. Commercial off-the-shelf components which can be used to make of a pico- to nanoliter cell deposition system; Tissue engineering. Proven cell-compatible support hydrogels and scaffolds can be modified to bind the deposited biomaterials of interest. Objectives: Feasibility and benefit analysis. Two mission contexts span the concept's scope (see below); Proof-of-concept demonstration. A simple grid of two proteins, fluorescent for easy detection, to validate the core technology concept; Proposed implementations for follow-on work. Avenues for future work on each core component (host cell, production control, material delivery, material binding, etc.); Complementary studies exploration. A survey of other emerging areas (in situ resource utilization, protein engineering, etc.) with the potential to multiply our technology's impact. Potential Impacts: This application could dramatically expand manufacturing capabilities on Earth and in space: In situ resource utilization. A far greater range of materials and products will be available from the limited palette offered by in situ resource extraction techniques; Reduced equipment and material upmass for off-Earth habitats. Ready- to-use highly specialized construction materials (radiation hardened, compressive/tensile, light or dense) from an extremely low starting mass; Structured biomaterial production. New ready-to-use macro, micro, and molecular manufacturing techniques for traditional biomaterials such as wood, bone and shell; New and novel biocomposite creation. The ability to create completely novel material composites from any base material that cells can be engineered to produce. Suggested Mission Contexts: ISS part manufacturing. A 'minimal working example' making a finished biomaterial part aboard the International Space Station; A long-term Mars habitat. 'Cradle-to-grave' use at a hypothetical Mars habitat, covering everything from tools to construction materials. Alternate Abstract: Imagine being able to print anything from tools and composite building materials to food and human tissues. Imagine being on Mars with the ability to replace any broken part, whether it's a part of your spacesuit, your habitat, or your own body. We propose a technique that would allow just that. By printing 3D arrays of cells engineered to secrete the necessary materials, the abundant in situ resources of atmosphere and regolith become organic, inorganic, or organic-inorganic composite materials. Such materials include novel, biologically derived materials not previously possible to fabricate.

habitats↗

Biomaterials Out of Thin Air: in Situ, On-demand Printing of Advanced Biocomposites

Upmass is the single most significant limitation of our current space mission capability. Although biomaterials and biocomposites have mass, strength, flexibility, and self-healing properties that could significantly reduce upmass, their use is limited by the following drawbacks: Expensive, specific production. Many biomaterials can only be produced as part of significant support ecosystem; Inaccessible functional customization. The grain of wood, the porosity of bone, and so on are an integral part of the materials' desired mechanical properties, but are not deterministic when the material is naturally grown; Limited compositions. Most biomaterials (unlike metal, plastic, etc.) cannot be easily combined or modified to produce new materials. This project builds on recent advances in: Synthetic biology. Libraries of standardized genetic parts which can be used for controlled cellular material production, delivery, and binding; 3D printing. Commercial off-the-shelf components which can be used to make of a pico- to nanoliter cell deposition system; Tissue engineering. Proven cell-compatible support hydrogels and scaffolds can be modified to bind the deposited biomaterials of interest. Objectives: Feasibility and benefit analysis. Two mission contexts span the concept's scope (see below); Proof-of-concept demonstration. A simple grid of two proteins, fluorescent for easy detection, to validate the core technology concept; Proposed implementations for follow-on work. Avenues for future work on each core component (host cell, production control, material delivery, material binding, etc.); Complementary studies exploration. A survey of other emerging areas (in situ resource utilization, protein engineering, etc.) with the potential to multiply our technology's impact. Potential Impacts: This application could dramatically expand manufacturing capabilities on Earth and in space: In situ resource utilization. A far greater range of materials and products will be available from the limited palette offered by in situ resource extraction techniques; Reduced equipment and material upmass for off-Earth habitats. Ready- to-use highly specialized construction materials (radiation hardened, compressive/tensile, light or dense) from an extremely low starting mass; Structured biomaterial production. New ready-to-use macro, micro, and molecular manufacturing techniques for traditional biomaterials such as wood, bone and shell; New and novel biocomposite creation. The ability to create completely novel material composites from any base material that cells can be engineered to produce. Suggested Mission Contexts: ISS part manufacturing. A 'minimal working example' making a finished biomaterial part aboard the International Space Station; A long-term Mars habitat. 'Cradle-to-grave' use at a hypothetical Mars habitat, covering everything from tools to construction materials. Alternate Abstract: Imagine being able to print anything from tools and composite building materials to food and human tissues. Imagine being on Mars with the ability to replace any broken part, whether it's a part of your spacesuit, your habitat, or your own body. We propose a technique that would allow just that. By printing 3D arrays of cells engineered to secrete the necessary materials, the abundant in situ resources of atmosphere and regolith become organic, inorganic, or organic-inorganic composite materials. Such materials include novel, biologically derived materials not previously possible to fabricate.

habitats↗

Can High-resolution WRF Simulations Be Used for Short-term Forecasting of Lightning?

A number of research teams have begun to make quasi-operational forecast simulations at high resolution with models such as the Weather Research and Forecast (WRF) model. These model runs have used horizontal meshes of 2-4 km grid spacing, and thus resolved convective storms explicitly. In the light of recent global satellite-based observational studies that reveal robust relationships between total lightning flash rates and integrated amounts of precipitation-size ice hydrometeors in storms, it is natural to inquire about the capabilities of these convection-resolving models in representing the ice hydrometeor fields faithfully. If they do, this might make operational short-term forecasts of lightning activity feasible. We examine high-resolution WRF simulations from several Southeastern cases for which either NLDN or LMA lightning data were available. All the WRF runs use a standard microphysics package that depicts only three ice species, cloud ice, snow and graupel. The realism of the WRF simulations is examined by comparisons with both lightning and radar observations and with additional even higher-resolution cloud-resolving model runs. Preliminary findings are encouraging in that they suggest that WRF often makes convective storms of the proper size in approximately the right location, but they also indicate that higher resolution and better hydrometeor microphysics would be helpful in improving the realism of the updraft strengths, reflectivity and ice hydrometeor fields.

Goodman, S. J.↗

Does Terrestrial Drought Explain Global CO2 Flux Anomalies Induced by El Nino?

The El Nino Southern Oscillation is the dominant year-to-year mode of global climate variability. El Nino effects on terrestrial carbon cycling are mediated by associated climate anomalies, primarily drought, influencing fire emissions and biotic net ecosystem exchange (NEE). Here we evaluate whether El Nino produces a consistent response from the global carbon cycle. We apply a novel bottom-up approach to estimating global NEE anomalies based on FLUXNET data using land cover maps and weather reanalysis. We analyze 13 years (1997-2009) of globally gridded observational NEE anomalies derived from eddy covariance flux data, remotely-sensed fire emissions at the monthly time step, and NEE estimated from an atmospheric transport inversion. We evaluate the overall consistency of biospheric response to El Nino and, more generally, the link between global CO2 flux anomalies and El Nino-induced drought. Our findings, which are robust relative to uncertainty in both methods and time-lags in response, indicate that each event has a different spatial signature with only limited spatial coherence in Amazonia, Australia and southern Africa. For most regions, the sign of response changed across El Nino events. Biotic NEE anomalies, across 5 El Nino events, ranged from -1.34 to +0.98 Pg Cyr(exp -1, whereas fire emissions anomalies were generally smaller in magnitude (ranging from -0.49 to +0.53 Pg C yr(exp -1). Overall drought does not appear to impose consistent terrestrial CO2 flux anomalies during El Ninos, finding large variation in globally integrated responses from 11.15 to +0.49 Pg Cyr(exp -1). Despite the significant correlation between the CO2 flux and El Nino indices, we find that El Nino events have, when globally integrated, both enhanced and weakened terrestrial sink strength, with no consistent response across events

Schwalm. C. R.↗

Three-Dimensional High-Order Spectral Volume Method for Solving Maxwell's Equations on Unstructured Grids

A three-dimensional, high-order, conservative, and efficient discontinuous spectral volume (SV) method for the solutions of Maxwell's equations on unstructured grids is presented. The concept of discontinuous 2nd high-order loca1 representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) method, but instead of using a Galerkin finite-element formulation, the SV method is based on a finite-volume approach to attain a simpler formulation. Conventional unstructured finite-volume methods require data reconstruction based on the least-squares formulation using neighboring cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In the SV method, one starts with a relatively coarse grid of triangles or tetrahedra, called spectral volumes (SVs), and partition each SV into a number of structured subcells, called control volumes (CVs), that support a polynomial expansion of a desired degree of precision. The unknowns are cell averages over CVs. If all the SVs are partitioned in a geometrically similar manner, the reconstruction becomes universal as a weighted sum of unknowns, and only a few universal coefficients need to be stored for the surface integrals over CV faces. Since the solution is discontinuous across the SV boundaries, a Riemann solver is thus necessary to maintain conservation. In the paper, multi-parameter and symmetric SV partitions, up to quartic for triangle and cubic for tetrahedron, are first presented. The corresponding weight coefficients for CV face integrals in terms of CV cell averages for each partition are analytically determined. These discretization formulas are then applied to the integral form of the Maxwell equations. All numerical procedures for outer boundary, material interface, zonal interface, and interior SV face are unified with a single characteristic formulation. The load balancing in a massive parallel computing environment is therefore easier to achieve. A parameter is introduced in the Riemann solver to control the strength of the smoothing term. Important aspects of the data structure and its effects to communication and the optimum use of cache memory are discussed. Results will be presented for plane TE and TM waves incident on a perfectly conducting cylinder for up to fifth order of accuracy, and a plane wave incident on a perfectly conducting sphere for up to fourth order of accuracy. Comparisons are made with exact solutions for these cases.

Liu, Yen↗

Estimating Evapotranspiration Using an Observation Based Terrestrial Water Budget

Evapotranspiration (ET) is difficult to measure at the scales of climate models and climate variability. While satellite retrieval algorithms do exist, their accuracy is limited by the sparseness of in situ observations available for calibration and validation, which themselves may be unrepresentative of 500m and larger scale satellite footprints and grid pixels. Here, we use a combination of satellite and ground-based observations to close the water budgets of seven continental scale river basins (Mackenzie, Fraser, Nelson, Mississippi, Tocantins, Danube, and Ubangi), estimating mean ET as a residual. For any river basin, ET must equal total precipitation minus net runoff minus the change in total terrestrial water storage (TWS), in order for mass to be conserved. We make use of precipitation from two global observation-based products, archived runoff data, and TWS changes from the Gravity Recovery and Climate Experiment satellite mission. We demonstrate that while uncertainty in the water budget-based estimates of monthly ET is often too large for those estimates to be useful, the uncertainty in the mean annual cycle is small enough that it is practical for evaluating other ET products. Here, we evaluate five land surface model simulations, two operational atmospheric analyses, and a recent global reanalysis product based on our results. An important outcome is that the water budget-based ET time series in two tropical river basins, one in Brazil and the other in central Africa, exhibit a weak annual cycle, which may help to resolve debate about the strength of the annual cycle of ET in such regions and how ET is constrained throughout the year. The methods described will be useful for water and energy budget studies, weather and climate model assessments, and satellite-based ET retrieval optimization.

Rodell, Matthew↗

Computational Algorithms for Unit Commitment with AC Power Flows (Final Report)

Security-constrained unit commitment (SCUC) is a key component in power system operations. When AC power flow constraints are considered in the SCUC model (AC-SCUC), the problem becomes extremely difficult due to its discrete and non-convex nature, as described in “Grid Optimization Competition Challenge 3 Problem Formulation (GOCC)”. There are four main challenges: (i) Discrete decisions regarding unit online/offline status and start-up/shut-down procedures for every single unit. The number of discrete decision variables increases considerably when a system integrates multiple generators; (ii) Configuration-based combined-cycle formulations, and multi-commodity models that include ramping products, spin/non-spin products, and regulation up/down products. The combined-cycle units introduce additional discrete decision variables and auxiliary service products further complicate the model by connecting multi-commodity products’ continuous and discrete variables; (iii) SCUC models with AC power flow constraints are far more complex due to massive bilinear terms in the large-scale nonlinear power balance equations. The nonlinear power balance equations are further complicated by the discrete step control variables of shunts; (iv) N − 1 contingency analysis. The size of the model increases linearly with the number of contingencies considered, greatly increasing the size of the optimization model. Accordingly, there is an emergent need to develop a robust algorithm capable of deriving a high-quality solution in a short time and passing through contingency tests simultaneously. In this project, we explore innovative techniques to address this challenging problem by integrating advanced polyhedral theory, approximation methods, relaxation strategies, decomposition techniques, and parallel computing. Each technique approaches the problem from a different perspective, leveraging its specific strengths to tackle distinct challenges. Each individual method has demonstrated its effectiveness in the PI’s previous research. Their integration is expected to significantly reduce the computational time required to solve the proposed complex problem. Successful completion of this project has the potential to transform the industry by enhancing optimization solvers capable of handling large-scale day-ahead energy market clearing models within strict time constraints, while incorporating AC power flow constraints. This advancement will lead to reduced overall generation costs and, consequently, increased social welfare.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A Hybrid Line List for CH4 and Hot Methane Continuum

Molecular line lists (a catalogue of transition frequencies and line strengths) are important for modelling absorption and emission processes in atmospheres of different astronomical objects, such as cool stars and exoplanets. In order to be applicable for high temperatures, line lists for molecules like methane must contain billions of transitions, which makes their direct (line-by-line) application in radiative transfer calculations impracticable. Here we suggest a new, hybrid line list format to mitigate this problem, based on the idea of temperature-dependent absorption continuum. Methods. The line list is partitioned into a large set of relatively weak lines and a small set of important, stronger lines. The weaker lines are then used either to construct a temperature-dependent (but pressure-independent) set of intensity cross sections or are blended into a greatly reduced set of super-lines. The strong lines are kept in the form of temperature independent Einstein A coefficients. Results. A line list for methane is constructed as a combination of 17 million strong absorption lines relative to the reference absorption spectra and a background methane continuum in two temperature-dependent forms, of cross sections and super-lines. This approach eases the use of large high temperature line lists significantly as the computationally expensive calculation of pressure dependent profiles only need to be performed for a relatively small number of lines. Both the line list and cross sections were generated using a new 34 billion methane line list (34 to10), which extends the 10to10 line list to higher temperatures (up to 2000 K). The new hybrid scheme can be applied to any large line lists containing billions of transitions. We recommend to use super-lines generated on a high resolution grid based on resolving power R = 1,000,000 to model the molecular continuum as a more flexible alternative to the temperature dependent cross sections.

Line Profiles↗