Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid singularity”

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.

117 records · Page 7

Contribution of Near Real Time MODIS-Based Forest Disturbance Detection Products to a National Forest Threat Early Warning System

This presentation discusses an effort to compute and post weekly MODIS forest change products for the conterminous US (CONUS), as part of national forest threat early warning system (EWS) known as the U.S. Forest Change Assessment Viewer (FCAV). The US Forest Service, NASA, USGS, and ORNL are working collaboratively to contribute weekly change products to this EWS. Large acreages of the nation's forests are being disturbed by a growing multitude of biotic and abiotic threats that can act either singularly or in combination. When common at regional scales, such disturbances can pose hazards and threats to floral and faunal bio-diversity, ecosystem sustainability, ecosystem services, and human settlements across the conterminous US. Regionally evident forest disturbances range from ephemeral periodic canopy defoliation to stand replacement mortality events due to insects, disease, fire, hurricanes, tornadoes, ice, hail, and drought. Mandated by the Healthy Forest Restoration Act of 2003, this forest threat EWS has been actively developed since 2006 and on-line since 2010. This FCAV system employs 250-meter MODIS NDVI-based forest change products as a key element of the system, providing regional and CONUS scale products in near real time every 8 days. Each forest change product in FCAV is based on current versus historical 24 day composite NDVI data gridded at 231.66 meter resolution. Current NDVI is derived from USGS eMODIS expedited products. MOD13 NDVI is used for constructing historical baselines. CONUS change products are computed for all forests as % change in the current versus historical NDVI. Change products are computed according to previous year, previous 3 years and previous 8 year historical baselines. The use of multiple baselines enables disturbance anomaly phenology to be more fully assessed. CONUS forest change products are posted each week on the FCAV, a web mapping service maintained by the National Environmental Modeling and Analysis Center. The FCAV EWS has been used to aid multiple Federal and State agency forest management activities, including aerial disturbance detection surveys, as well as rapid response preliminary assessments of timber loss due to tornadoes, regional drought studies, and fire damage assessments. The FCAV allows end-users to assess the context of apparent forest vegetation change with respect to ancillary data, such as land cover, topography, hydrology, climate variables, and administrative boundaries. Such change products are being evaluated through case studies involving comparison with higher spatial resolution satellite, aerial, and field data. The presentation will include multiple examples in which regionally evident forest disturbances were successfully detected and monitored with the MODIS-based change products, as part of the FCAV. FCAV's MODIS forest change products enable end-users (e.g., resource managers) to monitor forest hazards at regional scales throughout the year and across the nation.

Spruce, Joseph P.↗

A High-Order Finite Spectral Volume Method for Conservation Laws on Unstructured Grids

A time accurate, high-order, conservative, yet efficient method named Finite Spectral Volume (FSV) is developed for conservation laws on unstructured grids. The concept of a 'spectral volume' is introduced to achieve high-order accuracy in an efficient manner similar to spectral element and multi-domain spectral methods. In addition, each spectral volume is further sub-divided into control volumes (CVs), and cell-averaged data from these control volumes is used to reconstruct a high-order approximation in the spectral volume. Riemann solvers are used to compute the fluxes at spectral volume boundaries. Then cell-averaged state variables in the control volumes are updated independently. Furthermore, TVD (Total Variation Diminishing) and TVB (Total Variation Bounded) limiters are introduced in the FSV method to remove/reduce spurious oscillations near discontinuities. A very desirable feature of the FSV method is that the reconstruction is carried out only once, and analytically, and is the same for all cells of the same type, and that the reconstruction stencil is always non-singular, in contrast to the memory and CPU-intensive reconstruction in a high-order finite volume (FV) method. Discussions are made concerning why the FSV method is significantly more efficient than high-order finite volume and the Discontinuous Galerkin (DG) methods. Fundamental properties of the FSV method are studied and high-order accuracy is demonstrated for several model problems with and without discontinuities.

Wang, Z. J.↗

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗

The Global Modeling and Assimilation Office (GMAO) 4d-Var and its Adjoint-based Tools

The fifth generation of the Goddard Earth Observing System (GEOS-5) Data Assimilation System (DAS) is a 3d-var system that uses the Grid-point Statistical Interpolation (GSI) system developed in collaboration with NCEP, and a general circulation model developed at Goddard, that includes the finite-volume hydrodynamics of GEOS-4 wrapped in the Earth System Modeling Framework and physical packages tuned to provide a reliable hydrological cycle for the integration of the Modern Era Retrospective-analysis for Research and Applications (MERRA). This MERRA system is essentially complete and the next generation GEOS is under intense development. A prototype next generation system is now complete and has been producing preliminary results. This prototype system replaces the GSI-based Incremental Analysis Update procedure with a GSI-based 4d-var which uses the adjoint of the finite-volume hydrodynamics of GEOS-4 together with a vertical diffusing scheme for simplified physics. As part of this development we have kept the GEOS-5 IAU procedure as an option and have added the capability to experiment with a First Guess at the Appropriate Time (FGAT) procedure, thus allowing for at least three modes of running the data assimilation experiments. The prototype system is a large extension of GEOS-5 as it also includes various adjoint-based tools, namely, a forecast sensitivity tool, a singular vector tool, and an observation impact tool, that combines the model sensitivity tool with a GSI-based adjoint tool. These features bring the global data assimilation effort at Goddard up to date with technologies used in data assimilation systems at major meteorological centers elsewhere. Various aspects of the next generation GEOS will be discussed during the presentation at the Workshop, and preliminary results will illustrate the discussion.

Todling, Ricardo↗

Using Pilot Jobs and CernVM File System for Simplified Use of Containers and Software Distribution

High Energy Physics (HEP) experiments entail an abundance of computing resources, i.e. sites, to run simulations and analyses by processing data. This requirement is fulfilled by local batch farms, grid sites, private/commercial clouds, and supercomputing centers via High Throughput Computing (HTC). The growing needs of such experiments and resources being prone to trends of heterogeneity make it difficult for physicists to handle these resources directly. Additionally, HEP collaborations heavily rely on data and software releases, typically in the order of tens of gigabytes, while conducting simulations and analyses. Hence, aspects of scalability, reliability, and maintenance become crucial with regards to the distribution of the necessary data and software stack. The GlideinWMS [4] framework helps with the resource management problem by using pilot jobs, aka Glideins, to provision reliable elastic virtual clusters. Glideins are submitted to unreliable heterogeneous resources which are validated and customized by the Glideins to make the worker nodes available for end-user job execution. On the other hand, the CernVM File System (CernVM-FS or CVMFS) [1] helps with data distribution. It is a write-once, read-everywhere filesystem used to deploy scientific software to thousands of nodes on a worldwide distributed computing infrastructure. CVMFS is based on the Hyper Text Transfer Protocol and has been widely used within the particle physics community for (1) distributing experiment software and data such as calibrations, and (2) facilitating containerization by efficiently hosting container images along with providing containerization software, especially Singularity [3] GlideinWMS relies on CVMFS installed locally on the computing resources to satisfy the experiments' software needs. This requires system administrators' effort to install and maintain CVMFS at the sites and limits the use of sites, especially HPC resources, that do not have CVMFS installed. This poster presents a solution, taking advantage of Glideins to provide CVMFS at most sites without the need for a local installation. Doing so expands the pool of resources available for HEP experiments and reduces the effort of system administrators for current resources. Additionally, the proposed solution allows GlideinWMS to also start Singularity [3], a containerization software that can run unprivileged, on sites where neither CVMFS nor Singularity are available, including HPC sites. The benefits provided by this solution are: (1) lower overhead for site administrators in that they have less software to install, (2) an expanded pool of resources that run user jobs with easy access to software and data provided by CVMFS, thus making life easier for the scientists, and (3) improved flexibility to use HPC resources by enabling GlideinWMS pilot jobs to support HPC sites.

Urs, Namratha↗

Real-space analysis of radiation-induced specific changes with independent component analysis

A method of analysis is presented that allows for the separation of specific radiation-induced changes into distinct components in real space. The method relies on independent component analysis (ICA) and can be effectively applied to electron density maps and other types of maps, provided that they can be represented as sets of numbers on a grid. Here, for glucose isomerase crystals, ICA was used in a proof-of-concept analysis to separate temperature-dependent and temperature-independent components of specific radiation-induced changes for data sets acquired from multiple crystals across multiple temperatures. ICA identified two components, with the temperature-independent component being responsible for the majority of specific radiation-induced changes at temperatures below 130 K. The patterns of specific temperature-independent radiation-induced changes suggest a contribution from the tunnelling of electron holes as a possible explanation. In the second case, where a group of 22 data sets was collected on a single thaumatin crystal, ICA was used in another type of analysis to separate specific radiation-induced effects happening on different exposure-level scales. Here, ICA identified two components of specific radiation-induced changes that likely result from radiation-induced chemical reactions progressing with different rates at different locations in the structure. In addition, ICA unexpectedly identified the radiation-damage state corresponding to reduced disulfide bridges rather than the zero-dose extrapolated state as the highest contrast structure. The application of ICA to the analysis of specific radiation-induced changes in real space and the data pre-processing for ICA that relies on singular value decomposition, which was used previously in data space to validate a two-component physical model of X-ray radiation-induced changes, are discussed in detail. This work lays a foundation for a better understanding of protein-specific radiation chemistries and provides a framework for analysing effects of specific radiation damage in crystallographic and cryo-EM experiments.

Borek, Dominika (ORCID:0000000243216253)↗

Designing for Zero Energy and Zero Carbon on a Multi-Building Scale Using URBANopt: Preprint

Groundbreaking efforts are necessary to mitigate contributors increasing impacts of climate change. In parallel to inventing pioneering clean energy technologies it is even more fundamental to rethink designing energy systems within a singular facility and collectively to function as a district. Facilities should not be continuously passive by just consuming; there is a need to shift to perform more dynamically. Designing for zero energy and zero carbon on a multi-building scale can uncover opportunities for building energy efficiency, decarbonization, demand flexibility, and resiliency that are not accessible at an individual building scale. This approach can be challenging without innovative tools to evaluate the multitude of possibilities. As an investigated result, we highlight the use of a campus-scale energy modeling platform - URBANopt™ - for the expansion of the National Renewable Energy Laboratory's (NREL's) South Table Mountain campus in Golden, Colorado. Programmatic growth included the design of three new all-electric, zero-energy, and zero-carbon, mixed used buildings (a combination of research laboratories and office space). This investigation is critical to NREL reaching net-zero emissions for its operational footprint, which will occur in phases over the next decade. Leveraging URBANopt's capabilities, we evaluate 1) high-performance building energy efficiency and decarbonization measures, 2) 4th generation district heating and cooling (4th GDHC) systems, 3) optimized onsite generation and energy storage assets that meet zero-energy and zero-carbon targets at minimum life-cycle costs, and 4) cost-optimal distributed energy technology mixes, dispatch strategies, and associated capacities that increase resiliency to grid outages. This work demonstrates the use and capabilities of URBANopt through a real-world case study on a multi-building scale.

community energy model↗

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

The Reynolds numbers that characterize geophysical and astrophysical turbulence (Re approximately equals 10(exp 8) for the planetary boundary layer and Re approximately equals 10(exp 14) for the Sun's interior) are too large to allow a direct numerical simulation (DNS) of the fundamental Navier-Stokes and temperature equations. In fact, the spatial number of grid points N approximately Re(exp 9/4) exceeds the computational capability of today's supercomputers. Alternative treatments are the ensemble-time average approach, and/or the volume average approach. Since the first method (Reynolds stress approach) is largely analytical, the resulting turbulence equations entail manageable computational requirements and can thus be linked to a stellar evolutionary code or, in the geophysical case, to general circulation models. In the volume average approach, one carries out a large eddy simulation (LES) which resolves numerically the largest scales, while the unresolved scales must be treated theoretically with a subgrid scale model (SGS). Contrary to the ensemble average approach, the LES+SGS approach has considerable computational requirements. Even if this prevents (for the time being) a LES+SGS model to be linked to stellar or geophysical codes, it is still of the greatest relevance as an 'experimental tool' to be used, inter alia, to improve the parameterizations needed in the ensemble average approach. Such a methodology has been successfully adopted in studies of the convective planetary boundary layer. Experienc e with the LES+SGS approach from different fields has shown that its reliability depends on the healthiness of the SGS model for numerical stability as well as for physical completeness. At present, the most widely used SGS model, the Smagorinsky model, accounts for the effect of the shear induced by the large resolved scales on the unresolved scales but does not account for the effects of buoyancy, anisotropy, rotation, and stable stratification. The latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Canuto, V. M.↗

Restoration Hydro: A Watershed Approach to Standard Modular New Hydropower

The objectives of FOA DE- FOA-0001836- “Standard Modular Hydropower” included designing a standardized, modular, and environmentally compatible hydropower schematic for implementation in greenfield sites that generate up to 10 MW of capacity. Utilizing funds competitively awarded under DOE’s Water Power Technologies Office, the Natel Energy team developed a concept for modular new stream reach (NSR) hydropower that incorporates multi species upstream and downstream fish passage, improved river channel connectivity, and recreational modules. The in-stream design of the collective modules minimized site specific design and maximized the opportunities for modularity. Financial data was also presented using actual costs from regional suppliers, with figures provided in 2022 dollars. While the project team did not address potential permitting process improvements, the site selection criteria did consider established barriers to hydropower development such as tribal and preserved lands, interconnection proximity, and endangered species to exclude or deprioritize. The project’s design schematic met the objectives of the FOA, and presented a unique solution that targets alluvial pockets as natural features for sustainable development. Natel’s concept also incorporated the company's fish-safe Restoration Hydro Turbine for safe downstream passage, while featuring a rock arch that integrates fish passage, water, recreation, and grade control modules (including sediment). Alignment with the Department of Energy Office of Energy Efficiency and Renewable Energy (EERE) “Innovative Design Concepts for Standard Modular Hydropower and Pumped-Storage Hydropower” Program: According to the Hydropower Vision (DOE, 2016), approximately 16 GW of hydropower growth is possible with the development of technology solutions that balance efficiency, economics, and environmental sustainability. The desired outcome of the SMH program is transformational innovation specifically in the site identification, conceptual, and detailed design phases of technology development lifecycles (DOE, 2018). In developing the SMH design schematic, the team aimed to address the opportunities outlined in the Vision through an inverted design philosophy; rather than singularly prioritizing efficiency and power production, the team focused on integrating hydropower with restoration of degraded streams to optimal ecosystem function and provision of exceptional recreation value as design criteria. To achieve this, Restoration Hydro incorporates the principles of nature-based engineering (WWAP, 2018) and biomimicry (Biomimicry NL.) to strategically deploy complementary combinations of permanent, semi-permanent, and ephemeral low-head structures - such as natural and engineered log jams - that harness geomorphological and hydrological processes at the landscape-scale. Primary applications of Restoration Hydro include: 1) restoration of degraded watersheds’ natural ecological function and enhancement of hydrological connectivity; and 2) creation of associated co-benefits to hydro production, including increased groundwater recharge, improved sediment transport and management, improved water security and water quality. Restoration Hydro projects build upon proven watershed restoration engineering techniques by integrating hydropower turbines into low-head structures using innovative and evolving civil works concepts that facilitate fish and sediment passage, and in some cases create additional revenue-generating recreational opportunities. Powering low-head structures creates a directly monetizable layer of economic value in the form of flexible, reliable, renewable energy on top of the already high-value water, environmental and recreational benefits of watershed and river restoration. The approach aims to create a virtuous, self-reinforcing cycle whereby Restoration Hydro projects support the scaling of ecosystem restoration activities, creating a water-energy-carbon multiplier effect that, through the principles of adaptive change management: 1) improves the resilience of landscapes and downstream population centers for changing hydrological cycles; 2) creates a reliable energy resource that facilitates the integration of intermittent renewable power sources into grids; and 3) supports climate change mitigation through grid decarbonization and enhanced ecosystem carbon capture and retention.

13 HYDRO ENERGY↗