Search NASA⌕ Search

SEARCH · Search NASA

Results for “self-similar”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Optimization techniques in self-similar compressible flow

We investigate the one-dimensional (1D) inviscid compressible flow equations for an ideal gas through the lens of optimization techniques. It is the case that, to our knowledge, optimization analysis applied to the so-called “linear velocity” solutions of the Euler compressible flow equations has not been previously conducted. Through both gradient-based and variational techniques, new variants of well-studied flow scenarios, i.e., self-similar, 1D, linear velocity solution class to idealized inviscid compressible flow equations, are determined, as encoded in both the kinematic and thermodynamic properties of this self-similar solution class. With the kinematics of the said solutions being driven by a self-similar “scale radius” and the thermodynamics being driven separately through the appearance of an arbitrary function, a myriad of new solution classes is possible. Acting as a guide to more realistic physical circumstances as well as discovery, it is the hope that the presented cases serve as the framework for future investigations into the intersection of self-similarity and optimization techniques. Fields of study that may find this work to be of interest include aerodynamic design, flow control, inertial confinement fusion, physics-informed neural networks, and other related areas of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Self-similar Cosmic-Ray Transport in High-resolution Magnetohydrodynamic Turbulence

We study the propagation of cosmic rays (CRs) through a simulation of magnetohydrodynamic (MHD) turbulence at an unprecedented resolution of 10,240 3 . We drive turbulence that is subsonic and super-Alfvénic, characterized by δB rms /B 0 = 2. The high resolution enables an extended inertial range such that the Alfvén scale l A , where δB(l A ) ≈ B 0 , is well resolved. This allows us to properly capture how the cascade transitions from large amplitudes on large scales to small amplitudes on small scales. We find that sharp bends in the magnetic field are key mediators of particle transport even on small scales via resonant curvature scattering. We further find that particle scattering in the turbulence shows strong hints of self-similarity: (1) the diffusion has a weak energy dependence over almost two decades in particle energy, and (2) the particles’ random walk exhibits a broad power-law distribution of collision times such that the diffusion is dominated by the rarest, long-distance excursions. Our results suggest that large-amplitude MHD turbulence can provide efficient scattering over a wide range of CR energies and may help explain many CR observations above a ∼TeV: the flattening of the B/C spectrum, the hardening of CR primary spectra, and the weak dependence of arrival anisotropy on CR energy.

Kempski, Philipp [Princeton Univ., NJ (United Stat↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
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 OSTI2025

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY↗
DOE OSTI2026

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

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗
DOE OSTI2025

The effect on soot and its gas precursors of doping ethylene with 2,2,4,6,6-pentamethyl-heptane in the nitrogen-fuel stream of a laminar non-premixed Planar Mixing Layer Flame (PMLF)

Synthetic Aviation Turbine Fuels (SATFs) are promising for reducing soot emissions from the aviation sector and diversifying Jet Fuel (JF) sources. Accurately predicting the combustion and emissions behavior of SATFs (and other JFs) necessitates robust experimental databases to elucidate the chemistry of long-chain iso-paraffins, which can compose up to two-thirds of SATF blends and whose behavior is considered to be well-represented by that of iso-dodecane isomers. Here, this study characterizes two laminar non-premixed Planar Mixing Layer Flames (PMLFs) with mild soot loads fueled by nitrogen-diluted ethylene, pure and doped with 2,2,4,6,6-pentamethyl-heptane, respectively. The two PMLFs have the same stoichiometric mixture fraction and total hydrocarbon mole fraction in the fuel stream (X F,F =X C2H4,F +X C12H26,F = 0.260), resulting in nearly the same maximum temperature (T max ≈1800 K) and simple identification of the effects of doping. Importantly, any horizontal PMLF cross-section has a self-similar structure that can be modeled as an equivalent One-Dimensional Counterflow Flame (1D-CF) with vanishingly small strain rate (a). The cross-section at a Height Above the Burner (HAB) of 50 mm is characterized in terms of C 0 -C 18 gas species using capillary sampling followed by GC-MS analyses. Laser-Induced Emission Spectroscopy (LIES) quantifies the soot volume fraction (ƒ v ) profiles at HAB=25 and 50 mm where Elastic Laser Light Scattering (E-LLS) is performed to determine the a of the equivalent 1D-CFs and the profile of the E-LLS equivalent diameter ( d 6,3 ) of soot. The substitution of 1500 ppm of ethylene with 2,2,4,6,6-pentamethyl-heptane causes an increase of ≈1.5 in the concentrations of several polycyclic aromatic hydrocarbons and fv. Concurrently, the measured d 6,3 doubles in the oxidizer stream, yet remains the same in the fuel stream, at HAB=50 mm. Instead, at HAB= 25 mm, the iso-dodecane doping does not affect the d 6,3 profile in either stream. The experimental results partially validate the chemical reactions and soot formation kinetic model developed at Lawrence Livermore National Laboratory and provide directions to further improve its predictions.

Iso-dodecane (2,2,4,6,6-Pentamethyl-Heptane)↗
DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Scalable Hot-Water-Repellent Superhydrophobicity via Thermal Insulation

Superhydrophobic surfaces, which rely on a combination of surface texture and chemistry, often lose their repellent behavior when contacted by hot water (≳40 °C) because the impinging hot water replaces the requisite air layer within the surface texture via evaporation and recondensation. In contrast to previous approaches targeting this condensation-induced failure mode that rely on intricately tailored surface structures or complex chemical treatments, we present a scalable approach based on thermal design: the multilayered insulated superhydrophobic (MISH) coating mitigates condensation-induced failure by preventing heat transfer. Superhydrophobicity is retained at impinging water temperatures up to 90 °C, with durability demonstrated via long-term (>1 million impacts) droplet impingement experiments. We explain the mechanism for this approach with a detailed thermal model; the model reveals that the underlying physical behavior is self-similar across coating parameters and impinging fluid temperatures. Finally, the MISH coating accommodates curved geometries and large surfaces, and it is over 4 orders of magnitude less expensive than cleanroomnanofabricated alternatives, indicating promise for practical use in the energy sector, chemical processing, and the food and medical industries.

coating materials↗
DOE OSTI2026

Focusing of laser-accelerated ions with spherically curved targets

A theoretical model for laser-accelerated ion focusing from spherically curved targets using a Gaussian self-similar solution is presented. This model describes the evolution of the focal location and focal radius with final ion energy and the initial radial boundary of the ions, the latter being a function of both target geometry and the ion acceleration radius for a given ion energy. The theory is supported by particle-in-cell simulations of a variety of target shapes with varying radii of curvature and target opening angles, as well as variations in the injected electron beam radius and energy spectrum. The theory and simulations suggest that the focal location varies linearly with the radius of curvature, with the square root of the ratio of ion energy to effective electron temperature, and monotonically increases with the energy-dependent ratio of the initial ion radial boundary to the radius of curvature. Considering ponderomotive acceleration, this suggests that the focal length should scale inversely with intensity, $(I_Lλ_L^2)$ −1/4 , suggesting a 10× increase in intensity will reduce the focal distance by ̃1.8× for a given ion energy.

Physics - Plasma physics↗
DOE OSTI2026

Out-of-time-order correlators bridge classical transport and quantum dynamics

The out-of-time-order correlator (OTOC) has emerged as a central tool for quantifying decoherence across wide-ranging physical platforms. Here, we demonstrate its direct measurement in a classical ensemble using nuclear magnetic resonance with a modulated gradient spin echo sequence and extend the method into a multidimensional correlation to track exchange phenomena. Position is encoded through magnetic field gradients and momentum through the velocity autocorrelation function, enabling experimental access to OTOCs for proton motion confined within the self-similar lattice of the metal–organic framework MOF-808. Here, water confined to specified geometries within the MOF pores gives rise to spatially distinct diffusive eigenmodes with characteristic relative entropies. We demonstrate that periodic radio frequency driving combined with gradient modulation yields entropy evolution through the selection of distinct diffusion modes. Frequency-resolved diffusion spectra connect these entropy dynamics to classical heat exchange laws, revealing how operational features of quantum systems are mirrored in confined, macroscopic spin ensembles.

Fricke, Sophia N. [University of California, Berke↗
DOE OSTI2024

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science↗
DOE OSTI2026

Photon phase-space dynamics in a plasma wakefield accelerator

Frequency up-shifting of laser light in a beam-driven plasma wakefield has the potential to provide high-intensity sources of short wavelength radiation. Simulations have demonstrated that a laser pulse can undergo large frequency shifts, limited only by the drive beam energy, when the plasma density is tailored to match the accelerating phase of the wake to the group velocity of the pulse. Here, we study the dynamical evolution of photons in the 1D1P phase-space vicinity of the plasma wake-phase matching condition. Numerical calculations using a photon kinetic model are validated by direct comparison with 1D and quasi-3D particle-in-cell simulations. These calculations form the basis of a linear theory of the photon dynamics which reveals several important results, including scalings for the properties of the witness pulse and a self-similar solution for the photon phase-space dynamics. One prediction of the analytic theory is that the pulse can be compressed indefinitely, though the amount of compression would ultimately be limited by practical constraints. These results suggest that photon acceleration can provide a novel source of sub-femtosecond, short wavelength radiation.

XUV generation↗
DOE OSTI2026

Quantum critical collapse abhors a naked singularity

Classical critical collapse yields naked singularities from smooth initial data, challenging cosmic censorship, and shaping the spectrum of primordial black holes. We show that one-loop vacuum polarization near the threshold qualitatively changes this outcome by dressing the singularity with a horizon within a controlled semiclassical regime. In analytically tractable Einstein-scalar critical spacetimes, a one-loop $s$-wave treatment linearized around self-similar backgrounds shows that regularity uniquely selects an asymptotically Minkowskian, vacuum-polarization state. Its renormalized stress tensor carries a universal quantum growing mode that competes with the classical unstable mode, shifts the critical point, and generates a trapped surface along with a finite mass gap at the new threshold, thereby enforcing horizon formation even under arbitrary fine-tuning. In primordial collapse, the threshold shift enters exponentially into the formation fraction, while the mass gap truncates the low-mass tail, suggesting potentially important consequences for the predicted mass spectrum. Furthermore, these results provide a self-consistent semiclassical treatment of critical collapse and yield sharp predictions within the one-loop, near-critical, linearized regime.

Anomalies↗
DOE OSTI2026

Propagation of partially spatially coherent laser beams in instantaneous Kerr media

The propagation of intense, partially spatially coherent laser beams in a medium with instantaneous third-order susceptibility is studied analytically and numerically. For sufficiently high power relative to that required for nonlinear self-focusing, the propagation initially proceeds in two stages. In the first stage, spatial coherence builds up, and in the second stage, the number of speckles reduces. Once the degree of coherence is sufficiently high, whole-beam self-focusing occurs. The beam power is mostly confined within the initial spot radius. Two analytical approaches for describing the evolution of the beam are presented. The method of moments leads to an analytical solution for the rms spot radius that is in excellent agreement with simulations. This method does not require any knowledge of the field statistics beyond the initial conditions and provides no information about the evolution of the individual speckles. The other approach employs a self-similar solution for the second-order coherence function of the field and assumes that the fourth-order coherence function is factorizable and obeys complex circular Gaussian random statistics. The latter method also leads to an analytical expression for the spot radius, but its predictions for the qualitative evolution of the speckles disagree with wave-optics simulations.

lasers↗
DOE OSTI2026

Evaluating the Limits of QAOA Parameter Transfer at High-Rounds on Sparse Ising Models With Geometrically Local Cubic Terms

The emergent practical applicability of the Quantum Approximate Optimization Algorithm (QAOA) for approximate combinatorial optimization is a subject of considerable interest. One of the primary limitations of QAOA is the task of finding a set of good parameters, which is usually done using a variational optimization loop. Parameter transfer, or parameter concentration, is a phenomenon where QAOA angles trained on problem instances that are self-similar tend to perform well for other problem instances from that similar class. This suggests a potentially highly efficient and scalable non-variational learning method for QAOA angle finding. In this work, we systematically study QAOA parameter transferability from small problem sizes (16 and 27 decision variables) onto large problem instances (up to 156 qubits) for heavy-hex graph Ising models with geometrically local higher order terms using the Julia based QAOA simulation tool \texttt{JuliQAOA} to perform classical angle finding for up to $49$ QAOA layers ($p$). Parameter transfer of the fixed angles is validated using a combination of full statevector, Projected Entangled Pair States (PEPS), Matrix Product State (MPS), and LOWESA numerical simulations. We find that the QAOA parameter transfer from single instances applied to other (unseen) problem instances does not in general provide monotonically improving performance as a function of $p$ - there are many cases where the performance temporarily decreases as a function of $p$ - but despite this the transferred angles have a general trend of improved expectation value as the QAOA depth increases, in many cases converging close to the true ground-state energy of the $100+$ qubit instances. We also sample the hardware-compatible Ising models using the ensemble of transfer-learned QAOA parameters on several superconducting qubit IBM Quantum processors with 127, 133, and 156 qubits. We find continuous solution quality improvement of the hardware-compatible QAOA circuits run on the IBM NISQ processors up to $p=5$ on \texttt{ibm\_fez}, up to $p=9$ on \texttt{ibm\_torino}, and up to $p=10$ on \texttt{ibm\_pittsburgh}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗