Search NASA⌕ Search

SEARCH · Search NASA

Results for “parabolic”

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 OSTI2026

Multiscale Modeling and Experimental Insights into High-Temperature Soil Biodegradation Dynamics of Semi-Crystalline Poly(Lactic Acid) Nonwoven Fabrics

This study investigates the biodegradation of semi-crystalline poly(lactic acid) (PLA) nonwovens (NWs) in soil at 58 °C using both experimental and mathematical modeling approaches. The model utilizes a system of parabolic diffusion-reaction partial differential equations (PDEs) to elucidate chemical transformations over time and in space. It accounts for phenomena such as the diffusion of water and lactic acid monomers through the polymer matrix and into the surrounding soil, along with their microbial breakdown. It also accounts for the initial PLA crystallinity and predicts its evolution in time. The model is solved numerically for a single filament, and the results were used to shed light on PLA NW transformations observed in soil over a 180-day incubation period. Various characterization techniques, including scanning electron microscopy (SEM), differential scanning calorimetry (DSC), and Raman spectroscopy, were employed to assess morphological changes, crystallinity, and molecular changes in the NWs throughout the experiment. By comparing the experimental data with the model predictions, the hydrolysis rate coefficient was found to be 3.37 × 10 -7 s -1 , while the rate of microbial degradation of lactic acid monomers was faster, of the order of 9.63 × 10 -7 s -1 . The findings highlight the significant role of crystallinity in the biodegradation process. The PLA degradation ceases when no amorphous material remains, and the crystallinity reaches 0.8, as observed in the experiments by day 120. Furthermore, this research contributes to a deeper understanding of PLA biodegradation dynamics and offers insights for effectively managing biodegradable materials in environmental settings.

Diffusion−reaction modeling↗
DOE OSTI2024

Electromagnetic properties of indium isotopes illuminate the doubly magic character of 100 Sn

Understanding the nuclear properties in the vicinity of 100 Sn, which has been suggested to be the heaviest doubly magic nucleus with proton number Z equal to neutron number N, has been a long-standing challenge for experimental and theoretical nuclear physics. In particular, contradictory experimental evidence exists regarding the role of nuclear collectivity in this region of the nuclear chart. Here, we provide further evidence for the doubly magic character of 100 Sn by measuring the ground-state electromagnetic moments and nuclear charge radii of indium (Z = 49) isotopes as N approaches 50 from above using precision laser spectroscopy. Our results span almost the complete range between the two major closed neutron shells at N = 50 and N = 82 and reveal parabolic trends as a function of the neutron number, with a clear reduction towards these two closed neutron shells. Furthermore, a detailed comparison between our experimental results and numerical results from two complementary nuclear many-body frameworks (density functional theory and ab initio methods) exposes deficiencies in nuclear models and establishes a benchmark for future theoretical developments.

100Sn↗
DOE OSTI2024

Modification of wind turbine wakes by large-scale, convective atmospheric boundary layer structures

In this study, we consider the impact of large-scale, convective structures in an unstable atmospheric boundary layer on wind turbine wakes. Simulation data from a high-fidelity large-eddy simulation (LES) of the AWAKEN wind farm site matching unstable atmospheric conditions were analyzed, and both turbine performance and wake behavior were affected based on their location relative to the convective structures. Turbines located in updraft regions of the flow experienced lower inflow velocity and generated less power, but their wakes were observed to recover faster and saw greater turbulent kinetic energy mixing higher in the boundary layer. The opposite effect was found for turbines in the downdraft regions of the convective structures. A simplified model of this wake behavior was also developed based on a two-dimensional k–ε Reynolds-Averaged Navier–Stokes formulation. This simplified model included the effects of vertical transport, but could be efficiently solved as a parabolic system, and was found to capture similar wake modifications observed in the high-fidelity LES computations.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Laser-driven flash x-ray radiography of a shocked metallic foil

Characterizing hydrodynamic instability evolution in millimeter-scale, high-Z foils is crucial for understanding complex phenomena in high-energy-density physics. Here, we demonstrate a proof-of-concept, laser-driven flash x-ray radiography platform tailored for two-dimensional linear density mapping in shocked high-Z foils. Using chromium (Cr) foils with internal shockwaves (∼100 μm width), our platform achieves a spatial resolution of 59.8 ± 1.4 μm by employing a broadband x-ray source extending into the hundreds of keV range. The setup combines a compound parabolic concentrator cone with a tantalum wire target, a magnetic field to deflect residual transmitted electrons, and a copper casing to shield the sides and rear of the image plate pack. By varying the delay of the short-pulse beam driving the flash x-ray source, we resolve shockwave dynamics, specifically the velocity, position, width, and density profile, within the Cr foil. Reported experimental results are consistent with the corresponding hydrodynamics and radiation transport simulations, which accurately reproduce the measured electron and x-ray source terms. These developments enable the conversion of shockwave radiographs into two-dimensional density maps, enhancing interpretability for hydrodynamic instability evolution applications and validating the simulation approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Implementation of compound refractive lenses for large field-of-view x-ray phase-contrast imaging during hypervelocity impact experiments

Synchrotron x-ray phase-contrast imaging (XPCI) offers time-resolved visualization of dynamic compression phenomena, but its intrinsically small field-of-view (FOV) limits the time that key features remain in frame. A novel approach to enlarge the FOV is achieved by positioning a two-dimensional parabolic compound refractive lens (CRL) upstream of the sample to deliberately defocus the white beam. Ray-tracing simulations and XPCI measurements show that this CRL configuration can expand the beam by ∼50% vertically and ∼15% horizontally based on the full width at half-maximum of the beam. Implementing the CRL, however, attenuates the photon flux and lowers signal-to-noise ratio (SNR). Task-based analysis using a calibration grid (30 μm dots) showed that both setups fail to consistently meet the Rose criterion (SNR ≥ 5) for features of this size in single-bunch imaging. Extrapolating the measured SNR Rose values suggests that the minimum consistently detectable feature lies closer to 30–40 μm for the standard XPCI setup and above 40 μm for CRL-XPCI. Despite this limitation, the CRL configuration nearly doubles the illuminated area, enabling simultaneous tracking of front and rear observations of boron carbide targets subjected to rod and sphere impacts at 1.0–2.6 km/s. Image tracking algorithms and photonic Doppler velocimetry were used to measure penetration and rear-surface velocity histories. Together, these measurements capture crack fronts, penetration, and material breakout, offering new benchmark data for validating high-strain-rate constitutive models of ceramic materials.

Ceramic materials↗
DOE OSTI2025

Development and scaling of MeV x-ray radiography at NIF-ARC

We investigate the generation of MeV x-rays using the Advanced Radiography Capability laser system at the National Ignition Facility using 1, 10, and 38 ps pulse durations, with laser energies reaching up to 4 kJ and using compound parabolic concentrators. Hydrodynamic simulations using up-to-date measurements of the contrast of the ARC laser are conducted, allowing us to employ an electron scaling model that incorporates scale length and pulse duration, which aligns closely with the observed temperature distributions. Comparable x-ray sources, in terms of dose, are generated when using a 10 ps pulse duration with 2.4 kJ at ∼1.9 ± 0.4 × 10 18 W/cm 2 and when using a 38 ps pulse duration with 4 kJ at 9.9 ± 0.4 × 10 17 W/cm 2 , both achieving ∼16 rad in air at 1 m for x-ray energies > 0.5 MeV. Radiographs performed on the laser “line-of-sight” show significant improvements in image quality than radiographs performed at 65 degrees to the laser axis. We verify the radiography performance using Monte Carlo Simulations.

Physics - Plasma physics↗
DOE OSTI2026

Boundary-layer receptivity to oblique freestream vorticity waves for a high-enthalpy hypersonic flow

The receptivity of a Mach 15 straight-cone boundary layer to oblique freestream vorticity waves is investigated using direct numerical simulation (DNS) alongside linear stability theory and the linear parabolized stability equations. A thermochemical nonequilibrium gas model is used. Oblique freestream vorticity waves at frequencies of 400, 800, and 1200 kHz are considered, with incident angles ranging from 0° to 29.4° at 400 kHz, 0° to 15.7° at 800 kHz, and 0° to 20.6° at 1200 kHz. The 400 kHz case is of primary interest due to the strong second-mode amplification at this frequency. Although the underlying base flow is axisymmetric, the oblique vorticity waves lead to a fully three-dimensional boundary-layer disturbance whose characteristics vary depending on the azimuthal ray relative to the freestream wave. Moving downstream within the second-mode instability region, some clear trends emerge in terms of the boundary-layer disturbance amplitudes; that is, disturbance amplitudes are highest at the leeward ray (relative to the freestream wave), but weakest about halfway between the windward and leeward rays. Increasing the incident angle causes the amplitudes to increase on the leeward ray and decrease on the windward ray. Moreover, the boundary-layer disturbance throughout contains a wide spectrum of azimuthal wavenumbers in which the disturbance energy falls off at higher wavenumbers. Increasing the incident angle causes the azimuthal spectrum of the boundary-layer disturbance to broaden overall. Qualitatively similar results are found for the two higher frequencies leading up to the peak-amplitude locations corresponding to the second-mode instability.

42 ENGINEERING↗
DOE OSTI2026

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High-significance detection of correlation between the unresolved gamma-ray background and the large-scale cosmic structure

Our understanding of the γ-ray sky has improved dramatically in the past decade, however, the unresolved γ-ray background (UGRB) still has a potential wealth of information about the faintest γ-ray sources pervading the Universe. Statistical cross-correlations with tracers of cosmic structure can indirectly identify the populations that most characterize the γ-ray background. In this study, we analyze the angular correlation between the γ-ray background and the matter distribution in the Universe as traced by gravitational lensing, leveraging more than a decade of observations from the Fermi-Large Area Telescope (LAT) and 3 years of data from the Dark Energy Survey (DES). We detect a correlation at signal-to-noise ratio of 8.9. Most of the statistical significance comes from large scales, demonstrating, for the first time, that a substantial portion of the UGRB aligns with the mass clustering of the Universe as traced by weak lensing. Blazars provide a plausible explanation for this signal, especially if those contributing to the correlation reside in halos of large mass (∼ 10 14 M ⊙ ) and account for approximately 30–40% of the UGRB above 10 GeV. Additionally, we observe a preference for a curved γ-ray energy spectrum, with a log-parabolic shape being favored over a power-law. We also discuss the possibility of modifications to the blazar model and the inclusion of additional γ-ray sources, such as star-forming galaxies, misalinged active galactic nuclei, or particle dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Cosmological neutrino mass: a frequentist overview in light of DESI

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-α forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power σ of the different data combinations. We find that all profiles in the ΛCDM model are cut off by the ∑m ν ≥ 0 bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. The corresponding constraining power σ is 43 meV, which highlights the importance of the cut-off by negative values in the determination of the upper limit. Extending ΛCDM to non-zero curvature and w 0 w a CDM relaxes the constraints past 59 meV again, but only w 0 w a CDM exhibits profiles with a minimum at a positive value. Additionally, we extend the formalism to constrain the lightest neutrino mass. For DESI DR2 BAO, Planck PR4 and CMB lensing, we find confidence limits at 20 and 19 meV for normal and inverted ordering, respectively. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-α P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields m l ≤ 97 and 98 meV in the normal and inverted orderings, and total neutrino mass ∑m ν ≤ 285 meV (95% C.L.). The addition of DESI full-shape or Lyman-α P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-α free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

Frequentist statistics↗
DOE OSTI2025

Refining the two-band model for highly compensated semimetals using thermoelectric coefficients

In studying compensated semimetals, the two-band model has proven extremely useful in capturing electrical conductivity under magnetic field, as a function of density and mobility of electronlike and holelike carriers. However, it rarely offers practical insight into magnetothermoelectric properties. Here, in this work, we report the field dependence of thermoelectric (TE) coefficients in a highly compensated semimetal NbSb 2 , where we find the Seebeck (𝑆 𝑥⁢𝑥 ) and Nernst (𝑆 𝑥⁢𝑦 ) coefficients increase quadratically and linearly with applied magnetic field, respectively. Such field dependence was predicted in previous work that studied a system of two parabolic bands, within semiclassical Boltzmann transport theory when the following two conditions are simultaneously met: 𝜔 𝑐 ⁢𝜏 ≫ 1 and tan ⁡𝜃 𝐻 ≪ 1, where 𝜔 𝑐 and 𝜏 refer to the cyclotron frequency and relaxation time, respectively, and 𝜃 𝐻 is the Hall angle. Under these conditions, we find the field dependence of the TE coefficients directly provides a relation between the electronlike (𝑛 𝑒 ) and holelike (𝑛 ℎ ) carrier densities, which in turn can be used to refine two-band model fitting. With this, we find the compensation factor ($\frac{|Δ⁢𝑛|}{𝑛_{𝑒}}$, where Δ⁢𝑛 = 𝑛 𝑒 −𝑛 ℎ ) of NbSb 2 is two orders of magnitude smaller than what was found in unrestricted fitting, resulting in a larger saturation-field scale for magnetoresistance. Within the same framework of the semiclassical theory, we also deduce that the thermoelectric Hall angle tan⁡ 𝜃 𝛾 = $\frac{𝑆_{𝑥⁢𝑦}}{𝑆_{𝑥⁢𝑥}}$ can be expressed as ($\frac{Δ⁢𝑛}{𝑛_{𝑒}}$ ×𝜔 𝑐 ⁢𝜏) −1 , which serves as a parameter to predict the degree of compensation. Our findings offer crucial insights into identifying empirical conditions for field-induced enhancement of TE performance and into engineering-efficient thermoelectric devices based on semimetallic materials.

electrical conductivity↗
DOE OSTI2025

Franck-Condon electron emission from polar semiconductor photocathodes

An analytical formulation of (optical-)phonon-mediated and momentum-resonant Franck-Condon emission of photoexcited electrons from polar semiconductors is shown to be very consistent with (i) the observed emission properties of a cesiated Ga⁢As⁢(001) photocathode at 808 nm [J. Phys. D: Appl. Phys. 54, 205301 (2021)] and (ii) the measured spectral emission properties of a Ga⁢N(0001) photocathode from just below its bandgap energy to 5 eV. The theoretical analysis in the parabolic band approximation predicts the form of both the quantum efficiency and mean transverse energy of photoemission as a function of the photocathode’s electron affinity and the electron temperature in the vicinity of its emission face. The good agreement between theory and experimental data also suggests that sub-10-nm rms surface roughness effects are not significant for polar semiconductor photocathodes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Nonresonant central exclusive production of charged-hadron pairs in proton-proton collisions at s = 13 TeV

The central exclusive production of charged-hadron pairs in p p collisions at a center-of-mass energy of 13 TeV is examined, based on data collected in a special high- β * run of the LHC. The nonresonant continuum processes are studied with the invariant mass of the centrally produced two-pion system in the resonance-free region, m π + π − < 0.7 or m π + π − > 1.8 GeV . Differential cross sections as functions of the azimuthal angle between the surviving protons, squared exchanged four-momenta, and m π + π − are measured in a wide region of scattered proton transverse momenta, between 0.2 and 0.8 GeV, and for pion rapidities | y | < 2 . A rich structure of interactions related to double-pomeron exchange is observed. A parabolic minimum in the distribution of the two-proton azimuthal angle is observed for the first time. It can be interpreted as an effect of additional pomeron exchanges between the protons from the interference between the bare and the rescattered amplitudes. After model tuning, various physical quantities are determined that are related to the pomeron cross section, proton-pomeron and meson-pomeron form factors, pomeron trajectory and intercept, and coefficients of diffractive eigenstates of the proton. © 2024 CERN, for the CMS and TOTEMs Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Unconventional quantum oscillations and evidence of nonparabolic electronic states in quasi-two-dimensional electron system at complex oxide interfaces

The simultaneous occurrence of electric-field controlled superconductivity and spin-orbit interaction makes two-dimensional electron systems (2DES) constructed from perovskite transition metal oxides promising candidates for the next generation of spintronics and quantum computing. It is, however, essential to understand the electronic bands thoroughly and verify the predicted electronic states experimentally in these 2DES to advance technological applications. Here, we present insights into the electronic states of the 2DES at oxide interfaces through comprehensive investigations of Shubnikov–de Haas oscillations in three different systems: EuO/KTaO 3 , LaAlO 3 /SrTiO 3 , and amorphous-LaAlO 3 /KTaO 3 . To accurately resolve these oscillations, we conducted transport measurements in high magnetic fields up to 60 T and low temperatures down to 100 mK. For 2D confined electrons at these interfaces, we observed a progressive increase of oscillations frequency and cyclotron mass with the magnetic field. We interpret these universal and intriguing findings by considering the existence of nontrivial electronic bands, for which the 𝐸−𝑘 dispersion incorporates both linear and parabolic relations. In addition to providing experimental evidence for nonparabolic electronic states in KTaO 3 and SrTiO 3 2DES, the unconventional oscillations presented in this study establish a paradigm for quantum oscillations in 2DES based on perovskite transition metal oxides, where the oscillation frequencies in 1/𝐵 exhibit quadratic dependence on the magnetic field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Exotic edge states of C 3 high-fold fermions in honeycomb lattices

A generalization of the graphene honeycomb model to the case where each site in the honeycomb lattice contains a n -fold degenerate set of eigenstates of the C 3 symmetry has been recently proposed to describe several systems, including triangulene crystals and photonic lattices. These generalized honeycomb models are defined by ( n a , n b ) , the number of C 3 eigenstates in the a and b sites of the unit cell, resulting in n a + n b bands. Thus, the (1,1) case gives the coventional honeycomb model that describes the two low-energy bands in graphene. Generalizations, such as (2,1), (2,2), and (3,3) display several nontrivial features, such as coexisting graphenelike Dirac cones with flat bands, both at zero and finite energy, as well as robust degeneracy points where a flat band and a parabolic band meet at the Γ point. Here we explore the edge states of this class of crystals, using as reference triangulene crystals, and we find several types of edge states absent in the conventional (1,1) honeycomb case, associated to the nontrivial features of the two-dimensional bands of the high-fold case. First, we find dispersive edge states associated to the finite-energy flat bands, that occur both at the armchair and zigzag termination. Second, in the case of noncentrosymmetric triangulene crystals that lead to a S = 1 Dirac band, we have a bonding-antibonding pair of dispersive edge states, localized in the same edge so that their energy splitting is reduced as their localization increases, opposite to the conventional behavior of pairs of states localized in opposite edges. Third, for the (3,3) case, that hosts a gap separating a pair of flat conduction and valence bands, we find nondispersive edge states with E = 0 in all edge terminations. Published by the American Physical Society 2024

Madail, L. (ORCID:0000000208457748)↗
DOE OSTI2024

Mirror-centered representation of a focusing hyperbolic mirror for X-ray beamlines

Conic sections are commonly used in reflective X-ray optics. Hyperbolic mirrors can focus a converging light source and are frequently paired with elliptical or parabolic mirrors in Wolter type configurations. This paper derives the closed-form expression for a mirror-centered hyperbolic shape, with zero-slope at the origin. Combined with the slope and curvature, such an expression facilitates metrology, manufacturing and mirror-bending calculations. Previous works consider ellipses, parabolas, magnifying hyperbolas or employ lengthy approximations. Here, the exact shape function is given in terms of the mirror incidence angle and the source and image distances.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

3D Deep Learning Joint Inversion of Active Seismic Full Waveform and Passive Seismic Traveltime Data for Reservoir Imaging and Uncertainty Quantification

Here, we present deep learning (DL) networks for three-dimensional (3D) joint inversion of active seismic full waveform and passive seismic traveltime data to image reservoirs and their properties and quantify imaging uncertainties. Active seismic full-waveform data can provide high-resolution monitoring images but are collected only intermittently because of their high acquisition cost. In contrast, passive seismic data can be gathered at relatively low cost between regular active surveys, although their imaging quality can be compromised by factors such as low signal-to-noise ratios and limited ray coverage of the target. Although these datasets are routinely acquired together at CO 2 storage sites, their combined inversion within a 3D DL framework has not been previously demonstrated. To our knowledge, this is the first study to address this gap, combining the strength of both data types. For efficient data storage and DL training with large 3D seismic datasets, we use a 3D data matrix in which a random number of passive seismic traveltime data are stored as parabolic envelopes using one-hot encoding and a 3D full-waveform data matrix in which multiple shot gathers are summed. Two network architectures are evaluated: a single-encoder U-Net for single-data type inversion and a dual-encoder U-Net for joint inversion of active and passive seismic data. We also evaluate the single-encoder U-Net for joint inversion by concatenating full-waveform data and traveltime data. We propose a systematic approach for selecting an optimal dropout rate that balances regularization during training and Monte Carlo dropout-based uncertainty quantification during prediction by examining the correlation coefficient between standard deviation and prediction error, along with the training misfit, across a range of dropout rates. 3D DL inversion experiments include five different network configurations, with evaluations under ideal, noisy and dropout-enabled conditions. Both model and data uncertainties are assessed, as well as their combined effects. Across all conditions, the networks consistently predict accurate CO 2 saturation models with low prediction errors, such as a structural similarity index of 0.993 and CO 2 difference of 1.1%. Uncertainty estimates show strong spatial correlation with prediction errors, confirming the effectiveness of the proposed dropout selection approach. The results demonstrate that our DL approach, utilizing compact data representations and appropriate uncertainty quantification, yields accurate subsurface images under various inversion conditions and provides valuable insights into the reliability of predictions.

Um, Evan Schankee [Lawrence Berkeley National Labo↗
DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗