Search NASA⌕ Search

SEARCH · Search NASA

Results for “self similarity”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Total pressure loss in vortical solutions of the conical Euler equations

A technique for the solution of the conically self-similar form of the Euler equations is described. Solutions for the flow past a flat-plate delta wing at angle of attack are presented. These solutions show strong leading edge vortices with large total pressure losses in the cores. A study of the effects of various computational parameters on the total pressure loss is made. An explanation for the cause of the total pressure loss is presented. It is shown to be consistent with the results for both a quasi-one-dimensional model problem and the conically self-similar flow past the flat-plate delta wing.

Powell, K. G.↗
DOE OSTI2020

Nonlinear bubble competition of the multimode ablative Rayleigh–Taylor instability and applications to inertial confinement fusion

The self-similar nonlinear evolution of the multimode ablative Rayleigh–Taylor instability (RTI) and the ablation-generated vorticity effect are studied for a range of initial conditions. We show that, unlike classical RTI, the nonlinear multimode bubble-front evolution remains in the bubble competition regime due to ablation-generated vorticity, which accelerates the bubbles, thereby preventing a transition into the bubble-merger regime. We develop an analytical bubble competition model to describe the linear and nonlinear stages of ablative RTI. Here, we show that vorticity inside the multimode bubbles is most significant at small scales with large initial perturbation. Since these small scales persist in the bubble competition regime, the self-similar growth coefficient ab can be enhanced by up to 30% relative to ablative bubble competition without vorticity effects. We use the ablative bubble competition model to explain the hydrodynamic stability boundary observed in OMEGA low-adiabat implosion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Approximating Density Probability Distribution Functions Across Cosmologies

Using a suite of self-similar cosmological simulations, we measure the probability distribution functions (PDFs) of real-space density, redshift-space density, and their geometric mean. We find that the real-space density PDF is well-described by a function of two parameters: ns, the spectral slope, and σL, the linear rms density fluctuation. For redshift-space density and the geometric mean of real- and redshift-space densities, we introduce a third parameter, ${s}_{L}=\sqrt{\langle {\left({{dv}}_{\mathrm{pec}}^{L}/{dr}\right)}^{2}\rangle }/H$. We find that density PDFs for the LCDM cosmology is also well-parameterized by these three parameters. As a result, we are able to use a suite of self-similar cosmological simulations to approximate density PDFs for a range of cosmologies. We make the density PDFs publicly available and provide an analytical fitting formula for them.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Stability of a premixed flame in stagnation-point flow against general disturbances

Previously, the stability of a premixed flame in a stagnation flow was discussed for a restricted class of disturbances that are self-similar to the basic undisturbed flow; thus, flame fronts with corrugations only in the cross stream direction were considered. Here, we consider a more general class of three-dimensional flame front perturbations which also permits corrugations in the streamwise direction. It is shown that, because of the stretch experienced by the flame, the hydrodynamic instability is limited only to disturbances of short wavelength. If in addition diffusion effects have a stabilizing influence, as would be the case of mixtures with Lewis number greater than one, a stretched flame could be absolutely stable. Instabilities occur when the Lewis number is below some critical value less than one. Neutral stability boundaries are presented in terms of the Lewis number, the strain rate, and the appropriate wavenumbers. Beyond the stability threshold, the two-dimensional self-similar modes always grow first. However, if disturbances of long wavelength are excluded, it is possible for the three-dimensional modes to be the least stable one. Accordingly, the pattern that will be observed on the flame front, at the onset of instability, will consist of either ridges in the direction of stretch or the more common three-dimensional cellular structure.

Jackson, Thomas L.↗
NASA NTRS1993

A new look at rainfall fluctuations and scaling properties of spatial rainfall using orthogonal wavelets

It has been observed that the finite-dimensional distribution functions of rainfall cannot obey simple scaling laws due to rainfall intermittency (mixed distribution with an atom at zero) and the probability of rainfall being an increasing function of area. Although rainfall fluctuations do not suffer these limitations, it is interesting to note that very few attempts have been made to study them in terms of their self-similarity characteristics. This is due to the lack of unambiguous definition of fluctuations in multidimensions. This paper shows that wavelet transforms offer a convenient and consistent method for the decomposition of inhomogeneous and anisotropic rainfall fields in two dimensions and that the components of this decomposition can be looked at as fluctuations of the rainfall field. It is also shown that under some mild assumptions, the component fields can be treated as homogeneous and thus are amenable to second-order analysis, which can provide useful insight into the nature of the process. The fact that wavelet transforms are a space-scale method also provides a convenient tool to study scaling characteristics of the process. Orthogonal wavelets are used, and these properties are investigated for a squall-line storm to study the presence of self-similarity.

Kumar, Praveen↗
DOE OSTI2025

Spectroscopy of the fractal Hofstadter energy spectrum

Hofstadter’s butterfly, the predicted energy spectrum for non-interacting electrons confined to a two-dimensional lattice in a magnetic field, is one of the most remarkable fractal structures in nature. At rational ratios of magnetic flux quanta per lattice unit cell, this spectrum shows self-similar distributions of energy levels that reflect its recursive construction. For most materials, Hofstadter’s butterfly is predicted under experimental conditions that are unachievable using laboratory-scale magnetic fields. More recently, electrical transport studies have provided evidence for Hofstadter’s butterfly in materials engineered to have artificially large lattice constants, such as those with moiré superlattices. Yet, so far, direct spectroscopy of the fractal energy spectrum predicted by Hofstadter nearly 50 years ago has remained out of reach. Here, in this work, we use high-resolution scanning tunnelling microscopy/spectroscopy (STM/STS) to investigate the flat electronic bands in twisted bilayer graphene (TBG) near the predicted second magic angle, an ideal setting for spectroscopic studies of Hofstadter’s spectrum. Our study shows the fractionalization of flat moiré bands into discrete Hofstadter subbands and discerns experimental signatures of self-similarity of this spectrum. Moreover, our measurements uncover a spectrum that evolves dynamically with electron density, showing phenomena beyond that of Hofstadter’s original model owing to the combined effects of strong correlations, Coulomb interactions and the quantum degeneracy of electrons in TBG.

electronic properties and materials↗
NASA NTRS1989

Cross correlation and length scales in turbulent flows near surfaces

Two kinds of length scales are used in turbulent flows; 'functional length scales' such as mixing length, dissipation length L(sub epsilon), etc., and 'flow-field length scales' derived from cross correlations of velocity, pressure, etc. in the flow. Some connection between these scales are derived here. We first consider the cross correlation R(sub vv)(y,y(sub 1)) of the normal components u at two heights y, y(sub 1) above a rigid surface, normalized by the velocity y(sub 1) (greater than y). For shear-free boundary layers it is found theoretically, and in field and numerical experiments that R(sub vv) approximately equals y/y(sub 1). For shear layers it is also found that R(sub vv) approximately equals f(y/y(sub 1)) less than or equal to y,y(sub 1). This function f differs slightly between low Reynolds number numerical simulations and field experiments. The lateral structure defined by R(sub vv)(y,r(sub 3); y(sub 1),0) is also self similar and shows that the eddies centered at about y(sub 1) appear to have constant lateral width a(sub 3) above and below y(sub 1), where a(sub 3, sup +) approximately equals 7+1/(1.4dU(sup +)/dy(sup +)), when normalized on u(sub *) and v, where U is the mean velocity. Results for L(sub epsilon, sup -1) from direct numerical simulation are found to compare well with the formula L(sub epsilon, sup -1) = A(sub B)/y + A(sub S)dU/dy/v, for unidirectional and reversing turbulent boundary layers and channel flow, except near where dU/dy approximately equals 0. The conclusion is that the large-scale eddy structure and length scales in these flows are determined by a combination of shear and blocking, and that the vertical component of turbulence has a self-similar structure in both kinds of boundary layer.

Hunt, J. C. R.↗
NASA NTRS1993

Supersonic shock wave/vortex interaction

Although shock wave/vortex interaction is a basic and important fluid dynamics problem, very little research has been conducted on this topic. Therefore, a detailed experimental study of the interaction between a supersonic streamwise turbulent vortex and a shock wave was carried out at the Penn State Gas Dynamics Laboratory. A vortex is produced by replaceable swirl vanes located upstream of the throat of various converging-diverging nozzles. The supersonic vortex is then injected into either a coflowing supersonic stream or ambient air. The structure of the isolated vortex is investigated in a supersonic wind tunnel using miniature, fast-response, five-hole and total temperature probes and in a free jet using laser Doppler velocimetry. The cases tested have unit Reynolds numbers in excess of 25 million per meter, axial Mach numbers ranging from 2.5 to 4.0, and peak tangential Mach numbers from 0 (i.e., a pure jet) to about 0.7. The results show that the typical supersonic wake-like vortex consists of a non-isentropic, rotational core, where the reduced circulation distribution is self similar, and an outer isentropic, irrotational region. The vortex core is also a region of significant turbulent fluctuations. Radial profiles of turbulent kinetic energy and axial-tangential Reynolds stress are presented. The interactions between the vortex and both oblique and normal shock waves are investigated using nonintrusive optical diagnostics (i.e. schlieren, planar laser scattering, and laser Doppler velocimetry). Of the various types, two Mach 2.5 overexpanded-nozzle Mach disc interactions are examined in detail. Below a certain vortex strength, a 'weak' interaction exists in which the normal shock is perturbed locally into an unsteady 'bubble' shock near the vortex axis, but vortex breakdown (i.e., a stagnation point) does not occur. For stronger vortices, a random unsteady 'strong' interaction results that causes vortex breakdown. The vortex core reforms downstream of the rear stagnation point, and the reduced circulation distribution once again becomes self-similar in this region. A-new model of this interaction is proposed. Finally, a curve defining the approximate limits of supersonic vortex breakdown is presented.

Settles, G. S.↗
DOE OSTI2022

Combining Planck and SPT Cluster Catalogs: Cosmological Analysis and Impact on the Planck Scaling Relation Calibration

We provide the first combined cosmological analysis of the South Pole Telescope (SPT) and Planck cluster catalogs. The aim is to provide an independent calibration for Planck scaling relations, exploiting the cosmological constraining power of the SPT-SZ cluster catalog and its dedicated weak lensing (WL) and X-ray follow-up observations. We build a new version of the Planck cluster likelihood. In the νΛ CDM scenario, focusing on the mass slope and mass bias of Planck scaling relations, we find ${\alpha }_{\mathrm{SZ}}={1.49}_{-0.10}^{+0.07}$ and ${\left(1-b\right)}_{\mathrm{SZ}}={0.69}_{-0.14}^{+0.07}$, respectively. The results for the mass slope show a ~4 σ departure from the self-similar evolution, α SZ ~ 1.8. This shift is mainly driven by the matter density value preferred by SPT data, Ωm = 0.30 ± 0.03, lower than the one obtained by Planck data alone, ${{\rm{\Omega }}}_{m}={0.37}_{-0.06}^{+0.02}$. The mass bias constraints are consistent both with outcomes of hydrodynamical simulations and external WL calibrations, (1 – b) ~ 0.8, and with results required by the Planck cosmic microwave background cosmology, (1 – b) ~ 0.6. From this analysis, we obtain a new catalog of Planck cluster masses M 500 . We estimate the ratio between the published Planck M SZ masses and our derived masses M 500 , as a "measured mass bias," ${\left(1-b\right)}_{M}$. We analyze the mass, redshift, and detection noise dependence of ${\left(1-b\right)}_{M}$, finding an increasing trend toward high redshift and low mass. These results mimic the effect of departure from self-similarity in cluster evolution, showing different dependencies for the low-mass, high-mass, low-z, and high-z regimes.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

The Roles of Large- and Small-Scale Structures in Turbulent Free Shear Layers

The occurrence of large-scale coherent structures in turbulent free shear flows (especially the planar mixing layer) has been recognized for some time. Indeed, the observation of such structures in mixing layers did much to promote interest in the study of coherent structures in turbulence. It has been widely assumed that the large-scale structures in these flows are responsible for the entrainment of free-stream fluid and the overall growth of the layer, while the small-scale structures provide mixing and dissipation. A model of scalar mixing based on these ideas was proposed for these flows. However, recent experimental and computational evidence suggests that the dominance of the large-scale structures in turbulent mixing layers is not universal. In addition, there is a substantial variation among experiments in several statistical measures of self-similar mixing layers, for example growth rate and velocity variances. To investigate the importance of large-scale structures, several free shear flows (mixing layers and wakes) have been simulated via direct numerical simulation. The simulations are designed to mimic experimental mixing layers in which the splitter plate boundary layers are turbulent. Different levels of two-dimensional forcing are included resulting in large-scale structures of differing strength and importance. These simulations are used to investigate the role of large-scale coherent structures in free shear layers and the effect of these structures on relevant turbulence statistics and scalar mixing. It is found that the statistics and structures in several experiments involving turbulent mixing layers are in better agreement with simulations that do not exhibit dominant large-scale structures than those in which the common mixing layer structures do dominate. It is also found that the level of forcing can have a profound effect on the qualitative and quantitative features of these shear layer, even when they are nominally self-similar.

Moser, Robert D.↗
DOE OSTI2023

Rheo-electric measurements of carbon black suspensions containing polyvinylidene difluoride in N -methyl-2-pyrrolidone

Lithium-ion battery cathode slurries have a microstructure that depends sensitively on how they are processed due to carbon black's (CB) evolving structure when subjected coating flows. While polyvinylidene difluoride (PVDF), one of the main components of the cathode slurry, plays an important role in modifying the structure and rheology of CB, a quantitative understanding is lacking. In this work, we explore the role of PVDF in determining the structural evolution of Super C65 CB in N-methyl-2-pyrrolidinone (NMP) with rheo-electric measurements. We find that PVDF enhances the viscosity of NMP resulting in a more extensive structural erosion of CB agglomerates with increasing polymer concentration and molecular weight. We also show that the relative viscosity of all suspensions can be collapsed by the fluid Mason number (Mnf), which compares the hydrodynamic forces imposed by the medium to cohesive forces holding CB agglomerates together. Using simultaneous rheo-electric measurements, we find at high Mnf, the dielectric strength (Δε) scales with Mnf, and the power-law scaling can be quantitatively predicted by considering the self-similar break up of CB agglomerates. The collapse of the relative viscosity and scaling of Δε both suggest that PVDF increases the hydrodynamic force of the suspending medium without directly changing the CB agglomerate structure. These findings are valuable for optimizing the rheology of lithium ion battery cathode slurries. We also anticipate that these findings can be extended to understand the microstructure of similar systems under flow.

Mechanics↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗
NASA NTRS2009

Evolvable Neural Software System

The Evolvable Neural Software System (ENSS) is composed of sets of Neural Basis Functions (NBFs), which can be totally autonomously created and removed according to the changing needs and requirements of the software system. The resulting structure is both hierarchical and self-similar in that a given set of NBFs may have a ruler NBF, which in turn communicates with other sets of NBFs. These sets of NBFs may function as nodes to a ruler node, which are also NBF constructs. In this manner, the synthetic neural system can exhibit the complexity, three-dimensional connectivity, and adaptability of biological neural systems. An added advantage of ENSS over a natural neural system is its ability to modify its core genetic code in response to environmental changes as reflected in needs and requirements. The neural system is fully adaptive and evolvable and is trainable before release. It continues to rewire itself while on the job. The NBF is a unique, bilevel intelligence neural system composed of a higher-level heuristic neural system (HNS) and a lower-level, autonomic neural system (ANS). Taken together, the HNS and the ANS give each NBF the complete capabilities of a biological neural system to match sensory inputs to actions. Another feature of the NBF is the Evolvable Neural Interface (ENI), which links the HNS and ANS. The ENI solves the interface problem between these two systems by actively adapting and evolving from a primitive initial state (a Neural Thread) to a complicated, operational ENI and successfully adapting to a training sequence of sensory input. This simulates the adaptation of a biological neural system in a developmental phase. Within the greater multi-NBF and multi-node ENSS, self-similar ENI s provide the basis for inter-NBF and inter-node connectivity.

Curtis, Steven A.↗
NASA NTRS1999

Influence of Buoyant Convection on the Stability of Enclosed Laminar Flames

An investigation of the stability limits of Enclosed Laminar Flames (ELF) was conducted in the Middeck Glovebox (MGBX) facility on the STS-87 Space Shuttle mission (November to December 1997). The primary objective of the ELF glovebox investigation is to determine the effect of buoyancy on the stability of round, laminar, gas-jet diffusion flames in a co-flow air duct. Comparison tests were conducted in normal gravity to allow isolation and identification of the influence of buoyancy. The results were used to map the lift-off and blow-out stability limits as a function of the fuel and air velocities for the two buoyancy conditions. Approximately 50 tests were conducted during the Space Shuttle mission, using a 50/50 mixture (volume basis) of methane and nitrogen as the fuel. The experimental results verified the hypothesis that substantially greater velocities are required to destabilize the flame in microgravity. The increase in air velocity required to induce lift off in microgravity (compared to normal gravity) was nearly equal to the increase required to induce blow out. Furthermore, the air velocity increase was relatively independent of the fuel flow, except at low fuel flows. At high fuel flows, it was found that the microgravity flames tend to immediately blow out after lift off. This is in agreement with the free-jet theory which suggests that stable lifted flames are not possible for fuels with a Schmidt number of 0.5<Sc<l, such as the diluted methane used in this study. However, stable lifted flames were observed in microgravity at low fuel flows, and in normal gravity, despite the fuel's low Schmidt number. The discrepancy between the experimental data and the theory of Chung and Lee is presumably due to the differences in the flow condition; namely, the ducted flow versus a free jet, and a non-similar velocity profile versus a self-similar profile.

Brooker, John E.↗
DOE OSTI2020

Nanoscale morphological evolution of monocrystalline Pt surfaces during cathodic corrosion

This paper studies the cathodic corrosion of a spherical single crystal of platinum in an aqueous alkaline electrolyte, to map out the detailed facet dependence of the corrosion structures forming during this still largely unexplored electrochemical phenomenon. We find that anisotropic corrosion of the platinum electrode takes place in different stages. Initially, corrosion etch pits are formed, which reflect the local symmetry of the surface: square pits on (100) facets, triangular pits on (111) facets, and rectangular pits on (110) facets. We hypothesize that these etch pits are formed through a ternary metal hydride corrosion intermediate. In contrast to anodic corrosion, the (111) facet corrodes the fastest, and the (110) facet corrodes the slowest. For cathodic corrosion on the (100) facet and on higher-index surfaces close to the (100) plane, the etch pit destabilizes in a second growth stage, by etching faster in the (111) direction, leading to arms in the etch pit, yielding a concave octagon-shaped pit. In a third growth stage, these arms develop side arms, leading to a structure that strongly resembles a self-similar diffusion-limited growth pattern, with strongly preferred growth directions.

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

Nuclear response of additive manufactured GRCop-84 copper for use in Lower hybrid launchers in a fusion environment

Recent advances in selective laser melting (SLM) 3D printing technology allow additive manufacture of lower hybrid current drive (LHCD) RF launchers from a new material, Glenn Research Copper 84 (GRCop-84), a Cr 2 Nb (8 at. % Cr, 4 at. % Nb) precipitation hardened alloy, in configurations unachievable with conventional machining. Tensile strength and thermal conductivity are compared to other competing high conductivity copper alloys for fusion use. Swelling and conductivity at 100 DPA when exposed to the Hanford Fast Flux Test Facility neutron spectrum is predicted in the 1–2 % and 65-70% IACS range, respectively, by comparison to copper alloys with similar precipitate and grain size. Finally, self-similar ion irradiation of a GRCop-84 target at 430 °C to 20 DPA resulted in no noticeable void formation in FIB cross sections or TEM lamella.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗
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↗