Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numberical Integration”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Generation of Near-Equipartition Magnetic Fields in Turbulent Collisionless Plasmas

The mechanisms that generate “seed” magnetic fields in our Universe and that amplify them throughout cosmic time remain poorly understood. By means of fully kinetic particle-in-cell simulations of turbulent, initially unmagnetized plasmas, here we study the genesis of magnetic fields via the Weibel instability and follow their dynamo growth up to near-equipartition levels. In the kinematic stage of the dynamo, we find that the rms magnetic field strength grows exponentially with rate γ B ≃0.4⁢u rms /L, where L/2⁢π is the driving scale and u rms is the rms turbulent velocity. In the saturated stage, the magnetic field energy reaches about half of the turbulent kinetic energy. Here, magnetic field growth is balanced by dissipation via reconnection, as revealed by the appearance of plasmoid chains. At saturation, the integral-scale wave number of the magnetic spectrum approaches k int ≃12⁢π/L. Our results show that turbulence—induced by, e.g., the gravitational buildup of galaxies and galaxy clusters—can magnetize collisionless plasmas with large-scale near-equipartition fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Deep Generative Models that Solve PDEs: Distributed Computing for Training Large Data-Free Models

Recent progress in scientific machine learning (SciML) has opened up the possibility of training novel neural network architectures that solve complex partial differential equations (PDEs). Several (nearly data free) approaches have been recently reported that successfully solve PDEs, with examples including deep feed forward networks, generative networks, and deep encoder-decoder networks. However, practical adoption of these approaches is limited by the difficulty in training these models, especially to make predictions at large output resolutions (≥1024×1024). Here we report on a software framework for data parallel distributed deep learning that resolves the twin challenges of training these large SciML models - training in reasonable time as well as distributing the storage requirements. Our framework provides several out of the box functionality including (a) loss integrity independent of number of processes, (b) synchronized batch normalization, and (c) distributed higher-order optimization methods. We show excellent scalability of this framework on both cloud as well as HPC clusters, and report on the interplay between bandwidth, network topology and bare metal vs cloud. We deploy this approach to train generative models of sizes hitherto not possible, showing that neural PDE solvers can be viably trained for practical applications. We also demonstrate that distributed higher-order optimization methods are 2-3× faster than stochastic gradient-based methods and provide minimal convergence drift with higher batch-size.

PDEs↗
DOE OSTI2026

2013 California Vehicle Survey

The 2013 California Vehicle Survey (CVS) collected data on household and commercial vehicle usage, and on future vehicle purchases. ICF International conducted the survey on behalf of the California Energy Commission. Approximately 8,000 respondents, from California households and businesses, completed the survey. The household component of the CVS included a selection of households from the 2010-2012 California Household Travel Survey (CHTS), who had stated their intention to purchase a vehicle in the near future. Both household surveys used the same survey ID numbers enabling the integration of responses. The commercial vehicle component of the CVS—a stand-alone survey of commercial fleet owners in California—asked vehicle owners questions pertaining to economic and demographic attributes, current fleets, and preferences about planned vehicle purchases.

1Hz data↗
DOE OSTI2020

Support CEED-enabled ECP applications in their preparation for Aurora/Frontier

The goal of this milestone was to help CEED-enabled ECP applications (particularly ExaSMR, MARBL, ExaAM, ExaWind and E3SM) in their preparations for the Aurora and Frontier architectures. This work included collaboration with ECP vendors and porting and optimization of CEED’s benchmarks and miniapps to early access hardware. As part of this milestone, we also made the best bake-off problems and bake-off kernel implementation from Nek, MFEM, libParanumal and the external community available in the latest libCEED release, libCEED-0.7. During the milestone period we also organized, in virtual form, the fourth CEED Annual meeting (CEED4AM) which included representatives from ECP applications, vendors and software technology projects. The specific tasks addressed in this milestone were to: (1) Work with vendors to port and run CEED benchmarks on early access systems for Aurora and Frontier; (2) Make the best BP/BK implementations from Nek, MFEM, libParanumal and external community available in libCEED; (3) Organize the next CEED Annual meeting (CEED4AM); and (4) Optimize CEED applications and miniapps for Aurora and Frontier architectures. The artifacts delivered include the next libCEED release, libCEED-0.7, and a number of developments integrated within applications to improve their GPU performance and capabilities. See the CEED website, http://ceed.exascaleproject.org and the CEED GitHub organization, http://github.com/ceed for more details.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-Density, Actively Multiplexed µECoG Array on Reinforced Silicone Substrate

Simultaneous interrogation of electrical signals from wide areas of the brain is vital for neuroscience research and can aid in understanding the mechanisms of brain function and treatments for neurological disorders. There emerges a demand for development of devices with highly conformal interfaces that can span large cortical regions, have sufficient spatial resolution, and chronic recording capability while keeping a small implantation footprint. In this work, we have designed 61 channel and 48 channel high-density, cortical, micro-electrocorticographic electrode arrays with 400 µm pitch on an ultra-soft but durable substrate. We have also developed a custom multiplexing integrated circuit (IC), methods for packaging the IC in a water-tight liquid crystal polymer casing, and a micro-bonding method for attaching the electronics package to the electrode array. With the integrated multiplexer, the number of external wire connections can be reduced to 16 wires, thereby diminishing the invasive footprint of the device. Both the electrode array and IC were tested in vivo in a rat model to demonstrate the ability to sense finely-localized electrophysiological signals.

Rachinskiy, Iakov↗
DOE OSTI2020

Provision of Grid Services by PV Plants with Integrated Battery Energy Storage System: Preprint

Battery energy storage systems (BESS), due to their tremendous range of uses and configurations, may assist PV integration in any number of ways by increasing power system flexibility. In this paper we describe results of a research project conducted by NREL and First Solar for developing controls and demonstrating many use cases for PV-BESS systems including (1) matching generation to loads through time shifting, reducing PV curtailment,; (2) promoting higher levels of PV penetrations by balancing the grid through ancillary services; (3) use of PV-BESS systems to provide wide-area stability services in the form of oscillations damping controls; and (r) use of gird forming BESS for enabling black start and islanded applications for PV-BESS systems. In the paper we will also describe the unique multi-MW test and validation platform developed at NREL for conducting research on fielded utility-scale PV-BESS systems. The platform consists of First Solar 0.43 MW PV power plant, NREL’s 1MW/1MW BESS, 2MW PV power supplies and 7MVA power electronic grid simulator interfaced with industrial plant controllers and power-hardware-in-the-loop (PHIL) capability.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2023

CYDRES: CYber Defense and REsilient System for securing grid-interactive efficient buildings

Smart buildings, especially Grid-interactive Efficient Buildings (GEBs), suffer from cyber-attacks and physical faults due to the integration of a large number of sensors and controls, connected devices, and associated communication networks. This study demonstrated a real-time advanced building resilient platform, called CYber Defense and REsilient System (CYDRES), which is deployable for existing and emerging Building Automation Systems (BASs). CYDRES aims to empower GEBs with cyber-attack-immune capabilities through multi-layer prevention and adaptation mechanisms to monitor, detect, and respond to cyber-attacks and physical operational faults. CYDRES is demonstrated through real-time experiments in a Hardware-in-the-Loop (HIL) testbed.

Building automation system, Cyber-attacks, Physica↗
DOE OSTI2025

Integrated Distribution Planning

The contemporary distribution planning landscape is comprised of an increasing number of factors that require integration into the engineering of the modern electric grid. Expectations for electric utilities to accommodate heightened awareness of stakeholders' interest in things like decarbonization, resilience and equity are growing. As these interests are formed into objectives, many jurisdictions will experience increasing levels of load modifying technologies like DER, building and industrial electrification and electric vehicles which prove not only to challenge the capabilities of the grid; but the processes by which planning for it is traditionally done. Other related factors that strain the conventional distribution planning mold are the swelling amount and sources of data associated with these technologies and the need it creates for improved capabilities in the processes and tools that manage it. As the complexity of the distribution system expands, so will the distribution system's effects on the transmission and generation systems that it is a part of. Forecasting distribution system load and DER are examples of areas where this complexity will manifest, and harmonizing distribution forecasting with transmission and generation forecasting requires higher amounts of intentionality as these typically separate processes become a solitary one. Of course, core activities do not cease as a utility begins to integrate these other factors, and in this webinar we explore specifics of how distribution planning can be expected to evolve as progress towards Integrated Distribution System Planning is made.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

North American Renewables Integration Study (NARIS): Cooperative Research and Development Final Report, CRADA Number CRD-16-00633

The goal of the study is to help inform and assist power system planning entities, electricity system operators, government energy agencies, legislators and regulators to better understand the planning, engineering, operational, and economic implications of integrating large amounts of wind, solar and hydro energy into the Canadian, US and Mexican electrical systems.

14 SOLAR ENERGY↗
DOE OSTI2020

Parallel Time Integration for Constrained Optimization

The number of transistors in an average processor continues to increase, but individual clock speeds have plateaued. Those transistors are instead going into additional cores, increasing the number of different things that a processor can do at once and placing an emphasis on parallel computation. Many problems in scientific computing follow a time-evolution model, and it can be difficult to solve such problems in parallel across the temporal domain. The Multi-Grid Reduction In Time (MGRIT) algorithm, developed at Lawrence Livermore National Laboratory (LLNL), solves differential equations with a method designed specifically to take advantage of extreme numbers of processors by parallelizing across time. The Tri-diagonal MGRIT (TriMGRIT) algorithm, also developed at LLNL, is a generalization of MGRIT which enables parallel-in-time solving of a greater number of problems. Constrained optimization problems, in particular, may be solved in parallel using TriMGRIT. These consist of choosing a control function such that an objective functional is minimized, constrained by a differential-equation. We consider two such problems: applying torque to a pendulum to bring it to a gentle stop and moving a crowd of people from one distribution into another. We also perform some miscellaneous theoretical and practical research, including investigating the use of a line-search subroutine to refine intermediate TriMGRIT results and preliminary work on strategies for choosing operators for TriMGRIT to use.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Collective coordinate fix in the path integral

Collective coordinates are frequently employed in path integrals to manage divergences caused by fluctuations around saddle points that align with classical symmetries. These coordinates parametrize a manifold of zero modes and more broadly provide judicious coordinates on the space of fields. However, changing from local coordinates around a saddle point to more global collective coordinates is remarkably subtle. The main complication is that the mapping from local coordinates to collective coordinates is generically multivalued. Consequently one is forced to either restrict the domain of path integral in a delicate way, or otherwise correct for the multivaluedness by dividing the path integral by certain intersection numbers. We provide a careful treatment of how to fix collective coordinates while accounting for these intersection numbers, and then demonstrate the importance of the fix for free theories. We also provide a detailed study of the fix for interacting theories and show that the contributions of higher intersections to the path integral can be nonperturbatively suppressed. Using a variety of examples ranging from single-particle quantum mechanics to quantum field theory, we explain and resolve various pitfalls in the implementation of collective coordinates. Published by the American Physical Society 2024

Bhattacharya, Arindam (ORCID:0000000244578926)↗
DOE OSTI2025

A double copy from twisted (co)homology at genus g

We study a family of generalized hypergeometric integrals defined on punctured Riemann surfaces of genus g. These integrals are closely related to g-loop string amplitudes in chiral splitting, where one leaves the loop-momenta, moduli and all but one puncture un-integrated. We study the twisted homology groups associated to these integrals, and determine their intersection numbers. We make use of these homology intersection numbers to write a double-copy formula for the "complex" version of these integrals -- their closed-string analogues. To verify our findings, we develop numerical tools for the evaluation of the integrals in this work. This includes the recently introduced Enriquez kernels -- integration kernels for higher-genus polylogarithms.

Pokraka, Andrejz [Brown University; Amsterdam Univ↗
DOE OSTI2021

Synthesis of LaWN 3 nitride perovskite with polar symmetry

Perovskite structured materials have a variety of uses as photovoltaics, capacitors, and micromechanical actuators, along with other applications. Oxides, halides, and chalcogenides all have large numbers of perovskite structured materials. Examples of perovskite nitrides are conspicuously absent, but Talley et al managed to synthesize one (see the Perspective by Hong). Lanthanum tungsten nitride in the perovskite structure turns out to be piezoelectric, which is ideal for a variety of applications. Perovskite structured nitrides are very attractive because they could easily integrate with the large number of nitride-based semiconducting devices already in use.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Coherent Magnon–Photon Coupling in the Magnetic Semiconductor CrSBr

Magnon-based hybrid quantum systems are promising candidates for quantum interconnects and quantum sensors, and they offer a rich platform for exploring nonlinear magnonics and cavity–photon interactions. Two-dimensional (2D) van der Waals magnets provide a compact, atomically flat geometry that can be easily integrated into existing quantum circuits, such as superconducting resonators and qubits. Among various 2D magnets, the magnetic semiconductor CrSBr is particularly unique due to its strong spin–exciton, spin–lattice, and magnon–exciton interactions. In this work, we demonstrate coherent coupling between antiferromagnetic (AFM) magnons in CrSBr and microwave photons in a niobium-(Nb)-based-on-chip resonator. We tuned the magnon–photon coupling strength by changing the number of CrSBr flakes integrated into the Nb microwave photon resonators. Furthermore, this work demonstrates the first step toward integrating layered van der Waals 2D magnets into superconducting microwave circuits, with full access for microwave and optical probing.

Electromagnetic radiation↗
DOE OSTI2022

Level density within a micro-macroscopic approach

Statistical level density $ρ$($E, A$) is derived for nucleonic system with a given energy $E$, particle number $A$ and other integrals of motion in the micro-macroscopic approximation beyond the standard saddle-point method of the Fermi gas model. This level density reaches the two limits; the well-known Fermi gas grand-canonical ensemble limit for a large entropy $S$ related to large excitation energies, and the finite micro-canonical limit for a small combinatorical entropy $S$ at low excitation energies. In conclusion, the inverse level density parameter $K$ as function of the particle number $A$ in the semiclassical periodic orbit theory, taking into account the extended Thomas-Fermi and Strutinsky shell corrections, is calculated and compared with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Heterogeneous Composition of Oxygen-Evolving Complexes in Crystal Structures of Dark-Adapted Photosystem II

Photosystem II (PSII) is a homodimeric protein complex that catalyzes water oxidation at the oxygen-evolving complex (OEC), a heterocubanoid calcium-tetra-manganese cluster. Here, we analyze the omit electron-density peaks of the OEC’s metal ions in five X-ray free-electron laser PSII structures at a resolution between 2.15 Å and 1.95 Å. The omit peaks can be described by the total number of electrons and approximated by the variance of electron density distribution when the distributions are spherically symmetric. We show that the number of electrons of metal centers are different in the two OECs of PSII dimers, implying that either the oxidation states and/or occupancies of individual metal ions are different in the two monomers. In either case, we find that the two OECs of dark-adapted PSII dimers in crystals are not fully synchronized in the same expected S 1 state. Differences in redox states of the OEC in PSII only partially account for the observation that the electron densities integrate to a smaller number of electrons than expected. Differences between the determined and expected relative electron numbers are much larger than the estimated errors, indicating heterogeneity in the OEC composition. Here, the reported findings suggest that two OECs of dark-adapted PSII crystals are typically not fully synchronized in the S 1 state and likely generate heterogeneity in the crystallographical data of PSII intermediate structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Machine Learned Empirical Numerical Integrator from Simulated Data

Recently, a number of state-of-the-art surrogate machine learning (ML) models have been designed for global weather and climate prediction, which have been trained using reanalysis data products. Reanalysis data products are constructed using numerical model simulations that combine numerical integration of partial differential equations and parameterization schemes. These products are typically only archived and made available using coarsened spatial and temporal resolutions. This study explores the impact of the numerical generation methods used to produce the training datasets and the temporal resolution of those datasets on machine learning surrogate models. Using the nonlinear vector autoregression (NVAR) machine as an explainable ML technique, simple dynamical systems are emulated with ML models trained on data produced by three classical numerical integration schemes. NVAR is validated as a skillful ML method, capable of producing accurate predictions and, more importantly, reconstructing both the underlying dynamics and the numerical integration scheme used to generate the training data. However, the machine fails to generalize predictions on unseen test data generated by different numerical integration schemes, despite the underlying dynamical system being the same. This result provides a word of caution for the growing field of machine learning emulation of weather and climate dynamics. Furthermore, we illustrate using NVAR that training on temporally coarsened data may increase the required complexity of ML models and potentially introduce new numerical challenges. Finally, we discover that empirical integration schemes with arbitrary time-stepping sizes can be constructed directly from the data, which implies a potential for the development of empirical numerical integration schemes.

54 ENVIRONMENTAL SCIENCES↗