Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sundials”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

A Sundial-Atlas Precursor to the TIMED Mission: A Quick-Response Global Investigation into Coupled Lower Thermospheric, Ionospheric, and Mesospheric Physics

The SUNDIAL-ATLAS effort was a global-scale investigation which responded to the science priorities of the ITM Panel, the Integrated SPD Strategy Implementation Plan as a whole, and the need for potential cost-saving design criteria for the TIMED mission. The investigation focused on coupling processes in the ionospheric-thermospheric system, taking advantage of the timelines of the ATLAS-1 mission (March 1992), and the global-scale ground-based measurement and modeling activities of the SUNDIAL program. The collaborative SUNDIAL-ATLAS activity was the first opportunity for global measurements of the chemistry, kinetics, and electrodynamics which couple the E-, Fl-, and F2-regions into a single interactive system. As such, the program represented an important first step in studying global issues; and accordingly, was an important proof of concept experiment relevant to the strategic mission plans for the ITM community and the upcoming intermediate class satellite program called TIMED. To meet its projected goals, TIMED must perform a number of critical measurements and execute a number of correlations that were to be tried and tested for the first time in the SUNDIAL-ATLAS investigation. This was designed to include global correlations of thermospheric and ionospheric composition during quiet and disturbed conditions and the co-registration of global-scale ground-based measurements with along-track satellite diagnostics. The SUNDIAL component of the current investigation addressed this need by acquiring, reducing, and analyzing a multi-sensor database that complemented and extended that which was generated in the ATLAS mission (Atmospheric Laboratory for Applications and Science). The SUNDIAL data defined the state and condition of the global-scale ionosphere in the altitude range from 100 km to the F2-peak. These data specified the peak heights and densities of the E-, Fl-, and F2-regions, along with the global distributions of intermediate, descending, and sequential layers which play a critical role in the dynamo region of the lower ionospheric-thermospheric domain. The data were collected by the SUNDIAL network of more than 50 ground-based stations utilizing ionosondes, radars, photometers, Fabry-Perot interferometers, and total electron content measurements. The data were acquired during a three-week period centered on the eight-day ATLAS-1 mission, which provided image and photometric sensing of the altitude distributions of the major and minor ions and neutrals in the ITM system. This report focuses on the scientific contributions of the SUNDIAL component of the overall investigation. Specific findings are described in seven papers (attached) published in the Journal of Geophysical Research.

Szuszczewicz, E. P.↗
DOE OSTI2022

Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers

In recent years, the SUite of Nonlinear and DIfferential/ALgebraic equation Solvers (SUNDIALS) has been redesigned to better enable the use of application-specific and third-party algebraic solvers and data structures. Throughout this work, we have adhered to specific guiding principles that minimized the impact to current users while providing maximum flexibility for later evolution of solvers and data structures. The redesign was done through the addition of new linear and nonlinear solvers classes, enhancements to the vector class, and the creation of modern Fortran interfaces. The vast majority of this work has been performed “behind-the-scenes,” with minimal changes to the user interface and no reduction in solver capabilities or performance. These changes allow SUNDIALS users to more easily utilize external solver libraries and create highly customized solvers, enabling greater flexibility on extreme-scale, heterogeneous computational architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations

Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world’s largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach is used ubiquitously across scientific domains, and in many cases leads to a set of ordinary differential equations (ODEs) that need to be solved as part of a larger “outer-loop” time-stepping approach. The SUNDIALS library provides a plethora of robust time integration algorithms for solving ODEs, and the U.S. Department of Energy Exascale Computing Project (ECP) has supported its extension to applications on exascale-capable computing hardware. In this paper, we highlight some SUNDIALS capabilities and its deployment in combustion and cosmology application codes (Pele and Nyx, respectively) where operator splitting gives rise to numerous, small ODE systems that must be solved concurrently.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

SunDial – An Integrated SHINES System to Enable High-penetration Feeder-level PV

The Project Team of Fraunhofer USA, National Grid, and IPKeys developed and conducted a pilot deployment of the SunDial system, a virtual power plant platform that enables high-penetrations of solar PV to be integrated into the distribution grid. The pilot was conducted over a 15-month period from August 2018 through October 2019 on a National Grid distribution feeder in Shirley, MA. A vendor-agnostic control platform (the “Global Scheduler”) optimally shaped the net load for a virtual portfolio of non-co-located DERs based on user-defined policy objectives. The goal of the SunDial project was to simplify and reduce the risk associated with the deployment of solar in high-penetration environments by: (1) Developing an open-source, vendor-agnostic dispatch platform that can be readily adapted to optimize control of DERs over a variety of use cases; (2) Developing auto-calibrating load and solar prediction methodologies that can be readily implemented and scaled to new deployments; (3) Developing a methodology to use demand-side management with traditional electrochemical energy storage to provide “load shaping” services in high solar penetration environments; (4) Using grid-scale storage to minimize short-term intermittency association with PV production; and (5) Deploying on the National Grid distribution system to gain experience on the potential for (and limits of) integrated storage with demand-side management.

14 SOLAR ENERGY↗
DOE OSTI2026

New Time Integrators and Capabilities in SUNDIALS Versions 6.2.0-7.4.0

SUNDIALS is a well-established numerical library that provides robust and efficient time integrators and nonlinear solvers. This article overviews several significant improvements and new features added over the last 3 years to support scientific simulations run on high-performance computing systems. Notably, three new classes of one-step methods have been implemented: low storage Runge–Kutta, symplectic partitioned Runge–Kutta, and operator splitting. In addition, we describe new timestep adaptivity support for multirate methods, adjoint sensitivity analysis capabilities for explicit Runge–Kutta methods, additional options for Anderson acceleration in nonlinear solvers, and improved error handling and logging.

Computer science↗
NASA NTRS1990

Neutral winds derived from IRI parameters and from the HWM87 wind model for the SUNDIAL campaign of September, 1986

Meridional neutral winds derived from the height of the maximum ionization of the F2-layer are compared with values from results of the HWM87 empirical neutral wind model. The time period considered in this study is the SUNDIAL-2 campaign, September 21 through October 5, 1986. Winds have been derived from measurements by a global network of ionosondes, as well as from similar quantities generated by the International Reference Ionosphere. Global wind patterns from the three sources are similar. Differences tend to be the result of local or transient phenomena that are either too rapid to be described by the order of harmonics of the empirical models, or are the result of temporal changes not reproduced by models based on average conditions.

Miller, K. L.↗
NASA NTRS1990

Electric potential patterns deduced for the SUNDIAL period of September 23-26, 1986

High-latitude electric potential and ionospheric conductance patterns are presented and discussed for the Northern Hemisphere during the SUNDIAL period of September 23-26, 1986 using the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) technique of Richmond and Kamide (1988). Data sources used for the model are satellite and ground magnetometers, electron precipitation instruments, incoherent scatter radars, and ionospheric coherent radars. The period was characterized by many substorms, and a wide variety of instantaneous patterns of electrodynamic parameters on a hemispheric scale are derived throughout this period, of which only a few examples are displayed. The entire set of electric potentials and conductances are being made available through the NCAR CEDAR Data Base for further analysis and utilization in simulation models.

Emery, B. A.↗
NASA NTRS1990

Latitudinal variation of perturbation electric fields during magnetically disturbed periods - 1986 Sundial observations and model results

F-region incoherent scatter radar drift observations from Millstone Hill and Jicamarca, h-prime F observations from Huancayo, and high latitude ground-magnetometer measurements taken during the Sundial 1986 campaign are used to study the relationship between plasmaspheric electric field perturbations and high latitude currents during disturbed periods. The observations are in good agreement with numerical results from a Rice Covection Model run that involved a sharp increase in the polar cap potential drop followed by a subsequent decrease. The zonal disturbance electric field pattern is latitude independent, and the corresponding amplitudes change approximately as L exp n (where n is about 1.5). The meridional electric field patterns and amplitudes have larger latitudinal variations. The mid-, low, and equatorial electric fields from the Rice Convection Model are in good agreement with previous results from the semianalytic, Senior-Blanc (1987) model. Also discussed are three physical mechanisms (over-shielding, fossil winds, and magnetic reconfiguration) that contribute to the long lasting (1-2 h) equatorial zonal electric field perturbations associated with a sudden northward turning of the IMF. It is predicted that the penetration of high latitude electric fields to low latitudes should, in general, be closely related to the rate of motion of the shielding layer and the equatorward edge of the diffuse aurora.

Fejer, B. G.↗
NASA NTRS1990

Neutral winds derived from IRI parameters and from the HWM87 wind model for the sundial campaign of September, 1986

Meridional neutral winds derived from the height of the maximum ionization of the F2 layer are compared with values from results of the HWM87 empirical neutral wind model. The time period considered is the SUNDIAL-2 campaign, 21 Sept. through 5 Oct. 1986. Winds were derived from measurements by a global network of ionosondes, as well as from similar quantities generated by the International Reference Ionosphere. Global wind patterns from the three sources are similar. Differences tend to be the result of local or transient phenomena that are either too rapid to be described by the order of harmonics of the empirical models, or are the result of temporal changes not reproduced by models based on average conditions.

Miller, K. L.↗
NASA NTRS1988

Penetrating of high-latitude-electric-field effects to low latitudes during SUNDIAL 1984

Electric-field-penetration events have been identified using F-region vertical-drift measurements obtained in the October 6-13, 1984 period by the Jicamarcan incoherent-backscatter radar and corresponding h-prime F measurements from ionosondes at Fortaleza, Cachoeira Paulista, and Dakar. Predictions made using the Rice Convection Model for the pattern, strength, and duration of the low-latitude electric field occurring in response to an increasing high-latitude convection agree with observations. The observed 1-2 h duration of the low-latitude response to decreased convection can be explained by the fossil-wind theory of Richmond (1983).

Spiro, R. W.↗
DOE OSTI2025

Electromagnetic Transient Simulation of Photovoltaic Inverter Using Implicit-Explicit Solver

This paper introduces the implementation of electromagnetic transient (EMT) simulations of a photovoltaic (PV) inverter module using the Implicit-Explicit (ImEx) solver in the Suite of Nonlinear and Differential/Algebraic Equation Solvers (SUNDIALS). This study demonstrates the effectiveness of the ImEx solver in overcoming the challenges inherent in simulating the complex dynamics of PV inverter modules. Furthermore, using SUNDIALS’ ImEx solver module ARKODE for EMT simulation automates key aspects of the process, such as numerical integration, providing substantial benefits including enhanced consistency, faster implementation, reduced human error, and the capability to handle the complexities of advanced numerical integration. By conducting comparative simulations with an implicit method used in commercial software, the research showcases the ImEx solver’s capability in achieving high accuracy and reliability. Results indicate that leveraging the ImEx approach significantly enhances modeling fidelity and reduces simulation setup times, offering a promising tool for the EMT analysis of PV inverter systems in power electronics-dominated power grids.

Choi, Jongchan [ORNL] (ORCID:000000025952455X)↗
DOE OSTI2024

The Pele Simulation Suite for Reacting Flows at Exascale

In this work, we present the Pele suite of software tools for compressible and incompressible reacting flows. The Pele suite leverages several different libraries, notably AMReX and SUNDIALS, to achieve performance portability on heterogeneous computing architectures across the supercomputing landscape. The Pele suite is comprised of PeleC, a compressible reacting flow block-structured adaptive mesh refinement solver, PeleLMeX, a low-Mach number reacting flow block-structured adaptive mesh refinement solver, Pele-Physics, a library for transport, thermodynamics, finite rate chemistry, soot, spray and radiation physics. The objective of this paper is (i) to present the code development efforts necessary to achieve highly effective and scalable applications for exascale machines and (ii) to detail the performance results of the Combustion-Pele project applications on Oak Ridge National Laboratory's Frontier. We show good weak and strong scaling results for both PeleC and PeleLMeX up to more than 50 billion cells on more than 4096 Frontier graphics processing unit nodes. We also present a capability demonstration simulation of a dual-fuel pulse compression ignition engine (six adaptive mesh refinement levels, and 60 billion cells or 2.1 trillion degrees of freedom) on Frontier, to date one of the largest simulations performed on the first exascale-class supercomputer.

adaptive mesh refinement↗
DOE OSTI2023

ARKODE: A Flexible IVP Solver Infrastructure for One-step Methods

We describe the ARKODE library of one-step time integration methods for ordinary differential equation (ODE) initial-value problems (IVPs). In addition to providing standard explicit and diagonally implicit Runge–Kutta methods, ARKODE supports one-step methods designed to treat additive splittings of the IVP, including implicit-explicit (ImEx) additive Runge–Kutta methods and multirate infinitesimal (MRI) methods. We present the role of ARKODE within the SUNDIALS suite of time integration and nonlinear solver libraries, the core ARKODE infrastructure for utilities common to large classes of one-step methods, as well as its use of “time stepper” modules enabling easy incorporation of novel algorithms into the library. Numerical results show example problems of increasing complexity, highlighting the algorithmic flexibility afforded through this infrastructure, and include a larger multiphysics application leveraging multiple algorithmic features from ARKODE and SUNDIALS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Multiphysics Time-Integration for Turbulent Combustion at the Exascale

Turbulent reacting flow systems are often modeled with coupled time-dependent partial differential equations (PDEs). Solving such equations can easily tax the world's largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This generic operator-splitting strategy leads to a set of ordinary differential equations (ODEs) that need to be solved as part of an "outer-loop" time-stepping approach. In many combustion applications, the ODEs to be solved can be very stiff, exhibiting timescales that span many orders of magnitude. The SUNDIALS library provides a plethora of robust time integration algorithms for solving these ODEs on exascale-capable computing hardware, yet for many complex applications (such multicomponent fuels or emissions predictions), the chemical models remain too complex to solve using reasonable resources. The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of the simulations. In this talk, I will discuss the use of the SUDIALS library of ODE solvers together with automatic code generation tools to solve complex turbulent reacting flow problems using QSSA models.

chemistry↗
DOE OSTI2024

MCTrans++: a 0-D model for centrifugal mirrors

The centrifugal mirror confinement scheme incorporates supersonic rotation of a plasma into a magnetic mirror device. This concept has been shown experimentally to drastically decrease parallel losses and increase plasma stability as compared with prior axisymmetric mirrors. MCTrans++ is a dimensionless (0-D) scoping tool which rapidly models experimental operating points in the Centrifugal Mirror Fusion Experiment (CMFX) at the University of Maryland. In the low-collisionality regime, parallel losses can be modelled analytically. A confining potential is set up that is partially ambipolar and partially centrifugal. Due to the stabilizing effects of flow shear, the perpendicular losses can be modelled as classical. Radiation losses such as bremsstrahlung and cyclotron emission are taken into account. A neutrals model is included, and, in some circumstances, charge-exchange losses are found to exceed all other loss mechanisms. We use the SUNDIALS ARKODE library to solve the underlying equations of this model; the resulting software is suitable for scanning large parameter spaces, and can also be used to model time-dependent phenomena such as a capacitive discharge. MCTrans++ has been used to verify results from prior centrifugal mirrors, create an experimental plan for CMFX and find configurations for future reactor-scale fusion devices.

Schwartz, Nick R. (ORCID:0000000296846435)↗
DOE OSTI2024

scikit-SUNDAE ((SUN)DIALS Differential Algebraic Equations) [SWR-24-137]

Scikit-SUNDAE provides Python bindings to SUNDIALS integrators. The implicit differential algebraic (IDA) solver and C-based variable-coefficient ordinary differential equations (CVODE) solver are both included. The name SUNDAE combines (SUN)DIALS and DAE, which stands for differential algebraic equations. Solvers specific to DAE problems are not frequently available in Python. An ordinary differential equation (ODE) solver is also included for completeness. ODEs can be categorized as a subset of DAEs (i.e., DAEs with no algebraic constraints). https://pypi.org/project/scikit-sundae

Randall, Corey↗
DOE OSTI2025

ORNL/PFiSM

Simple phase-field model (PFM) based on SAMRAI, Sundials and hypre

Fattebert, Jean-Luc [Oak Ridge National Laboratory↗