Search NASA⌕ Search

SEARCH · Search NASA

Results for “test cases”

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 343 records · Page 19

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

Linear solvers for power grid optimization problems: A review of GPU-accelerated linear solvers

The linear equations that arise in interior methods for constrained optimization are sparse symmetric indefinite, and they become extremely ill-conditioned as the interior method converges. These linear systems present a challenge for existing solver frameworks based on sparse LU or LDL T decompositions. Here, we benchmark five well known direct linear solver packages on CPU- and GPU-based hardware, using matrices extracted from power grid optimization problems. The achieved solution accuracy varies greatly among the packages. None of the tested packages delivers significant GPU acceleration for our test cases. For completeness of the comparison we include results for MA57, which is one of the most efficient and reliable CPU solvers for this class of problem.

97 MATHEMATICS AND COMPUTING↗

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗

Time-series learning of latent-space dynamics for reduced-order model closure

In this work, we study the performance of long short-term memory networks (LSTMs) and neural ordinary differential equations (NODEs) in learning latent-space representations of dynamical equations for an advection-dominated problem given by the viscous Burgers equation. Our formulation is devised in a nonintrusive manner with an equation-free evolution of dynamics in a reduced space with the latter being obtained through a proper orthogonal decomposition. In addition, we leverage the sequential nature of learning for both LSTMs and NODEs to demonstrate their capability for closure in systems that are not completely resolved in the reduced space. We assess our hypothesis for two advection-dominated problems given by the viscous Burgers equation. We observe that both LSTMs and NODEs are able to reproduce the effects of the absent scales for our test cases more effectively than does intrusive dynamics evolution through a Galerkin projection. This result empirically suggests that time-series learning techniques implicitly leverage a memory kernel for coarse-grained system closure as is suggested through the Mori–Zwanzig formalism.

97 MATHEMATICS AND COMPUTING↗

Spatiotemporally dynamic implicit large eddy simulation using machine learning classifiers

In this article, we utilize machine learning to dynamically determine if a point on the computational grid requires implicit numerical dissipation for large eddy simulation (LES). The decision making process is learnt through a priori training on quantities derived from direct numerical simulation (DNS) data. In particular, we compute eddy-viscosities obtained through the coarse-graining of DNS quantities and utilize their projection onto a Gaussian distribution to categorize areas that may require dissipation. If our learning determines that closure is necessary, an upwinded scheme is utilized for computing the non-linear Jacobian. In contrast, if it is determined that closure is unnecessary, a symmetric and second-order accurate energy and enstrophy preserving Arakawa scheme is utilized instead. This results in a closure framework that precludes the specification of any model-form for the small scale contributions of turbulence but deploys an appropriate numerical dissipation from explicit closure driven hypotheses. Here, this methodology is deployed for the Kraichnan turbulence test-case and assessed through various statistical quantities such as angle-averaged kinetic energy spectra and vorticity structure functions. Our framework thus establishes a link between the use of explicit LES ideologies for closure and numerical dissipation-based modeling of turbulence leading to improved statistical fidelity of a posteriori simulations.

42 ENGINEERING↗

Influence of fertilization on the dynamics of energy use in wheat

Plant energy use is fundamental to plant survival and growth. However, we still lack effective means to quantify plant energy use strategies. This study introduced a concept quantifying the light level at which photochemical and non-photochemical energy use in plants are in equilibrium — the photochemical compensation point (PCCP) which can be determined with chlorophyll fluorescence measurements. We used winter wheat as a test case to explore the dynamics of PCCP and its physiological and biochemical regulations. Winter wheat PCCP decreased significantly across growth stages from jointing to grain filling. Long-term nitrogen and phosphate (NP) fertilization significantly increased PCCP, whereas potassium (K) and manure (M) fertilizer supplementation had negligible effects. PCCP exhibited significant positive correlations with leaf thickness, leaf P and sulfur (S), and stomatal conductance (gs) across all growth stages. All manure-amended treatments exhibited positive correlations of PCCP with leaf N, P, K and gs, and negative correlations with leaf calcium (Ca). Random forest analysis revealed that gs was the most significant predictor of PCCP variation, followed by leaf P, iWUE, and leaf thickness across all treatments. We suggest that plant energy use strategies are strongly coupled with plant water use strategies and nutrient availability through a complex interplay of effects on physiological and biochemical traits.

energy allocation↗

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗

Impacts of the morphology of new neighborhoods on microclimate and building energy

In anticipation of emerging global urbanization and consequent increases in energy use and carbon dioxide emissions, better understanding and quantification of climate effects on energy use in cities are needed, requiring coordinated research into large-scale, regional, and microclimate impacts to and from the city structure. The methodology described here addresses this need by (1) demonstrating a process for creating and testing example morphologies for new neighborhoods for their impact on local and regional meteorology within a two-way-coupled four-domain nested mesoscale weather model (6 km horizontal resolution outer domain, 90 m horizontal innermost domain) and (2) allocating resulting building-level meteorological profiles to each building in a neighborhood for parallel computation of building-by-building energy use. Our Chicago Loop test case shows that the morphology of even a small new added development to a neighborhood affects not only its own microclimate, but also the microclimate of the original neighborhood to which the development was added, and that these changes in microclimate affect both neighborhoods’ building energy use. Furthermore, this method represents an important step toward quantifying and analyzing the relationships among climatic conditions, urban morphology, and energy use and using these relationships to inform energy-efficient urban development and planning.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Urban building energy modeling (UBEM) tools: A state-of-the-art review of bottom-up physics-based approaches

Regulations corroborate the importance of retrofitting existing building stocks or constructing new energy-efficient districts. There is, thus, a need for modeling tools to evaluate energy scenarios to better manage and design cities, and numerous methodologies and tools have been developed. Among them, Urban Building Energy Modelling (UBEM) tools allow the energy simulation of buildings at large scales. Choosing an appropriate UBEM tool, balancing the level of complexity, accuracy, usability, and computing needs, remains a challenge for users. The review focuses on the main bottom-up physics-based UBEM tools, comparing them from a user-oriented perspective. Five categories are used: (i) the required inputs, (ii) the reported outputs, (iii) the exploited workflow, (iv) the applicability of each tool, and (v) the potential users. Moreover, a critical discussion is proposed, focusing on interests and trends in research and development. The results highlighted major differences between UBEM tools that must be considered to choose the proper one for an application. Finally, barriers of adoption of UBEM tools include the needs of a standardized ontology, a common three-dimensional city model, a standard procedure to collect data, and a standard set of test cases. This feeds into future development of UBEM tools to support cities’ sustainability goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Renewable energy analysis in indigenous communities using bottom-up demand prediction

This paper provides a methodology for the holistic analysis of hybrid renewable energy systems in rural communities. Electric demand is an important component for modeling and analysis of renewable energy systems. Typically, electric demand data is not available due to the internal privacy policies of utility providers. Therefore, this study proposes the use of bottom-up approaches for the development of the electric demand profile, considering the general homogeneity of residential and commercial buildings in rural communities. As a test case, this study develops the electric demand profile and investigates the technical and environmental feasibility of a hybrid renewable energy system for the New Town community on the Fort Berthold Indian Reservation (FBIR) in North Dakota. This study conducts the hybrid renewable energy system’s analysis by developing scripts in the LK scripting language and integrating System Advisor Model software’s open-source modules for modeling of renewable energy systems. Here, the results for the validation testbed of this study show that hybrid renewable resources have higher ratios of energy used for self-consumption to the total energy generated compared to stand-alone wind and PV farms.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Deep reinforcement learning assisted co-optimization of Volt-VAR grid service in distribution networks

With the increasing penetration of distributed energy resources in distribution networks, Volt-VAR control and optimization (VVC/VVO) have become very important to ensure an acceptable quality of service to all customers. System operators can rely on slow-responding utility devices, including capacitor banks and on-load tap changing transformers, along with fast-responding battery and photovoltaic (PV) inverters for the VVC/VVO implementation. Because of variations in response time of these two classes of devices, and different control actions (discrete versus continuous), coordinated and optimal scheduling and operation have become of utmost importance. Here, this paper develops a look-ahead deep reinforcement learning (DRL)-based multi-objective VVO technique to improve the voltage profile of active distribution networks, decrease network and inverter power loss, and save the operational cost of the grid. It proposes a deep deterministic policy gradient (DDPG)-based approach to schedule the optimal reactive and/or active power set-points of fast-responding inverters, and a deep Q-network (DQN)-based DRL agent to schedule the discrete decisions variables of slow-responding assets. The reactive power output of PV and battery smart inverters are scheduled at 30-minute intervals and the capacitors’ commitment status is scheduled with several hour intervals. The proposed framework is validated on the modified IEEE 34-bus and 123-bus test cases with embedded PV and PV-plus-storage. To validate the efficacy of the proposed VVO, it is compared with several scenarios, including the base case without VVO, localized droop control of DERs, DDPG-only, and twin delayed DDPG (TD3) agent-based DRL techniques. The results justify the superior performance of the proposed method to improve the voltage profile, reduce network power loss, and minimize the look-ahead grid operational cost while minimizing the undesirable power losses in inverters as a result of power factor adjustments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

How equation of state selection impacts accuracy near the critical point: Forced convection supercritical CO 2 flow over a cylinder

Increasing application of supercritical fluids leads to increased use of numerical simulations to guide designs. However, the impact of Equations of State (EoS) selection on computational cost and accuracy is not well understood. Laminar flow around a heated cylinder is used as a test case. A robust matrix of 99-simulations is used to analyze the flow of carbon dioxide just above the critical point at the isobaric condition of 8 MPa and temperatures between 305 and 390 K. Four EoS were studied: Ideal Gas Law, cubic EoS’s Peng-Robinson and Soave-Redlich-Kwong, and the Span-Wagner EoS. Results were identical for isothermal cases, but for heated or cooled cases, the high-fidelity Span-Wagner EoS simulations presented drastically different wake structures. Compared to the high-fidelity SW EoS, cEoS’s underestimated Nu values by up to 38% and cEoS’s underestimated CD values by 18–42%. In conclusion, the Span-Wagner EoS required only 17% more computation time then the cubic EoS simulations, making the Span-Wagner EoS recommended for engineering design.

42 ENGINEERING↗

A Bayesian framework to investigate radiation reaction in strong fields

Recent experiments aiming to measure phenomena predicted by strong-field quantum electrodynamics (SFQED) have done so by colliding relativistic electron beams and high-power lasers. In such experiments, measurements of collision parameters are not always feasible. However, precise knowledge of these parameters is required to accurately test SFQED. Here, we present a novel Bayesian inference procedure that infers collision parameters that could not be measured on-shot. This procedure is applicable to all-optical non-linear Compton scattering experiments investigating radiation reaction. The framework allows multiple diagnostics to be combined self-consistently and facilitates the inclusion of known information pertaining to the collision parameters. Using this Bayesian analysis, the relative validity of the classical, quantum-continuous and quantum-stochastic models of radiation reaction was compared for several test cases, which demonstrates the accuracy and model selection capability of the framework and highlight its robustness if the experimental values of fixed parameters differ from their values in the models.

47 OTHER INSTRUMENTATION↗

Improved performance of stellarator coil design optimization

Following up on earlier work which demonstrated an improved numerical stellarator coil design optimization performance by the use of stochastic optimization (Lobsien et al., Nucl. Fusion, vol. 58 (10), 2018, 106013), it is demonstrated here that significant further improvements can be made – lower field errors and improved robustness – for a Wendelstein 7-X test case. This is done by increasing the sample size and applying fully three-dimensional perturbations, but most importantly, by changing the design sequence in which the optimization targets are applied: optimization for field error is conducted first, with coil shape penalties only added to the objective function at a later step in the design process. A robust, feasible coil configuration with a local maximum field error of 3.66 % and an average field error of 0.95 % is achieved here, as compared to a maximum local field error of 6.08 % and average field error of 1.56 % found in our earlier work. These new results are compared to those found without stochastic optimization using the FOCUS and ONSET suites. The relationship between local minima in the optimization space and coil shape penalties is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Strong gradient neoclassical transport in the plateau regime

Strong gradient regions in tokamaks such as the pedestal or internal transport barriers are regions of reduced turbulence where neoclassical transport can play a dominant role. In pedestals, gradient lengths comparable to the ion poloidal gyroradius have been measured. Standard neoclassical theory can miss important strong gradient effects in these regions because it assumes that the gradient length scales of density, temperature and potential are larger than the ion poloidal gyroradius. We extend plateau regime neoclassical theory into regions of gradients of the order of the ion poloidal gyroradius to capture strong gradient effects on transport processes in the pedestal and internal transport barriers. The fundamental idea behind our new framework is to keep a scale separation between the orbit widths and the gradient length scales by performing a large aspect ratio expansion. In the plateau regime, strong gradients cause poloidal variation that is in–out as well as up–down asymmetric. We study two different test cases assuming either radial force balance or the absence of turbulence and show that strong gradient effects can enhance or reduce standard neoclassical theory predictions in the plateau regime in strong gradient regions.

fusion plasma↗

Community-Driven Methods for Open and Reproducible Software Tools for Analyzing Datasets from Atom Probe Microscopy

Atom probe tomography, and related methods, probe the three-dimensional architecture of a material. The software tools that microscopists use, and how these tools are connected into workflows, makes a substantial contribution to the accuracy and precision of such a material characterization experiment. Typically, we adapt methods from other communities like mathematics, data science, computational geometry, artificial intelligence, or scientific computing. We also realize that improving on research data management is a challenge when it comes to align with the FAIR data stewardship principles. Faced with this global challenge, we are convinced that collaborating is useful. Here, we report the results and challenges with an inter-laboratory call for developing test cases for several types of atom probe software tools. The results support why defining detailed recipes of software workflows and sharing these recipes is necessary and rewarding: Open source tools and (meta)data exchange can help to make our day-to-day data processing tasks become more efficient, the training of new users and knowledge transfer become easier, and assist us with automated quantification of uncertainties to gain access to substantiated results.

36 MATERIALS SCIENCE↗