SEARCH · Search NASA
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.
This report details an end of summer internship. The report lists the objective, "Add more multirate time integration methods to Amanzi" and concludes that, "Multirate methods can be used to speed up simulations and to get higher orders of accuracy", and "Coupled surface and subsurface simulations could benefit from using multirate schemes".
New developments describing the theoretical basis towards an effective virtual-pulse time integral procedure are presented for general nonlinear transient heat transfer problems. To validate the proposed methodology of computation, numerical test examples are given and comparisons are drawn with the implicit second-order accurate Crank-Nicolson technique. The proposed method shows superior or improved accuracy and stability characteristics for the models tested.
The paper describes the development and applicability of a generalized VIrtual-Pulse (VIP) time integral method of computation for thermal problems. Unlike past approaches for general heat transfer computations, and with the advent of high speed computing technology and the importance of parallel computations for efficient use of computing environments, a major motivation via the developments described in this paper is the need for developing explicit computational procedures with improved accuracy and stability characteristics. As a consequence, a new and effective VIP methodology is described which inherits these improved characteristics. Numerical illustrative examples are provided to demonstrate the developments and validate the results obtained for thermal problems.
The use of the area-time-integral (ATI) technique, based only on satellite data, to estimate convective rain volume over a moving target is examined. The technique is based on the correlation between the radar echo area coverage integrated over the lifetime of the storm and the radar estimated rain volume. The processing of the GOES and radar data collected in 1981 is described. The radar and satellite parameters for six convective clusters from storm events occurring on June 12 and July 2, 1981 are analyzed and compared in terms of time steps and cluster lifetimes. Rain volume is calculated by first using the regression analysis to generate the regression equation used to obtain the ATI; the ATI versus rain volume relation is then employed to compute rain volume. The data reveal that the ATI technique using satellite data is applicable to the calculation of rain volume.
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.
We measure the time-integrated C P asymmetry in D 0 → K S 0 K S 0 decays reconstructed in e + e − → c c ¯ events collected by the Belle and Belle II experiments. The corresponding data samples have integrated luminosities of 980 and 428 fb − 1 , respectively. The D 0 decays are required to originate from the D * + → D 0 π + decay, which determines the charm flavor at production time. A control sample of D 0 → K + K − decays is used to correct for production and detection asymmetries. The result, ( − 1.4 ± 1.3 ( stat ) ± 0.1 ( syst ) ) % , is consistent with previous determinations and with C P symmetry. Published by the American Physical Society 2025
The research described here was performed as part of the DOE SciDAC project Coupling Approaches for Next Generation Architectures (CANGA). A framework was developed for the derivation of novel algorithms for the multirate time integration of two-component systems coupled across an interface between spatial domains. The multirate aspect means that different time steps are allowed by each component integrator. The framework provides a way to construct multirate integrators with desirable properties related to stability, accuracy and preservation of system invariants. This report describes the framework and summarizes the major results, examples and research products.
The principal goal of this project is to establish relationships that would allow application of area-time integral (ATI) calculations based upon satellite data to estimate rainfall volumes. The research has been pursued using two different approaches, which for convenience can be designated as the 'fixed-threshold approach' and the 'variable-threshold approach'. In the former approach, an attempt is made to determine a single temperature threshold in the satellite infrared data that would yield ATI values for identifiable cloud clusters which are most closely related to the corresponding rainfall amounts. Results thus far have indicated that a strong correlation exists between the rain volumes and the satellite ATI values, but the optimum threshold for this relationship seems to differ from one geographic location to another. The difference is probably related to differences in the basic precipitation mechanisms that dominate in the different regions. The average rainfall rate associated with each cloudy pixel is also found to vary across the spectrum of ATI values. Work on the second, or 'variable-threshold', approach for determining the satellite ATI values was essentially suspended during this period due to exhaustion of project funds. Most of the ATI work thus far has dealt with cloud clusters from the Lagrangian or 'floating-target' point of view. For many purposes, however, the Eulerian or 'fixed-target' perspective is more appropriate. For a very large target area encompassing entire cluster life histories, the rain volume-ATI relationship would obviously be the same in either case. The important question for the Eulerian perspective is how small the fixed area can be made while maintaining consistency in that relationship.
We present what are to our knowledge first-time calculations from vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional dielectric waveguides. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and the nonlinear convolution accounts for two quantum effects, the Kerr and Raman interactions. By retaining the optical carrier, the new method solves for fundamental quantities - optical electric and magnetic fields in space and time - rather than a nonphysical envelope function. It has the potential to provide an unprecedented two- and three-dimensional modeling capability for millimeter-scale integrated-optical circuits with submicrometer engineered inhomogeneities.
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.
The principal goal of this project is to establish relationships that would allow application of area-time integral (ATI) calculations based upon satellite data to estimate rainfall volumes. The research is being carried out as a collaborative effort between the two participating organizations, with the satellite data analysis to determine values for the ATIs being done primarily by the STC-METSAT scientists and the associated radar data analysis to determine the 'ground-truth' rainfall estimates being done primarily at the South Dakota School of Mines and Technology (SDSM&T). Synthesis of the two separate kinds of data and investigation of the resulting rainfall-versus-ATI relationships is then carried out jointly. The research has been pursued using two different approaches, which for convenience can be designated as the 'fixed-threshold approach' and the 'adaptive-threshold approach'. In the former, an attempt is made to determine a single temperature threshold in the satellite infrared data that would yield ATI values for identifiable cloud clusters which are closely related to the corresponding rainfall amounts as determined by radar. Work on the second, or 'adaptive-threshold', approach for determining the satellite ATI values has explored two avenues: (1) attempt involved choosing IR thresholds to match the satellite ATI values with ones separately calculated from the radar data on a case basis; and (2) an attempt involved a striaghtforward screening analysis to determine the (fixed) offset that would lead to the strongest correlation and lowest standard error of estimate in the relationship between the satellite ATI values and the corresponding rainfall volumes.
We measure the time-integrated 𝐶𝑃 asymmetry in 𝐷 0 → 𝐾$^{0}_{S}$𝐾$^{0}_{S}$ decays reconstructed in 𝑒 + 𝑒 − → $c\bar{c}$ events collected by the Belle and Belle II experiments. The corresponding data samples have integrated luminosities of 980 and 428 fb −1 , respectively. To infer the flavor of the 𝐷 0 meson, we exploit the correlation between the flavor of the reconstructed decay and the electric charges of particles reconstructed in the rest of the 𝑒 + 𝑒 − → $c\bar{c}$ event. This results in a sample which is independent from any other previously used at Belle or Belle II. The result, 𝐴 𝐶𝑃 (𝐷 0 → 𝐾$^{0}_{S}$𝐾$^{0}_{S}$)=(1.3±2.0±0.2)%, where the first uncertainty is statistical and the second systematic, is consistent with previous determinations and with 𝐶𝑃 symmetry.
Discrete dislocation dynamics (DDD) simulations are being increasingly employed to investigate metal plasticity at the mesoscale. However, in spite of its ability to access beyond the length and time limits of atomistic methods, the DDD model is still limited by its high computational cost, with ranges of achievable strains too low and strain rates too high by several orders of magnitude compared with typical experimental conditions. By combining the efficiency of the recently developed subcycling time-integrator with the highly parallel architecture of graphics processing unit (GPU) devices, we developed a DDD model that provides significant acceleration compared to existing implementations. Our GPU-accelerated implementation enables large-scale DDD simulations that can reach relevant levels of strain using a moderate amount of computational resources.
The initial results for femtosecond pulse propagation and scattering interactions for a Lorentz medium obtained by a direct time integration of Maxwell's equations are reported. The computational approach provides reflection coefficients accurate to better than 6 parts in 10,000 over the frequency range of dc to 3 x 10 to the 16th Hz for a single 0.2-fs Gaussian pulse incident upon a Lorentz-medium half-space. New results for Sommerfeld and Brillouin precursors are shown and compared with previous analyses. The present approach is robust and permits 2D and 3D electromagnetic pulse propagation directly from the full-vector Maxwell's equations.
An implicitly Coulomb-gauged A-ϕ formulation has previously been proposed and validated for finite ele- ment simulations of low-frequency and multiscale electromag- netic problems in the frequency domain. This formulation has demonstrated numerical stability across all frequencies, with its accuracy, efficiency, and iterative convergence established in various frequency-domain scenarios. However, direct time- domain computation is often preferable for wideband electro- magnetic problems and is typically indispensable in nonlinear and multiphysics simulations. In this work, the A-ϕ formulation is extended to the time domain. By incorporating the well-known Newmark-β time integration scheme, the proposed formulation is validated through capacitive and inductive test cases. The results confirm the solution’s accuracy and demonstrate the formulation’s stability in the time domain.
Plasmas with varying collisionalities occur in many applications, such as tokamak edge regions, where the flows are characterized by significant variations in density and temperature. While a kinetic model is necessary for weakly-collisional high-temperature plasmas, high collisionality in colder regions render the equations numerically stiff due to disparate time scales. In this study, we propose an implicit–explicit algorithm for such cases, where the collisional term is integrated implicitly in time, while the advective term is integrated explicitly in time, thus allowing time step sizes that are comparable to the advective time scales. This partitioning results in a more efficient algorithm than those using explicit time integrators, where the time step sizes are constrained by the stiff collisional time scales. Finally, we implement semi-implicit additive Runge–Kutta methods in COGENT, a high-order finite-volume gyrokinetic code and test the accuracy, convergence, and computational cost of these semi-implicit methods for test cases with highly-collisional plasmas.
A stable explicit time-scale splitting algorithm for stiff chemical Langevin equations (CLEs) is developed, based on the concept of computational singular perturbation. The drift term of the CLE is projected onto basis vectors that span the fast and slow subdomains. The corresponding fast modes exhaust quickly, in the mean sense, and the system state then evolves, with a mean drift controlled by slow modes, on a random manifold. The drift-driven time evolution of the state due to fast exhausted modes is modeled algebraically as an exponential decay process, while that due to slow drift modes and diffusional processes is integrated explicitly. This enables time integration step sizes much larger than those required by typical explicit numerical methods for stiff stochastic differential equations. The algorithm is motivated and reported on, and extensive numerical experiments are conducted to illustrate its accuracy and stability with a number of model systems.
A variable order method of integrating the structural dynamics equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. When the time variation of the system can be modeled exactly by a polynomial it produces nearly exact solutions for a wide range of time step sizes. Solutions of a model nonlinear dynamic response exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with solutions obtained by established methods.