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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2021

Self-Similarity, Fractality and Entropy Principle in Collisions of Hadrons and Nuclei at Tevatron, RHIC and LHC

$z$-Scaling of inclusive spectra as a manifestation of self-similarity and fractality of hadron interactions is illustrated.The scaling for negative particle production in $Au+Au$ collisionsfrom BES-I program at RHIC is demonstrated. The scaling variable $z$ depends on the momentum fractions of the colliding objects carried by the interacting constituents, and momentum fractions of the produced objects in the scattered and recoil directions, which carries the inclusive particle and its counterpart, respectively. Structures of the colliding objects and fragmentation processes in final stateare expressed by fractal dimensions. Medium produced in the collisions is described by a specific heat. The scaling function $\psi(z)$ reveals energy, angular, multiplicity, and flavor independence.It has a power behavior at high $z$ (high $p_T$). Based on entropy principle and $z$-scaling, energy loss as a function of the collision energy, centrality and transverse momentum of inclusive particle is estimated.New conservation law including fractal dimensions is found. Quantization of fractal dimensions is discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Scalar flux transport models for self-similar turbulent mixing

A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.

Reynolds-averaged Navier Stokes↗
DOE OSTI2021

Self-similar mesocrystals form via interface-driven nucleation and assembly

Crystallization by particle attachment (CPA) is a frequently occurring mechanism of colloidal crystallization that results in hierarchical morphologies. CPA has been exploited to create nanomaterials with unusual properties and is implicated in the development of complex mineral textures. Oriented attachment–a form of CPA in which particles align along specific crystallographic directions–produces mesocrystals that diffract as single crystals do, although the constituent particles are still discernible. The conventional view of CPA is that nucleation provides a supply of particles that aggregate via Brownian motion biased by attractive interparticle potentials. However, mesocrystals often exhibit regular morphologies and uniform sizes. Although many crystal systems form mesocrystals and individual attachment events have been directly visualized10, how random attachment events lead to well defined, self-similar morphologies remains unknown, as does the role of surface-bound ligands, which are ubiquitous in nanoparticle systems. Attempts to understand mesocrystal formation are further complicated in many systems by the presence of precursor nanoparticles with a phase distinct from that of the bulk. Some studies propose that such particles convert before attachment, whereas others attribute conversion to the attachment process itself and yet others conclude that transformation occurs after the mesocrystals exceed a characteristic size. Here we investigate mesocrystal formation by iron oxides, which are important colloidal phases in natural environments and classic examples of systems forming ubiquitous precursor phases and undergoing CPA accompanied by phase transformations. Combining in situ transmission electron microscopy (TEM) at 80 degrees Celsius with ‘freeze-and-look’ TEM, we tracked the formation of haematite (Hm) mesocrystals in the presence of oxalate (Ox), which is abundant in soils, where iron oxides are common. In this work, we find that isolated Hm particles rarely appear, but once formed, interfacial gradients at the Ox-covered surfaces drive Hm particles to nucleate repeatedly about two nanometres from the surfaces, to which they then attach, thereby generating mesocrystals. Comparison to natural and synthetic systems suggests that interface-driven pathways are widespread.

36 MATERIALS SCIENCE↗
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 OSTI2022

Self-similarity of scalar isosurface area density in a temporal mixing layer

Understanding the behaviour of scalar isosurfaces in a turbulent flow is of particular interest for many problems in turbulent mixing that contain sharp interfaces between regions of the flow. Common examples include combustion, where the chemical reactions occur in thin regions within the flow, and the turbulent/non-turbulent interface in shear flows, where a thin region separates the rotational, turbulent motions from the irrotational, non-turbulent background. Recent advances in computing technology allow for in-depth analysis of these interface problems that are difficult to quantify in a laboratory setting. In this paper, the results of a direct numerical simulation of a passive scalar $\varPhi$ evolving on a turbulent, temporally developing mixing layer are described. A novel approach has been taken to calculate the surface area of individual scalar isosurfaces, $A_{iso}$ , throughout the simulation, as well as the mean isosurface area density, $\varSigma$ , as a function of the cross-stream direction and time. A notable finding is that the profiles of $\varSigma$ develop in a self-similar manner when scaled by the Taylor scale of the scalar field, $\lambda _\phi$ . Remarkably, the scaling appears to hold for a wide range of isovalues. A rough scaling argument based on the formal definition of $\varSigma$ and properties of a temporal mixing layer is presented which also exposes a dependence on $\lambda _\phi$ . Based on these results, a possible scaling for the isosurface area is presented as $A_{iso}/A_0 \sim (Re \, Sc)^{1/2}$ , where $Re$ and $Sc$ are local Reynolds and Schmidt numbers, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Gaps labeling theorem for the bubble-diamond self-similar graphs

Abstract Motivated by the appearance of fractals in several areas of physics, especially in solid state physics and the physics of aperiodic order, and in other sciences, including the quantum information theory, we present a detailed spectral analysis for a new class of fractal-type diamond graphs, referred to as bubble-diamond graphs, and provide a gap-labeling theorem in the sense of Bellissard for the corresponding probabilistic graph Laplacians using the technique of spectral decimation. Labeling the gaps in the Cantor set by the normalized eigenvalue counting function, also known as the integrated density of states, we describe the gap labels as orbits of a second dynamical system that reflects the branching parameter of the bubble construction and the decimation structure. The spectrum of the natural Laplacian on limit graphs is shown generically to be pure point supported on a Cantor set, though one particular graph has a mixture of pure point and singularly continuous components.

Physics↗
DOE OSTI2021

Self-similarity of k -nearest neighbour distributions in scale-free simulations

ABSTRACT We use the k-nearest neighbour probability distribution function (kNN-PDF; Banerjee & Abel 2021a) to assess convergence in a scale-free N-body simulation. Compared to our previous two-point analysis, the kNN-PDF allows us to quantify our results in the language of haloes and numbers of particles, while also incorporating non-Gaussian information. We find good convergence for 32 particles and greater at densities typical of haloes, while 16 particles and fewer appear unconverged. Halving the softening length extends convergence to higher densities, but not to fewer particles. Our analysis is less sensitive to voids, but we analyse a limited range of underdensities and find evidence for convergence at 16 particles and greater even in sparse voids.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2023

On the evolution of scalar iso-surface area density in a turbulent mixing layer

In this paper, the transport equation for scalar iso-surface area density (Σ) in a turbulent, temporally developing mixing layer is examined. Exploring the spatial and temporal evolution of the terms in the Σ transport equation is vital to improving our understanding of turbulent flows characterized by distinct interfaces, e.g. the flame surface or the turbulent/non-turbulent interface. Previous work reported by the authors identified that Σ exhibits self-similar behaviour consistent with the development of the temporal mixing layer. Accordingly, each of the terms in the Σ transport equation is found to behave in a self-similar manner, though there are notable differences in the self-similar behaviours for each term. Based on the results presented herein, it is suggested that the rate of change of Σ (∂Σ/∂t) and the advection term scale with hλ Φ /ΔU, where h is the width of the mixing layer, λ Φ is the scalar Taylor length scale and ΔU is the velocity difference. The production and destruction terms are found to scale with an additional factor (ReSc) 1/2 . In contrast, the molecular diffusion term is found to scale with a factor (ReSc) –1/2 compared to ∂Σ/∂t. Importantly, it is found that the difference between the production and destruction terms, or net surface ‘stretch’, scales with the same factor as ∂Σ/∂t and the advection term, which may have a significant impact on how the evolution of Σ is understood and modelled in turbulent flows.

42 ENGINEERING↗
DOE OSTI2020

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Optimal white-noise stochastic forcing for linear models of turbulent channel flow

In the present study an optimisation problem is formulated to determine the forcing of an eddy-viscosity-based linearised Navier–Stokes model in channel flow at $Re_\tau \approx 5200$ ( $Re_\tau$ is the friction Reynolds number), where the forcing is white-in-time and spatially decorrelated. The objective functional is prescribed such that the forcing drives a response to best match a set of velocity spectra from direct numerical simulation (DNS), as well as remaining sufficiently smooth. Strong quantitative agreement is obtained between the velocity spectra from the linear model with optimal forcing and from DNS, but only qualitative agreement between the Reynolds shear stress co-spectra from the model and DNS. The forcing spectra exhibit a level of self-similarity, associated with the primary peak in the velocity spectra, but they also reveal a non-negligible amount of energy spent in phenomenologically mimicking the non-self-similar part of the velocity spectra associated with energy cascade. By exploiting linearity, the effect of the individual forcing components is assessed and the contributions from the Orr mechanism and the lift-up effect are also identified. Finally, the effect of the strength of the eddy viscosity on the optimisation performance is investigated. The inclusion of the eddy viscosity diffusion operator is shown to be essential in modelling of the near-wall features, while still allowing the forcing of the self-similar primary peak. In particular, reducing the strength of the eddy viscosity results in a considerable increase in the near-wall forcing of wall-parallel components.

Mechanics↗