Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local time stepping”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Parallel algorithm for dominant points correspondences in robot binocular stereo vision

This paper presents an algorithm to find the correspondences of points representing dominant feature in robot stereo vision. The algorithm consists of two main steps: dominant point extraction and dominant point matching. In the feature extraction phase, the algorithm utilizes the widely used Moravec Interest Operator and two other operators: the Prewitt Operator and a new operator called Gradient Angle Variance Operator. The Interest Operator in the Moravec algorithm was used to exclude featureless areas and simple edges which are oriented in the vertical, horizontal, and two diagonals. It was incorrectly detecting points on edges which are not on the four main directions (vertical, horizontal, and two diagonals). The new algorithm uses the Prewitt operator to exclude featureless areas, so that the Interest Operator is applied only on the edges to exclude simple edges and to leave interesting points. This modification speeds-up the extraction process by approximately 5 times. The Gradient Angle Variance (GAV), an operator which calculates the variance of the gradient angle in a window around the point under concern, is then applied on the interesting points to exclude the redundant ones and leave the actual dominant ones. The matching phase is performed after the extraction of the dominant points in both stereo images. The matching starts with dominant points in the left image and does a local search, looking for corresponding dominant points in the right image. The search is geometrically constrained the epipolar line of the parallel-axes stereo geometry and the maximum disparity of the application environment. If one dominant point in the right image lies in the search areas, then it is the corresponding point of the reference dominant point in the left image. A parameter provided by the GAV is thresholded and used as a rough similarity measure to select the corresponding dominant point if there is more than one point the search area. The correlation is used as a final decision tool when there is still more than one point in the search area. If there is no dominant point in the search area of if the points in the search area are below a correlation threshold, then the dominant point in the reference image is occluded and can not be corresponded. The algorithm has been modeled, implemented and shown to be fast, robust and parallel. The parallelism is created from three main features: locality of the operators; a memory optimization scheme; and the ability to fully parallelize the extraction phase which is the most computational intensive task in the algorithm. The last feature is achieved by performing the extraction phase on the two images simultaneously.

Al-Tammami, A.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗
DOE OSTI2025

HPC Campaign Management: Remote data access with user-defined error bound using ADIOS and ZFP

Remote access to large-scale scientific datasets, like those generated by combustion simulations or other high-performance computing (HPC) applications, presents a significant challenge. Downloading entire datasets is often impractical due to their size and the bandwidth limitations of typical networks. To address this challenge, we propose a novel approach that enables efficient remote access to large datasets distributed across multiple facilities. Our method enables technologies to download only the data values of a select variable, in a select region of interest, to a user-defined accuracy. For this purpose, we extended the ADIOS IO library to provide read functions with user-defined accuracy, a remote data server that understands multidimensional selections of specific variables, steps and accuracy from an ADIOS dataset, and which uses lossy compression on the remote site to reduce the data to be transferred back to the client. In addition, our extension of the ADIOS library collects metadata from multiple datasets in small files called Campaign Archives, which can be shared among project participants on any HPC, cloud or laptop, and which can easily facilitate the discovery of content and pointers to the data location as well as remote access to the data by local tools as if data was local. This feature called Campaign Management, enables a group of scientists to manage related datasets stored in multiple files, across multiple facilities as if it was in a single file/database. We demonstrate the effectiveness of our approach using a 1.5 TB dataset from the S3D combustion simulation on Frontier at the Oak Ridge Leadership Facility. Even a single variable from this dataset, at 64 GB, is too large to be processed on a standard laptop. We show two different reading patterns for 2D plots and 3D visualization, with careful settings that a scientist studying combustion data would do and show that running the same Python scripts on Frontier directly takes comparable time than running them on the local laptop with remote access to the data on Frontier.

Podhorszki, Norbert [ORNL] (ORCID:000000019647542X↗
NASA NTRS2011

Diagnosing the Sensitivity of Local Land-Atmosphere Coupling via the Soil Moisture-Boundary Layer Interaction

The inherent coupled nature of earth s energy and water cycles places significant importance on the proper representation and diagnosis of land atmosphere (LA) interactions in hydrometeorological prediction models. However, the precise nature of the soil moisture precipitation relationship at the local scale is largely determined by a series of nonlinear processes and feedbacks that are difficult to quantify. To quantify the strength of the local LA coupling (LoCo), this process chain must be considered both in full and as individual components through their relationships and sensitivities. To address this, recent modeling and diagnostic studies have been extended to 1) quantify the processes governing LoCo utilizing the thermodynamic properties of mixing diagrams, and 2) diagnose the sensitivity of coupled systems, including clouds and moist processes, to perturbations in soil moisture. This work employs NASA s Land Information System (LIS) coupled to the Weather Research and Forecasting (WRF) mesoscale model and simulations performed over the U.S. Southern Great Plains. The behavior of different planetary boundary layers (PBL) and land surface scheme couplings in LIS WRF are examined in the context of the evolution of thermodynamic quantities that link the surface soil moisture condition to the PBL regime, clouds, and precipitation. Specifically, the tendency toward saturation in the PBL is quantified by the lifting condensation level (LCL) deficit and addressed as a function of time and space. The sensitivity of the LCL deficit to the soil moisture condition is indicative of the strength of LoCo, where both positive and negative feedbacks can be identified. Overall, this methodology can be applied to any model or observations and is a crucial step toward improved evaluation and quantification of LoCo within models, particularly given the advent of next-generation satellite measurements of PBL and land surface properties along with advances in data assimilation schemes.

Santanello, Joseph A., Jr.↗
DOE OSTI2026

Ultrafast Formation of Jahn–Teller Polarons Revealed by State-Selective Excitation in Correlated Spinel Co 3 O 4

Jahn–Teller polarons are quasiparticles that stem from symmetry breaking and strong local electron–phonon coupling. They originate from an excess charge carrier being dressed by a local lattice distortion, caused by the Jahn–Teller effect, and they critically impact electrical, structural, and magnetic properties in transition metal oxides. The observation of the microscopic steps involved in their formation is essential for enabling control over material properties through the targeted activation of local, site-specific modifications with light pulses. While Jahn–Teller distortion associated with polaron formation was predicted to contribute significantly to changes in electronic band gap and optical properties in Co 3 O 4 , its experimental observation remains elusive, requiring signatures of local symmetry reduction. In this work, we demonstrate Jahn–Teller polaron formation in spinel Co 3 O 4 . By exciting electronic transitions at 3.10 eV and 1.55 eV, we target either the Oh Cobalt(III) or the Td Cobalt(II) ions, and drive the subsequent coherent responses of the system through two different pathways. For the former, we demonstrate that ligand-to-metal charge transfer leads to Jahn–Teller polaron formation, which is linked to the deformation potential and magnetoelastic coupling. For the latter, we identify the coherent excitation of a T 2g phonon mode launched by on-site d-d electronic transitions. Key to our observations is the ability to target site-specific electronic excitations in spinel Co 3 O 4 using ultrafast optical pulses, while monitoring the ensuing low-energy collective modes through the coherent time-domain response of the material. Our approach, which combines the comparative analysis of experimental fingerprints with the support from density functional theory calculations, is broadly applicable to systems in which Jahn–Teller polaron physics has been theoretically predicted but remains experimentally unverified, and it underscores the potential of electronic-state targeting as a route to selectively excite and probe quasiparticle dynamics in solids.

Lattices↗
NASA NTRS2017

Nano Icy Moons Propellant Harvester

As one of just a few bodies identified in the solar system with a liquid ocean, Europa has become a top priority in the search for life outside of Earth. However, cost estimates for exploring Europa have been prohibitively expensive, with estimates of a NASA Flagship class orbiter and lander approaching $5 billion. ExoTerra's NIMPH offers an affordable solution that can not only land, but return a sample from the surface to Earth. NIMPH combines solar electric propulsion (SEP) technologies being developed for the asteroid redirect mission and microsatellite electronics to reduce the cost of a full sample return mission below $500 million. A key to achieving this order-of-magnitude cost reduction is minimizing the initial mass of the system. The cost of any mission is directly proportional to its mass. By keeping the mission within the constraints of an Atlas V 551 launch vehicle versus an SLS, we can significantly reduce launch costs. To achieve this we reduce the landed mass of the sample return lander, which is the largest multiplier of mission mass, and shrink propellant mass through high efficiency SEP and gravity assists. The NIMPH projects first step in reducing landed mass focuses on development of a micro-In Situ Resource Utilization (micro-ISRU) system. ISRU allows us to minimize landed mass of a sample return mission by converting local ice into propellants. The project reduces the ISRU system to a CubeSat-scale package that weighs just 1.74 kg and consumes just 242 W of power. We estimate that use of this ISRU vs. an identical micro-lander without ISRU reduces fuel mass by 45 kg. As the dry mass of the lander grows for larger missions, these savings scale exponentially. Taking full advantage of the micro-ISRU system requires the development of a micro-liquid oxygen-liquid hydrogen engine. The micro-liquid oxygen-liquid hydrogen engine is tailored for the mission by scaling it to match the scale of the micro-lander and the low gravity of the target moon. We also tailor the engine for a near stoichiometric mixture ratio of 7.5. Most high-performance liquid oxygen-liquid hydrogen engines inject extra liquid hydrogen to lower the average molecular weight of the exhaust, which improves specific impulse. However, this extra liquid hydroden requires additional power and processing time on the surface for the ISRU to create. This increases mission cost, and on missions within high radiation environments such as Europa, increases radiation shielding mass. The resulting engine weighs just 1.36 kg and produces 71.5 newton of thrust at 364 s specific impulse. Finally, the mission reduces landed mass by taking advantage of the SEP modules solar power to beam energy to the surface using a collimated laser. This allows us to replace an 45 kg MMRTG with a 2.5 kg resonant array. By using the combination of ISRU, a liquid oxygen-liquid hydrogen engine, and beamed power, we reduce the initial mass of the lander to just 51.5 kg. When combined with an SEP module to ferry the lander to Europa the initial mission mass is just 6397 kg - low enough to be placed on an Earth escape trajectory using an Atlas V 551 launch vehicle. By comparison, we estimate a duplicate lander using an MMRTG and semi-storable propellants such as liquid oxygen-methane would result in an order of magnitude increase in initial lander mass to 445 kg. Attempting to perform the trajectory with a 450 s liquid oxygen-liquid hydrogen engine would increase initial mass to approximately 135,000 kg. Using an Atlas V 1 U.S. Dollar per kg rate to Earth escape value of $27.7k per kg, just the launch savings are over $3.5 billion.

liquid↗
NASA NTRS2021

GRACE Fast Mascons from Spherical Harmonics and a Regularization Design Trade Study

Mass concentration (mascon) solutions have become a prominent medium for investigating time-variable gravity recovered by GRACE and GRACE Follow-On. While GRACE Level-2 spherical harmonic products require various post-processing techniques to eliminate correlated noise, mascon formulations employ spatial regularization strategies during the estimation step to improve signal recovery. However, mascon estimation has traditionally required large computing resources and GRACE Level-1B processing capabilities. In this study, we show that a typical mascon estimation system can be reformulated to allow for the estimation of regularized mascons from Level-2 spherical harmonics. Provided that spherical harmonic solution covariances are available, the computed mascons will be mathematically equivalent to similar mascons estimated from Level-1B observations. This method is computationally efficient, better leverages GRACE spherical harmonics than past methods, and matches the performance of typical mascon solutions without locking scientists into pre-determined regularization designs. We develop a proof-of-concept solution using ITSG-Grace2018 and compare results with traditional mascons from JPL and NASA GSFC. We then assess the effects of spherical harmonic truncation and use of regularization correlations on basin signal recovery. We find that spherical expansions to degree and order 60 provide the minimum expansion necessary to study most basins, while larger expansions help further localize signals. We also find that diagonal regularizations (i.e., regularizations that do not contain inter-mascon correlations) are adversely affected by leakage, especially across boundaries such as coastlines where signals are not highly correlated, whereas including inter-mascon correlations and regional boundaries in the regularization greatly improves signal recovery.

M. J. Croteau↗
DOE OSTI2023

Statistical Learning for Nonlinear Model Reduction from Local Simulations of Stochastic and Particle- and Agent-Based Systems

Stochastic physical systems across the sciences that have very high-dimensional state spaces, with a large number of fast degrees of freedom that force direct simulators to proceed by integration steps that are orders of magnitude smaller than events of interests (e.g., particle collisions). Examples range from molecular motion to dynamics of large populations of cells. A grand challenge in the simulation and understanding of such systems is the systematic construction of accurate, interpretable, reduced models, enabling faster simulations, revealing fundamental properties of the dynamics, and predicting phenomena of interest that the original simulator could not reached with sufficient accuracy or within a given computational budget. In this projected we developed novel statistical estimation/machine learning techniques for analyzing and building empirical reduced models for important families of high-dimensional stochastic systems, in particular: - we developed techniques for estimating interaction kernels in interacting particle- and agent-based systems, which are ubiquitous in Physics, Biology and many other sciences, given observed trajectories of the system; - we developed techniques for nonlinear model reduction for high-dimensional stochastic systems that have a small number of unknown, nonlinear slow variables, and a large number of fast modes, that are possibly of large magnitude, given observed short trajectories of the system in the form of bursts of trajectories from different initial conditions; - we developed novel techniques for estimating linear dynamical systems on graphs when both the dynamics and the underlying graph are unknown, and we have a sparse set of space-time observations; - we considered the problem of estimating an unknown nonlinear observation function of a standard process (e.g. Brownian motion), so that we can recognized if an observed dynamics is "just" a nonlinear version of a known dynamics; we also developed benchmarks for learning algorithms aimed at learning and classifying diffusion processes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1997

Reduced Navier Stokes Relaxation Procedures for Internal Flows

In spite of significant advancement in the field of high speed computing, flow calculations involving complex geometries and/or flow behavior still require large amounts of CPU time and memory. In order to predict such flows without sacrificing grid convergence and accuracy, adaptive gridding techniques that provide optimal resolution are highly desirable. The present work combines multigrid techniques and domain decomposition concepts to provide local, solution adaptive, grid refinement. Several viscous compressible and incompressible, two and three-dimensional, flows with strong inviscid interaction and/or axial flow reversal, are considered with a segmented multigrid domain decomposition (SMGDD) procedure for which uniform meshes result in each domain. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable non-staggered grid computation. The SMGDD technique presented herein has previously been applied for incompressible two dimensional flows. The present work offers improvement in the gridding strategy, by allowing for disjoint subdomains that provide optimal resolution of disparate flow features. It also extends the SMGDD technique to three dimensional compressible flows. Laminar and turbulent flow in a backward facing step channel is considered; although the procedure is applicable to more severe geometries. The standard K-epsilon model is applied for turbulence closure. For Re greater than 400, differences between two-dimensional theory and experiment are resolved through a three dimensional simulation, which confirms the experimentally observed three dimensionality of the recirculation patterns on the upper and lower surfaces.

Rubin, Stanley G.↗
NASA NTRS2011

Understanding Geospace on a Grand Scale: The Global Ionoshphere/Thermosphere Constellation

We present the concept of a constellation of polar orbiting satellites equally spaced in longitude (local time) to systematically sample both the neutral and ionzed gas components of the Earth in circular orbits near 350 km, including their density, temperature, and velocities. The instrumentation would include techniques to measure the height of the ionospheric "F-peak" and its variations along the orbit. The number of satellites (l2? 24? 48?) and their configuration would be determined from modeling analysis and expected geophysical phenomena, including their drivers and characteristic time scales. Together with imaging data from separate satellites, the array of satellites with in situ probes would be expected to provide a new picture of (1) high latitude electrodynamics and atmospheric processes and associated coupling with magnetospheric mass and momentum input, (2) the response of the global ionosphere and thermosphere to magnetic storms, and (3) global neutral wind circulation patterns, neutral density structure, tides, planetary waves, and gravity waves. The comprehensive measurements gathered by the IT-Constellation envisioned here would provide a major leap forward in each of these areas, addressing global physical processes and providing fundamental, new knowledge of Geospace. In particular, by its very nature, the constellation addresses "system science", revealing how the ionosphere-thermosphere connects globally to the magnetosphere above and the troposphere below. We present this concept as the next logical step in observing the "whole" space environment using in situ probes in conjunction with imagers. We invite modelers to not only comment on this concept but also to become actively engaged in helping to define it.

Pfaff, Robert↗
DOE OSTI2026

Adaptive Cybersecurity for Distributed Energy Resources (AdCyDER): Online Reinforcement Learning with Stackelberg-Optimized Defenses — Pipeline Architecture, Evaluation Methodology, and Findings from a Synthetic-Data Evaluation

This report documents the design and evaluation of an integrated online-learning pipeline developed within the AdCyDER project for Distributed Energy Resource (DER) cybersecurity. The pipeline couples a Reinforcement Learning (RL) attack classifier — which produces an attack-type probability distribution — with a Stackelberg game-theoretic (GT) defense selector that consumes those distributions alongside SME-encoded priors over (defense, attack) effectiveness pairings and perdefense costs to choose grid-health-preserving defenses. The objective is not attack classification per se but production of distributions that drive effective defense selection through the Stackelberg layer, learned from delayed grid-health feedback rather than labeled attack data. AdCyDER as a whole is broader than the work presented here; this report covers the specific RL/GT loop integration and its evaluation. We present the integrated pipeline (SCADA telemetry with Fronius inverter physics, Suricata IDS, time-windowed aggregation, per-facility LSTM classifier, Stackelberg optimizer, OpenC2 actuators), an experimental campaign of 28 eight-hour iterations across three baseline modes, and a pipeline-ordered diagnostic protocol. The protocol identifies two distinct failure modes within the loop: paired supervised ceilings on the same features establish that the deployed online RL classifier (macro F1 ≈ 0.07) sits at least 4.7× below a same-architecture supervised LSTM (≈ 0.34) and 10–11× below a linear feature-signal ceiling (≈ 0.70–0.79 depending on per-facility isolation), localizing the dominant failure to the training procedure; and the reward signal driving online updates carries weak directional coupling with classifier correctness in the methodology-expected direction (multi-lens convergent: top-decile P(true) records produce more frequent state changes and slightly larger improvements, top-vs-bot Cohen’s 𝑑 ≈ −0.19), but at effect magnitudes too small to drive gradient-based learning at the campaign sample size. The original learning hypothesis is not supported by the data. The primary contributions are the diagnostic methodology — proposed as a transferable falsification protocol for online RL/GT defense pipelines learning from delayed environmental reward — and the open, reproducible experimental infrastructure. We outline reward reformulation as the highest-priority aspirational next step given the underpowered-but-aligned Q6 reading, with hardware-in-the-loop evaluation as the broadest scope-expansion option.

Blakely, Benjamin [Argonne National Laboratory (AN↗
NASA NTRS2000

High Order Discontinuous Gelerkin Methods for Convection Dominated Problems with Application to Aeroacoustics

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. On the analysis side, we have studied the efficient and stable discontinuous Galerkin framework for small second derivative terms, for example in Navier-Stokes equations, and also for related equations such as the Hamilton-Jacobi equations. This is a truly local discontinuous formulation where derivatives are considered as new variables. On the applied side, we have implemented and tested the efficiency of different approaches numerically. Related issues in high order ENO and WENO finite difference methods and spectral methods have also been investigated. Jointly with Hu, we have presented a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the RungeKutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method. Jointly with Hu, we have constructed third and fourth order WENO schemes on two dimensional unstructured meshes (triangles) in the finite volume formulation. The third order schemes are based on a combination of linear polynomials with nonlinear weights, and the fourth order schemes are based on combination of quadratic polynomials with nonlinear weights. We have addressed several difficult issues associated with high order WENO schemes on unstructured mesh, including the choice of linear and nonlinear weights, what to do with negative weights, etc. Numerical examples are shown to demonstrate the accuracies and robustness of the methods for shock calculations. Jointly with P. Montarnal, we have used a recently developed energy relaxation theory by Coquel and Perthame and high order weighted essentially non-oscillatory (WENO) schemes to simulate the Euler equations of real gas. The main idea is an energy decomposition under the form epsilon = epsilon(sub 1) + epsilon(sub 2), where epsilon(sub 1) is associated with a simpler pressure law (gamma)-law in this paper) and the nonlinear deviation epsilon(sub 2) is convected with the flow. A relaxation process is performed for each time step to ensure that the original pressure law is satisfied. The necessary characteristic decomposition for the high order WENO schemes is performed on the characteristic fields based on the epsilon(sub l) gamma-law. The algorithm only calls for the original pressure law once per grid point per time step, without the need to compute its derivatives or any Riemann solvers. Both one and two dimensional numerical examples are shown to illustrate the effectiveness of this approach.

Shu, Chi-Wang↗
NASA NTRS2015

On the Feasibility of Monitoring Carbon Monoxide in the Lower Troposphere from a Constellation of Northern Hemisphere Geostationary Satellites (PART 1)

By the end of the current decade, there are plans to deploy several geostationary Earth orbit (GEO) satellite missions for atmospheric composition over North America, East Asia and Europe with additional missions proposed. Together, these present the possibility of a constellation of geostationary platforms to achieve continuous time-resolved high-density observations over continental domains for mapping pollutant sources and variability at diurnal and local scales. In this paper, we use a novel approach to sample a very high global resolution model (GEOS-5 at 7 km horizontal resolution) to produce a dataset of synthetic carbon monoxide pollution observations representative of those potentially obtainable from a GEO satellite constellation with predicted measurement sensitivities based on current remote sensing capabilities. Part 1 of this study focuses on the production of simulated synthetic measurements for air quality OSSEs (Observing System Simulation Experiments). We simulate carbon monoxide nadir retrievals using a technique that provides realistic measurements with very low computational cost. We discuss the sampling methodology: the projection of footprints and areas of regard for geostationary geometries over each of the North America, East Asia and Europe regions; the regression method to simulate measurement sensitivity; and the measurement error simulation. A detailed analysis of the simulated observation sensitivity is performed, and limitations of the method are discussed. We also describe impacts from clouds, showing that the efficiency of an instrument making atmospheric composition measurements on a geostationary platform is dependent on the dominant weather regime over a given region and the pixel size resolution. These results demonstrate the viability of the "instrument simulator" step for an OSSE to assess the performance of a constellation of geostationary satellites for air quality measurements.

GEOS-5↗
DOE OSTI2020

Solar Resource and Infrastructure Assessment for the Town of Blandford

This report is a solar resource and infrastructure assessment for the town of Blandford, Massachusetts. The assessment was funded through the National Renewable Energy Laboratory, Solar Energy Innovation Network (NREL SEIN) Solar in Rural Communities Program, as part of a project to develop a Community-Informed Proactive Solar Siting and Financing Model. As a first step, the project lead organization, UMass Clean Energy Extension prepared an assessment of existing infrastructure, resources, and potential solar development opportunities in participating municipalities, including Blandford. This assessment was designed to describe relevant bylaws and infrastructure within the town, identify the types of solar facilities that could be developed, and quantify the total space available for each type of facility. In this report, we reviewed existing electricity grid infrastructure, and the potential to interconnect additional solar facilities. At the present time, most distribution lines providing electricity to Blandford are over-saturated with authorized and proposed solar projects, and cannot accommodate additional solar projects to interconnect to the grid. We can expect significant upgrades to these circuits, if any of the in-process projects are to proceed, which might then free up additional capacity for new projects. There is one circuit which serves a section of North Blandford Road that is not over-saturated; for the immediate future, this is the most cost-effective location for new large-scale projects to be sited. Meanwhile, most three-phase lines could likely accommodate additional small-to-medium scale projects (under 200 kW), and most single-phase lines could likely accommodate additional projects under 50 kW in size. This description represents the local grid infrastructure as it is - planning for future scenarios of development could include recommendations for areas of grid infrastructure improvement to allow siting of distributed generation in preferred locations. Future scenarios may also include the addition energy storage and other "non-wires alternatives."

14 SOLAR ENERGY↗
DOE OSTI2021

Deployment of Community Scale Solar Photovoltaic Energy Systems in Forest County Potawatomi Community

FCPC is a federally recognized tribe with a membership of more than 1400. The Tribe has a reservation in Forest County, Wisconsin, and tribal trust lands in Milwaukee County, Wisconsin, with significant governmental and business operations in both areas. Until relatively recently, FCPC was one of the poorest tribes in the Midwest. Due to its successful gaming operations, it now has improved its financial situation and is using its resources to benefit both its members and the general public through environmental and other community efforts, many of which are long-delayed because of the extended period of abject poverty. The Tribe has developed the long-term goal of becoming energy independent using renewable resources. In order to meet this goal, the Tribe has taken several important steps. As an interim step, FCPC is purchasing nearly 55 million kilowatt hours of renewable energy certificates (“RECs”) annually from certified renewable energy facilities in the United States. These RECs offset more than 100% of the present electricity use at all of the Tribe’s facilities. The Tribe is currently ranked number 11 on the Environmental Protection Agency’s (EPA) Top 20 Local Government list of green power purchasers and is a member of the EPA’s Green Power Partnership. The original project suggested small scale installations on tribal homes and certain government buildings. The Tribe suggested that the project would result in 730 kw installed. By revising the project, the Tribe was able to install 938 kw and owns the largest installation in Northern Wisconsin with the 800 kW installation at the Potawatomi Carter Casino Hotel (“PCCH”). The PCCH installation, in summer months, can provide near 90% of the use during daylight hours. Although tracking was not immediately available FCPC has since confirmed significant production with certain quarterly savings nearing $210,000. These installations provide a reduction in facility costs providing more revenue for Tribal services including cultural, health, sustainable food sources, environmental preservation and other Tribal initiatives. Further, the installations have resulted in positive steps towards the Tribe’s goals of energy sovereignty and the use of 100% renewable clean energy. Finally, the PCCH project, specifically, has had the additional benefit of providing a likely area for a microgrid. The size of the installation and offset achieved have pushed the Tribe toward planning for storage and a microgrid system that would both further the Tribe’s energy goals and provide a safe haven for the Community during times of need.

14 SOLAR ENERGY↗
DOE OSTI2022

Title Coherent X-ray Studies of Surface Growth and Patterning Processes

X-ray Photon Correlation Spectroscopy (XPCS) is being developed as a tool to study nanoscale dynamics of fluctuations during thin film growth and surface patterning. XPCS examines the evolution of the X-ray scattering speckle pattern in reciprocal space to reveal dynamics information not accessible through any other means. Unlike low-coherence conventional X-ray scattering, which incoherently averages over different regions of a sample, coherent X-ray scattering is sensitive to the detailed structure of a given sample at that moment in time. In thin film growth processes, heterodyning, which occurs due to coherent mixing of two scattered signals, can be used to investigate the relationship between surface growth velocity and defect propagation. And, for polycrystalline thin film growth, the spatial coherence of the X-ray beam can substitute for the missing spatial coherence of the growth process to track layered growth in detail, even in a growth regime in which there are no conventional growth oscillations of the X-ray intensity. This has opened the door for detailed dynamics studies of individual atomic layers during real-world growth and patterning, not just in perfect single-crystal growth cases, which greatly expands the applicability of in-situ X-ray scattering methods. Step-flow dynamics in mounded polycrystalline growth has been studied separately for two thin film organic semiconductors deposited by thermal deposition in a vacuum environment, C60 and diindenoperylene (DIP). Highly oriented polycrystalline thin films are readily obtained in both systems, where mounds are composed of crystalline monolayer-height steps and terraces in a so-called wedding cake morphology. The formation of mounds is understood to be due to significant Ehrlich-Schwoebel step edge barriers that inhibit molecules from hopping down from one layer to the one below. The mounds exhibit local step flow, which can be monitored using coherent X-ray scattering. This is made possible due to heterodyning between scattering from the average mounds and the moving steps, which becomes visible in XPCS analysis. The effect of desorption, i.e. re-evaporation of deposited molecules is found to be important for understanding these processes. The impact of this work is to enable testing of models of step dynamics, the shape of mounds, and merging of mounds to form continuous thin films. Highly ordered polycrystalline thin films deposited on inexpensive substrates have applications in thin film solar cells and other organic electronic devices. In a separate set of experiments, speckle analysis during self-organized ion beam nanopatterning reveals memory stretching back to the beginning of patterning in the early stages and enables measurement of the velocity of self-organized patterns across surfaces providing the possibility of stringent new tests of the theory of pattern formation and motion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Guide to Updating Interconnection Rules and Incorporating IEEE Standard 1547

The intent of this paper is to support local authorities in adopting the most current technical requirements for interconnecting distributed energy resources (DERs) to the electric distribution system as specified in the Institute of Electrical and Electronics Engineers Standard for Interconnection and Interoperability of Distributed Energy Resources with Associated Electric Power Systems Interfaces (IEEE Std 1547-2018). To create interconnection rules that successfully address policy goals, market trends, and technical requirements, state regulators need to be well informed of DER integration considerations. Interconnection processes and technical issues often overlap; the two are not easily segregated, and they are often considered and evaluated concurrently within the context of interconnection rulemaking. Accordingly, this paper addresses both the technical issues and related process considerations. A lack of uniform and transparent procedures for addressing interconnection rule changes can result in implementation issues or inefficiencies, such as unclear, lengthy, and complicated interconnection rules that can increase distributed generation "soft costs" (i.e., non-hardware costs). This can further delay the deployment of DERs and jeopardize time-constrained national and state policy goals, such as renewable energy targets, or limit participation in emerging techno-economic trends. Additionally, misconstructed interconnection rules are at risk of becoming mired in stakeholder conflicts. This guide addresses the concerns of electric service regulators from both the process and technical standpoints by presenting a structured, step-by-step approach to developing and updating existing interconnection rules. In this document, the process of developing and updating interconnection rules is subdivided into three steps: (1) determining the context (stakeholders and major drivers); (2) developing the rule, including updating technical requirements; and (3) maintaining and revising the rule over time.

14 SOLAR ENERGY↗